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
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
Comment