Question: how can I align several fastq files with hisat2 and get htseq-count
gravatar for Learner
14 months ago by
Learner 160
Learner 160 wrote:


I have been strugelling to run Hisat2 code for alignment of several fastq files from human I have read as many post as I could online and I went through the protocol here

At the moment, i have the latest binary version of hisat2 in $HOME/bin When I do hisat2 in the terminal and shows that I don't import any input but it indicates that it is there and can be invoked

No index, query, or output file specified! HISAT2 version 2.1.0 by Daehwan Kim (, Usage:

I have 10 fasta files on a folder named data on my desktop. I have put the three scripts attached to the protocl in that folder and the files are shown below

SRR1_1.fastq.gz SRR9_1.fastq.gz SRR1_2.fastq.gz SRR9_2.fastq.gz SRR2_1.fastq.gz SRR10_1.fastq.gz SRR2_2.fastq.gz SRR10_2.fastq.gz SRR3_1.fastq.gz SRR3_2.fastq.gz SRR4_1.fastq.gz SRR4_2.fastq.gz SRR5_1.fastq.gz SRR5_2.fastq.gz SRR6_1.fastq.gz SRR6_2.fastq.gz SRR7_1.fastq.gz SRR7_2.fastq.gz rnaseq_ballgown.R SRR8_1.fastq.gz SRR8_2.fastq.gz

Can you please let me know how I can run the code to align these fastq files with a Human gene reference ?

I did not change anything on any of the script which I believe that I must change the one in

At this moment, it is

HISAT2=$(which hisat2)
STRINGTIE=$(which stringtie)
SAMTOOLS=$(which samtools)
BASEDIR=$(pwd -P)/chrX_data
TEMPLOC="./tmp" #this will be relative to the output directory

All what I want is to get the htseq-count. txt files for each sample

based on this post Run mutltiple Fastq files on Hisat2 If I do the following

for f in `ls -1 *_1.fastq.gz | sed 's/_1.fastq.gz//' `
echo hisat2 -x mm10idx -1 ${f}_1.fastq.gz -2 ${f}_2.fastq.gz -S ${f}.bam

It will generate all the bam files but all empty. if I echo the do , it will be the output as an example

hisat2 -x mm10idx -1 SRR1_1.fastq.gz -2 SRR1_2.fastq.gz -S SRR1.bam
hisat2 -x mm10idx -1 SRR2_1.fastq.gz -2 SRR42_2.fastq.gz -S SRR2.bam
hisat2 -x mm10idx -1 SRR3_1.fastq.gz -2 SRR3_2.fastq.gz -S SRR3.bam
hisat2 -x mm10idx -1 SRR4_1.fastq.gz -2 SRR_2.fastq.gz -S SRR4.bam
hisat2 -x mm10idx -1 SRR5_1.fastq.gz -2 SRR5_2.fastq.gz -S SRR5.bam
hisat2 -x mm10idx -1 SRR6_1.fastq.gz -2 SRR6_2.fastq.gz -S SRR6.bam
hisat2 -x mm10idx -1 SRR7_1.fastq.gz -2 SRR7_2.fastq.gz -S SRR7.bam
hisat2 -x mm10idx -1 SRR8_1.fastq.gz -2 SRR8_2.fastq.gz -S SRR8.bam
hisat2 -x mm10idx -1 SRR9_1.fastq.gz -2 SRR9_2.fastq.gz -S SRR9.bam
hisat2 -x mm10idx -1 SRR10_1.fastq.gz -2 SRR10_2.fastq.gz -S SRR10.bam

On the other hand, where should I involve the reference genome? for example, I downloaded the hg19 from here

rna-seq • 1.1k views
ADD COMMENTlink modified 14 months ago • written 14 months ago by Learner 160

Maybe DEWE (Differential Workflow Executor for RNA-Seq) can be helpful to you: It provides a user friendly GUI to run HISAT2, convert aligments to bam and run htseq-count. Regards.

ADD REPLYlink written 14 months ago by Hugo150

Better to use Galaxy, it's widely used and provides most of the common tools.

ADD REPLYlink written 14 months ago by Devon Ryan89k

@Devon Ryan can you please tell me how to use galaxy for it? I know how to upload the data and how to run the hisat but some parameters I don't know. For example, single end or paired reads? I guess pair read? or input read as individual or paird? and the rest default values. but the question is how then get the htseq-count? because when I try to use htseq-count function, it needs gff file which I don't have

ADD REPLYlink written 14 months ago by Learner 160

You have paired-end data, so you should select paired end. Run featureCounts on the resulting BAM files. It's in Galaxy too, as is DESeq2. In fact, there's training material available online (it's what we made to use when we train people here in Freiburg).

ADD REPLYlink modified 14 months ago • written 14 months ago by Devon Ryan89k

@Devon Ryan thank you . actually many people think of htseq-count rather than feature counts. I do't know which one is the best but I will try to follow your treating protocol and contact you if needed. Thank you again I believe Galaxy is amazing and easier to work with

ADD REPLYlink written 14 months ago by Learner 160

@Devon Ryan can you please tell me how to get this gff for human ? can I use and then send it to galaxy or do you know a better way?

ADD REPLYlink written 14 months ago by Learner 160

Never use GFF files or annotations from UCSC. Download a GTF file from Gencode or Ensembl.

ADD REPLYlink written 14 months ago by Devon Ryan89k

@Devon Ryan can you please give me a link to the right gft? why would UCSC be a problem?

ADD REPLYlink written 14 months ago by Learner 160

Here is the link to the Gencode, release 27 primary assembly GTF. UCSC annotations have historically both lagged behind those of Ensembl/Gencode and been of lower quality. They have also historically had a variety of odd issues, such as having multiple instances of the same gene on multiple strands/chromosomes all sharing the same ID...which makes no biological sense (they can have the same name, but they need different IDs since they're different loci). Ensembl/Gencode GTFs won't break tools, UCSC annotations have been known to do exactly that.

ADD REPLYlink written 14 months ago by Devon Ryan89k
gravatar for Devon Ryan
14 months ago by
Devon Ryan89k
Freiburg, Germany
Devon Ryan89k wrote:

You're trying to do too many things at once, which is just going to cause confusion for someone new to the field. Start out with a single sample that manually type the command for. Once that works properly you can move on to multiple samples in a loop.

Since you apparently downloaded the hg19 prebuilt reference, you need to specify where it is with the -x option. In your examples, that was -x mm10idx, but presumably something like -x /some/path/to/the/downloaded/files/hg19 would make more sense.

ADD COMMENTlink written 14 months ago by Devon Ryan89k

@Devon Ryan I am trying to run this, I will accept your answer once I get it to run

ADD REPLYlink written 14 months ago by Learner 160

@Devon Ryan that helped me thanks although I had to do a lot of work to get it to run :-D

ADD REPLYlink modified 14 months ago • written 14 months ago by Learner 160
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1133 users visited in the last hour