- Any question or suggest or feedback is welcome!
- Feel free to write to us! Your emails encourage us to improve the package.
- We keep on updating selected and well-organized Q&A on this page.
- We welcome your evolvment/submission if you are interested to improve the code throught github.
- CGmap format
- 1.1 Generate CGMapfiles from BAM files
- 1.2 CGmap format by BS-Seeker2 and CGmapTools
- 1.3 merge multiple CGmap files
- 2.1 R program error “optparse”
- Differientially methylated region
- 3.1 DMR analysis when serveral replicates
- 3.1 DMR analysis when lots of replicates
- 3.3 annotate DMRs
Question (by Zhao, Sep. 13th, 2017)
Can I use
cgmaptools convert bam2cgmapto generate CGmap files from the BAM files generated by BSMAP?
BS-Seeker2, Bismark and BSMAP are the mainstream aligners for BS-seq data. We designed CGmapTools to be more general as a downstream analysis pipeline.
From our test, CGmapTools works in right-way for the BAM files generated by both BS-Seeker2 and Bismark.
However, it shows strange result for the BAM file BSMAP when generating ATCGmap files. It is not suggest to directly feed the BAM of BSMAP to CGmapTools currently.
Let us know if you have good solution.
Question (by Fan, Mar. 5th, 2017)
Here I know the BS-Seeker2 and CGmapTools both can get the ATCGmap file and CGmap file. Are they the same file?
The file format is consistent between two softwares.
- BS-Seeker2 is implemented in python (depend on pysam package).
- CGmapTools implemented generating of CGmap from BAM in C language, which would be significantly faster for large genome. Also, CGmapTools supported the BAM files being generated by Bismark.
Question (by Jingyuan, Sep. 18th, 2017)
How could I merge multiple CGmap files using CGmapTools?
You can try two ways using CGmapTools:
The 1st way:
cgmaptools merge2 cgmap -1 A.CGmap.gz -2 B.CGmap.gz | cgmaptools merge2 cgmap -1 C.CGmap.gz | cgmaptools merge2 cgmap -1 D.CGmap.gz -o ABCD_merged.CGmap.gz
The 2nd way:
cgmaptools mergelist tosingle -i A.CGmap.gz,B.CGmap.gz,C.CGmap.gz,D.CGmap.gz -f cgmap -o ABCG_merged.CGmap.gz
Both works well. For two samples,
cgmaptools merge2 cgmaphas better computational efficiency.
Question (by Zhao, Sep. 13th, 2017)
As the package optparse is required when running
cgmaptools lollipop, I use
install.packages("optparse)in R enviroment, but encounter the following error message
package \'optparse\' is not available (for R version x.x.x)
You can try to download the optparse to your local computer, and use
install.packagesto install the library locally.
Question (by Christelle, Aug. 1st, 2017)
How could I used the DMR analysis function in CGmapTools if I have multiple replicates in both case and control? For example, 6 vs. 6.
cgmaptools dmris designed to work with two CGmap files.
If you have multiple replicates, you can first merge the replicates to a single file (to increase coverage and detected region) use the command
cgmaptools mergelist tosingle.
cgmaptools dmrfor the two merged files.
Question (by Yu-Ling, Sep. 5th, 2017)
I worked on RRBS libraries, high-coverage for specific regions, 33 samples vs. 33 samples. How could I use CGmapTools?
When there are lots of replicates, and they are not sufferring from low-coverage, you can
(1) use the command
cgmaptools findCCGGto find regions for RRBS liraries.
(2) use the command
cgmaptools mtrto caculate the average mC levels in region region for each sample.
(3) merge result by (2) into a matrix
(4) write your own program to do t-test.
Question (by Christelle, Sep. 16th, 2017)
How could I annotate the DMR find by CGmapTools.
You may want to annotate DMRs to see the overlapping with genes, or specific functional regions. The software bedtools would be helpful for investigating overlaps.
If you find specific genes are intersting, you can visuazlise the DNA methylations levels, DMR or DMS in local region, using the command