View Single Post
Old 06-08-2020, 12:40 PM   #1
drdna
Member
 
Location: Kentucky

Join Date: May 2012
Posts: 76
Default GATK calling false SNPs in duplicated regions

Is there a way to prevent GATK from calling SNPs in regions containing sequences that are duplicated in the reference genome? I'm having a problem with it calling false SNPs based on alignments of multi-mapped reads.

Tried running GenotypeGVCFs with the --read-filter NotSecondaryAlignmentReadFilter flag but this made no difference.
drdna is offline   Reply With Quote