To celebrate the release of GATK 4.0, we are giving away free credits for running the GATK4 Best Practices pipelines in FireCloud, our secure online analysis portal. It’s first come first serve, so sign up now to claim your free credits worth $250. Sponsored by Google Cloud. Learn more at https://software.broadinstitute.org/firecloud/documentation/freecredits

MNP detection using ReadBackedPhasing

GATK team,

I have used HaplotypeCaller to call variations in the exome data, performing joint calling on the samples. Then, ReadBackedPhasing was running on each sample individually with "-enableMergeToMNP" option. I was expecting to get the phased SNPs merged into MNPs. In the attached screenshot of one of the samples, I expect TG>CA, since I have not changed the default option of -maxDistMNP. This is the screen shot of phased VCF file output from ReadBackedPhasing.

19      1037715 rs200177867     T       C       68058.90        PASS    .       GT:AD:DP:GQ:HP:PGT:PID:PL       0/1:23,7:30:99:1037715-1,1037715-2:0|1:1037715_T_C:198,0,1019
19      1037716 rs201532581     G       A       68349.90        PASS    .       GT:AD:DP:GQ:HP:PGT:PID:PL:PQ    0/1:24,6:30:99:1037715-1,1037715-2:0|1:1037715_T_C:198,0,1019:1092.06
19      1037718 rs199741851     G       T       67457.30        PASS    .       GT:AD:DP:GQ:HP:PGT:PID:PL:PQ    0/1:24,6:30:99:1037715-1,1037715-2:0|1:1037715_T_C:177,0,1022:1100.75

Command used :

java -jar GATK \
    -T ReadBackedPhasing \
     -R EnsemblhgGRCh37_71.karyoorder.gatk.fa \
      -I file.bam \
      --variant f.recode.vcf \
      -o phased_f.vcf \
     -enableMergeToMNP \
     --phaseQualityThresh 20.0

Please let me know if the command needs additional parameters. Thanks.

Answers

Sign In or Register to comment.