/oak/stanford/groups/akundaje/marinovg/Symbiodinium/2023-09-18-Fugacium_kawagutii-Hi-C/juicer-Fugacium_kawagutii_V3-Hi-C-3D_DNA-r20/run-asm-pipeline.sh -r 20 --sort-output --build-gapped-map /oak/stanford/groups/akundaje/marinovg/genomes/Dinoflagellata/Fugacium_kawagutii_V3/bwa-indexes/Fugacium_kawagutii_V3_genome_scaffold.fa /oak/stanford/groups/akundaje/marinovg/Symbiodinium/2023-09-18-Fugacium_kawagutii-Hi-C/juicer-Fugacium_kawagutii_V3-Hi-C/aligned/merged_nodups.txt
version: 180922
 -r|--rounds flag was triggered, will run 20 round(s) of misjoin correction.
 --sort-output was triggered, will sort output scaffolds by size.
 --build-gapped-map was triggered, will build an additional hic file corresponding to final assembly with gaps added between draft sequences.
###############
Starting iterating scaffolding with editing:
...starting round 0 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.0.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.0.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
:) DONE!
...visualizing round 0 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 0 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34621
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.1.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 1 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.1.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.1.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
:) DONE!
...visualizing round 1 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 1 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.3488
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.2.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 2 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.2.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.2.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
:) DONE!
...visualizing round 2 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 2 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34706
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.3.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 3 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.3.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.3.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
:) DONE!
...visualizing round 3 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 3 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34719
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.4.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 4 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.4.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.4.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
...Starting iteration # 44
...Starting iteration # 45
...Starting iteration # 46
:) DONE!
...visualizing round 4 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 4 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34659
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.5.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 5 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.5.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.5.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
...Starting iteration # 44
...Starting iteration # 45
...Starting iteration # 46
:) DONE!
...visualizing round 5 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 5 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34669
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.6.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 6 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.6.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.6.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
...Starting iteration # 44
...Starting iteration # 45
...Starting iteration # 46
...Starting iteration # 47
:) DONE!
...visualizing round 6 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 6 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34563
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.7.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 7 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.7.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.7.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
...Starting iteration # 44
...Starting iteration # 45
:) DONE!
...visualizing round 7 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 7 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34066
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.8.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 8 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.8.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.8.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
:) DONE!
...visualizing round 8 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 8 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34473
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.9.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 9 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.9.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.9.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
...Starting iteration # 44
...Starting iteration # 45
:) DONE!
...visualizing round 9 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 9 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34348
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.10.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 10 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.10.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.10.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
:) DONE!
...visualizing round 10 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 10 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.3426
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.11.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 11 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.11.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.11.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
...Starting iteration # 44
...Starting iteration # 45
...Starting iteration # 46
:) DONE!
...visualizing round 11 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 11 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34626
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.12.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 12 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.12.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.12.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
:) DONE!
...visualizing round 12 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 12 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34503
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.13.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 13 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.13.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.13.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
...Starting iteration # 44
:) DONE!
...visualizing round 13 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 13 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34426
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.14.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 14 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.14.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.14.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
...Starting iteration # 44
...Starting iteration # 45
...Starting iteration # 46
...Starting iteration # 47
:) DONE!
...visualizing round 14 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 14 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34435
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.15.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 15 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.15.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.15.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
...Starting iteration # 44
...Starting iteration # 45
...Starting iteration # 46
...Starting iteration # 47
...Starting iteration # 48
:) DONE!
...visualizing round 15 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 15 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34358
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.16.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 16 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.16.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.16.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
...Starting iteration # 44
...Starting iteration # 45
...Starting iteration # 46
...Starting iteration # 47
:) DONE!
...visualizing round 16 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 16 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34149
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.17.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 17 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.17.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.17.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
...Starting iteration # 44
...Starting iteration # 45
...Starting iteration # 46
...Starting iteration # 47
...Starting iteration # 48
...Starting iteration # 49
...Starting iteration # 50
...Starting iteration # 51
...Starting iteration # 52
...Starting iteration # 53
...Starting iteration # 54
...Starting iteration # 55
...Starting iteration # 56
...Starting iteration # 57
:) DONE!
...visualizing round 17 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 17 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34331
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.18.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 18 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.18.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.18.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
...Starting iteration # 44
...Starting iteration # 45
...Starting iteration # 46
...Starting iteration # 47
...Starting iteration # 48
...Starting iteration # 49
...Starting iteration # 50
...Starting iteration # 51
...Starting iteration # 52
...Starting iteration # 53
...Starting iteration # 54
...Starting iteration # 55
...Starting iteration # 56
...Starting iteration # 57
:) DONE!
...visualizing round 18 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 18 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.34289
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.19.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 19 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.19.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.19.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
...Starting iteration # 44
...Starting iteration # 45
...Starting iteration # 46
...Starting iteration # 47
...Starting iteration # 48
...Starting iteration # 49
...Starting iteration # 50
...Starting iteration # 51
...Starting iteration # 52
:) DONE!
...visualizing round 19 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
...detecting misjoins in round 19 assembly:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 25000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 125000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 25000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.3428
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
:) -w flag was triggered, performing cursory search for repeat at 25000 resolution
:) -t flag was triggered, flagging regions with coverage higher than 2
...Dumping 25000 resolution coverage track
HiC file version: 8
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.20.*
...applying edits to cprops file
...applying edits to mnd file
...starting round 20 of scaffolding:
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -s flag was triggered, starting calculations with 15000 threshold starting contig/scaffold size
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.20.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.mnd.20.txt
...Scaffolding all scaffolds and contigs greater or equal to 15000 bp.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
...Starting iteration # 13
...Starting iteration # 14
...Starting iteration # 15
...Starting iteration # 16
...Starting iteration # 17
...Starting iteration # 18
...Starting iteration # 19
...Starting iteration # 20
...Starting iteration # 21
...Starting iteration # 22
...Starting iteration # 23
...Starting iteration # 24
...Starting iteration # 25
...Starting iteration # 26
...Starting iteration # 27
...Starting iteration # 28
...Starting iteration # 29
...Starting iteration # 30
...Starting iteration # 31
...Starting iteration # 32
...Starting iteration # 33
...Starting iteration # 34
...Starting iteration # 35
...Starting iteration # 36
...Starting iteration # 37
...Starting iteration # 38
...Starting iteration # 39
...Starting iteration # 40
...Starting iteration # 41
...Starting iteration # 42
...Starting iteration # 43
...Starting iteration # 44
...Starting iteration # 45
...Starting iteration # 46
...Starting iteration # 47
...Starting iteration # 48
...Starting iteration # 49
...Starting iteration # 50
...Starting iteration # 51
...Starting iteration # 52
...Starting iteration # 53
...Starting iteration # 54
...Starting iteration # 55
...Starting iteration # 56
:) DONE!
...visualizing round 20 results:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
###############
Starting polish:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, performing polishing taking into account signal with minimum 1 mapping quality
...-j flag was triggered, will use Juicebox map Fugacium_kawagutii_V3_genome_scaffold.resolved.hic
...-a flag was triggered, will use scaffold annotation file Fugacium_kawagutii_V3_genome_scaffold.resolved_asm.scaffold_track.txt
...-b flag was triggered, will use superscaffold annotation file Fugacium_kawagutii_V3_genome_scaffold.resolved_asm.superscaf_track.txt
...-s flag was triggered, will ignore all scaffolds shorter than 1000000 for polishing
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 100000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 3000000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 100000 resolution
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -d flag was triggered, depletion score will be averaged across a region bounded by 3000000 superdiagonal
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 100000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.32672
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.polish.*
...applying edits to cprops file
...applying edits to mnd file
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
:) -s flag was triggered, starting calculations with 1000000 threshold starting contig/scaffold size
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.resolved.polish.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.resolved.polish.mnd.txt
...Explicit scaffold set has been listed as input. Using set as a first iteration.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
...Starting iteration # 6
...Starting iteration # 7
...Starting iteration # 8
...Starting iteration # 9
...Starting iteration # 10
...Starting iteration # 11
...Starting iteration # 12
:) DONE!
:) -p flag was triggered. Running LIGer with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations with 1 threshold mapping quality
:) -s flag was triggered, starting calculations with 1000000 threshold starting contig/scaffold size
...Using cprops file: Fugacium_kawagutii_V3_genome_scaffold.resolved.polish.cprops
...Using merged_nodups file: Fugacium_kawagutii_V3_genome_scaffold.resolved.polish.mnd.txt
...Explicit scaffold set has been listed as input. Using set as a first iteration.
...Starting iteration # 1
...Starting iteration # 2
...Starting iteration # 3
...Starting iteration # 4
...Starting iteration # 5
:) DONE!
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -q flag was triggered, starting calculations for 1 threshold mapping quality
:) -i flag was triggered, building mapq without
:) -c flag was triggered, will remove temporary files after completion
...Remapping contact data from the original contig set to assembly
...Building track files
...Building the hic file
Not including fragment map
Start preprocess
Writing header
Writing body
..
Writing footer

Finished preprocess
HiC file version: 8

Calculating norms for zoom BP_2500000
Calculating norms for zoom BP_1000000
Calculating norms for zoom BP_500000
Calculating norms for zoom BP_250000
Calculating norms for zoom BP_100000
Calculating norms for zoom BP_50000
Calculating norms for zoom BP_25000
Calculating norms for zoom BP_10000
Calculating norms for zoom BP_5000
Calculating norms for zoom BP_1000
Writing expected
Writing norms
Finished writing norms
###############
Starting split:
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
...-q flag was triggered, will ignore all reads with mapping quality less then 1 for polishing
...-j flag was triggered, will use Juicebox map Fugacium_kawagutii_V3_genome_scaffold.polished.hic
...-a flag was triggered, will use scaffold annotation file Fugacium_kawagutii_V3_genome_scaffold.polished_asm.scaffold_track.txt
...-b flag was triggered, will use superscaffold annotation file Fugacium_kawagutii_V3_genome_scaffold.polished_asm.superscaf_track.txt
...-s flag was triggered, will ignore all scaffolds shorter than 100000 for polishing
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 100000 resolution
:) -d flag was triggered, depletion score will be averaged across a region bounded by 3000000 superdiagonal
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
:) -p flag was triggered. Running with GNU Parallel support parameter set to true.
:) -c flag was triggered, starting calculations with 5% saturation level
:) -w flag was triggered, performing cursory search for mismatches at 100000 resolution
:) -k flag was triggered, starting calculations with 55% depletion as mismatch threshold
:) -d flag was triggered, depletion score will be averaged across a region bounded by 3000000 superdiagonal
:) -n flag was triggered, performing mismatch region thinning at 1000 resolution
...Dumping 100000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Coarse resolution saturation level = 4.32691
...Analyzing near-diagonal mismatches
...Filtering mismatch calls without better alternative: Not functional yet, skipping
...Thinning mismatch region boundaries.
...Dumping 1000 resolution matrix
HiC file version: 8
...Estimating necessary saturation level for requested misassembly sensitivity
...Fine resolution saturation level = 1
...Analyzing near-diagonal mismatches
...-p flag was triggered. Running with GNU Parallel support parameter set to true.
...-r flag was triggered, output will be labeled as *.split.*
...applying edits to cprops file
...applying edits to mnd file