Question: cuffdiff does not accept BAM files produced by HISAT2
0
gravatar for abryants
21 months ago by
abryants0
abryants0 wrote:

I ran HISAT2 on my several  RNAseq samples and then wanted to process the resulting BAMs with cuffdiff for splicing differences. The cuffdiff does not accept these BAM files and states the need for the files to be in the same order. I sorted the files with BAMtools, but it didn't fix the issue. I will appreciate a solution to this problem or any useful comment.

Thanks!

 

AB

rna-seq cufflinks bam • 1.3k views
ADD COMMENTlink modified 21 months ago • written 21 months ago by abryants0

UPDATE:

Well, I guess i don't do it right. TopHat BAMs are giving the same error:

Fatal error: Exit code 1 () Error: sort order of reads in BAMs must be the same

ADD REPLYlink written 21 months ago by abryants0
0
gravatar for Jennifer Hillman Jackson
20 months ago by
United States
Jennifer Hillman Jackson22k wrote:

Hello,

As a quick guess, the sort order of the reference annotation used could potentially be a factor. Sort the GTF dataset by chromosome and start coordinate and re-run to see if that makes a difference (use the tool "Sort").

If that does not work, would you like to submit the problematic run as a bug report? This will allow us to look over the data and provide the best feedback. 

Should you be executing the job on a local/cloud, if some portion of it can be replicated on http://usegalaxy.org and submitted as a bug, that will be enough. It the problem cannot be replicated on Main, then that is also useful information to know about.

This question was a bit difficult to track because of the update submitted as a reply to the original post (indicating a solution) instead of as a comment. I'll adjust that now.

Thanks, Jen, Galaxy team

 

ADD COMMENTlink written 20 months ago by Jennifer Hillman Jackson22k
Please log in to add an answer.

Help
Access

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