SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
Compare "-GTF" and "-GTF-guide" of cufflink options jiyan Bioinformatics 1 01-23-2014 08:42 AM
Some "wrong" XS:A in Tophat output for strand specific pair-end RNA-Seq data ct586 Bioinformatics 4 05-08-2013 05:15 PM
"allele balance ratio" and "quality by depth" in VCF files efoss Bioinformatics 2 10-25-2011 11:13 AM
The position file formats ".clocs" and "_pos.txt"? Ist there any difference? elgor Illumina/Solexa 0 06-27-2011 07:55 AM
"Systems biology and administration" & "Genome generation: no engineering allowed" seb567 Bioinformatics 0 05-25-2010 12:19 PM

Reply
 
Thread Tools
Old 06-20-2013, 07:21 AM   #1
EA01
Junior Member
 
Location: Lisbon

Join Date: Feb 2013
Posts: 5
Default how to compare tophat output files with and without "_random" sequences

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
EA01 is offline   Reply With Quote
Old 06-20-2013, 09:29 AM   #2
bruce01
Senior Member
 
Location: .

Join Date: Mar 2011
Posts: 157
Default

You can just remove the lines with the different chromosome_ID using samtools:

samtools view aligned.bam chr1 chr2 ... chrN > aligned.chr1-N.bam

BAM needs to be coordinate sorted (samtools sort), indexed (samtools index).

Did you use the same reference to align to? You could remove these 'random' chromosomes/contigs from the reference so you just have the 'consensus' chromosomes (this presumes you are working on a well enough annotated species).
bruce01 is offline   Reply With Quote
Old 06-21-2013, 12:05 AM   #3
EA01
Junior Member
 
Location: Lisbon

Join Date: Feb 2013
Posts: 5
Default

Thanks! I'll try this.
Yes, I'm using the mm9, so it should be well annotated.
EA01 is offline   Reply With Quote
Reply

Tags
cufflinks, tophat, _random sequences

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 05:11 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO