getting the coordinates of mapped regions from a bam file

0

Entering edit mode

2.5 years ago

sciencercur • 0

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 • 997 views

ADD COMMENT • link updated 2.5 years ago by seidel 11k • written 2.5 years ago by sciencercur • 0

1

Entering edit mode

What have you tried?

ADD REPLY • link 2.5 years ago by Ram 44k

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 • link 2.5 years ago by Medhat 9.8k

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 • link 2.5 years ago by seidel 11k

Login before adding your answer.