Question: how to get a read sequence from reference genome
0
gravatar for zhangdengwei
17 months ago by
zhangdengwei130
zhangdengwei130 wrote:

hi, If I have the coordinate and length of a read, how can I fastly get its sequence from the reference genome? Thanks!!!

sequence • 254 views
ADD COMMENTlink written 17 months ago by zhangdengwei130

How To Use Coordinates In Order To Extract Sequences In Fasta File?

ADD REPLYlink written 17 months ago by Nicolas Rosewick9.0k

Yes! I want to use the coordinates of each reads in a bam or sam file to extract its corresponding sequence on reference genome.

ADD REPLYlink written 17 months ago by zhangdengwei130

I want to use the coordinates of each reads in a bam or sam file

You actually want coordinates of the hit that read has in the reference. This would require some manipulation of your alignment file (e.g. extract the reference chromosome name and the start of the hit). Then you could use that information to extract the relevant sequence from the reference.

ADD REPLYlink modified 17 months ago • written 17 months ago by genomax89k

yes, thanks! That's what I want to ask. And is there any python module which can handle this problem quickly?

ADD REPLYlink written 17 months ago by zhangdengwei130
Please log in to add an answer.

Help
Access

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