Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Quantification in Alignment mode for Nanopore Data #903

Open
Harshangda opened this issue Dec 8, 2023 · 2 comments
Open

Quantification in Alignment mode for Nanopore Data #903

Harshangda opened this issue Dec 8, 2023 · 2 comments

Comments

@Harshangda
Copy link

I am using Salmon (Salmon 1.3 installed through homebrew on mac ventura 13.5.2) for quantification in alignment mode. I have bam files (aligned to transcriptome) for nanopore data from Minimap2, which has primary and secondary alignments both (50% of total mapped reads are primary and other 50% are secondary alignments). Salmon is only performing quantification for primary alignments, but I also want to include the secondary alignments.

  1. Is there any argument I could use to include secondary alignments for quantification in Salmon Alignment mode?
  2. Is it a good practice to include secondary alignments? 3. Are there any other important flags I should use which might be specific to "Nanopore" or "Splicing Analysis".

This is the code I am using at the moment:

salmon quant -t gencode.vM33.transcripts.fa -l A -g gencode.vM32.annotation.gtf -a barcode07.bam -o barcode_07_salmon_quant

@Harshangda Harshangda changed the title Quantifyication in Alignment mode for Nanopore Data Quantification in Alignment mode for Nanopore Data Dec 8, 2023
@alexyfyf
Copy link

alexyfyf commented Jan 8, 2024

also interested to know how Salmon uses secondary alignment. Because I found this tutorial https://combine-lab.github.io/salmon-tutorials/2021/ont-long-read-quantification/ actually includes secondary alignments.
And based on my experience, secondary alignments are used by Salmon, because when I give a BAM before and after removing secondary (-F 256 flag in samtools), the results are different.

@Harshangda
Copy link
Author

Harshangda commented Jan 10, 2024 via email

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants