Genome Analysis Toolkit

Variant Discovery in High-Throughput Sequencing Data

GATK process banner

Need Help?

Search our documentation

Community Forum

Hi, How can we help?

Developed in the Data Sciences Platform at the Broad Institute, the toolkit offers a wide variety of tools with a primary focus on variant discovery and genotyping. Its powerful processing engine and high-performance computing features make it capable of taking on projects of any size. Learn more

Using GATK ASE read counter with strelka

0

2 comments

  • Avatar
    Gökalp Çelik

    Hi Ramiro Barrantes

    Strelka VCF does not contain a GT field therefore each variant site is ignored. ASEReadCounter tool requires heterozygous marked sites for proper function. Mutect2 adds GT field as a compatibility measure and does not warrant any proper genotype information to be present at FORMAT/GT field. You may want to modify strelka vcf by adding FORMAT/GT field to each variant and add 0/1 to each biallelic site. 

    Regards. 

    1
    Comment actions Permalink
  • Avatar
    Ramiro Barrantes

    Thank you, this is exactly what I needed!!

    0
    Comment actions Permalink

Please sign in to leave a comment.

Powered by Zendesk