1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
Set up
$ mkdir -p $OUTDIR
Test dataset.xml as input
$ $BLASR_EXE $DATDIR/test_dataset/chunking.subreadset.xml $DATDIR/ecoli_reference.fasta --unaligned $OUTDIR/unaligned.txt --noPrintUnalignedSeqs --concordant 1>/dev/null && echo $?
[INFO]* (glob)
[INFO]* (glob)
0
$ head -5 $OUTDIR/unaligned.txt
m150404_101626_42267_c100807920800000001823174110291514_s1_p0/55/0_1380
m150404_101626_42267_c100807920800000001823174110291514_s1_p0/55/1432_3136
m150404_101626_42267_c100807920800000001823174110291514_s1_p0/480/11699_11988
m150404_101626_42267_c100807920800000001823174110291514_s1_p0/480/12033_13456
m150404_101626_42267_c100807920800000001823174110291514_s1_p0/480/13519_14067
$ $BLASR_EXE $DATDIR/ecoli_subset.fasta $DATDIR/ecoli_reference.fasta --unaligned $OUTDIR/unaligned.txt --noPrintUnalignedSeqs 1>/dev/null && echo $?
[INFO]* (glob)
[INFO]* (glob)
0
$ head -5 $OUTDIR/unaligned.txt
m121004_000921_42130_c100440700060000001523060402151341_s1_p0/8/2724_3021
m121004_000921_42130_c100440700060000001523060402151341_s1_p0/13/0_278
m121004_000921_42130_c100440700060000001523060402151341_s1_p0/13/327_954
m121004_000921_42130_c100440700060000001523060402151341_s1_p0/13/1004_1580
m121004_000921_42130_c100440700060000001523060402151341_s1_p0/13/1625_2202
|
