Determining alignment coordinates from BAM file
0
0
Entering edit mode
anikajohn • 0
@e0dc22d0
Last seen 2.3 years ago
Switzerland

Hi,

I apologize in advance as this question might be rather basic. I have an alignment BAM file and I'd like to find out where base x on read y maps to on the reference sequence. In the end I would like to have a list of reference sequence coordinates which correspond to the location where every base of read y maps to. For example:

aligned read: AAGT-> 1,2,6,7 (mapping position on reference genome)

I saw that the GenomicAlignments package has a function called "mapFromAlignments" which sort of does what I want but I am not certain I have the right input to use it. Plus the problem seems easy in principle so I was wondering if there is a more straight forward way to do it.

Thank you already!

GenomicAlignments GenomicA • 1.4k views
ADD COMMENT

Login before adding your answer.

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