This is a test version of Biostars. For the public version, visit https://www.biostars.org.
cuffmerge output question

Hello,

I used cuffmerge to merge all assemblies to create a single merged transcriptome annotation.

cufflinks version 2.2.1

cuffmerge -g Homo_sapiens.GRCh37.82.chr_patch_hapl_scaff.gtf -s human_g1k_v37_decoy_bowtie2 -p 8 assemblies.txt

Output:

[Thu Jun  1 17:56:10 2017] Beginning transcriptome assembly merge
-------------------------------------------

[Thu Jun  1 17:56:10 2017] Preparing output location ./merged_asm/
[Thu Jun  1 17:56:22 2017] Converting GTF files to SAM
[17:56:22] Loading reference annotation.
[17:56:25] Loading reference annotation.
[17:56:27] Loading reference annotation.
[17:56:30] Loading reference annotation.
[17:56:33] Loading reference annotation.
[17:56:36] Loading reference annotation.
[Thu Jun  1 17:56:39 2017] Quantitating transcripts
Warning: Could not connect to update server to verify current version. Please check at the Cufflinks website (http://cufflinks.cbcb.umd.edu).
Command line:
cufflinks -o ./merged_asm/ -F 0.05 -g Homo_sapiens.GRCh37.82.chr_patch_hapl_scaff.gtf -q --overhang-tolerance 200   --library-type=transfrags -A 0.0 --min-frags-per-transfrag 0 --no-5-extend -p 8 ./merged_asm/tmp/mergeSam_fileckgq08 
[bam_header_read] EOF marker is absent. The input is probably truncated.
[bam_header_read] invalid BAM binary header (this is not a BAM file).
File ./merged_asm/tmp/mergeSam_fileckgq08 doesn't appear to be a valid BAM file, trying SAM...
[17:56:39] Loading reference annotation.
[17:57:07] Inspecting reads and determining fragment length distribution.
Processed 57942 loci.                       
> Map Properties:
>       Normalized Map Mass: 405038.00
>       Raw Map Mass: 405038.00
>       Fragment Length Distribution: Truncated Gaussian (default)
>                     Default Mean: 200
>                  Default Std Dev: 80
[17:57:11] Assembling transcripts and estimating abundances.
Processed 57942 loci.                       
[Thu Jun  1 18:02:31 2017] Comparing against reference file Homo_sapiens.GRCh37.82.chr_patch_hapl_scaff.gtf
Warning: Could not connect to update server to verify current version. Please check at the Cufflinks website (http://cufflinks.cbcb.umd.edu).
Error: file/directory human_g1k_v37_decoy_bowtie2 does not exist!
    [FAILED]
Error: could not execute cuffcompare

There is no merged.gtf file in merged_asm. It has genes.fpkm_tracking, isoforms.fpkm_tracking, skipped.gtf, transcripts.gtf and logs, tmp folders. Why there is no merged.gtf file and what about cuff compare?

cuffmerge cufflinks rna-seq

Have you looked into this error: file/directory human_g1k_v37_decoy_bowtie2 does not exist!

Yes, I looked into it. It exists. I used the same for tophat also. Actually In b37 folder human_g1k_v37_decoy_bowtie2.1.bt2, human_g1k_v37_decoy_bowtie2.2.bt2, human_g1k_v37_decoy_bowtie2.3.bt2, human_g1k_v37_decoy_bowtie2.4.bt2, human_g1k_v37_decoy_bowtie2.rev.1.bt2, human_g1k_v37_decoy_bowtie2.rev.2.bt2 files were there. I gave the path b37/human_g1k_v37_decoy_bowtie2

It showed same cuffcompare error and there is no merged.gtf output.

Does cuffmerge comes with cufflinks installation?

Are the genome fasta files in that same directory?

-s/–ref-sequence <seq_dir>/<seq_fasta>

This argument should point to the genomic DNA sequences for the reference. If a directory, it should contain one fasta file per contig. If a multifasta file, all contigs should be present.

No, there are no genome fast files. I have only .bt2 files. How to generate genome fast file?

Where did you get the index files from? The fasta genome files should be available from the same source since they need to match what you used for alignments.

Nope I didn't have any fasta genome files with that. Is it compulsory to use fasta genome files for cuffmerge or index files are fine which are used for tophat?

0 answers

No answers yet.

Log in to answer this question.