Hi,
I am trying to compare two RNAseq datasets obtained using Ilumina in two different dates.
I have the aligned files using Tophat but when I try to run Cufflinks I get the error message "Error: sort order of reads in BAMs must be the same".
I have checked the header of the aligned.bam files (with samtools view) and realized that one of the datasets had sequences aligned to a genome containing sequences which could not be specifically located on a chromosome (these are sequences called chr6_random, chr9_random etc. etc.) and the other did not. This means that one file has some sequences aligned to a location which is not present at all in the other and this is why Cufflinks is not able to work.
Does anyone know what is the problem with my files and have any suggestion on how to fix this without remapping?
I have actually tried to remap some of the files but I am still getting the same type of aligned sequences in each case...
Thanks!
EA
I am trying to compare two RNAseq datasets obtained using Ilumina in two different dates.
I have the aligned files using Tophat but when I try to run Cufflinks I get the error message "Error: sort order of reads in BAMs must be the same".
I have checked the header of the aligned.bam files (with samtools view) and realized that one of the datasets had sequences aligned to a genome containing sequences which could not be specifically located on a chromosome (these are sequences called chr6_random, chr9_random etc. etc.) and the other did not. This means that one file has some sequences aligned to a location which is not present at all in the other and this is why Cufflinks is not able to work.
Does anyone know what is the problem with my files and have any suggestion on how to fix this without remapping?
I have actually tried to remap some of the files but I am still getting the same type of aligned sequences in each case...
Thanks!
EA
Comment