*** Bismark methylation extractor version v0.19.0 ***

Setting option '--no_overlap' since this is (normally) the right thing to do for paired-end data

Summarising Bismark methylation extractor parameters:
===============================================================
Bismark paired-end SAM format specified (default)
Number of cores to be used: 8
Output will be written to the current directory ('/Users/yaamini/Documents/project-virginica-oa/analyses/2018-10-03-Bismark-Parameter-Testing/2018-10-03-L012-Output')


Summarising bedGraph parameters:
===============================================================
Generating additional output in bedGraph and coverage format
bedGraph format:	<Chromosome> <Start Position> <End Position> <Methylation Percentage>
coverage format:	<Chromosome> <Start Position> <End Position> <Methylation Percentage> <count methylated> <count non-methylated>

Using a cutoff of 1 read(s) to report cytosine positions
Reporting and sorting cytosine methylation information in CpG context only (default)
The bedGraph UNIX sort command will use the following memory setting:	'2G'. Temporary directory used for sorting is the output directory

Checking file >>zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam<< for signs of file truncation...
...passed!
Writing result file containing methylation information for C in CpG context from the original top strand to CpG_OT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original top strand to CpG_CTOT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original bottom strand to CpG_CTOB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the original bottom strand to CpG_OB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHG context from the original top strand to CHG_OT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original top strand to CHG_CTOT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original bottom strand to CHG_CTOB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the original bottom strand to CHG_OB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHH context from the original top strand to CHH_OT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original top strand to CHH_CTOT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original bottom strand to CHH_CTOB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the original bottom strand to CHH_OB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt


Now reading in Bismark result file zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam


Now reading in Bismark result file zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam


Now reading in Bismark result file zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam


Now reading in Bismark result file zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam


Now reading in Bismark result file zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam


Now reading in Bismark result file zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Now reading in Bismark result file zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Now reading in Bismark result file zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

skipping SAM header line:	@HD	VN:1.0	SO:unsorted
skipping SAM header line:	@SQ	SN:NC_035780.1	LN:65668440
skipping SAM header line:	@SQ	SN:NC_035781.1	LN:61752955
skipping SAM header line:	@SQ	SN:NC_035782.1	LN:77061148
skipping SAM header line:	@SQ	SN:NC_035783.1	LN:59691872
skipping SAM header line:	@SQ	SN:NC_035784.1	LN:98698416
skipping SAM header line:	@SQ	SN:NC_035785.1	LN:51258098
skipping SAM header line:	@SQ	SN:NC_035786.1	LN:57830854
skipping SAM header line:	@SQ	SN:NC_035787.1	LN:75944018
skipping SAM header line:	@SQ	SN:NC_035788.1	LN:104168038
skipping SAM header line:	@SQ	SN:NC_035789.1	LN:32650045
skipping SAM header line:	@SQ	SN:NC_007175.2	LN:17244
skipping SAM header line:	@PG	ID:Bismark	VN:v0.19.0	CL:"bismark --non_directional -p 4 -u 10000 -score_min L,0,-1.2 --genome /Users/yaamini/Documents/project-virginica-oa/analyses/2018-04-27-Bismark/2018-04-27-Bismark-Inputs/ -1 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_10_s1_R1_val_1.fq.gz -2 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_10_s1_R2_val_2.fq.gz"
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Now waiting for all child processes to complete

Merging individual splitting reports into overall report: 'zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated_splitting_report.txt'
Merging from these individual files:

Processed 5299 lines in total
Total number of methylation call strings processed: 10598

Final Cytosine Methylation Report
=================================
Total number of C's analysed:	96314

Total methylated C's in CpG context:	9242
Total methylated C's in CHG context:	552
Total methylated C's in CHH context:	2490

Total C to T conversions in CpG context:	3980
Total C to T conversions in CHG context:	20569
Total C to T conversions in CHH context:	59481

C methylated in CpG context:	69.9%
C methylated in CHG context:	2.6%
C methylated in CHH context:	4.0%



Merging individual M-bias reports into overall M-bias statistics from these 8 individual files:
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Deleting unused files ...

CpG_OT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_OB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept


Using these input files: CpG_OT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_OB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Summary of parameters for bismark2bedGraph conversion:
======================================================
bedGraph output:		zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
output directory:		><
remove whitespaces:		no
CX context:			no (CpG context only, default)
No-header selected:		no
Sorting method:			Unix sort-based (smaller memory footprint, but slower)
Sort buffer size:		2G
Coverage threshold:		1
=============================================================================
Methylation information will now be written into a bedGraph and coverage file
=============================================================================

Using the following files as Input:
Writing bedGraph to file: zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
Also writing out a coverage file including counts methylated and unmethylated residues to file: zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz

The genome of interest was specified to contain gazillions of chromosomes or scaffolds. Merging all input files and sorting everything in memory instead of writing out individual chromosome files...
Writing all merged methylation calls to temp file zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged

Finished writing methylation calls from CpG_OT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOT_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_OB_zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Sorting input file zr2096_10_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged by positions (using -S of 2G)
sort: stray character in field spec: invalid field specification `3,3V'
Died at /Users/Shared/Apps/Bismark_v0.19.0/bismark2bedGraph line 486.
Finished BedGraph conversion ...

Checking file >>zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam<< for signs of file truncation...
...passed!
Writing result file containing methylation information for C in CpG context from the original top strand to CpG_OT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original top strand to CpG_CTOT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original bottom strand to CpG_CTOB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the original bottom strand to CpG_OB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHG context from the original top strand to CHG_OT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original top strand to CHG_CTOT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original bottom strand to CHG_CTOB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the original bottom strand to CHG_OB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHH context from the original top strand to CHH_OT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original top strand to CHH_CTOT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original bottom strand to CHH_CTOB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the original bottom strand to CHH_OB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt


Now reading in Bismark result file zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
Now reading in Bismark result file zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Now reading in Bismark result file zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
skipping SAM header line:	@HD	VN:1.0	SO:unsorted
skipping SAM header line:	@SQ	SN:NC_035780.1	LN:65668440
skipping SAM header line:	@SQ	SN:NC_035781.1	LN:61752955
skipping SAM header line:	@SQ	SN:NC_035782.1	LN:77061148
skipping SAM header line:	@SQ	SN:NC_035783.1	LN:59691872
skipping SAM header line:	@SQ	SN:NC_035784.1	LN:98698416
skipping SAM header line:	@SQ	SN:NC_035785.1	LN:51258098
skipping SAM header line:	@SQ	SN:NC_035786.1	LN:57830854
skipping SAM header line:	@SQ	SN:NC_035787.1	LN:75944018
skipping SAM header line:	@SQ	SN:NC_035788.1	LN:104168038
skipping SAM header line:	@SQ	SN:NC_035789.1	LN:32650045
skipping SAM header line:	@SQ	SN:NC_007175.2	LN:17244
skipping SAM header line:	@PG	ID:Bismark	VN:v0.19.0	CL:"bismark --non_directional -p 4 -u 10000 -score_min L,0,-1.2 --genome /Users/yaamini/Documents/project-virginica-oa/analyses/2018-04-27-Bismark/2018-04-27-Bismark-Inputs/ -1 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_1_s1_R1_val_1.fq.gz -2 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_1_s1_R2_val_2.fq.gz"
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Now waiting for all child processes to complete

Merging individual splitting reports into overall report: 'zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated_splitting_report.txt'
Merging from these individual files:

Processed 2872 lines in total
Total number of methylation call strings processed: 5744

Final Cytosine Methylation Report
=================================
Total number of C's analysed:	46083

Total methylated C's in CpG context:	3937
Total methylated C's in CHG context:	313
Total methylated C's in CHH context:	2987

Total C to T conversions in CpG context:	1686
Total C to T conversions in CHG context:	8808
Total C to T conversions in CHH context:	28352

C methylated in CpG context:	70.0%
C methylated in CHG context:	3.4%
C methylated in CHH context:	9.5%



Merging individual M-bias reports into overall M-bias statistics from these 8 individual files:
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Deleting unused files ...

CpG_OT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_OB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept


Using these input files: CpG_OT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_OB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Summary of parameters for bismark2bedGraph conversion:
======================================================
bedGraph output:		zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
output directory:		><
remove whitespaces:		no
CX context:			no (CpG context only, default)
No-header selected:		no
Sorting method:			Unix sort-based (smaller memory footprint, but slower)
Sort buffer size:		2G
Coverage threshold:		1
=============================================================================
Methylation information will now be written into a bedGraph and coverage file
=============================================================================

Using the following files as Input:
Writing bedGraph to file: zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
Also writing out a coverage file including counts methylated and unmethylated residues to file: zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz

The genome of interest was specified to contain gazillions of chromosomes or scaffolds. Merging all input files and sorting everything in memory instead of writing out individual chromosome files...
Writing all merged methylation calls to temp file zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged

Finished writing methylation calls from CpG_OT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOT_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_OB_zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Sorting input file zr2096_1_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged by positions (using -S of 2G)
sort: stray character in field spec: invalid field specification `3,3V'
Died at /Users/Shared/Apps/Bismark_v0.19.0/bismark2bedGraph line 486.
Finished BedGraph conversion ...

Checking file >>zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam<< for signs of file truncation...
...passed!
Writing result file containing methylation information for C in CpG context from the original top strand to CpG_OT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original top strand to CpG_CTOT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original bottom strand to CpG_CTOB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the original bottom strand to CpG_OB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHG context from the original top strand to CHG_OT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original top strand to CHG_CTOT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original bottom strand to CHG_CTOB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the original bottom strand to CHG_OB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHH context from the original top strand to CHH_OT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original top strand to CHH_CTOT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original bottom strand to CHH_CTOB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the original bottom strand to CHH_OB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt


Now reading in Bismark result file zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam


Now reading in Bismark result file zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
Now reading in Bismark result file zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Now reading in Bismark result file zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
skipping SAM header line:	@HD	VN:1.0	SO:unsorted
skipping SAM header line:	@SQ	SN:NC_035780.1	LN:65668440
skipping SAM header line:	@SQ	SN:NC_035781.1	LN:61752955
skipping SAM header line:	@SQ	SN:NC_035782.1	LN:77061148
skipping SAM header line:	@SQ	SN:NC_035783.1	LN:59691872
skipping SAM header line:	@SQ	SN:NC_035784.1	LN:98698416
skipping SAM header line:	@SQ	SN:NC_035785.1	LN:51258098
skipping SAM header line:	@SQ	SN:NC_035786.1	LN:57830854
skipping SAM header line:	@SQ	SN:NC_035787.1	LN:75944018
skipping SAM header line:	@SQ	SN:NC_035788.1	LN:104168038
skipping SAM header line:	@SQ	SN:NC_035789.1	LN:32650045
skipping SAM header line:	@SQ	SN:NC_007175.2	LN:17244
skipping SAM header line:	@PG	ID:Bismark	VN:v0.19.0	CL:"bismark --non_directional -p 4 -u 10000 -score_min L,0,-1.2 --genome /Users/yaamini/Documents/project-virginica-oa/analyses/2018-04-27-Bismark/2018-04-27-Bismark-Inputs/ -1 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_2_s1_R1_val_1.fq.gz -2 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_2_s1_R2_val_2.fq.gz"
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Now waiting for all child processes to complete

Merging individual splitting reports into overall report: 'zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated_splitting_report.txt'
Merging from these individual files:

Processed 4976 lines in total
Total number of methylation call strings processed: 9952

Final Cytosine Methylation Report
=================================
Total number of C's analysed:	75465

Total methylated C's in CpG context:	6975
Total methylated C's in CHG context:	519
Total methylated C's in CHH context:	3320

Total C to T conversions in CpG context:	2934
Total C to T conversions in CHG context:	15916
Total C to T conversions in CHH context:	45801

C methylated in CpG context:	70.4%
C methylated in CHG context:	3.2%
C methylated in CHH context:	6.8%



Merging individual M-bias reports into overall M-bias statistics from these 8 individual files:
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Deleting unused files ...

CpG_OT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_OB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept


Using these input files: CpG_OT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_OB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Summary of parameters for bismark2bedGraph conversion:
======================================================
bedGraph output:		zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
output directory:		><
remove whitespaces:		no
CX context:			no (CpG context only, default)
No-header selected:		no
Sorting method:			Unix sort-based (smaller memory footprint, but slower)
Sort buffer size:		2G
Coverage threshold:		1
=============================================================================
Methylation information will now be written into a bedGraph and coverage file
=============================================================================

Using the following files as Input:
Writing bedGraph to file: zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
Also writing out a coverage file including counts methylated and unmethylated residues to file: zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz

The genome of interest was specified to contain gazillions of chromosomes or scaffolds. Merging all input files and sorting everything in memory instead of writing out individual chromosome files...
Writing all merged methylation calls to temp file zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged

Finished writing methylation calls from CpG_OT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOT_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_OB_zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Sorting input file zr2096_2_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged by positions (using -S of 2G)
sort: stray character in field spec: invalid field specification `3,3V'
Died at /Users/Shared/Apps/Bismark_v0.19.0/bismark2bedGraph line 486.
Finished BedGraph conversion ...

Checking file >>zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam<< for signs of file truncation...
...passed!
Writing result file containing methylation information for C in CpG context from the original top strand to CpG_OT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original top strand to CpG_CTOT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original bottom strand to CpG_CTOB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the original bottom strand to CpG_OB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHG context from the original top strand to CHG_OT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original top strand to CHG_CTOT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original bottom strand to CHG_CTOB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the original bottom strand to CHG_OB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHH context from the original top strand to CHH_OT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original top strand to CHH_CTOT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original bottom strand to CHH_CTOB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the original bottom strand to CHH_OB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt


Now reading in Bismark result file zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
Now reading in Bismark result file zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Now reading in Bismark result file zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
skipping SAM header line:	@HD	VN:1.0	SO:unsorted
skipping SAM header line:	@SQ	SN:NC_035780.1	LN:65668440
skipping SAM header line:	@SQ	SN:NC_035781.1	LN:61752955
skipping SAM header line:	@SQ	SN:NC_035782.1	LN:77061148
skipping SAM header line:	@SQ	SN:NC_035783.1	LN:59691872
skipping SAM header line:	@SQ	SN:NC_035784.1	LN:98698416
skipping SAM header line:	@SQ	SN:NC_035785.1	LN:51258098
skipping SAM header line:	@SQ	SN:NC_035786.1	LN:57830854
skipping SAM header line:	@SQ	SN:NC_035787.1	LN:75944018
skipping SAM header line:	@SQ	SN:NC_035788.1	LN:104168038
skipping SAM header line:	@SQ	SN:NC_035789.1	LN:32650045
skipping SAM header line:	@SQ	SN:NC_007175.2	LN:17244
skipping SAM header line:	@PG	ID:Bismark	VN:v0.19.0	CL:"bismark --non_directional -p 4 -u 10000 -score_min L,0,-1.2 --genome /Users/yaamini/Documents/project-virginica-oa/analyses/2018-04-27-Bismark/2018-04-27-Bismark-Inputs/ -1 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_3_s1_R1_val_1.fq.gz -2 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_3_s1_R2_val_2.fq.gz"
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Now waiting for all child processes to complete
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..

Merging individual splitting reports into overall report: 'zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated_splitting_report.txt'
Merging from these individual files:

Processed 5358 lines in total
Total number of methylation call strings processed: 10716

Final Cytosine Methylation Report
=================================
Total number of C's analysed:	96336

Total methylated C's in CpG context:	10719
Total methylated C's in CHG context:	553
Total methylated C's in CHH context:	2242

Total C to T conversions in CpG context:	3168
Total C to T conversions in CHG context:	20508
Total C to T conversions in CHH context:	59146

C methylated in CpG context:	77.2%
C methylated in CHG context:	2.6%
C methylated in CHH context:	3.7%



Merging individual M-bias reports into overall M-bias statistics from these 8 individual files:
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Deleting unused files ...

CpG_OT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_OB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept


Using these input files: CpG_OT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_OB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Summary of parameters for bismark2bedGraph conversion:
======================================================
bedGraph output:		zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
output directory:		><
remove whitespaces:		no
CX context:			no (CpG context only, default)
No-header selected:		no
Sorting method:			Unix sort-based (smaller memory footprint, but slower)
Sort buffer size:		2G
Coverage threshold:		1
=============================================================================
Methylation information will now be written into a bedGraph and coverage file
=============================================================================

Using the following files as Input:
Writing bedGraph to file: zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
Also writing out a coverage file including counts methylated and unmethylated residues to file: zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz

The genome of interest was specified to contain gazillions of chromosomes or scaffolds. Merging all input files and sorting everything in memory instead of writing out individual chromosome files...
Writing all merged methylation calls to temp file zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged

Finished writing methylation calls from CpG_OT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOT_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_OB_zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Sorting input file zr2096_3_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged by positions (using -S of 2G)
sort: stray character in field spec: invalid field specification `3,3V'
Died at /Users/Shared/Apps/Bismark_v0.19.0/bismark2bedGraph line 486.
Finished BedGraph conversion ...

Checking file >>zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam<< for signs of file truncation...
...passed!
Writing result file containing methylation information for C in CpG context from the original top strand to CpG_OT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original top strand to CpG_CTOT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original bottom strand to CpG_CTOB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the original bottom strand to CpG_OB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHG context from the original top strand to CHG_OT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original top strand to CHG_CTOT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original bottom strand to CHG_CTOB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the original bottom strand to CHG_OB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHH context from the original top strand to CHH_OT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original top strand to CHH_CTOT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original bottom strand to CHH_CTOB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the original bottom strand to CHH_OB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt


Now reading in Bismark result file zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam


Now reading in Bismark result file zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
Now reading in Bismark result file zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Now reading in Bismark result file zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
skipping SAM header line:	@HD	VN:1.0	SO:unsorted
skipping SAM header line:	@SQ	SN:NC_035780.1	LN:65668440
skipping SAM header line:	@SQ	SN:NC_035781.1	LN:61752955
skipping SAM header line:	@SQ	SN:NC_035782.1	LN:77061148
skipping SAM header line:	@SQ	SN:NC_035783.1	LN:59691872
skipping SAM header line:	@SQ	SN:NC_035784.1	LN:98698416
skipping SAM header line:	@SQ	SN:NC_035785.1	LN:51258098
skipping SAM header line:	@SQ	SN:NC_035786.1	LN:57830854
skipping SAM header line:	@SQ	SN:NC_035787.1	LN:75944018
skipping SAM header line:	@SQ	SN:NC_035788.1	LN:104168038
skipping SAM header line:	@SQ	SN:NC_035789.1	LN:32650045
skipping SAM header line:	@SQ	SN:NC_007175.2	LN:17244
skipping SAM header line:	@PG	ID:Bismark	VN:v0.19.0	CL:"bismark --non_directional -p 4 -u 10000 -score_min L,0,-1.2 --genome /Users/yaamini/Documents/project-virginica-oa/analyses/2018-04-27-Bismark/2018-04-27-Bismark-Inputs/ -1 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_4_s1_R1_val_1.fq.gz -2 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_4_s1_R2_val_2.fq.gz"
Finished processing child process. Exiting..
Now waiting for all child processes to complete
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..

Merging individual splitting reports into overall report: 'zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated_splitting_report.txt'
Merging from these individual files:

Processed 5293 lines in total
Total number of methylation call strings processed: 10586

Final Cytosine Methylation Report
=================================
Total number of C's analysed:	85446

Total methylated C's in CpG context:	7365
Total methylated C's in CHG context:	454
Total methylated C's in CHH context:	2475

Total C to T conversions in CpG context:	3497
Total C to T conversions in CHG context:	17953
Total C to T conversions in CHH context:	53702

C methylated in CpG context:	67.8%
C methylated in CHG context:	2.5%
C methylated in CHH context:	4.4%



Merging individual M-bias reports into overall M-bias statistics from these 8 individual files:
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Deleting unused files ...

CpG_OT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_OB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept


Using these input files: CpG_OT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_OB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Summary of parameters for bismark2bedGraph conversion:
======================================================
bedGraph output:		zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
output directory:		><
remove whitespaces:		no
CX context:			no (CpG context only, default)
No-header selected:		no
Sorting method:			Unix sort-based (smaller memory footprint, but slower)
Sort buffer size:		2G
Coverage threshold:		1
=============================================================================
Methylation information will now be written into a bedGraph and coverage file
=============================================================================

Using the following files as Input:
Writing bedGraph to file: zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
Also writing out a coverage file including counts methylated and unmethylated residues to file: zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz

The genome of interest was specified to contain gazillions of chromosomes or scaffolds. Merging all input files and sorting everything in memory instead of writing out individual chromosome files...
Writing all merged methylation calls to temp file zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged

Finished writing methylation calls from CpG_OT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOT_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_OB_zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Sorting input file zr2096_4_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged by positions (using -S of 2G)
sort: stray character in field spec: invalid field specification `3,3V'
Died at /Users/Shared/Apps/Bismark_v0.19.0/bismark2bedGraph line 486.
Finished BedGraph conversion ...

Checking file >>zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam<< for signs of file truncation...
...passed!
Writing result file containing methylation information for C in CpG context from the original top strand to CpG_OT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original top strand to CpG_CTOT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original bottom strand to CpG_CTOB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the original bottom strand to CpG_OB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHG context from the original top strand to CHG_OT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original top strand to CHG_CTOT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original bottom strand to CHG_CTOB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the original bottom strand to CHG_OB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHH context from the original top strand to CHH_OT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original top strand to CHH_CTOT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original bottom strand to CHH_CTOB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the original bottom strand to CHH_OB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt


Now reading in Bismark result file zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam


Now reading in Bismark result file zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
Now reading in Bismark result file zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Now reading in Bismark result file zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
skipping SAM header line:	@HD	VN:1.0	SO:unsorted
skipping SAM header line:	@SQ	SN:NC_035780.1	LN:65668440
skipping SAM header line:	@SQ	SN:NC_035781.1	LN:61752955
skipping SAM header line:	@SQ	SN:NC_035782.1	LN:77061148
skipping SAM header line:	@SQ	SN:NC_035783.1	LN:59691872
skipping SAM header line:	@SQ	SN:NC_035784.1	LN:98698416
skipping SAM header line:	@SQ	SN:NC_035785.1	LN:51258098
skipping SAM header line:	@SQ	SN:NC_035786.1	LN:57830854
skipping SAM header line:	@SQ	SN:NC_035787.1	LN:75944018
skipping SAM header line:	@SQ	SN:NC_035788.1	LN:104168038
skipping SAM header line:	@SQ	SN:NC_035789.1	LN:32650045
skipping SAM header line:	@SQ	SN:NC_007175.2	LN:17244
skipping SAM header line:	@PG	ID:Bismark	VN:v0.19.0	CL:"bismark --non_directional -p 4 -u 10000 -score_min L,0,-1.2 --genome /Users/yaamini/Documents/project-virginica-oa/analyses/2018-04-27-Bismark/2018-04-27-Bismark-Inputs/ -1 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_5_s1_R1_val_1.fq.gz -2 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_5_s1_R2_val_2.fq.gz"
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Now waiting for all child processes to complete
Finished processing child process. Exiting..

Merging individual splitting reports into overall report: 'zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated_splitting_report.txt'
Merging from these individual files:

Processed 5173 lines in total
Total number of methylation call strings processed: 10346

Final Cytosine Methylation Report
=================================
Total number of C's analysed:	88261

Total methylated C's in CpG context:	7407
Total methylated C's in CHG context:	500
Total methylated C's in CHH context:	2394

Total C to T conversions in CpG context:	3462
Total C to T conversions in CHG context:	17842
Total C to T conversions in CHH context:	56656

C methylated in CpG context:	68.1%
C methylated in CHG context:	2.7%
C methylated in CHH context:	4.1%



Merging individual M-bias reports into overall M-bias statistics from these 8 individual files:
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Deleting unused files ...

CpG_OT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_OB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept


Using these input files: CpG_OT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_OB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Summary of parameters for bismark2bedGraph conversion:
======================================================
bedGraph output:		zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
output directory:		><
remove whitespaces:		no
CX context:			no (CpG context only, default)
No-header selected:		no
Sorting method:			Unix sort-based (smaller memory footprint, but slower)
Sort buffer size:		2G
Coverage threshold:		1
=============================================================================
Methylation information will now be written into a bedGraph and coverage file
=============================================================================

Using the following files as Input:
Writing bedGraph to file: zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
Also writing out a coverage file including counts methylated and unmethylated residues to file: zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz

The genome of interest was specified to contain gazillions of chromosomes or scaffolds. Merging all input files and sorting everything in memory instead of writing out individual chromosome files...
Writing all merged methylation calls to temp file zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged

Finished writing methylation calls from CpG_OT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOT_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_OB_zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Sorting input file zr2096_5_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged by positions (using -S of 2G)
sort: stray character in field spec: invalid field specification `3,3V'
Died at /Users/Shared/Apps/Bismark_v0.19.0/bismark2bedGraph line 486.
Finished BedGraph conversion ...

Checking file >>zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam<< for signs of file truncation...
...passed!
Writing result file containing methylation information for C in CpG context from the original top strand to CpG_OT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original top strand to CpG_CTOT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original bottom strand to CpG_CTOB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the original bottom strand to CpG_OB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHG context from the original top strand to CHG_OT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original top strand to CHG_CTOT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original bottom strand to CHG_CTOB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the original bottom strand to CHG_OB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHH context from the original top strand to CHH_OT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original top strand to CHH_CTOT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original bottom strand to CHH_CTOB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the original bottom strand to CHH_OB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt


Now reading in Bismark result file zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
Now reading in Bismark result file zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Now reading in Bismark result file zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
skipping SAM header line:	@HD	VN:1.0	SO:unsorted
skipping SAM header line:	@SQ	SN:NC_035780.1	LN:65668440
skipping SAM header line:	@SQ	SN:NC_035781.1	LN:61752955
skipping SAM header line:	@SQ	SN:NC_035782.1	LN:77061148
skipping SAM header line:	@SQ	SN:NC_035783.1	LN:59691872
skipping SAM header line:	@SQ	SN:NC_035784.1	LN:98698416
skipping SAM header line:	@SQ	SN:NC_035785.1	LN:51258098
skipping SAM header line:	@SQ	SN:NC_035786.1	LN:57830854
skipping SAM header line:	@SQ	SN:NC_035787.1	LN:75944018
skipping SAM header line:	@SQ	SN:NC_035788.1	LN:104168038
skipping SAM header line:	@SQ	SN:NC_035789.1	LN:32650045
skipping SAM header line:	@SQ	SN:NC_007175.2	LN:17244
skipping SAM header line:	@PG	ID:Bismark	VN:v0.19.0	CL:"bismark --non_directional -p 4 -u 10000 -score_min L,0,-1.2 --genome /Users/yaamini/Documents/project-virginica-oa/analyses/2018-04-27-Bismark/2018-04-27-Bismark-Inputs/ -1 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_6_s1_R1_val_1.fq.gz -2 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_6_s1_R2_val_2.fq.gz"
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Now waiting for all child processes to complete

Merging individual splitting reports into overall report: 'zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated_splitting_report.txt'
Merging from these individual files:

Processed 5375 lines in total
Total number of methylation call strings processed: 10750

Final Cytosine Methylation Report
=================================
Total number of C's analysed:	86395

Total methylated C's in CpG context:	8372
Total methylated C's in CHG context:	429
Total methylated C's in CHH context:	2295

Total C to T conversions in CpG context:	3055
Total C to T conversions in CHG context:	18194
Total C to T conversions in CHH context:	54050

C methylated in CpG context:	73.3%
C methylated in CHG context:	2.3%
C methylated in CHH context:	4.1%



Merging individual M-bias reports into overall M-bias statistics from these 8 individual files:
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Deleting unused files ...

CpG_OT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_OB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept


Using these input files: CpG_OT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_OB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Summary of parameters for bismark2bedGraph conversion:
======================================================
bedGraph output:		zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
output directory:		><
remove whitespaces:		no
CX context:			no (CpG context only, default)
No-header selected:		no
Sorting method:			Unix sort-based (smaller memory footprint, but slower)
Sort buffer size:		2G
Coverage threshold:		1
=============================================================================
Methylation information will now be written into a bedGraph and coverage file
=============================================================================

Using the following files as Input:
Writing bedGraph to file: zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
Also writing out a coverage file including counts methylated and unmethylated residues to file: zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz

The genome of interest was specified to contain gazillions of chromosomes or scaffolds. Merging all input files and sorting everything in memory instead of writing out individual chromosome files...
Writing all merged methylation calls to temp file zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged

Finished writing methylation calls from CpG_OT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOT_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_OB_zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Sorting input file zr2096_6_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged by positions (using -S of 2G)
sort: stray character in field spec: invalid field specification `3,3V'
Died at /Users/Shared/Apps/Bismark_v0.19.0/bismark2bedGraph line 486.
Finished BedGraph conversion ...

Checking file >>zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam<< for signs of file truncation...
...passed!
Writing result file containing methylation information for C in CpG context from the original top strand to CpG_OT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original top strand to CpG_CTOT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original bottom strand to CpG_CTOB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the original bottom strand to CpG_OB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHG context from the original top strand to CHG_OT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original top strand to CHG_CTOT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original bottom strand to CHG_CTOB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the original bottom strand to CHG_OB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHH context from the original top strand to CHH_OT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original top strand to CHH_CTOT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original bottom strand to CHH_CTOB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the original bottom strand to CHH_OB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt


Now reading in Bismark result file zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
Now reading in Bismark result file zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Now reading in Bismark result file zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam


Now reading in Bismark result file zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
Warning: unable to close filehandle IN properly.
skipping SAM header line:	@HD	VN:1.0	SO:unsorted
skipping SAM header line:	@SQ	SN:NC_035780.1	LN:65668440
skipping SAM header line:	@SQ	SN:NC_035781.1	LN:61752955
skipping SAM header line:	@SQ	SN:NC_035782.1	LN:77061148
skipping SAM header line:	@SQ	SN:NC_035783.1	LN:59691872
skipping SAM header line:	@SQ	SN:NC_035784.1	LN:98698416
skipping SAM header line:	@SQ	SN:NC_035785.1	LN:51258098
skipping SAM header line:	@SQ	SN:NC_035786.1	LN:57830854
skipping SAM header line:	@SQ	SN:NC_035787.1	LN:75944018
skipping SAM header line:	@SQ	SN:NC_035788.1	LN:104168038
skipping SAM header line:	@SQ	SN:NC_035789.1	LN:32650045
skipping SAM header line:	@SQ	SN:NC_007175.2	LN:17244
skipping SAM header line:	@PG	ID:Bismark	VN:v0.19.0	CL:"bismark --non_directional -p 4 -u 10000 -score_min L,0,-1.2 --genome /Users/yaamini/Documents/project-virginica-oa/analyses/2018-04-27-Bismark/2018-04-27-Bismark-Inputs/ -1 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_7_s1_R1_val_1.fq.gz -2 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_7_s1_R2_val_2.fq.gz"
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Now waiting for all child processes to complete

Merging individual splitting reports into overall report: 'zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated_splitting_report.txt'
Merging from these individual files:

Processed 5144 lines in total
Total number of methylation call strings processed: 10288

Final Cytosine Methylation Report
=================================
Total number of C's analysed:	91819

Total methylated C's in CpG context:	8022
Total methylated C's in CHG context:	528
Total methylated C's in CHH context:	2444

Total C to T conversions in CpG context:	3243
Total C to T conversions in CHG context:	19168
Total C to T conversions in CHH context:	58414

C methylated in CpG context:	71.2%
C methylated in CHG context:	2.7%
C methylated in CHH context:	4.0%



Merging individual M-bias reports into overall M-bias statistics from these 8 individual files:
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Deleting unused files ...

CpG_OT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_OB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept


Using these input files: CpG_OT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_OB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Summary of parameters for bismark2bedGraph conversion:
======================================================
bedGraph output:		zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
output directory:		><
remove whitespaces:		no
CX context:			no (CpG context only, default)
No-header selected:		no
Sorting method:			Unix sort-based (smaller memory footprint, but slower)
Sort buffer size:		2G
Coverage threshold:		1
=============================================================================
Methylation information will now be written into a bedGraph and coverage file
=============================================================================

Using the following files as Input:
Writing bedGraph to file: zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
Also writing out a coverage file including counts methylated and unmethylated residues to file: zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz

The genome of interest was specified to contain gazillions of chromosomes or scaffolds. Merging all input files and sorting everything in memory instead of writing out individual chromosome files...
Writing all merged methylation calls to temp file zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged

Finished writing methylation calls from CpG_OT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOT_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_OB_zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Sorting input file zr2096_7_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged by positions (using -S of 2G)
sort: stray character in field spec: invalid field specification `3,3V'
Died at /Users/Shared/Apps/Bismark_v0.19.0/bismark2bedGraph line 486.
Finished BedGraph conversion ...

Checking file >>zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam<< for signs of file truncation...
...passed!
Writing result file containing methylation information for C in CpG context from the original top strand to CpG_OT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original top strand to CpG_CTOT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original bottom strand to CpG_CTOB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the original bottom strand to CpG_OB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHG context from the original top strand to CHG_OT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original top strand to CHG_CTOT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original bottom strand to CHG_CTOB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the original bottom strand to CHG_OB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHH context from the original top strand to CHH_OT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original top strand to CHH_CTOT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original bottom strand to CHH_CTOB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the original bottom strand to CHH_OB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt


Now reading in Bismark result file zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
Now reading in Bismark result file zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Now reading in Bismark result file zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
skipping SAM header line:	@HD	VN:1.0	SO:unsorted
skipping SAM header line:	@SQ	SN:NC_035780.1	LN:65668440
skipping SAM header line:	@SQ	SN:NC_035781.1	LN:61752955
skipping SAM header line:	@SQ	SN:NC_035782.1	LN:77061148
skipping SAM header line:	@SQ	SN:NC_035783.1	LN:59691872
skipping SAM header line:	@SQ	SN:NC_035784.1	LN:98698416
skipping SAM header line:	@SQ	SN:NC_035785.1	LN:51258098
skipping SAM header line:	@SQ	SN:NC_035786.1	LN:57830854
skipping SAM header line:	@SQ	SN:NC_035787.1	LN:75944018
skipping SAM header line:	@SQ	SN:NC_035788.1	LN:104168038
skipping SAM header line:	@SQ	SN:NC_035789.1	LN:32650045
skipping SAM header line:	@SQ	SN:NC_007175.2	LN:17244
skipping SAM header line:	@PG	ID:Bismark	VN:v0.19.0	CL:"bismark --non_directional -p 4 -u 10000 -score_min L,0,-1.2 --genome /Users/yaamini/Documents/project-virginica-oa/analyses/2018-04-27-Bismark/2018-04-27-Bismark-Inputs/ -1 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_8_s1_R1_val_1.fq.gz -2 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_8_s1_R2_val_2.fq.gz"
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Now waiting for all child processes to complete
Finished processing child process. Exiting..

Merging individual splitting reports into overall report: 'zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated_splitting_report.txt'
Merging from these individual files:

Processed 4760 lines in total
Total number of methylation call strings processed: 9520

Final Cytosine Methylation Report
=================================
Total number of C's analysed:	79951

Total methylated C's in CpG context:	6529
Total methylated C's in CHG context:	440
Total methylated C's in CHH context:	2567

Total C to T conversions in CpG context:	3223
Total C to T conversions in CHG context:	16761
Total C to T conversions in CHH context:	50431

C methylated in CpG context:	67.0%
C methylated in CHG context:	2.6%
C methylated in CHH context:	4.8%



Merging individual M-bias reports into overall M-bias statistics from these 8 individual files:
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Deleting unused files ...

CpG_OT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_OB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept


Using these input files: CpG_OT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_OB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Summary of parameters for bismark2bedGraph conversion:
======================================================
bedGraph output:		zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
output directory:		><
remove whitespaces:		no
CX context:			no (CpG context only, default)
No-header selected:		no
Sorting method:			Unix sort-based (smaller memory footprint, but slower)
Sort buffer size:		2G
Coverage threshold:		1
=============================================================================
Methylation information will now be written into a bedGraph and coverage file
=============================================================================

Using the following files as Input:
Writing bedGraph to file: zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
Also writing out a coverage file including counts methylated and unmethylated residues to file: zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz

The genome of interest was specified to contain gazillions of chromosomes or scaffolds. Merging all input files and sorting everything in memory instead of writing out individual chromosome files...
Writing all merged methylation calls to temp file zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged

Finished writing methylation calls from CpG_OT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOT_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_OB_zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Sorting input file zr2096_8_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged by positions (using -S of 2G)
sort: stray character in field spec: invalid field specification `3,3V'
Died at /Users/Shared/Apps/Bismark_v0.19.0/bismark2bedGraph line 486.
Finished BedGraph conversion ...

Checking file >>zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam<< for signs of file truncation...
...passed!
Writing result file containing methylation information for C in CpG context from the original top strand to CpG_OT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original top strand to CpG_CTOT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the complementary to original bottom strand to CpG_CTOB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CpG context from the original bottom strand to CpG_OB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHG context from the original top strand to CHG_OT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original top strand to CHG_CTOT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the complementary to original bottom strand to CHG_CTOB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHG context from the original bottom strand to CHG_OB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Writing result file containing methylation information for C in CHH context from the original top strand to CHH_OT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original top strand to CHH_CTOT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the complementary to original bottom strand to CHH_CTOB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt
Writing result file containing methylation information for C in CHH context from the original bottom strand to CHH_OB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt


Now reading in Bismark result file zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Now reading in Bismark result file zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam
Warning: unable to close filehandle IN properly.

Now reading in Bismark result file zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.bam

Warning: unable to close filehandle IN properly.
skipping SAM header line:	@HD	VN:1.0	SO:unsorted
skipping SAM header line:	@SQ	SN:NC_035780.1	LN:65668440
skipping SAM header line:	@SQ	SN:NC_035781.1	LN:61752955
skipping SAM header line:	@SQ	SN:NC_035782.1	LN:77061148
skipping SAM header line:	@SQ	SN:NC_035783.1	LN:59691872
skipping SAM header line:	@SQ	SN:NC_035784.1	LN:98698416
skipping SAM header line:	@SQ	SN:NC_035785.1	LN:51258098
skipping SAM header line:	@SQ	SN:NC_035786.1	LN:57830854
skipping SAM header line:	@SQ	SN:NC_035787.1	LN:75944018
skipping SAM header line:	@SQ	SN:NC_035788.1	LN:104168038
skipping SAM header line:	@SQ	SN:NC_035789.1	LN:32650045
skipping SAM header line:	@SQ	SN:NC_007175.2	LN:17244
skipping SAM header line:	@PG	ID:Bismark	VN:v0.19.0	CL:"bismark --non_directional -p 4 -u 10000 -score_min L,0,-1.2 --genome /Users/yaamini/Documents/project-virginica-oa/analyses/2018-04-27-Bismark/2018-04-27-Bismark-Inputs/ -1 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_9_s1_R1_val_1.fq.gz -2 /Volumes/web/Athaliana/20180411_trimgalore_10bp_Cvirginica_MBD/zr2096_9_s1_R2_val_2.fq.gz"
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Now waiting for all child processes to complete
Finished processing child process. Exiting..
Finished processing child process. Exiting..
Finished processing child process. Exiting..

Merging individual splitting reports into overall report: 'zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated_splitting_report.txt'
Merging from these individual files:

Processed 4994 lines in total
Total number of methylation call strings processed: 9988

Final Cytosine Methylation Report
=================================
Total number of C's analysed:	90373

Total methylated C's in CpG context:	9629
Total methylated C's in CHG context:	636
Total methylated C's in CHH context:	2423

Total C to T conversions in CpG context:	3041
Total C to T conversions in CHG context:	19213
Total C to T conversions in CHH context:	55431

C methylated in CpG context:	76.0%
C methylated in CHG context:	3.2%
C methylated in CHH context:	4.2%



Merging individual M-bias reports into overall M-bias statistics from these 8 individual files:
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Determining maximum read lengths for M-Bias plots
Maximum read length of Read 1: 80
Maximum read length of Read 2: 80

Perl module GD::Graph::lines is not installed, skipping drawing M-bias plots (only writing out M-bias plot table)
Deleting unused files ...

CpG_OT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_CTOB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CpG_OB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_CTOB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHG_OB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_CTOB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept
CHH_OB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt contains data ->	kept


Using these input files: CpG_OT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_CTOB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CpG_OB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_CTOB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHG_OB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_CTOB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt CHH_OB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt

Summary of parameters for bismark2bedGraph conversion:
======================================================
bedGraph output:		zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
output directory:		><
remove whitespaces:		no
CX context:			no (CpG context only, default)
No-header selected:		no
Sorting method:			Unix sort-based (smaller memory footprint, but slower)
Sort buffer size:		2G
Coverage threshold:		1
=============================================================================
Methylation information will now be written into a bedGraph and coverage file
=============================================================================

Using the following files as Input:
Writing bedGraph to file: zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz
Also writing out a coverage file including counts methylated and unmethylated residues to file: zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.bismark.cov.gz

The genome of interest was specified to contain gazillions of chromosomes or scaffolds. Merging all input files and sorting everything in memory instead of writing out individual chromosome files...
Writing all merged methylation calls to temp file zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged

Finished writing methylation calls from CpG_OT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOT_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_CTOB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Finished writing methylation calls from CpG_OB_zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.txt to merged temp file
Sorting input file zr2096_9_s1_R1_val_1_bismark_bt2_pe.deduplicated.bedGraph.gz.methylation_calls.merged by positions (using -S of 2G)
sort: stray character in field spec: invalid field specification `3,3V'
Died at /Users/Shared/Apps/Bismark_v0.19.0/bismark2bedGraph line 486.
Finished BedGraph conversion ...