How to extract unique mapped reads from your bam files?
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.
- 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?
- 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.
- 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.
• 5,077 views
•
link
updated
by
Ram
•
written
by
343924049 •
•
written
by
343924049 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
RNA-seq analysis
written by sk_24 •
-
whole exome sequence fragment size analysis
written by adR
-
featureCounts for WGS instead of RNA-seq
written by BioDH •
-
Why there are low assigned reads in featureCounts output and different outputs for quality checks?
written by newbie
-
IGV how to calculate the "expected insert size"?
written by 864606804 •
-
samtools - Finding Average Fragment size for pair?
written by oars
-
mapping rate of small rna-seq with different references
written by woongjaej •
-
Need Clarity/ Suggestions on Medip Seq Data analysis
written by gskbioinfo143 •
-
number of Unique reads in bowtie2 alignment
written by onkar •
-
Gene expression level with RNA-seq
written by bharata1803
I think these [simple] questions were asked before and you can find answers easily. Here is the one about uniquely mapped reads and about the insert size estimation.
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)
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?