@@ -17,7 +17,7 @@ Nextflow pipeline to perform BAM realignment or fastq alignment and QC, with/wit
17
17
1 . Nextflow : for common installation procedures see the [ IARC-nf] ( https://github.com/IARCbioinfo/IARC-nf ) repository.
18
18
19
19
### Basic fastq alignment
20
- 2 . [ * bwa* ] ( https://github.com/lh3/bwa )
20
+ 2 . [ * bwa2 * ] ( https://github.com/bwa-mem2/bwa-mem2 ) (default) or [ * bwa* ] ( https://github.com/lh3/bwa )
21
21
3 . [ * samblaster* ] ( https://github.com/GregoryFaust/samblaster )
22
22
4 . [ * sambamba* ] ( https://github.com/lomereiter/sambamba )
23
23
@@ -75,6 +75,7 @@ Nextflow pipeline to perform BAM realignment or fastq alignment and QC, with/wit
75
75
| --feature_file | null | Path to feature file for qualimap |
76
76
| --multiqc_config | null | config yaml file for multiqc |
77
77
| --adapterremoval_opt | null | Command line options for AdapterRemoval |
78
+ | --bwa_mem | bwa-mem2 mem | bwa-mem command; use "bwa mem" to switch to regular bwa-mem (both are in the docker and singularity containers) |
78
79
79
80
* #### Flags
80
81
@@ -89,20 +90,26 @@ Flags are special parameters without value.
89
90
| --bwa_option_M | Trigger the -M option in bwa and the corresponding compatibility option in samblaster (marks shorter split hits as secondary) |
90
91
91
92
## Usage
92
-
93
-
94
93
To run the pipeline on a series of fastq or BAM files in folder * input* and a fasta reference file hg19.fasta, one can type:
95
94
96
95
``` bash
97
- nextflow run iarcbioinfo/alignment-nf -r v1.1 -profile singularity --input_folder input/ --ref hg19.fasta --output_folder output
96
+ nextflow run iarcbioinfo/alignment-nf -r v1.2 -profile singularity --input_folder input/ --ref hg19.fasta --output_folder output
98
97
```
99
98
100
99
To run the pipeline without singularity just remove "-profile singularity". Alternatively, one can run the pipeline using a docker container (-profile docker) the conda receipe containing all required dependencies (-profile conda).
101
100
101
+
102
+ ### Use bwa-mem instead of bwa-mem2
103
+ To use bwa-mem, one can type:
104
+
105
+ ``` bash
106
+ nextflow run iarcbioinfo/alignment-nf -r v1.2 -profile singularity --input_folder input/ --ref hg19.fasta --output_folder output --bwa_mem " bwa mem"
107
+ ```
108
+
102
109
### Enable adapter trimming
103
110
To use the adapter trimming step, you must add the *** --trim* option** , as well as satisfy the requirements above mentionned. For example:
104
111
``` bash
105
- nextflow run iarcbioinfo/alignment-nf -r v1.1 -profile singularity --input_folder input/ --ref hg19.fasta --output_folder output --trim
112
+ nextflow run iarcbioinfo/alignment-nf -r v1.2 -profile singularity --input_folder input/ --ref hg19.fasta --output_folder output --trim
106
113
```
107
114
108
115
### Enable ALT mode
0 commit comments