SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
HISAT2 / hisat2_extract_snps_haplotypes_VCF.py: only first variant in output kevinrue RNA Sequencing 1 07-18-2016 04:25 AM
hisat2 output more reads than are in the file frymor Bioinformatics 1 04-12-2016 11:01 AM
pblat empty output Bachbioinfo SOLiD 6 07-07-2015 05:30 AM
Pindel- empty output? Seq_student Bioinformatics 12 08-15-2013 12:30 PM
Bfast output and "Empty Sequence Dictionary" in .sam output aiden Bioinformatics 1 05-28-2010 06:50 PM

Reply
 
Thread Tools
Old 08-07-2021, 04:37 AM   #1
Niloofar
Junior Member
 
Location: Australia

Join Date: May 2012
Posts: 9
Default Hisat2 output is empty

I run this command:

hisat2 -p 40 -q -x Ene_index --sra-acc SRR13696274 --al alignedRNA_1 -S alignment1

and the output says:
4249510 reads; of these:
4249510 (100.00%) were unpaired; of these:
3194713 (75.18%) aligned 0 times
1045690 (24.61%) aligned exactly 1 time
9107 (0.21%) aligned >1 times
24.82% overall alignment rate

However, the produced alignedRNA_1 file is empty

Does anyone know what the problem may be?

Thanks
Niloofar is offline   Reply With Quote
Old 08-14-2021, 07:48 PM   #2
Niloofar
Junior Member
 
Location: Australia

Join Date: May 2012
Posts: 9
Default

Ok - here's the reply I got from the developers, which solved the problem

Hi Niloofar,

Thank you for your reporting.
Could you download the SRA file (use fastq-dump program from SRAToolkit) and run HISAT2 with the file again?
There is a problem when you use --sra-acc and --al (--al*, --al-conc* , --un-conc* options) together, HISAT2 doesn't create read files properly. We'll fix it.

Thanks
Chanhee
Niloofar is offline   Reply With Quote
Reply

Tags
--all, empty output, hisat2

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:09 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO