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

Option -P leads to signal 8 #200

Open
subwaystation opened this issue Sep 14, 2023 · 1 comment
Open

Option -P leads to signal 8 #200

subwaystation opened this issue Sep 14, 2023 · 1 comment

Comments

@subwaystation
Copy link
Contributor

I was trying to use wfmash in mapping mode to generate all versus all mappings of 2k e. coli. @AndreaGuarracino you already have the FASTA.

However, when no minmers are found, I get a signal 8:

[mashmap] MashMap v3.1.1
[mashmap] Reference = [/home-link/afaaj01/data/ecoli_2146/ecoli_2146.fa.gz]
[mashmap] Query = [/home-link/afaaj01/data/ecoli_2146/ecoli_2146.fa.gz]
[mashmap] Kmer size = 19
[mashmap] Sketch size = 298
[mashmap] Segment length = 5000 (read split allowed)
[mashmap] Block length min = 25000
[mashmap] Chaining gap max = 20000
[mashmap] Mappings per segment = 2145
[mashmap] Percentage identity threshold = 90%
[mashmap] Skip self mappings
[mashmap] Hypergeometric filter w/ delta = 0 and confidence 0.999
[mashmap] Mapping output file = /dev/stdout
[mashmap] Filter mode = 1 (1 = map, 2 = one-to-one, 3 = none)
[mashmap] Execution threads  = 20
[mashmap::skch::Sketch::build] minmer windows picked from reference = 0
[mashmap::skch::Sketch::index] unique minmers = 0
[mashmap::skch::Sketch::computeFreqHist] No minmers.
[wfmash::map] time spent computing the reference index: 0.0131097 sec
[mashmap::skch::Map::mapQuery] mapped  0.00% @ 0.00e+00 bp/s elapsed: 00:00:00:00 remain: 00:00:00:00
Command terminated by signal 8

Is there a way to circumvent this? Is this intended?

I am still wondering, if running wfmash in all versus all mapping mode on all sample at once will give the equivalent results to running wfmash in all versus all mapping mode for each prefix in the given sequences.

Thanks for any tips!

@subwaystation
Copy link
Contributor Author

In total, I can observe this for 3 sequences:

wfmash_map_0110_ecoli_2146_s5k_p90.BL21-Gold(DE3)pLysS_AG_f399789#1#NC_012947.1.time:[mashmap::skch::Sketch::computeFreqHist] No minmers.
wfmash_map_0110_ecoli_2146_s5k_p90.clone_D_i14_e789781#1#NC_017652.1.time:[mashmap::skch::Sketch::computeFreqHist] No minmers.
wfmash_map_0110_ecoli_2146_s5k_p90.clone_D_i2_32d971c#1#NC_017651.1.time:[mashmap::skch::Sketch::computeFreqHist] No minmers.

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