This is a test version of Biostars. For the public version, visit https://www.biostars.org.
HISAT2 running pair end

Hi is my first time to use HISAT2 but it does not want to run. I just get:

Warning: Could not open read file "S1_F_paired.fq" for reading; skipping...

I understand that by adding -q it should read the .fq files. Do I need to change my files names to 1.fq and 2.fq to each of my sample reads? I have 48 samples in pairs.

/space/home/aguilar/Programs/hisat2-2.1.0/hisat2 -x /space/home/aguilar/Ofav_temp/HiSat/Ofav -q -1 S1_F_paired.fq,S10_F_paired.fq,S11_F_paired.fq,S12_F_paired.fq -2 S1_R_paired.fq,S10_R_paired.fq,S11_R_paired.fq,S12_R_paired.fq -S /space/home/aguilar/Ofav_temp/HiSat/Ofav.sam

Thanks,

Catalina

hisat2 rna-seq
• 2,953 views • link
written by catagui • catagui

1 answer

I would try changing the file name to 'S1_F_paired_R1.fq' for R1 read and 'S1_F_paired_R2.fq' for R2 read.

• 0 views • link
written by biplab biplab

Thanks, it seems that the program only reads the 1.fq or 2.fq for pair-reads. I wanted to make sure this is the case before editing all my files.

• 0 views • link
written by catagui • catagui

Log in to answer this question.

More posts like this