Running Hisat2 but it seems not working at all
1
0
Entering edit mode
3.1 years ago
chansik ▴ 10

I'm running Hisat2 with the code below:

hisat2 -p 9 -f ../0.Reference/CHO-PICRH -1 ../2.Trimmomatic/ERR2593198_1_trimmed.fastq.gz -2 ../2.Trimmomatic/ERR2593198_2_trimmed.fastq.gz 2> ../3.HISAT/ERR25593198.log| samtools view -Sbo ../3.HISAT/ERR2593198.bam

It's been running over 24 hours but it doesn't seem work well.

Can you please check the code above and help me solve the problem? (I also attached a screen capture)

Thanks,

enter image description here

HISAT2 Mapping RNA-seq • 1.4k views
ADD COMMENT
0
Entering edit mode

Have you opened another terminal to 1) check that log file you are creating for errors or progress messages and 2) use top or htop to take a look at running processes?

ADD REPLY
0
Entering edit mode

yes, hisat2+ is keep running for several hours on 'top' but as shown in the screen capture, the .log is empty

ADD REPLY
0
Entering edit mode
3.1 years ago

as far as I know hisat2 takes the index via the -x flag

how come your command line does not list that flag

then you are passing the -f flag (FASTA query) but your input is in (FASTQ)

ADD COMMENT
0
Entering edit mode

Thanks for your comment Istvan.

I didn't realize the typo. it should be -x, not -f.

thanks again!

ADD REPLY

Login before adding your answer.

Traffic: 2400 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