Commit 0fce02b9 authored by Ryan Herbert's avatar Ryan Herbert
Browse files

Merge branch 'dev' of gitlab.inria.fr:vidjil/vidjil into dev

parents 2e0dbfcf 132c251c
......@@ -8,21 +8,20 @@
!LAUNCH: (file1=$(mktemp); \
(for i in {1..10}; do \
echo '>seq'\\$i; echo ggctggagtgggtttcatacattagtagtaatagtggtgccatatactacgcagactctgtgaagggccgattcaccatctccagaaacaatgccaaggactcactgtatctgcaaatgaacagcctgagagccgaggacacggctgtgtattactgtgcgagagcgatcccccggtattactatgatact\\$(cat /dev/urandom | tr -dc 'acgt' | head -c 4)ggcccaaacgactactggggccagggaaccctggtcaccgtctcctcag; \
echo '>seq'\\$i; echo ggctggagtgggtttcatacattagtagtaatagtggtgccatatactacgcagactctgtgaagggccgattcaccatctccagaaacaatgccaaggactcactgtatctgcaaatgaacagcctgagagccgaggacacggctgtgtattactgtgcgagagcgatcccccggtattactatgatact\\$(cat /dev/urandom | LC_CTYPE=C tr -dc 'acgt' | head -c 4)ggcccaaacgactactggggccagggaaccctggtcaccgtctcctcag; \
done;\
cat $VIDJIL_DIR/data/buggy-D.fa) > \\$file1; \
$VIDJIL_DIR/vidjil $VIDJIL_DEFAULT_OPTIONS -d -r 1 -w 60 -z 100 -G $VIDJIL_DIR/germline/IGH \\$file1; \
file2=$(mktemp); \
(for i in {1..5}; do \
echo '>seq'\\$i; echo ggctggagtgggtttcatacattagtagtaatagtggtgccatatactacgcagactctgtgaagggccgattcaccatctccagaaacaatgccaaggactcactgtatctgcaaatgaacagcctgagagccgaggacacggctgtgtattactgtgcgagagcgatcccccggtattactatgatact\\$(cat /dev/urandom | tr -dc 'acgt' | head -c 4)ggcccaaacgactactggggccagggaaccctggtcaccgtctcctcag; \
echo '>seq'\\$i; echo ggctggagtgggtttcatacattagtagtaatagtggtgccatatactacgcagactctgtgaagggccgattcaccatctccagaaacaatgccaaggactcactgtatctgcaaatgaacagcctgagagccgaggacacggctgtgtattactgtgcgagagcgatcccccggtattactatgatact\\$(cat /dev/urandom | LC_CTYPE=C tr -dc 'acgt' | head -c 4)ggcccaaacgactactggggccagggaaccctggtcaccgtctcctcag; \
done;\
cat $VIDJIL_DIR/data/buggy-D.fa) > \\$file2; \
$VIDJIL_DIR/vidjil $VIDJIL_DEFAULT_OPTIONS -d -r 1 -w 60 -z 100 -G $VIDJIL_DIR/germline/IGH \\$file2;\
$VIDJIL_DIR/vidjil $VIDJIL_DEFAULT_OPTIONS -d -r 1 -w 60 -z 100 -G $VIDJIL_DIR/germline/IGH $VIDJIL_DIR/data/buggy-D.fa) | \
sort | uniq -c
$VIDJIL_DIR/vidjil $VIDJIL_DEFAULT_OPTIONS -d -r 1 -w 60 -z 100 -G $VIDJIL_DIR/germline/IGH $VIDJIL_DIR/data/buggy-D.fa)
$ Three times the same window
1: 3 TGTGCGGGATCTTCGTCCTCTTATCATAATAATGGTTTTTTGGCGGGGGAGTCATGGGGC
3: TGTGCGGGATCTTCGTCCTCTTATCATAATAATGGTTTTTTGGCGGGGGAGTCATGGGGC
$ Three times the same segmentation
1: 3 IGHV4-34.*IGHD.*IGHJ4
3: IGHV4-34.*IGHD.*IGHJ4
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