Skip to content
GitLab
Menu
Projects
Groups
Snippets
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
genotoul-bioinfo
metagWGS
Commits
c1241934
Commit
c1241934
authored
Jan 21, 2022
by
VIENNE MAINA
Browse files
Update host_filter.nf
parent
f454f4e4
Pipeline
#48509
skipped with stage
Changes
1
Pipelines
1
Hide whitespace changes
Inline
Side-by-side
modules/host_filter.nf
View file @
c1241934
...
...
@@ -25,10 +25,12 @@ process HOST_FILTER {
mkdir host_filter_flagstat
samtools flagstat ${sampleId}.bam > ${sampleId}.no_filter.flagstat
samtools flagstat ${sampleId}.without_host.bam >> host_filter_flagstat/${sampleId}.host_filter.flagstat
bamToFastq -i ${sampleId}.without_host.bam -fq cleaned_${sampleId}_R1.fastq -fq2 cleaned_${sampleId}_R2.fastq
samtools sort -n -o ${sampleId}.without_host_sort.bam ${sampleId}.without_host.bam
bamToFastq -i ${sampleId}.without_host_sort.bam -fq cleaned_${sampleId}_R1.fastq -fq2 cleaned_${sampleId}_R2.fastq
gzip cleaned_${sampleId}_R1.fastq
gzip cleaned_${sampleId}_R2.fastq
rm ${sampleId}.bam
rm ${sampleId}.without_host.bam
rm ${sampleId}.without_host_sort.bam
"""
}
\ No newline at end of file
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment