Analyze single-end reads in MMDiff2
0
0
Entering edit mode
cattapre • 0
@cattapre-11119
Last seen 4.8 years ago

I am trying to test the MMDiff2 package on some of my single-end read ChIP-seq data, but getPeakReadsalways returns me error messages (warnings). I cannot get the full package to work. The messages are:

In getReads(P, chroms[j], bam.file, pairedEnd = pairedEnd) :

  reads have different read lengths:

The command I used was

MMD <- getPeakReads(MMD, PeakBoundary = 200, pairedEnd = FALSE, run.parallel = FALSE)

I looked at the getPeakReads code ( https://github.com/Bioconductor-mirror/MMDiff2/blob/master/R/getPeakReads.R ) to try to figure it out why I could not get through. If I understood it correctly, I think that pairedEnds is forced to be on all the time.

Could you please help?

Thanks!

Rinaldo

 

software error mmdiff2 • 462 views
ADD COMMENT

Login before adding your answer.

Traffic: 250 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6