How to get MD tag sequence of a read from a sam file given a bed file?
0
0
Entering edit mode
8 months ago
n.a.tekkey • 0

Hello! I am attempting to get a MD tag sequence from a read given a bed file with start and stop positions. For example, my sam file has a read like:

chr1    0   chr1    78142606    60  78278871S73M2D...

and I have a bed file with a line like:

chr1    78142736    78143236

I want to get the sequence from indexed position 78142736-78143236 so that I can calculate a percent identity score for that specific bin. I've attempted using samtools view but I can't seem to get it to do what I want it to. Please help!

samtools bed sam cigar • 377 views
ADD COMMENT
0
Entering edit mode

What was the exact options you used with samtools view and what were the results?

ADD REPLY

Login before adding your answer.

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