large-r.should 520 Bytes
Newer Older
1
!LAUNCH: (i=1; while [ $i -le 100000 ]; do echo '>read' ; echo ccgtgtattactgtgcgagagagctgaatacttccagcactg ; i=$((i+1)); done ;) > same-igh-100k.fa ; $LAUNCHER $VIDJIL_DIR/$EXEC $EXTRA -g $VIDJIL_DIR/germline/homo-sapiens.g:IGH -r 5000  -w 15  same-igh-100k.fa; rm -f same-igh-100k.fa
2 3

$ Find a unique clone with all reads
4
1:>clone-001--IGH--0100000--100.--window
5 6

$ Find a representative
7
1:>clone-001--IGH--0100000--100.--read-.0,41.
8

9 10 11 12
$ Compute the diversity. No diversity here.
1: 1   0.000
1:  E = 0.000
1: Ds = 0.000