Description

Align reads to a reference genome using hisat2 then sort with samtools

Input

name
description
pattern

reads

List of input FastQ files of size 1 and 2 for single-end and paired-end data,
respectively.

Structure: [ val(meta), [ path(reads) ] ]

index

HISAT2 genome index file.

Structure: [ val(meta), path(index) ]

*.ht2

splicesites

Splice sites in gtf file.

Structure: [ val(meta), path(splicesites) ]

*.txt

ch_fasta_fai

Reference genome fasta file and its index.

Structure: [ val(meta), path(fasta), path(fai) ]

*.{fasta,fa,fai}

save_unaligned

Save unaligned reads to FastQ files

Output

name
description
pattern

orig_bam

Unsorted BAM file containing read alignments straight from HISAT2.

Structure: [ val(meta), path(bam) ]

*.bam

summary

Alignment log.

Structure: [ val(meta), path(log) ]

*.log

fastq

Optional output FASTQ files containing unaligned reads.

Structure: [ val(meta), [ path(fastq) ] ]

*.fastq.gz

bam

Sorted BAM/CRAM/SAM file.

Structure: [ val(meta), path(bam) ]

*.{bam,cram,sam}

index

BAM/CRAM/SAM index file.

Structure: [ val(meta), path(index) ]

*.{bai,crai,csi}

stats

File containing samtools stats output.

Structure: [ val(meta), path(stats) ]

*.stats

flagstat

File containing samtools flagstat output.

Structure: [ val(meta), path(flagstat) ]

*.flagstat

idxstats

File containing samtools idxstats output.

Structure: [ val(meta), path(idxstats) ]

*.idxstats