This is a test version of Biostars. For the public version, visit https://www.biostars.org.
bsmap /methratio/ paired end mode and unpaired aligned reads

Hi everyone. Anyone knows if, following paired-end alignement with bsmap, methratio.py will use also the unpaired mapped reads to extract the ratio or just the properly paired?

Thanks

Francesca

paired methylation bisulfite bsmap alignment
• 1,770 views • link
written by francesca3 • francesca3

1 answer

The default is to use all mapped data, if you want methratio.py to use only properly paired reads, use --pair parameter.

 $ ./methratio.py -h
Usage: methratio.py [options] BSMAP_MAPPING_FILES
[..]
  -u, --unique          process only unique mappings/pairs.
  -p, --pair            process only properly paired mappings.
[..]
• 0 views • link
written by sklages sklages

Thanks, one more question. Since I noticed that in my experiment (but I think it happens all the times), there are some "mates" that mapped on different chromosomes, do you know how these reads are treated?

Thanks a lot

Francesca

• 0 views • link
written by francesca3 • francesca3

Log in to answer this question.

More posts like this