Merge branch 'fq2fa' into 'dev'
Skip fastq>a conversion See merge request epi2melabs/workflow-containers/wf-isoforms!58
This commit is contained in:
commit
63d1d894ee
@ -4,6 +4,10 @@ All notable changes to this project will be documented in this file.
|
||||
The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
|
||||
and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
|
||||
|
||||
## [unreleased]
|
||||
### Changed
|
||||
- Skip unnecessary conversion to fasta from fastq.
|
||||
|
||||
## [v0.1.4]
|
||||
### Changed
|
||||
- Args parser for fastqingress
|
||||
|
||||
@ -15,7 +15,7 @@
|
||||
"out_dir": {
|
||||
"type": "string",
|
||||
"default": "output",
|
||||
"format": "path",
|
||||
"format": "directory-path",
|
||||
"description": "Directory for output of all user-facing files."
|
||||
},
|
||||
"fastq": {
|
||||
|
||||
@ -20,8 +20,7 @@ process map_reads{
|
||||
RightShift: ${params.poly_context}, RegexEnd: "[Aa]{${params.max_poly_run},}",
|
||||
Stranded: True,Invert: True, Tsv: "internal_priming_fail.tsv"} """
|
||||
"""
|
||||
seqkit fq2fa ${fastq_reads} -o "reads.fa";
|
||||
minimap2 -t ${params.threads} -ax splice ${params.minimap2_opts} ${index} "reads.fa"\
|
||||
minimap2 -t ${params.threads} -ax splice ${params.minimap2_opts} ${index} ${fastq_reads}\
|
||||
| samtools view -q ${params.minimum_mapping_quality} -F 2304 -Sb -\
|
||||
| seqkit bam -j ${params.threads} -x -T '${ContextFilter}' -\
|
||||
| samtools sort -@ ${params.threads} -o "${sample_id}_reads_aln_sorted.bam" - ;
|
||||
@ -45,4 +44,4 @@ workflow reference_assembly {
|
||||
emit:
|
||||
bam = map_reads.out.bam
|
||||
stats = map_reads.out.stats
|
||||
}
|
||||
}
|
||||
|
||||
Loading…
Reference in New Issue
Block a user