Skip to content

Commit d4ad75f

Browse files
authored
Merge pull request #27 from B-UMMI/fix_gap_plot
[Patch fix] - Gap and SPN plots
2 parents 7d2aba7 + 2ab52d1 commit d4ad75f

12 files changed

+10
-58
lines changed

README.md

-4
Original file line numberDiff line numberDiff line change
@@ -118,10 +118,6 @@ To use LMAS the following options are available:
118118
--abyssBloomSize Bloom filter size for the ABySS assembler.
119119
It must be a sting with a value and an unit.
120120
(default: 2G)
121-
--bcalm Boolean controling the execution of the BCALM2 assembler.
122-
(default: true)
123-
--bcalmKmerSize K-mer size for the BCALM2 assembler, as an intiger.
124-
(default 31)
125121
--gatb_minia Boolean controling the execution of the GATB Minia Pipeline assembler.
126122
(default: true)
127123
--gatbKmerSize K-mer sizes for the GATB Minia Pipeline assembler.

docs/dev/add_process.rst

+2-5
Original file line numberDiff line numberDiff line change
@@ -12,7 +12,6 @@ The `assembly.nf <https://github.com/B-UMMI/LMAS/tree/main/modules/assembly/asse
1212
The current available assemblers are:
1313

1414
* `ABySS <https://github.com/bcgsc/abyss>`_
15-
* `BCALM 2 <https://github.com/GATB/bcalm>`_
1615
* `GATB-Minia Pipeline <https://github.com/GATB/gatb-minia-pipeline>`_
1716
* `IDBA <https://github.com/loneknightpy/idba>`_
1817
* `MEGAHIT <https://github.com/voutcn/megahit>`_
@@ -152,8 +151,7 @@ It should look like:
152151
153152
.. code-block:: bash
154153
155-
all_assemblies = ABYSS.out.assembly | mix(BCALM2.out.assembly,
156-
GATBMINIAPIPELINE.out.assembly,
154+
all_assemblies = ABYSS.out.assembly | mix(GATBMINIAPIPELINE.out.assembly,
157155
IDBA.out.assembly,
158156
MEGAHIT.out.assembly,
159157
METAHIPMER2.out.assembly,
@@ -175,8 +173,7 @@ It should look like:
175173
176174
.. code-block:: bash
177175
178-
all_versions = ABYSS.out.version | mix(BCALM2.out.version,
179-
GATBMINIAPIPELINE.out.version,
176+
all_versions = ABYSS.out.version | mix(GATBMINIAPIPELINE.out.version,
180177
IDBA.out.version,
181178
MEGAHIT.out.version,
182179
METAHIPMER2.out.version,

docs/getting_started/installation.rst

-4
Original file line numberDiff line numberDiff line change
@@ -146,10 +146,6 @@ To use LMAS the following options are available:
146146
--abyssBloomSize Bloom filter size for the ABySS assembler.
147147
It must be a sting with a value and an unit.
148148
(default: 2G)
149-
--bcalm Boolean controling the execution of the BCALM2 assembler.
150-
(default: true)
151-
--bcalmKmerSize K-mer size for the BCALM2 assembler, as an intiger.
152-
(default 31)
153149
--gatb_minia Boolean controling the execution of the GATB Minia Pipeline assembler.
154150
(default: true)
155151
--gatbKmerSize K-mer sizes for the GATB Minia Pipeline assembler.

docs/user/assemblers.rst

-12
Original file line numberDiff line numberDiff line change
@@ -35,18 +35,6 @@ graph, reducing memory requirements for de novo assembly. **It's a traditional s
3535
* **Date of last release:** 22/04/2021
3636
* **Container:** `cimendes/abyss:2.3.1-1 <https://hub.docker.com/repository/docker/cimendes/abyss>`_
3737

38-
BCALM2
39-
^^^^^^
40-
41-
This assembler, published by `Chikhi et al, 2016 <https://academic.oup.com/bioinformatics/article/32/12/i201/2289008>`_ in
42-
*Bioinformatics*, is a fast and low memory algorithm for graph compaction, consisting of three stages: careful distribution
43-
of input k-mers into buckets, parallel compaction of the buckets, and a parallel reunification step to glue together
44-
the compacted strings into unitigs. **It's a traditional single k-mer value De Bruijn assembler.**
45-
46-
* **Source code:** https://github.com/GATB/bcalm
47-
* **Date of last release:** 22/05/2020
48-
* **Container:** `cimendes/bcalm:2.2.3-1 <https://hub.docker.com/repository/docker/cimendes/bcalm>`_
49-
5038
GATB-Minia Pipeline
5139
^^^^^^^^^^^^^^^^^^^
5240

docs/user/basic_usage.rst

-4
Original file line numberDiff line numberDiff line change
@@ -105,10 +105,6 @@ To use LMAS the following options are available:
105105
--abyssBloomSize Bloom filter size for the ABySS assembler.
106106
It must be a sting with a value and an unit.
107107
(default: 2G)
108-
--bcalm Boolean controling the execution of the BCALM2 assembler.
109-
(default: true)
110-
--bcalmKmerSize K-mer size for the BCALM2 assembler, as an intiger.
111-
(default 31)
112108
--gatb_minia Boolean controling the execution of the GATB Minia Pipeline assembler.
113109
(default: true)
114110
--gatbKmerSize K-mer sizes for the GATB Minia Pipeline assembler.

docs/user/parameters.rst

-18
Original file line numberDiff line numberDiff line change
@@ -120,24 +120,6 @@ ABySS
120120

121121
* **Default:** 2G
122122

123-
BCALM2
124-
^^^^^^
125-
126-
* **Param:** :code:`--bcalm`
127-
128-
* **Definition:** Boolean controling the execution of the BCALM2 assembler.
129-
130-
* **Default:** true
131-
132-
----------------
133-
134-
* **Param:** :code:`--bcalmKmerSize`
135-
136-
* **Definition:** K-mer size for the BCALM2 assembler, as an intiger.
137-
138-
* **Default:** 31
139-
140-
141123
GATB Minia Pipeline
142124
^^^^^^^^^^^^^^^^^^^
143125

modules/assembly/README.md

+1-3
Original file line numberDiff line numberDiff line change
@@ -5,9 +5,7 @@ This module contains the following processes:
55
- REFORMAT
66
- merges the paired-end read files with [BBtools reformat.sh](https://sourceforge.net/projects/bbmap/)
77
- ABYSS
8-
- Assemble the input fastq data with [ABySS](https://github.com/bcgsc/abyss) assembler. Returns it's assembly and the software version.
9-
- BCALM2
10-
- Assemble the input fastq data with [BCALM 2](https://github.com/GATB/bcalm) assembler. Returns it's assembly and the software version.
8+
- Assemble the input fastq data with [ABySS](https://github.com/bcgsc/abyss) assembler. Returns it's assembly and the software version.
119
- GATBMINIAPIPELINE
1210
- Assemble the input fastq data with [GATB-Minia Pipeline](https://github.com/GATB/gatb-minia-pipeline) assembler. Returns it's assembly and the software version.
1311
- IDBA

modules/postprocessing/postprocessing.nf

+3-3
Original file line numberDiff line numberDiff line change
@@ -273,10 +273,10 @@ workflow postprocessing_wf {
273273
PROCESS_SHRIMP_PLOT(phred_csv)
274274
PLOT_CONTIG_DISTRIBUTION(mapping_df_csv)
275275
GAP_ASSESSMENT(paf, triple_reference)
276-
PLOT_GAP_BOXPLOT(GAP_ASSESSMENT.out.json)
277-
PLOT_GAP_REFERENCE(GAP_ASSESSMENT.out.csv)
276+
PLOT_GAP_BOXPLOT(GAP_ASSESSMENT.out.json | collect)
277+
PLOT_GAP_REFERENCE(GAP_ASSESSMENT.out.csv | collect)
278278
SNP_ASSESSMENT(paf, triple_reference)
279-
PLOT_SNP_REFERENCE(SNP_ASSESSMENT.out.csv)
279+
PLOT_SNP_REFERENCE(SNP_ASSESSMENT.out.csv | collect)
280280
MISASSEMBLY(paf)
281281
PROCESS_MISASSEMBLY(MISASSEMBLY.out.trace_pkl | collect, MISASSEMBLY.out.contig_length_pkl | collect, MISASSEMBLY.out.misassembly_json | collect, MISASSEMBLY.out.misassembled_reference_json | collect)
282282
PLOT_MISASSEMBLY(MISASSEMBLY.out.csv | collect)

templates/process_assembly_stats_global.py

+1-1
Original file line numberDiff line numberDiff line change
@@ -36,7 +36,7 @@
3636
ASSEMBLY_STATS_GLOBAL_FILE = '$assembly_stats_global_files'.split()
3737
ASSEMBLY_STATS_GLOBAL_FILE_JSON = "$json_report".split()
3838
N_TARGET = float("$params.n_target")
39-
ASSEMBLER_SKIP = {"ABySS":json.loads("$params.abyss"), "BCALM2": json.loads("$params.bcalm"), "GATBMiniaPipeline": json.loads("$params.gatb_minia"),
39+
ASSEMBLER_SKIP = {"ABySS":json.loads("$params.abyss"), "GATBMiniaPipeline": json.loads("$params.gatb_minia"),
4040
"MetaHipMer2": json.loads("$params.metahipmer2"), "MINIA": json.loads("$params.minia"), "MEGAHIT": json.loads("$params.megahit"),
4141
"metaSPAdes": json.loads("$params.metaspades"), "Unicycler": json.loads("$params.unicycler"), "SPAdes": json.loads("$params.spades"),
4242
"SKESA": json.loads("$params.skesa"), "VelvetOptimiser": json.loads("$params.velvetoptimiser"), "IDBA-UD": json.loads("$params.idba")}

templates/process_assembly_stats_mapping.py

+1-1
Original file line numberDiff line numberDiff line change
@@ -33,7 +33,7 @@
3333

3434
if __file__.endswith(".command.sh"):
3535
ASSEMBLY_STATS_GLOBAL_FILE_JSON = "$json_report".split()
36-
ASSEMBLER_SKIP = {"ABySS":json.loads("$params.abyss"), "BCALM2": json.loads("$params.bcalm"), "GATBMiniaPipeline": json.loads("$params.gatb_minia"),
36+
ASSEMBLER_SKIP = {"ABySS":json.loads("$params.abyss"), "GATBMiniaPipeline": json.loads("$params.gatb_minia"),
3737
"MetaHipMer2": json.loads("$params.metahipmer2"), "MINIA": json.loads("$params.minia"), "MEGAHIT": json.loads("$params.megahit"),
3838
"metaSPAdes": json.loads("$params.metaspades"), "Unicycler": json.loads("$params.unicycler"), "SPAdes": json.loads("$params.spades"),
3939
"SKESA": json.loads("$params.skesa"), "VelvetOptimiser": json.loads("$params.velvetoptimiser"), "IDBA-UD": json.loads("$params.idba")}

templates/process_versions.py

-1
Original file line numberDiff line numberDiff line change
@@ -37,7 +37,6 @@
3737

3838
ASSEMBLER_PROCESS_DICT = {
3939
"ABySS": "ABYSS",
40-
"BCALM2": "BCALM2",
4140
"GATBMiniaPipeline": "GATBMINIAPIPELINE",
4241
"MetaHipMer2": "METAHIPMER2",
4342
"MINIA": "MINIA",

templates/utils.py

+2-2
Original file line numberDiff line numberDiff line change
@@ -14,10 +14,10 @@
1414
COLOURS = ['#a6cee3', '#1f78b4', '#b2df8a', '#33a02c', '#fb9a99', '#e31a1c',
1515
'#fdbf6f', '#ff7f00', '#cab2d6', '#6a3d9a', '#ebdb75', '#b15928']
1616

17-
ASSEMBLER_NAMES = ["ABySS", "BCALM2", "GATBMiniaPipeline", "MetaHipMer2", "MINIA", "MEGAHIT", "metaSPAdes", "Unicycler", "SPAdes",
17+
ASSEMBLER_NAMES = ["ABySS", "GATBMiniaPipeline", "MetaHipMer2", "MINIA", "MEGAHIT", "metaSPAdes", "Unicycler", "SPAdes",
1818
"SKESA", "VelvetOptimiser", "IDBA-UD"]
1919

20-
ASSEMBLER_PROCESS_LIST = ["ABYSS", "BCALM2", "GATBMINIAPIPELINE", "MINIA", "METAHIPMER2", "MEGAHIT", "METASPADES", "UNICYCLER", "SPADES",
20+
ASSEMBLER_PROCESS_LIST = ["ABYSS", "GATBMINIAPIPELINE", "MINIA", "METAHIPMER2", "MEGAHIT", "METASPADES", "UNICYCLER", "SPADES",
2121
"SKESA", "VELVETOPTIMISER", "IDBA"]
2222

2323

0 commit comments

Comments
 (0)