Attention une mise à jour du serveur va être effectuée le lundi 17 mai entre 13h et 13h30. Cette mise à jour va générer une interruption du service de quelques minutes.

Commit d3efe5e2 authored by Mikaël Salson's avatar Mikaël Salson

revcomp.should_get: Test that the same e-value appears twice

We prefix the e-values with a string "e-value:" to make sure that we are matching e-values.
parent 88c3213c
!LAUNCH: $LAUNCHER ../../vidjil -k 9 -G ../../germline/IGH -K -c clones ../../data/revcomp.fa ; grep 'X.X.X' out/revcomp.affects | sed 's/[^X]//g' | sort -u ; grep '#>' out/revcomp.affects | sed 's/.*SEG...//' | cut -f 1 -d' ' | nl -bn -s 'evalue: ' | uniq -c
!LAUNCH: $LAUNCHER ../../vidjil -k 9 -G ../../germline/IGH -K -c clones ../../data/revcomp.fa ; grep 'X.X.X' out/revcomp.affects | sed 's/[^X]//g' | sort -u ; grep '#>' out/revcomp.affects | sed 's/.*SEG.../e-value:/' | cut -f 1 -d' '
$ Segments both reads, normal and reverse
1:junction detected in 2 reads
......@@ -10,4 +10,4 @@ $ Only one window
1:==> found 1 ..-window
$ Same e-value for both reads
1: 2 .* evalue: .[.].*e-..
2: e-value:[0-9][.][0-9]+e-..
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