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

Error in RG argument does not stop minimap2 #541

Closed
dariober opened this issue Jan 6, 2020 · 1 comment
Closed

Error in RG argument does not stop minimap2 #541

dariober opened this issue Jan 6, 2020 · 1 comment

Comments

@dariober
Copy link

dariober commented Jan 6, 2020

If an invalid RG line is passed to option -R, minimap detects the error but keeps going without adding RG tags. Maybe better to make it quit? E.g:

minimap2 -V
2.17-r941

minimap2 -R 'foobar' -a some.fasta some.fastq > /dev/null

[M::mm_idx_gen::2.132*1.00] collected minimizers
[M::mm_idx_gen::2.469*1.25] sorted minimizers
[ERROR] the read group line is not started with @RG             <============
[M::main::2.469*1.25] loaded/built the index for 36 target sequence(s)
[M::mm_mapopt_update::2.633*1.24] mid_occ = 32
[M::mm_idx_stat] kmer size: 15; skip: 10; is_hpc: 0; #seq: 36
[M::mm_idx_stat::2.762*1.23] distinct minimizers: 5072250 (89.73% are singletons); average occurrences: 1.216; average spacing: 5.327
[M::worker_pipeline::8.651*2.27] mapped 4000 sequences
[M::main] Version: 2.17-r941
[M::main] CMD: minimap2 -R foobar -a ref/TriTrypDB-45_LmajorFriedlin_Genome.fasta /export/projects/III-data/wcmp_bioinformatics/db291g/data/20190906_catarina_nanopore/20190902_1335_MN23371_FAK87904_27b810d7/fastq_pass/FAK87904_a513394c0ef96eff721553a5edf39959f27553c9_44.fastq.gz
[M::main] Real time: 8.718 sec; CPU: 19.723 sec; Peak RSS: 0.641 GB
@lh3
Copy link
Owner

lh3 commented Jan 6, 2020

Maybe better to make it quit?

Fair enough. I will consider this in the next release.

@lh3 lh3 closed this as completed in 3c0d05d Jan 7, 2020
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants