I have 2 questions:
- Regarding nf4-science (RNA-seq) in this process of HiSat2
modules/hisat2_align.nf
#!/usr/bin/env nextflow
process HISAT2_ALIGN {
container "community.wave.seqera.io/library/hisat2_samtools:5e49f68a37dc010e"
publishDir "results/align", mode: 'copy'
input:
path reads
path index_zip
output:
path "${reads.simpleName}.bam", emit: bam
path "${reads.simpleName}.hisat2.log", emit: log
script:
"""
tar -xzvf $index_zip
hisat2 -x ${index_zip.simpleName} -U $reads \
--new-summary --summary-file ${reads.simpleName}.hisat2.log | \
samtools view -bS -o ${reads.simpleName}.bam
"""
}
and this main pipeline (rnaseq.nf)
/*
* Pipeline parameters
*/
params.hisat2_index_zip = "data/genome_index.tar.gz"
The code import index_zip
as an input, however, I cannot do so. I have chack GitHub repo, i can’t find genome_index.tar.gz
=========================
2. I entry the training because I want to modify nf-core/chipseq for my own analysis (start from .bam instead of .fastq). After finish nf4-science (rna-seq) and open up the nfcore/rnaseq and nfcore/chipseq code, I still feel that it is quite a lot complicate. The question is which training modules or where should I take the skill to fulfill my purpose?
Thank you for your answer,
Namthip