Somatic Snp Mutation Igv Codon
1
1
Entering edit mode
11.5 years ago
Hmm ▴ 500

I am looking at a tumor/normal pair illumina HiSeq data paired end. When i open my normal and tumor bam files in IGV, i expect to see a mutation at a specific location(which i do see). But the base there is A. The mutation is suppose to change G -> T but what the reads show me only is A -> C and on the Reference sequence in IGV i see a "C" (so that is fine). How can I explain this to someone?

The codon i am looking at is "ACC" so the complement of that is "TGG"..right...why don't i see any T and G in my reads? Why are there only A and C?

somatic snp • 2.9k views
ADD COMMENT
2
Entering edit mode
11.5 years ago
Mitch Bekritsky ★ 1.3k

This may not generally be the case, but when I use BWA to align reads, it reports aligned reads relative to the original genome you used for alignment (i.e. the forward strands of all the chromosomes). Therefore, if a read maps to the reverse strand of a chromosome, BWA will report the reverse complement of the read you originally submitted as the sequence after alignment in the SAM/BAM file.

If your coding sequence is on the reverse strand of the original genome you used for alignment, then it will only report G->T, not A->C, since that is what the sequence was aligned to.

To double check, look for your gene in the genome browser. If it is transcribed from the reverse strand, then that could be why it's reporting the G->T substitution.

ADD COMMENT
0
Entering edit mode

Thank you..Youe explanation was really helpful!

ADD REPLY
0
Entering edit mode

It was my pleasure!

ADD REPLY

Login before adding your answer.

Traffic: 1993 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6