SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Tophat: generate junctions.bed file from BAM file Julien Roux Bioinformatics 5 01-14-2016 03:19 PM
How to deal with multiple accepted_hits.bam file as input into cuffdiff? canhu RNA Sequencing 1 06-07-2015 06:17 AM
cuffdiff uses bam file or Sam file wilson90 Bioinformatics 5 08-19-2013 01:38 AM
Cuffdiff - 'the file doesn't appear to be a valid BAM file ekimmike Bioinformatics 2 10-17-2012 01:49 AM
Cuffdiff won't open TopHat accepted_hits.bam Richard Barker RNA Sequencing 12 08-14-2012 10:14 AM

Reply
 
Thread Tools
Old 09-10-2013, 09:41 AM   #1
gwilymh
Member
 
Location: Milwaukee

Join Date: Dec 2011
Posts: 72
Default CuffDiff cannot opn BAM file made in TopHat

I am running a number of Tuxedo pipelines, starting with TopHat and ending in CuffDiff. Most of these pipelines work fine. However, I encountered the following error in CuffDiff for one of my pipelines:

open: No such file or directory
File /gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exonCDS/AR_DM1005_Female/accepted_hits.bam doesn't appear to be a valid BAM file, trying SAM...
Error: cannot open alignment file /gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exonCDS/AR_DM1005_Female/accepted_hits.bam for reading

I have checked over the input file carefully, and cannot find any errors in the file paths.

The accepted_hits.bam file was made in TopHat. I have tried remaking it in TopHat and get the same result. I have looked over it in SamTools and can see no obvious errors.

The only thing that separates this pipelines from my others is that I am using a gff file (made in TopHat from a Flybase GFF3 file) that contains exon and CDS data. Other GFF files with different data combinations seem to work (e.g. exons only, CDS + UTR). Has anyone encountered a similar problem? Is this a bug in CuffDiff?

My full CuffDiff and CuffMerge pipeline is below.



module load tophat/2.0.9
module load bowtie/2.1.0
module load cufflinks/2.1.1
cuffmerge \
-o /gpfs/group1/f/flyinv/Outputs_CuffMerge/exonCDS/Test_ARDM1005 \
-g /gpfs/group1/f/flyinv/working_index/ExonCDS_.gff \
-s /gpfs/group1/f/flyinv/RNASeq/Dpse3_0.fasta \
/gpfs/group1/f/flyinv/Outputs_CuffMerge/exonCDS/CuffLinks_Output_GTF_List_ARDM1005.txt
cuffdiff \
-o "/gpfs/group1/f/flyinv/Outputs_CuffDiff/exonCDS/Test_ARDM1005" \
-L AR_DM1005_Male,AR_DM1005_Female \
--total-hits-norm \
--frag-bias-correct /gpfs/group1/f/flyinv/working_index/Dpse3_0_1.fa \
--multi-read-correct \
--library-norm-method classic-fpkm \
/gpfs/group1/f/flyinv/Outputs_CuffMerge/exonCDS/Test_ARDM1005/merged.gtf \
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_ExonCDS/AR_DM1005_Male/accepted_hits.bam \
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exonCDS/AR_DM1005_Female/accepted_hits.bam
gwilymh is offline   Reply With Quote
Old 09-11-2013, 04:12 PM   #2
tmarria
Junior Member
 
Location: Kansas

Join Date: Jun 2012
Posts: 1
Default Help!! Cuffdiff cannot open alignment file for reading

Hi gwilymh!

I am also getting the exact same error message that you've gotten when trying to run cuffdiff

open: No such file or directory
File /Volumes/Data/2013-08-20_Gp_Cell_Cycle_Transcriptome/Olson_Samples_Run2/SampleXp1/accepted_hits.bam doesn't appear to be a valid BAM file, trying SAM...
Error: cannot open alignment file /Volumes/Data/2013-08-20_Gp_Cell_Cycle_Transcriptome/Olson_Samples_Run2/SampleXp1/accepted_hits.bam for reading

My accepted_hits.bam file was generated in Tophat v2.0.8. I used Cufflinks v2.1.1 but got an error when I tried to use cuffmerge, so I used a slightly older version of Cuffmerge (v2.0.2) to merge my transcript.gtf files generated by Cufflinks (I can provide this error message later if anyone is interested).

I have successfully used previous versions of the Tuxedo software to analyze my RNA-seq data (Tophat v2.0.8, Bowtie 2.1.0.0, and Cufflinks v2.0.2), so maybe it is a bug with the new version of Cufflinks?

I did notice that the version of samtools has changed to 0.1.19.0, where the analyses that I've done with the previous Tuxedo software used samtools 0.1.18.0. Could it also be a samtools issue?


If anyone else is getting this error message, please help us!!!
tmarria is offline   Reply With Quote
Old 09-12-2013, 03:04 AM   #3
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

The first thing to check is if you can "samtools view" the accepted_hits.bam file. If not, it was likely corrupted at some point. gwilymh's problem might be a bug in cuffdiff, it'd be good to make stripped down BAM and annotation files and see if the error persists.

I should note that the initial "open: No such file or directory" error would suggest that the path for the accepted_hits.bam file is simply being misspecified.

Last edited by dpryan; 09-12-2013 at 03:05 AM. Reason: I should really proof read before posting!
dpryan is offline   Reply With Quote
Reply

Tags
bam file, cuffdiff, cufflinks, cuffmerge, tophat

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 10:51 PM.


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