Hello Everyone,
I have made bam files using hisat2 and then I am trying to use these bam files for salmon alignment.
hisat2-build genome.fa hisat
hisat2 -x hisat -1 R1_paired.fastq.gz -2 R2_paired.fastq.gz -S out.sam
samtools view -bS out.sam > out.bam
java -Xms1g -Xmx3g -jar picard.jar MarkDuplicates \
I=out.bam \
O=marked_duplicates.bam \
M=metrics.txt
salmon quant -t transcripts.fasta -l A -a marked_duplicates.bam -o out.sf
But it gives an error:
Please provide a reference FASTA file that includes all targets present in the BAM header
If you have access to the genome FASTA and GTF used for alignment
consider generating a transcriptome fasta using a command like:
gffread -w output.fa -g genome.fa genome.gtf
Thank you for any help!
1 answer
You aligned your reads to the genome with HISAT2, and you are trying to quantify the transcriptome with salmon. The references are not the same, and salmon does not (currently) support consuming genome-aligned BAM files for transcriptome quantification. You either need to align the reads directly against the transcriptome (which could be done with HISAT2, or with a tool like Bowtie2), or use the STAR aligner and take advantage of it's built-in option to project genomic alignments onto transcriptomic coordinates using the --quantMode TranscriptomeSAM option.
No, this will simply prevent HISAT2 from aligning reads outside of annotated transcripts. These will still be spliced alignments in genomic coordinates. As far as I know, HISAT2 has no option to project genomic alignments to transcriptomic coordinates, so the only way to accomplish this with HISAT2 would be to align the reads directly against the transcriptome rather than the genome.
Log in to answer this question.
-
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 •
-
Salmon in alignment mode - what do they mean in the documentation, when they say to use transcripto…
written by Shraddha •
-
RNAseq - GBFF cannot be converted to GTF
written by sacdefruits •
-
3' RNA-seq (Lexogen Quant-seq) Mapping Problems with Salmon
written by bdy8 •
-
error in MarkDuplicates Picard
written by evelyn
-
hisat2 error index out of range
written by evelyn
-
salmon light weight alignment mode
written by evelyn
-
Sanity Check on RNA-Seq featureCounts / DESeq2
written by gwebste7 •
