Commit 17056853 authored by Mikaël Salson's avatar Mikaël Salson
Browse files

chimera.should-get: We segment again on IGL by chance

This is due to the Aho-Corasick automaton where we now have an index load for
the V gene and for the J gene. There is only one matching kmer from IGLJ but
that's enough to have an evalue < 1 because the index load for IGLJ is very
low and because there are only two reads in the dataset.

We could put a -e 1e-1 on the command line to make the test pass. Would this
be hacky or not?
parent 005881f8
!LAUNCH: $VIDJIL_DIR/vidjil $VIDJIL_DEFAULT_OPTIONS -A -uU -g $VIDJIL_DIR/germline $VIDJIL_DIR/data/chimera-trg.fa
$ Do not segment on IGL by chance
1:IGL .* -> .* 0
f1:IGL .* -> .* 0
$ Report the sequence as AMBIGUOUS
1:UNSEG ambiguous -> 1
f1:UNSEG ambiguous -> 1
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment