Question: Cuffdiff - Error in running cuffdiff by using the merged transcriptome assembly along with the BAM files from TopHat
gravatar for evashan96
15 months ago by
evashan960 wrote:

Hi I'm new to RNA-Seq analysis. I ran the following command with cuffdiff: cuffdiff -o /home/silicon/mouse/diff_out -b /home/silicon/mouse/GRCm38/Sequence/Bowtie2Index/genome.fa -p 7 -L g1,g2 -u /home/silicon/mouse/assembly/merged_asm/merged.gtf, \ .//home/silicon/mouse/alignment/Tophat_S2/accepted_hits.bam, /home/silicon/mouse/alignment/samp12_alignment/accepted_hits.bam

and I got the following error: Warning: Could not connect to update server to verify current version. Please check at the Cufflinks website ( Error: cannot open reference GTF file /home/silicon/mouse/assembly/merged_asm/merged.gtf, for reading

Can anyone help me out?

rna-seq • 366 views
ADD COMMENTlink written 15 months ago by evashan960

See below link regarding deprecation of Tophat2-cufflink-cuffdiff

Diverging read alignment and counts in TopHat vs. HISAT2

RNA-Seq Cufflinks normalization methods

RNA Seq Pipeline Question (CuffMerge)

ADD REPLYlink modified 15 months ago • written 15 months ago by Prakash2.0k

Try this if you have one control and one treated sample cufflinks-2.2.1.Linux_x86_64/cuffdiff -u -p 25 --library-type fr-firststrand -o /Cuffdiff/Control_VS_Treated /Cuffmerge/merged.gtf Control/A_Sorted.bam Treated/C_Sorted.bam, -L Control,Treated

In case of replicates try the following command (replicates should be comma separated) cufflinks-2.2.1.Linux_x86_64/cuffdiff -u -p 25 --library-type fr-firststrand -o /Cuffdiff/Control_VS_Treated /Cuffmerge/merged.gtf Control/A_Sorted.bam,Control_2/A1_Sorted.bam Treated/C_Sorted.bam,Treated_2/C1_Sorted.bam -L Control,Treated Also check the location of your merged gtf file (cuffmerge output)

ADD REPLYlink written 15 months ago by bioinfo_ga20
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: 1254 users visited in the last hour