Question: SNP calling in Identical exons of two different genes
gravatar for Sudhir Jadhao
8 days ago by
Sudhir Jadhao60 wrote:

Hello Everyone,

I am calling germline mutation using BWA and GATK pipeline.

RHD and RHCE gene share 97% homology.

The exon2 and exon8 of these genes are identical. So there is a high chance that the reads from these regions may map to each other and call false negative SNPs

How can I differentiate the true mapping read for each gene?

How to deal with identical exon region of Homologs gene for variant calling?


Are these things taken care of by BWA and GATK pipeline

Thank you in advance

ADD COMMENTlink modified 8 days ago • written 8 days ago by Sudhir Jadhao60

GATK should normally use only primary alignment for SNP calling (or am I wrong ?)

ADD REPLYlink written 8 days ago by Nicolas Rosewick7.3k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1786 users visited in the last hour