How to align multiple fastq files with Hisat2 ?
0
0
Entering edit mode
3.2 years ago
sarai • 0

Hi everyone! I'm trying to align nine fastq.gz files (3 samples in triplicate) with a mouse index using Hisat2. I am using the following loop but it seem that nothing happens when I hit run. Can someone check and correct the commands I'm using? Thanks in advance.

for sample in `ls *.fastq.gz | sed 's/.fastq.gz//g'`
do 
    filename=$(basename $sample)
    echo "file path:" $sample
    echo "file name:" $filename
    ~/tools/hisat2-2.2.1/hisat2 \
    -x Hisat2Index/Mus_musculus.GRCm38.dna.primary_assembly.hisat2index \
    -U $sample.fastq.gz \
    --no-unal \
    --summary-file hisat2_output/$filename.log \
    | samtools view -bS | samtools sort -o hisat2_output/$filename.bam 
done
rna-seq • 619 views
ADD COMMENT

Login before adding your answer.

Traffic: 3221 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6