Resolving Salmon Errors?
I'm running Salmon with this line:
*./salmon-1.8.0_linux_x86_64/bin/salmon quant -p 12 -t microbial_all_cds.fasta -l A -a Sample1_megahit.annotation_bwa.sam -o Sample1_salmon*
And I encounter these errors:
*" Transcript appears twice in the transcript FASTA file" and "Transcript appears in the reference but did not appear in the BAM."*
Is there a way to, in command line,
- process my transcript FASTA file so that there are no duplicates and
- process my reference file so that it's filtered, containing only transcripts that are also in my SAM file?
• 1,940 views
•
link
written
by
pearl2070 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Low mapping rate using SALMON
written by Pegasus •
-
Salmon Quantification in Alignment based-mode
written by Patadu94 •
-
Salmon error (ex. transcript id appears in the BAM head, but was not in the provided FASTA file)
written by tonyya21 •
-
Meaningful Names in BWA Alignment For Salmon?
written by pearl2070 •
-
Running salmon in alignment mode?
written by pearl2070 •
-
Kallisto and Salmon give different results
written by mailard •
-
Trouble with choosing the right reference transcriptome file for mm10
written by wyoo083 •
-
Salmon quant in Galaxy vs Command Line
written by NewtoBioinfo •
-
Salmon Counts differ in Alignment and Quasi-Mapping mode
written by MaxF
-
loaded a sorted bam file in IGV but I can't see anything
written by saadleeshehreen
When you generated your BAM did you align your reads to the genome or transcriptome? Salmon requires alignment to the transcriptome. You may also want to include your code for that part.
Ah, I aligned the reads to a reference database (microbial_all_cds.fasta). For BWA, I ran:
Followed by:
Should I instead be doing:
And then:
If
Sample1_megahit.contigs.fais your assembled transcriptome then yes.Got it, thank you!
Prior reference: Running salmon in alignment mode?