1
0
Fork 0
mirror of https://github.com/MillironX/taxprofiler.git synced 2024-11-26 03:59:54 +00:00

Include the mapped reads in samtools stats

This commit is contained in:
sofstam 2022-11-22 16:50:50 +01:00
parent 10c72f1265
commit 30af26a8b3
2 changed files with 5 additions and 5 deletions

View file

@ -44,8 +44,8 @@ workflow LONGREAD_HOSTREMOVAL {
SAMTOOLS_INDEX ( SAMTOOLS_VIEW.out.bam ) SAMTOOLS_INDEX ( SAMTOOLS_VIEW.out.bam )
ch_versions = ch_versions.mix( SAMTOOLS_INDEX.out.versions.first() ) ch_versions = ch_versions.mix( SAMTOOLS_INDEX.out.versions.first() )
bam_bai = SAMTOOLS_VIEW.out.bam bam_bai = MINIMAP2_ALIGN.out.bam
.join(SAMTOOLS_INDEX.out.bai, remainder: true) .join(SAMTOOLS_INDEX.out.bai,by:[0], remainder: true)
SAMTOOLS_STATS ( bam_bai, reference ) SAMTOOLS_STATS ( bam_bai, reference )
ch_versions = ch_versions.mix(SAMTOOLS_STATS.out.versions.first()) ch_versions = ch_versions.mix(SAMTOOLS_STATS.out.versions.first())

View file

@ -25,7 +25,7 @@ workflow SHORTREAD_HOSTREMOVAL {
ch_bowtie2_index = index.first() ch_bowtie2_index = index.first()
} }
BOWTIE2_ALIGN ( reads, ch_bowtie2_index, true, false ) BOWTIE2_ALIGN ( reads, ch_bowtie2_index, true, true)
ch_versions = ch_versions.mix( BOWTIE2_ALIGN.out.versions.first() ) ch_versions = ch_versions.mix( BOWTIE2_ALIGN.out.versions.first() )
ch_multiqc_files = ch_multiqc_files.mix( BOWTIE2_ALIGN.out.log ) ch_multiqc_files = ch_multiqc_files.mix( BOWTIE2_ALIGN.out.log )
@ -41,8 +41,8 @@ workflow SHORTREAD_HOSTREMOVAL {
SAMTOOLS_INDEX ( SAMTOOLS_VIEW.out.bam ) SAMTOOLS_INDEX ( SAMTOOLS_VIEW.out.bam )
ch_versions = ch_versions.mix( SAMTOOLS_INDEX.out.versions.first() ) ch_versions = ch_versions.mix( SAMTOOLS_INDEX.out.versions.first() )
bam_bai = SAMTOOLS_VIEW.out.bam bam_bai = BOWTIE2_ALIGN.out.bam
.join(SAMTOOLS_INDEX.out.bai, remainder: true) .join(SAMTOOLS_INDEX.out.bai, by:[0], remainder: true)
SAMTOOLS_STATS ( bam_bai, reference ) SAMTOOLS_STATS ( bam_bai, reference )
ch_versions = ch_versions.mix(SAMTOOLS_STATS.out.versions.first()) ch_versions = ch_versions.mix(SAMTOOLS_STATS.out.versions.first())