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

No -I option #40

Open
soisa001 opened this issue Aug 19, 2023 · 3 comments
Open

No -I option #40

soisa001 opened this issue Aug 19, 2023 · 3 comments

Comments

@soisa001
Copy link

When running winnowmap, the -I option is not recognized.
e.g. after generating the repetitive_k15.txt with meryl:

winnowmap -W repetitive_k15.txt -a -x map-pb -Y -L --eqx --cs -I 32G ref.fa.gz reads.fastq.gz | samtools view -hb | samtools sort -@8 > alignment_sorted.bam

Yields the following error:

[ERROR] unknown option in "-I"

The -I option is needed for a multi-part index.
Thanks.

@cjain7
Copy link
Contributor

cjain7 commented Aug 21, 2023

Sorry, multi-part indexing is not supported yet.

@diego-rt
Copy link

Hi @cjain7

Does this mean that it's not possible to map to genomes larger than 4G while getting accurate mapQs? For minimap2 this would be the case in the absence of the -I flag.

Thanks!

@skoren
Copy link
Member

skoren commented May 23, 2024

I saw this change was added post the last v2.0.3 release version so the condo-installed versions allow using the -I option. I do see slight differences in alignments when increasing -I on genomes w/>4gb genome size. I wanted to confirm if it is safe to use this option assuming no saved index is used or was it removed because it was not working correctly in v2.0.3 as well?

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

4 participants