w-100-short-shift.should 1.02 KB
Newer Older
1
!LAUNCH: $VIDJIL_DIR/$EXEC -K -e 10 -r 1 -y 0 -s '#####-#####' -w 100 -V $VIDJIL_DIR/germline/homo-sapiens/IGHV.fa -D $VIDJIL_DIR/germline/homo-sapiens/IGHD.fa -J $VIDJIL_DIR/germline/homo-sapiens/IGHJ.fa -b stanford-w100 $VIDJIL_DATA/Stanford_S22.fasta ; grep 'w.*/.* seed' out/stanford-w100.affects; grep -v 'UNSEG strand' out/stanford-w100.vidjil
Mikaël Salson's avatar
Mikaël Salson committed
2

Mathieu Giraud's avatar
Mathieu Giraud committed
3 4 5
$ Look for windows up to 100bp
1:windows up to 100bp

Mikaël Salson's avatar
Mikaël Salson committed
6
$ Find the good number of "too short sequences" for windows of size 100
7 8 9
1: UNSEG too short w ->      0

$ Some reads have shortened or shifted windows
10
1: SEG changed w -> 1370
Mikaël Salson's avatar
Mikaël Salson committed
11

Mathieu Giraud's avatar
Mathieu Giraud committed
12
$ Most changed windows are just shifted
13
915: w100/-5
Mathieu Giraud's avatar
Mathieu Giraud committed
14 15 16 17 18
360: w100/-10

$ Some changed windows are lighlty shortened
63: w95/-10
29: w90/-10
19

Mikaël Salson's avatar
Mikaël Salson committed
20
$ Find the good number of segmented sequences (including "too short sequences")
21
1: junction detected in 13153 reads .100%.
Mikaël Salson's avatar
Mikaël Salson committed
22

Mathieu Giraud's avatar
Mathieu Giraud committed
23
$ Find the good number of windows in Stanford S22
24
1: found 11836 windows in 13153 reads
25

Mathieu Giraud's avatar
Mathieu Giraud committed
26
$ Find the correct number of clones with shifted of shortened windows
27 28
1245: "Short or shifted window"
1245: "W50"
Mathieu Giraud's avatar
Mathieu Giraud committed
29

30