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

Mapping 2 wheat genomes #1115

Open
mfruzan opened this issue Sep 8, 2023 · 1 comment
Open

Mapping 2 wheat genomes #1115

mfruzan opened this issue Sep 8, 2023 · 1 comment

Comments

@mfruzan
Copy link

mfruzan commented Sep 8, 2023

I used below command to map one wheat genome against reference in a 32 cpu and 1TByte machine.

minimap2 -x asm5 -f100 -K4g --cs=long -t 32 --secondary=no iwgsc_refseqv2.1_assembly.fa arinalr_sm.fa | paftools.js view -f maf - >cs_arina.maf;

It only uses 6 CPUs, after 3 hours it generated 7Gbyte of maf file (up to chromosome 2) and remained at following line:
[M::main] Version: 2.22-r1105-dirty

Is there any other way to speed up, such as increasing k-mer size or changing -f, or even align each chromosome of query genome one by one?

Thanks ,
Mario

@xbarvazx
Copy link

try -f 0.05

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