Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-103_S27_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-103_S27_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1084861 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 838867 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 675348 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 771355 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 694256 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 676942 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 495336 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 695303 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 438486 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 575576 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 570579 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 606443 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 489725 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 626011 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 499547 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 395422 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 420528 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 304048 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-103_S27_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-103_S27_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-104_S28_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-104_S28_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1245587 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 972582 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 779148 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 885398 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 806823 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 781504 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 574871 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 805014 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 508231 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 670718 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 654979 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 701313 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 568484 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 707519 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 574717 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 454207 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 481518 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 346412 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-104_S28_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-104_S28_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-111_S29_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-111_S29_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1191107 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 919576 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 734661 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 845147 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 756515 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 738196 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 545035 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 758039 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 480139 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 634255 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 620582 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 664082 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 538729 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 676868 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 545617 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 432366 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 462079 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 328702 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-111_S29_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-111_S29_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-113_S30_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-113_S30_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1144955 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 887443 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 713156 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 820426 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 742515 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 722827 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 533945 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 737135 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 466646 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 616450 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 600686 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 641027 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 521484 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 651840 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 530329 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 417711 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 446720 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 324079 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-113_S30_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-113_S30_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-119_S31_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-119_S31_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1258095 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 981838 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 783744 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 903172 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 814296 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 793418 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 581495 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 816384 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 515848 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 677051 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 672248 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 707857 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 578501 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 726275 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 589938 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 466175 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 495298 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 357030 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-119_S31_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-119_S31_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-120_S32_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-120_S32_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1149296 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 891788 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 709909 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 825490 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 734685 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 713899 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 529135 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 735663 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 466720 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 618795 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 604403 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 644350 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 526388 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 660854 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 532091 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 418026 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 452442 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 318699 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-120_S32_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-120_S32_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-127_S33_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-127_S33_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1080767 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 835656 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 669892 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 769248 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 695838 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 674370 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 497398 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 693054 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 441134 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 576963 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 571834 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 603117 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 492545 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 624704 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 499341 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 393815 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 418916 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 298019 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-127_S33_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-127_S33_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-128_S34_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-128_S34_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1118299 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 868235 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 697636 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 794861 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 718042 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 698463 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 516544 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 718157 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 454963 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 598576 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 587298 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 625777 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 511044 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 646813 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 519831 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 408189 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 436717 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 309215 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-128_S34_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-128_S34_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-135_S35_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-135_S35_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1111554 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 861610 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 696231 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 791967 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 710008 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 693184 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 512744 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 717754 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 453428 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 598852 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 583083 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 623176 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 508375 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 640525 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 514190 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 405937 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 432890 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 307670 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-135_S35_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-135_S35_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-136_S36_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-136_S36_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1138865 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 879482 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 710165 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 808271 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 725471 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 714051 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 521796 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 728520 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 460988 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 609333 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 598730 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 634410 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 519292 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 654163 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 521834 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 415638 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 443714 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 314317 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-136_S36_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-136_S36_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-143_S37_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-143_S37_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1031704 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 802134 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 645547 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 733058 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 658967 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 644780 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 476265 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 665879 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 423164 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 553777 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 544486 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 579307 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 471307 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 598154 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 478612 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 375442 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 405105 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 287266 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-143_S37_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-143_S37_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-145_S38_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-145_S38_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1123355 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 869315 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 694748 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 795166 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 718802 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 703924 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 517083 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 720656 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 456431 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 601856 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 595335 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 628128 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 514060 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 640967 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 519907 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 406366 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 438128 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 313701 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-145_S38_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-145_S38_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-151_S2_L002_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-151_S2_L002_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 917488 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 710836 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 572506 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 651179 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 586577 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 571817 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 421052 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 589920 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 373543 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 493551 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 477687 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 512479 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 417542 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 534622 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 426534 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 331754 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 359835 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 260112 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-151_S2_L002_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-151_S2_L002_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-152_S3_L002_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-152_S3_L002_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 985275 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 761753 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 612427 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 692270 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 628856 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 611275 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 456942 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 635587 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 399396 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 531151 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 510000 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 552500 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 451232 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 569697 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 455064 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 355309 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 384776 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 274441 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-152_S3_L002_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-152_S3_L002_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-153_S4_L002_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-153_S4_L002_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 823252 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 634980 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 510329 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 585349 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 522095 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 513299 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 378483 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 528072 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 332036 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 441562 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 429418 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 459135 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 375203 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 472706 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 379642 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 297137 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 321938 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 232067 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-153_S4_L002_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-153_S4_L002_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-154_S5_L002_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-154_S5_L002_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 888035 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 690016 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 550880 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 625712 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 567151 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 554398 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 410273 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 568429 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 361707 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 480747 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 464292 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 494758 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 404268 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 507274 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 409208 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 322155 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 346467 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 253569 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-154_S5_L002_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-154_S5_L002_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-159_S6_L002_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-159_S6_L002_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 696991 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 538435 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 435628 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 493156 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 444835 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 437174 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 322082 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 447375 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 281339 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 376613 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 362371 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 387040 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 315456 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 407522 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 320771 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 249595 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 271640 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 197802 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-159_S6_L002_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-159_S6_L002_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-160_S7_L002_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-160_S7_L002_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 983464 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 762039 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 612215 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 695096 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 633032 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 619295 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 452038 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 630025 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 399644 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 528350 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 514303 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 549748 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 448936 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 568370 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 453690 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 356275 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 386758 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 276622 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-160_S7_L002_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-160_S7_L002_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-161_S8_L002_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-161_S8_L002_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 839695 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 653856 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 521097 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 600273 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 535499 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 527366 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 386794 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 539081 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 343247 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 448148 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 441994 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 469099 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 381538 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 484678 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 387699 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 303897 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 327884 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 234600 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-161_S8_L002_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-161_S8_L002_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-162_S9_L002_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-162_S9_L002_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 824811 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 636674 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 511422 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 588916 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 523360 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 513804 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 381151 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 528499 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 333536 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 439322 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 433963 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 460897 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 373712 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 475525 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 380650 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 296799 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 322633 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 231606 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-162_S9_L002_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-162_S9_L002_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-167_S10_L002_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-167_S10_L002_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1018116 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 785397 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 626447 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 711392 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 640091 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 628486 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 460858 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 643157 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 405664 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 546922 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 523529 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 568951 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 458151 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 575423 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 458729 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 365934 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 389376 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 281952 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-167_S10_L002_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-167_S10_L002_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-168_S11_L002_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-168_S11_L002_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 851833 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 658637 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 532651 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 607722 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 541913 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 533405 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 395443 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 550033 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 348188 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 457276 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 445413 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 475085 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 388593 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 497461 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 395878 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 306998 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 335908 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 241347 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-168_S11_L002_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-168_S11_L002_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-169_S12_L002_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-169_S12_L002_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 813911 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 625817 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 501452 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 563569 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 514149 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 500743 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 370520 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 515777 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 326970 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 436873 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 418853 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 449813 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 363957 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 466678 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 368345 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 291196 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 313484 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 227693 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-169_S12_L002_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-169_S12_L002_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-170_S13_L002_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-170_S13_L002_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 922019 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 718685 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 579664 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 667234 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 598437 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 575309 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 428688 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 595151 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 374691 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 496883 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 490605 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 517070 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 424042 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 539467 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 429896 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 335247 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 366582 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 267188 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-170_S13_L002_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-170_S13_L002_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-175_S14_L003_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-175_S14_L003_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 988407 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 764177 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 611384 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 690495 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 630755 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 615274 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 451643 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 633856 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 399555 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 533542 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 512287 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 553228 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 449652 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 567070 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 454946 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 358216 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 381508 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 279672 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-175_S14_L003_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-175_S14_L003_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-176_S15_L003_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-176_S15_L003_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1159032 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 900348 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 718766 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 828066 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 741726 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 717732 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 538566 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 744939 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 469032 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 621815 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 611769 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 651006 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 527137 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 661544 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 538992 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 423563 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 451365 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 322183 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-176_S15_L003_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-176_S15_L003_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-181_S16_L003_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-181_S16_L003_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 998161 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 768701 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 619104 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 699129 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 637275 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 615659 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 457724 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 638718 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 402123 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 536706 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 513651 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 555522 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 454680 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 574243 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 456918 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 360109 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 387428 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 279820 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-181_S16_L003_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-181_S16_L003_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-182_S17_L003_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-182_S17_L003_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1151279 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 888521 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 711485 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 820744 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 734933 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 718122 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 529592 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 737582 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 463641 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 616317 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 599860 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 642837 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 520472 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 652934 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 530515 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 414750 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 445036 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 319418 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-182_S17_L003_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-182_S17_L003_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-184_S18_L003_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-184_S18_L003_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 848319 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 655019 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 517902 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 584845 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 533780 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 516401 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 381508 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 535872 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 337659 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 458185 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 432639 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 474384 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 381630 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 492486 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 379007 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 304134 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 322647 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 237367 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-184_S18_L003_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-184_S18_L003_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-185_S19_L003_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-185_S19_L003_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 484179 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 373957 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 300673 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 342166 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 305901 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 302013 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 223663 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 310818 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 196060 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 260328 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 251300 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 268719 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 218854 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 287116 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 223166 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 170762 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 190139 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 137016 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-185_S19_L003_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-185_S19_L003_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-187_S20_L003_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-187_S20_L003_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 881972 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 676395 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 537813 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 606694 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 553812 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 536520 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 398885 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 561750 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 348321 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 480620 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 446585 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 490597 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 397847 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 501399 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 393666 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 310045 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 333037 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 245945 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-187_S20_L003_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-187_S20_L003_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-188_S21_L003_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-188_S21_L003_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 798150 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 615466 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 495784 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 567021 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 511133 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 497142 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 363280 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 511105 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 322414 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 427587 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 412612 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 442809 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 360304 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 459543 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 363610 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 288904 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 311855 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 224087 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-188_S21_L003_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-188_S21_L003_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-193_S22_L003_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-193_S22_L003_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 940810 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 725522 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 581439 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 659857 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 597631 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 583141 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 431107 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 601439 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 379198 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 504234 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 488112 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 524874 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 426238 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 543391 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 430378 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 339210 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 365463 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 262264 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-193_S22_L003_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-193_S22_L003_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-194_S23_L003_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-194_S23_L003_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1062314 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 821406 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 659527 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 749875 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 678861 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 660722 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 485774 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 682495 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 428876 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 576099 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 553837 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 593795 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 485470 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 609850 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 485152 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 384247 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 411437 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 297624 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-194_S23_L003_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-194_S23_L003_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-199_S24_L003_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-199_S24_L003_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 638642 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 489136 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 397405 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 447812 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 402771 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 391983 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 293167 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 410189 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 258109 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 343796 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 328259 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 353616 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 290825 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 371946 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 290202 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 227350 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 246994 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 179607 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-199_S24_L003_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-199_S24_L003_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-200_S25_L003_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-200_S25_L003_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 844120 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 649218 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 524205 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 593654 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 537636 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 521971 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 388534 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 540152 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 337515 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 456011 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 433925 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 470310 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 383574 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 486062 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 384431 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 301035 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 324842 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 234051 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-200_S25_L003_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-200_S25_L003_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-205_S26_L004_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-205_S26_L004_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 633868 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 489365 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 392582 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 447283 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 405765 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 389583 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 288420 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 405302 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 255438 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 343287 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 326515 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 353084 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 287684 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 367743 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 289319 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 227685 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 246330 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 178113 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-205_S26_L004_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-205_S26_L004_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-206_S27_L004_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-206_S27_L004_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 819271 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 633748 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 508285 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 573892 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 523880 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 502290 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 369655 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 524936 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 329216 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 446368 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 418828 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 460621 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 370383 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 477303 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 369033 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 292288 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 315950 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 228844 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-206_S27_L004_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-206_S27_L004_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-208_S28_L004_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-208_S28_L004_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 959798 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 739451 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 593722 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 682325 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 613105 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 592973 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 438640 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 613937 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 386665 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 513649 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 501028 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 535057 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 434593 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 553280 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 438881 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 348769 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 372849 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 275459 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-208_S28_L004_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-208_S28_L004_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-209_S29_L004_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-209_S29_L004_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 991820 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 765645 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 617532 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 698975 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 636316 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 613630 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 453739 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 635411 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 402034 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 535118 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 517565 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 552793 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 454112 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 568648 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 453748 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 356583 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 384548 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 278130 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-209_S29_L004_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-209_S29_L004_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-214_S30_L004_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-214_S30_L004_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 1040752 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 804070 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 644840 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 732741 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 667086 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 647371 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 471837 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 665764 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 417522 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 557644 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 541195 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 582086 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 470372 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 592504 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 475296 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 375326 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 402364 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 290665 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-214_S30_L004_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-214_S30_L004_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-215_S31_L004_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-215_S31_L004_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 787457 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 605843 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 488162 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 548026 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 498521 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 483934 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 360961 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 504654 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 317439 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 425813 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 405567 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 433779 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 356256 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 452628 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 360234 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 278919 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 303620 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 219955 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-215_S31_L004_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-215_S31_L004_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-220_S32_L004_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-220_S32_L004_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 651275 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 502417 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 404651 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 455894 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 410467 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 403915 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 296192 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 413495 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 262974 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 351401 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 335396 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 363128 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 295032 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 377557 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 294735 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 233828 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 252190 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 181839 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-220_S32_L004_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-220_S32_L004_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-221_S33_L004_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-221_S33_L004_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 694072 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 530785 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 427521 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 488137 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 439980 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 428529 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 317419 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 443343 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 276962 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 372833 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 358735 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 386741 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 313721 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 404708 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 315587 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 247440 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 268934 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 193081 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-221_S33_L004_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-221_S33_L004_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-226_S34_L004_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-226_S34_L004_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 557259 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 430644 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 345154 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 389519 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 353048 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 342963 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 254345 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 355825 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 225260 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 299356 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 287375 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 309429 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 253934 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 321961 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 251222 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 197802 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 217183 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 154555 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-226_S34_L004_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-226_S34_L004_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-227_S35_L004_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-227_S35_L004_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 639892 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 494526 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 398276 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 450713 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 404972 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 401346 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 292955 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 411968 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 258569 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 345000 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 334987 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 356896 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 291728 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 376146 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 290620 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 228970 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 250136 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 179827 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-227_S35_L004_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-227_S35_L004_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-229_S36_L004_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-229_S36_L004_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 793018 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 612412 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 491510 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 554188 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 503590 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 489180 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 359190 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 507185 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 316121 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 428924 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 408830 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 441795 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 358005 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 459699 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 360378 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 284460 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 306547 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 223890 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-229_S36_L004_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-229_S36_L004_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-230_S37_L004_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-230_S37_L004_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 884635 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 682411 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 552029 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 629842 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 565750 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 551130 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 407491 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 570448 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 357899 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 475053 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 466137 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 495323 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 402587 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 514149 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 409282 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 321136 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 346471 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 247340 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-230_S37_L004_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-230_S37_L004_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-41_S38_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-41_S38_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 818143 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 632253 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 509202 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 578985 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 522764 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 511518 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 377164 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 528289 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 330752 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 442308 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 427749 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 456164 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 370653 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 476452 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 374465 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 293334 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 320271 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 232542 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-41_S38_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-41_S38_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-42_S39_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-42_S39_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 884520 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 681912 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 549373 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 624791 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 570838 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 550131 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 405945 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 571687 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 359550 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 478814 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 461898 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 493727 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 402130 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 507384 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 405569 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 316393 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 342993 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 248931 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-42_S39_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-42_S39_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-43_S40_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-43_S40_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 814989 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 627668 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 507139 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 575366 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 518695 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 503899 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 372483 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 523135 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 328660 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 435901 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 424382 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 452889 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 370802 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 470545 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 368403 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 292391 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 318579 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 226235 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-43_S40_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-43_S40_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq! Summary of parameters for genome-wide cytosine report: ============================================================================== Coverage infile: /gscratch/srlab/sr320/data/geoduck/cov_files/EPI-44_S41_L005_R1_001_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz Output directory: >< Parent directory: >/gscratch/scrubbed/sr320/1104c/< Genome directory: >/gscratch/srlab/sr320/data/geoduck/v074/< CX context: no (CpG context only, default) Pooling CpG top/bottom strand evidence: yes Genome coordinates used: 1-based (default) GZIP compression: no Split by chromosome: no Now reading in and storing sequence information of the genome specified in: /gscratch/srlab/sr320/data/geoduck/v074/ chr PGA_scaffold1__77_contigs__length_89643857 (89643857 bp) chr PGA_scaffold2__36_contigs__length_69596280 (69596280 bp) chr PGA_scaffold3__111_contigs__length_57743597 (57743597 bp) chr PGA_scaffold4__129_contigs__length_65288255 (65288255 bp) chr PGA_scaffold5__109_contigs__length_67248332 (67248332 bp) chr PGA_scaffold6__104_contigs__length_61759565 (61759565 bp) chr PGA_scaffold7__69_contigs__length_43120122 (43120122 bp) chr PGA_scaffold8__63_contigs__length_61151155 (61151155 bp) chr PGA_scaffold9__45_contigs__length_38581958 (38581958 bp) chr PGA_scaffold10__49_contigs__length_53961475 (53961475 bp) chr PGA_scaffold11__79_contigs__length_51449921 (51449921 bp) chr PGA_scaffold12__71_contigs__length_50438331 (50438331 bp) chr PGA_scaffold13__52_contigs__length_44396874 (44396874 bp) chr PGA_scaffold14__91_contigs__length_45393038 (45393038 bp) chr PGA_scaffold15__101_contigs__length_47938513 (47938513 bp) chr PGA_scaffold16__33_contigs__length_31980953 (31980953 bp) chr PGA_scaffold17__51_contigs__length_34923512 (34923512 bp) chr PGA_scaffold18__69_contigs__length_27737463 (27737463 bp) Stored sequence information of 18 chromosomes/scaffolds in total ============================================================================== Methylation information will now be written into a genome-wide cytosine report ============================================================================== >>> Writing genome-wide cytosine report to: EPI-44_S41_L005_.CpG_report.txt <<< Storing all covered cytosine positions for chromosome: PGA_scaffold1__77_contigs__length_89643857 Writing cytosine report for chromosome PGA_scaffold1__77_contigs__length_89643857 (stored 588975 different covered positions) Writing cytosine report for chromosome PGA_scaffold2__36_contigs__length_69596280 (stored 453991 different covered positions) Writing cytosine report for chromosome PGA_scaffold3__111_contigs__length_57743597 (stored 365994 different covered positions) Writing cytosine report for chromosome PGA_scaffold4__129_contigs__length_65288255 (stored 416373 different covered positions) Writing cytosine report for chromosome PGA_scaffold5__109_contigs__length_67248332 (stored 373321 different covered positions) Writing cytosine report for chromosome PGA_scaffold6__104_contigs__length_61759565 (stored 363497 different covered positions) Writing cytosine report for chromosome PGA_scaffold7__69_contigs__length_43120122 (stored 273122 different covered positions) Writing cytosine report for chromosome PGA_scaffold8__63_contigs__length_61151155 (stored 377813 different covered positions) Writing cytosine report for chromosome PGA_scaffold9__45_contigs__length_38581958 (stored 236744 different covered positions) Writing cytosine report for chromosome PGA_scaffold10__49_contigs__length_53961475 (stored 316503 different covered positions) Writing cytosine report for chromosome PGA_scaffold11__79_contigs__length_51449921 (stored 305652 different covered positions) Writing cytosine report for chromosome PGA_scaffold12__71_contigs__length_50438331 (stored 327312 different covered positions) Writing cytosine report for chromosome PGA_scaffold13__52_contigs__length_44396874 (stored 268339 different covered positions) Writing cytosine report for chromosome PGA_scaffold14__91_contigs__length_45393038 (stored 341289 different covered positions) Writing cytosine report for chromosome PGA_scaffold15__101_contigs__length_47938513 (stored 267283 different covered positions) Writing cytosine report for chromosome PGA_scaffold16__33_contigs__length_31980953 (stored 210410 different covered positions) Writing cytosine report for chromosome PGA_scaffold17__51_contigs__length_34923512 (stored 229577 different covered positions) Writing cytosine report for last chromosome PGA_scaffold18__69_contigs__length_27737463 (stored 164350 different covered positions) Finished writing out cytosine report for covered chromosomes (processed 18 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 >>EPI-44_S41_L005_.CpG_report.txt<<, in output directory '') >>> Writing a new coverage file with top and bottom strand CpG methylation evidence merged to EPI-44_S41_L005_.CpG_report.merged_CpG_evidence.cov <<< CpG merging complete. Good luck finding DMRs with bsseq!