SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   TopHat - Cufflinks problems in a multifasta file (http://seqanswers.com/forums/showthread.php?t=17252)

ccstaats 01-31-2012 04:35 AM

TopHat - Cufflinks problems in a transcript multifasta file whitout coordinates
 
Hello all
I'm new to rnaseq analysis, so I wonder if you guys can help me with some issues.
So, I'm trying to evaluate differential expression from the WT and a mutant strain of a yeast, whose genome was sequenced at Broad Institute.
In a first attempt, I'm trying only to evaluate the genes that have been annotated.
The format that comes from Broad is a fasta file with up to 10.000 genes. The genome is not annotated.
In this way, I've downloaded the fasta file with the genes and started to wor on it.
First, I've build a database using bowtie-build.
Then I mapped the reads to the db using tophat
(tophat -r 30 -o /Users/Analises/db/tophat/delta genesmod4.fas 101.fastq) It runs fine
Bou when I tried to run cufflins (cufflinks -o /Users/Analises/db/cufflinks_delta/ /Users/Analises/db/tophat/delta/accepted_hits.bam), I've got this message:

You are using Cufflinks v1.3.0, which is the most recent release.
[10:06:58] Inspecting reads and determining fragment length distribution.
> Processing Locus ABC0001 [ ] 0%
Error: this SAM file doesn't appear to be correctly sorted!
current hit is at ABC_0002:3, last one was at ABC_0001:5679
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.

Any clues?
Thank you in advance

Charley

jbrwn 01-31-2012 10:08 AM

there may be problems in other areas, but it's difficult to troubleshoot with the given information. things that would help are fragment and read length, whether the reads are actually paired-end or not.

for the final problem though, it's likely as simple as it suggests. just sort the output using:
samtools sort aln.bam aln.sorted

see: http://samtools.sourceforge.net/samtools.shtml

ccstaats 01-31-2012 10:51 AM

Hello jbrwn
Thank you for your answer
I think I've found a way to escape away from this problem
I was only considering the transcript fasta file from the Broad Institute.
Know, I got the whole genome fasta file and the gft.
Alll the pathway from tophat to cuffdiff was simple.
Thank you very much
Charley


All times are GMT -8. The time now is 02:21 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.