Question: How to extract unique mapped reads from your bam files?
gravatar for 343924049
4.5 years ago by
3439240490 wrote:

Hi, guys!

Recently I aligned pair end MeDIP-seq reads using bowtie to hg19, I was confused about three problems although I had searched a lot of related posts.

1.I want to get the fragment of MeDIP-seq from the bam files, for the single end reads, after I aligned to ref genome, I extend my reads to the average fragment size for the following analysis. However, for paired end reads, the fragment could be accurate estimated by insert size. I want to know if existed some kinds of software for me to accomplish the job?

2.For my bam files, I want to get the uniquely mapped reads (fragment) for the following analysis, and this step is before or after the step of "get the fragment"? how do you achieve the goals? Please give me some advice.

3. Would you give me some pipeline or parameters when you process from the fastq files to  unique mapped reads?

This is my first post and looking forward for your help! Many thanks.

ADD COMMENTlink modified 4.5 years ago • written 4.5 years ago by 3439240490

I think these [simple] questions were asked before and you can find answers easily. Here is the one about Bwa-Mem: Discriminate Between Reads Mapping Uniquely And Those Mapping In Multiple Positions and about the Estimate Insert Size In Paired-End/Mate-Pair.

ADD REPLYlink written 4.5 years ago by Pavel Senin1.9k

Thank you for your rely!

I have carefully read the post,  the following fomat is what I want from my bam file(pair end reads)

name start end

chr1 100 300

chr1 120 330


that is each row represent a fragment. I know how to convert bam to bed format, the main confusion is how to get each fragment? Need I to refer the TLEN of sam files?

ADD REPLYlink written 4.5 years ago by 3439240490
Please log in to add an answer.


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