getting the coordinates of mapped regions from a bam file
0
0
Entering edit mode
5 weeks ago

Once I mapped a genome to a reference genome, how can I get the mapped region coordinates and chromosome?

In addition, if possible, how can I get a bam file with the regions of the mapped + a padding (so maybe 2000kb above and below)?

Reads Mapping SAM BAM BWA • 299 views
ADD COMMENT
1
Entering edit mode

What have you tried?

ADD REPLY
0
Entering edit mode

It looks like a homework; but anyways, I did not get the first part of the question, for the second part, did you try to look at bedtools?

ADD REPLY
0
Entering edit mode

I mapped a genome to a reference genome

What does this mean? Please explain in more detail.

how can I get a bam file with the regions of the mapped + a padding

You get a BAM file from an alignment program. The BAM file will not give you any padding. Rather, you would read the BAM file - giving you an alignment position for each read. From this you can collapse (or reduce) overlapping reads to regions, and add padding using things like bedtools, R and the GenomicRanges library, etc. etc.

ADD REPLY

Login before adding your answer.

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