Commit aaeea2c6 authored by Mathieu Giraud's avatar Mathieu Giraud

tests: testing 'unable to segment' on some germlines when using -c -v -g

parent a0d13502
!LAUNCH: ../../vidjil -G ../../germline/IGH -d -c segment ../../data/segment_simul.fa | grep '^>'
!LAUNCH: ../../vidjil -v -g ../../germline -d -c segment ../../data/segment_simul.fa | grep '^[>#]'
$ First sequence, easy segmentation (no error, few deletions at the windows, small N)
......@@ -12,5 +12,8 @@ $ First sequence, easy segmentation (no error, few deletions at the windows, sma
# -2/CTTC/-3 -1/GTTA/5
1:^>seq1 \\+ VDJ 0 72 77 89 94 139 IGHV5-10-1.0[1-4] 2/CTTC/3 IGHD1-14.01 1/GTTA/5 IGHJ5.01
$ Do not manage to segment on TRG
2:TRG: unable to segment
$ seq1_polyT is the same as seq1 but starts and ends with 10Ts, the recombination is in the middle. start position should not be 0 and end position should not be length-1
f1:^>seq1_polyT \\+ VDJ 10 82 87 99 104 149 IGHV5-10-1.0[1-4] 2/CTTC/3 IGHD1-14.01 1/GTTA/5 IGHJ5.01
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