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

Adherence to the "-G" parameter #1265

Open
sanyalab opened this issue Nov 28, 2024 · 0 comments
Open

Adherence to the "-G" parameter #1265

sanyalab opened this issue Nov 28, 2024 · 0 comments

Comments

@sanyalab
Copy link

Hello,

In minimap the -G parameter decides the maximum gap on the reference during spliced alignments. So I take this as either intron length or intergenic length. Is my assumption correct? This is set to 200000 by default. I set the -G value to 50000 for my species, as the 90th percentile value of intron length was 48000. It turns out that the gene I was interested in had a single 60000 bp intron. Therefore, FLNC reads that should have aligned did not, as I set -G 50000. When I ran minimap with default -G the reads aligned fine. Therefore, a few questions.

  1. Is it advisable to reduce this size closer to the anticipated intron size of a species?
  2. How strictly does the tool adhere to the -G parameter value? You mention that this has an effect on the chaining and alignment band width. Which are these parameters? and what do they control?

Thanks
Abhijit

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