cuffmerge badly broken
Cufflinks v1.0.3, I get
__________________________________________________________
You are using Cufflinks v1.0.3, which is the most recent release.
[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 allMerged_refseq/tmp/mergeSam_filenInW01 doesn't appear to be a valid BAM file, trying SAM...
[16:10:58] Loading reference annotation.
[16:11:00] Inspecting reads and determining fragment length distribution.
Error: this SAM file doesn't appear to be correctly sorted!
current hit is at Zv9_NA110:3369, last one was at chrM:11038
Cufflinks requires that if your file has SQ records in
the SAM header that they appear in the same order as the chromosomes names
in the alignments.
If there are no SQ records in the header, or if the header is missing,
the alignments must be sorted lexicographically by chromsome
name and by position.
[FAILED]
_______________________________________________________________
I had to edit the merged_blabla and re-arrange the chromosome order and run
the commands from the log file.
Any fix?
Cufflinks v1.0.3, I get
__________________________________________________________
You are using Cufflinks v1.0.3, which is the most recent release.
[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 allMerged_refseq/tmp/mergeSam_filenInW01 doesn't appear to be a valid BAM file, trying SAM...
[16:10:58] Loading reference annotation.
[16:11:00] Inspecting reads and determining fragment length distribution.
Error: this SAM file doesn't appear to be correctly sorted!
current hit is at Zv9_NA110:3369, last one was at chrM:11038
Cufflinks requires that if your file has SQ records in
the SAM header that they appear in the same order as the chromosomes names
in the alignments.
If there are no SQ records in the header, or if the header is missing,
the alignments must be sorted lexicographically by chromsome
name and by position.
[FAILED]
_______________________________________________________________
I had to edit the merged_blabla and re-arrange the chromosome order and run
the commands from the log file.
Any fix?
Comment