Extract nucleotide from read at SNP position
0
1
Entering edit mode
7.2 years ago
bdeonovic ▴ 210

I have read alignments in SAM format, I have read fastqs, I have genome fasta, and I have a VCF file with locations of SNPs. Is there an easy way to pull out the sequenced nucleotide for each read at the SNPs?

I'm not interested in calling SNPs (this is already done, I have my VCF) but I would like to know what nucleotide a read had at particular SNPs

RNA-Seq sam SNP • 2.9k views
ADD COMMENT
1
Entering edit mode

The following posts are relevant. There is also a python code in the first post that works with the pileup format.

ADD REPLY
0
Entering edit mode

SAM has read sequences for each alignment as well as the alignment position -- you can calculate offset into the read that corresponds to your SNP's location and look at the base

ADD REPLY
0
Entering edit mode

I am trying to do the same. Did you find a solution?

ADD REPLY

Login before adding your answer.

Traffic: 1613 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