Hello. everyone
I am trying to count aligned reads using Salmon in alignment-based mode.
My command is:
salmon quant -t transcript.fa -l A -a aligned_result.bam -o salmon_quant
I am encountering the following error:
Transcript chr8 appeared in the BAM header, but was not in the provided FASTA file.
The preprocessing steps I performed are as follows:
I obtained the transcript fasta file from https://www.gencodegenes.org/human/ and used it for analysis. I also tried creating the transcript fasta file using gffread, but the same error persists.
Could this be a problem with the transcript fasta file or the BAM file?
Thank you in advance for your help.
•
written
by
tonyya21 •
You are using different sources of annotation. NCBI genome (which might use chromosome identifiers such as 1,2,3) and GENCODE for the fasta which uses chr1,2,3... Stick with one source.
By the way, steps 4-8 are not necessary to run salmon, and actually you could skip STAR entirely if the only goal is to run salmon, and use salmon for the mapping itself.
Oh, that was my mistake. I accidentally used gencode instead of ncbi in step 1. then, I'll proceed by running salmon for the mapping itself. Thank you
Log in to answer this question.
