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

Empty alignment gaf output #93

Open
abanerjee10 opened this issue Mar 3, 2023 · 4 comments
Open

Empty alignment gaf output #93

abanerjee10 opened this issue Mar 3, 2023 · 4 comments

Comments

@abanerjee10
Copy link

I get empty output gaf files when aligning reads to genome graphs. How do I resolve this issue?

Data:
vg_toy_x.zip
Command:
minigraph -t 24 -cx lr vg_toy_x.gfa vg_toy_x.fastq

@colindaven
Copy link

You've set parameters for long reads
-cx lr

but the toy example has one short read, so you need 
-cx sr

@abanerjee10
Copy link
Author

abanerjee10 commented Mar 7, 2023

Sorry, I tried both arguments and pasted the command with long read mode. Both give me empty output gafs.

Data:
vg_toy_x.zip

Commands:
minigraph -t 24 -cx sr vg_toy_x.gfa vg_toy_x.fastq
minigraph -t 24 -cx lr vg_toy_x.gfa vg_toy_x.fastq

Output:
[M::main::0.001*3.34] loaded the graph from "vg_toy_x.gfa" [M::mg_index::0.003*4.91] indexed the graph [M::mg_opt_update::0.004*4.82] occ_max1=50; lc_max_occ=2 [M::worker_pipeline::0.005*3.85] mapped 2 sequences [M::main] Version: 0.20-r559 [M::main] CMD: minigraph -t 24 -cx lr vg_toy_x.gfa vg_toy_x.fastq [M::main] Real time: 0.005 sec; CPU: 0.019 sec; Peak RSS: 0.002 GB

I have the same error for larger genome graphs with short and long reads that I can't upload due to file sizes.

@Szhiha
Copy link

Szhiha commented Nov 13, 2023

I have the same issue.

@aleixcanalda
Copy link

@Szhiha I had the same problem but apparently the problem is with the v20 release (@lh3 not sure if this was intentional or not?). If you use minigraph-v19, -xsr does work

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