Skip to content

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.

License

Notifications You must be signed in to change notification settings

uhkniazi/dismiss

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

14 Commits
 
 
 
 
 
 
 
 
 
 

Repository files navigation

dismiss

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 xls results file for methylated regions for peak calling performed by MACS2.
  • Alignment file in .bam format (the same file used to perform peak calling)
  • Bam index file with .bai extension (if bam file is called file.bam then index file should be called file.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

Example data of a completed experiment is available to download and view in IGV here;

http://bioinformatics.ibers.aber.ac.uk/IGVdownloads.zip

About

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.

Topics

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages