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

BWA and short read mapping #424

Open
alessandromacc opened this issue Jul 10, 2024 · 0 comments
Open

BWA and short read mapping #424

alessandromacc opened this issue Jul 10, 2024 · 0 comments

Comments

@alessandromacc
Copy link

alessandromacc commented Jul 10, 2024

Dear developers,

I am currently working with miRNA-seq reads from Illumina platforms, which are very short reads, ranging from 16 to roughly 30 nucleotides, and I want to map them onto a reference database for miRNAs, much shorter than the reference genome.

Is there a specific reason, in your opinion, to claim that BWA can outperform other popular mappers like Bowtie, when it comes to this kind of short reads?

Are the default parameters suitable for miRNA sequencing analysis in your opinion?

Also, I am currently interested in counting all the alternative hits of each mapped read (those listed in the XA tag in the SAM records) for a differential expression analysis. After doing some research, I don't believe it to be possible to instruct BWA to write a separate SAM record for an alternative mapping of a read (say like Bowtie does by default).
Can you suggest how to count all the alignments reported by BWA in the SAM file (primary and alternative)?

Thank you kindly for your help!

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

1 participant