No called peaks from macse 2 and the alignment rate 90%
Hi I am facing a problem in analysing my chip seq data. I ran the bowtie 2 to generate aligment read SAM, then I used the samtools view, sort,rm to generate cleam bam file ready for pack calling. The aligment rate was 90% , 80% in some of my data. But the peack called was 100-800/ and in one case0. I do not know if my chip seq failled and I must reapeat the experiment or there is a technical issue.
• 1,308 views
•
link
updated
by
h.mon
•
written
by
sahhomaha •
•
written
by
sahhomaha •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
CHIP-Seq data analysis
written by prs •
-
High Redundancy Rate in MACS2 Peak Calling for TF ChIP-Seq
written by ilovesuperheroes1993 •
-
RepeatMasker: alignment sometimes doesn't fit
written by Ondřej Moravčík •
-
Aligment rate vs called peaks
written by sahhomaha •
-
All my reads fall in intergenic regions ?
written by debitboro
-
empty SAM file
written by indrani.sarkar89 •
-
Chip-Seq analysis in Biostar handbook
written by s_bio •
-
Validating ChIP-seq peak-calling output across replicates
written by bioinfc37 •
-
Does MNase increase Multi mapper fraction?
written by mkrause •
-
Junctions.Bed File Produced From Tophat
written by Varun Gupta
There could be possibly two issue I can see here ..either there is very low sequencing depth, if not then there is ChIP efficiency problem.. also did you use input sample for peak calling?
load the bam file or bigwig (make it from bam) to IGV, and visually check if there any enrichment in the genome or not. if not, the experiment failed most likely.