DISMISS is an R script, which as an additional step in MeDIP-Seq data analysis workflow, enables the allocation of strands to methylated DNA regions. It does this by analyzing the proportions of first mate reads aligning to the methylated locus from the plus and minus strands.
##Requirements
Dismiss is a collection of functions implemented in R and require the following to be present on the system:
- I. R (version 3.1.2)
- II. Bioconductor , (version 3.0) with packages
- A. GenomicAlignments
- B. GenomicFeatures
- C. rtracklayer
##Running dismiss
To run dismiss, the following files should be in the current path or directory:
dismiss_header.R(the script with all the required functions)dismiss_macs2_extractor.R(utility script to process MACS2 results)- Tab separated
xlsresults file for methylated regions for peak calling performed by MACS2. - Alignment file in
.bamformat (the same file used to perform peak calling) - Bam index file with
.baiextension (if bam file is calledfile.bamthen index file should be calledfile.bam.bai)
From the GNU/Linux command line the following command will run dismiss assuming that the files above are in your current working directory;
Rscript dismiss_macs2_extractor.R MACS2_result.xls File.bam [p or s]
where the last argument is essential and should be set to either;
p = for paired end data OR
s = for single end data
##Additional Information
- dismiss on Galaxy http://share-galaxy.ibers.aber.ac.uk/dismiss
- Manual http://uhkniazi.github.io/dismiss
Example data of a completed experiment is available to download and view in IGV here;