call peaks of single-end mapping bam files
Hi, for certain reasons, I need to map double-end sequencing result seperately, for example, I mapped A_R1.fastq to A_R1.bam, and A_R2.fastq to A_R2.bam, now I want to call peaks of this file, how can I call peaks of these two files? Do I need to merge them and call peaks?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to merge samples of the same cell type to do differential peak calling?
written by zau saa •
-
How to analyze the generated VCF file, what to do if you have multiple VCF file for the same gene?
written by catalyst •
-
Best way to handle samples with both single end and paired end?
written by ichbinlynn93 •
-
How to treat the pair-end files as single-end when RNA-Seq align?
written by Grace_G •
-
Best practices for variant calling on multiple sequencing runs of the same sample
written by Damian Kao
-
question about calling peak by macs2
written by ch8316f5eyu •
-
BWA: Why paired reads mapped to different chromosome?
written by lghust2011
-
How to merge multiple bam files into single bam file
written by Bioinfonext
-
normalize on chip-seq data
written by catherine
-
Chip-Seq Peak Calling Using Macs In Galaxy- Overlapping Peaks For Replicates
written by anon111 •
What are the reasons to align R1 and R2 seperately??
And if you are aligning them seperately and has some reason to process them as a single end sequencing reads then why would you merge them?
I hope this will help you to think in a right way.
Thanks, Nitin
Now I am analysing hichip data, a software needs to map the reads seperately and then use the two mapped bam files to generate a matrix, this is for the software. But I also need to call peaks from the bam files, so I don't know how to call peak for the seperately mapping files.