Entering edit mode
7.9 years ago
Satyajeet Khare
★
1.6k
Biostars,
I am trying cuffmerge using following command
cuffmerge -o cuffmerge_out -g gene.GTF -s genome.fa -p 10 assembly_list.txt
The command ends with following error
[17:56:54] Inspecting reads and determining fragment length distribution.
Processed 40099 loci.
> Map Properties:
> Normalized Map Mass: 1975849.00
> Raw Map Mass: 1975849.00
> Fragment Length Distribution: Truncated Gaussian (default)
> Default Mean: 200
> Default Std Dev: 80
[17:57:04] Assembling transcripts and estimating abundances.
Processed 40099 loci.
[Sat Dec 24 17:59:32 2016] Comparing against reference file gene.GTF
You are using Cufflinks v2.2.1, which is the most recent release.
Error (GFaSeqGet): subsequence cannot be larger than 1346
Error getting subseq for CUFF.1401.1 (1..1350)!
[FAILED]
Error: could not execute cuffcompare
But if run cuffmerge without the fasta file, it works.
cuffmerge -o cuffmerge_out -g gene.GTF -p 10 assembly_list.txt
Any idea whats going on?
P.S. I am merging 48 samples.
Best,