Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: ../output/02.20-bismark-methylation-extraction/CF03-CM05-Larvae_R1_001.fastp-trim.20220827_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/home/shared/8TB_HDD_03/sr320/github/ceasmallr/code/< Genome directory: >/home/shared/8TB_HDD_03/sr320/github/ceasmallr/data/genome/< Coverage threshold 0 (all positions will be reported, default) CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 0-based (user specified) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /home/shared/8TB_HDD_03/sr320/github/ceasmallr/data/genome/ chr NC_035780.1 (65668440 bp) chr NC_035781.1 (61752955 bp) chr NC_035782.1 (77061148 bp) chr NC_035783.1 (59691872 bp) chr NC_035784.1 (98698416 bp) chr NC_035785.1 (51258098 bp) chr NC_035786.1 (57830854 bp) chr NC_035787.1 (75944018 bp) chr NC_035788.1 (104168038 bp) chr NC_035789.1 (32650045 bp) chr NC_007175.2 (17244 bp) Stored sequence information of 11 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== Adding context-specific methylation summaries >>> Writing genome-wide cytosine report to: ../output/03-methylation/CF03-CM05-Larvae.CpG_report.txt <<< >>> Writing all cytosine context summary file to: ../output/03-methylation/CF03-CM05-Larvae.cytosine_context_summary.txt <<< Storing all covered cytosine positions for chromosome: NC_007175.2 Writing cytosine report for chromosome NC_007175.2 (stored 1126 different covered positions) Writing cytosine report for chromosome NC_035780.1 (stored 1709148 different covered positions) Writing cytosine report for chromosome NC_035781.1 (stored 1543127 different covered positions) Writing cytosine report for chromosome NC_035782.1 (stored 1844277 different covered positions) Writing cytosine report for chromosome NC_035783.1 (stored 1535522 different covered positions) Writing cytosine report for chromosome NC_035784.1 (stored 2591979 different covered positions) Writing cytosine report for chromosome NC_035785.1 (stored 1060633 different covered positions) Writing cytosine report for chromosome NC_035786.1 (stored 1223008 different covered positions) Writing cytosine report for chromosome NC_035787.1 (stored 1637973 different covered positions) Writing cytosine report for chromosome NC_035788.1 (stored 2172332 different covered positions) Writing cytosine report for last chromosome NC_035789.1 (stored 747271 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 11 chromosomes/scaffolds in total) Now processing chromosomes that were not covered by any methylation calls in the coverage file... All chromosomes in the genome were covered by at least some reads. coverage2cytosine processing complete. Now merging top and bottom strand CpGs into a single CG dinucleotide entity (reading from file >>../output/03-methylation/CF03-CM05-Larvae.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to ../output/03-methylation/CF03-CM05-Larvae.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq!