How does MEDIPS handle multiple mapping of a read

"Allen [guest]" <[email protected]> Tue, 26 Aug 2014 02:49:38 -0700 (PDT)
Newsgroups gmane.science.biology.informatics.conductor
Message-ID <[email protected]>
Hi,
I was wondering if someone can tell me what happens in MEDIPS a particular read maps to multiple locations on the genome? I have just aligned my MeDIPS-seq data using BWA and only about 40% have a mapping quality of above 30. Does MEDIPS only take into account unique reads and throws out the rest? 

I looked at the MEDIPS manual (http://www.bioconductor.org/packages/release/bioc/vignettes/MEDIPS/inst/doc/MEDIPS.pdf) and can't find any information on this.

I read on various forums and many people suggest just using unique reads. However, for ChIP-seq, it is suggested that one could use CSEM so that information from multiple mapping is not completely lost.

Thanks in advance for your help.

Allen



 -- output of sessionInfo(): 

error 

--
Sent via the guest posting facility at bioconductor.org.

_______________________________________________
Bioconductor mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor