Entering edit mode
4.4 years ago
eholleman
•
0
I am aligning a consensus sequence to a reference genome using Bowtie2 and I am interested in the nucleotide sequence of the reference where the consensus sequence aligned. Is it possible to infer the reference sequence using only information in the SAM file such as the CIGAR and MD strings? Thanks!
I think OP only want the sequence from the BAM, not the REF.
I am looking for the sequence in the reference that the query aligned to
Thank you, this is what I am looking for. Would this be faster than building an indexed BLAST database and then using blastdbcmd to accomplish the same task?
Yes I would assume so