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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
|
@HD VN:1.4
@SQ SN:insert LN:599
@SQ SN:ref1 LN:45
@SQ SN:ref2 LN:40
@SQ SN:ref3 LN:4
@RG ID:test_input_1_a PG:donkey
@RG ID:test_input_1_b PG:llama
@RG ID:test_input_1_c
@CO
@CO Do you know?
@CO Do you know?
@CO Another comment from test_input_1_c
r000 99 insert 50 30 10M = 80 30 ATTTAGCTAC AAAAAAAAAA PG:Z:bull RG:Z:test_input_1_a
r000 211 insert 80 30 10M = 50 -30 CCCAATCATT AAAAAAAAAA PG:Z:bull RG:Z:test_input_1_a
r001 163 ref1 7 30 8M4I4M1D3M = 37 39 TTAGATAAAGAGGATACTG * XX:B:S,12561,2,20,112 YY:i:100 RG:Z:test_input_1_a
r005 163 ref1 7 30 8M4I4M1D3M = 37 39 TTAGATAAAGAGGATACTG * XX:B:S,12561,2,20,112 YY:i:100 PG:Z:donkey-3A2CCEF5 RG:Z:test_input_1_b
r008 163 ref1 7 30 8M4I4M1D3M = 37 39 TTAGATAAAGAGGATACTG * XX:B:S,12561,2,20,112 YY:i:100 PG:Z:donkey-2B019719 RG:Z:test_input_1_c
r002 0 ref1 9 30 1S2I6M1P1I1P1I4M2I * 0 0 AAAAGATAAGGGATAAA * XA:Z:abc XB:i:-10 RG:Z:test_input_1_a
r003 0 ref1 9 30 5H6M * 0 0 AGCTAA * RG:Z:test_input_1_a
r006 0 ref1 9 30 1S2I6M1P1I1P1I4M2I * 0 0 AAAAGATAAGGGATAAA * XA:Z:abc XB:i:-10 PG:Z:donkey-3A2CCEF5 RG:Z:test_input_1_b
r007 0 ref1 9 30 5H6M * 0 0 AGCTAA * PG:Z:donkey-3A2CCEF5 RG:Z:test_input_1_b
r009 0 ref1 9 30 1S2I6M1P1I1P1I4M2I * 0 0 AAAAGATAAGGGATAAA * XA:Z:abc XB:i:-10 PG:Z:donkey-2B019719 RG:Z:test_input_1_c
r010 0 ref1 9 30 5H6M * 0 0 AGCTAA * PG:Z:donkey-2B019719 RG:Z:test_input_1_c
r004 0 ref1 16 30 6M14N1I5M * 0 0 ATAGCTCTCAGC * RG:Z:test_input_1_a
r007 0 ref1 16 30 6M14N1I5M * 0 0 ATAGCTCTCAGC * PG:Z:donkey-3A2CCEF5 RG:Z:test_input_1_b
r010 0 ref1 16 30 6M14N1I5M * 0 0 ATAGCTCTCAGC * PG:Z:donkey-2B019719 RG:Z:test_input_1_c
r003 16 ref1 29 30 6H5M * 0 0 TAGGC * RG:Z:test_input_1_a
r006 16 ref1 29 30 6H5M * 0 0 TAGGC * PG:Z:donkey-3A2CCEF5 RG:Z:test_input_1_b
r009 16 ref1 29 30 6H5M * 0 0 TAGGC * PG:Z:donkey-2B019719 RG:Z:test_input_1_c
r001 83 ref1 37 30 9M = 7 -39 CAGCGCCAT * RG:Z:test_input_1_a
r005 83 ref1 37 30 9M = 7 -39 CAGCGCCAT * PG:Z:donkey-3A2CCEF5 RG:Z:test_input_1_b
r008 83 ref1 37 30 9M = 7 -39 CAGCGCCAT * PG:Z:donkey-2B019719 RG:Z:test_input_1_c
x1 0 ref2 1 30 20M * 0 0 AGGTTTTATAAAACAAATAA * PG:Z:bull RG:Z:test_input_1_a
x7 0 ref2 1 30 20M * 0 0 AGGTTTTATAAAACAAATAA * PG:Z:bull-055424A4 RG:Z:test_input_1_b
x10 0 ref2 1 30 20M * 0 0 AGGTTTTATAAAACAAATAA * PG:Z:bull-6ADB4A65 RG:Z:test_input_1_c
x2 0 ref2 2 30 21M * 0 0 GGTTTTATAAAACAAATAATT ????????????????????? PG:Z:bull RG:Z:test_input_1_a
x8 0 ref2 2 30 21M * 0 0 GGTTTTATAAAACAAATAATT ????????????????????? PG:Z:bull-055424A4 RG:Z:test_input_1_b
x11 0 ref2 2 30 21M * 0 0 GGTTTTATAAAACAAATAATT ????????????????????? PG:Z:bull-6ADB4A65 RG:Z:test_input_1_c
x3 0 ref2 6 30 9M4I13M * 0 0 TTATAAAACAAATAATTAAGTCTACA ?????????????????????????? PG:Z:bull RG:Z:test_input_1_a
x9 0 ref2 6 30 9M4I13M * 0 0 TTATAAAACAAATAATTAAGTCTACA ?????????????????????????? PG:Z:bull-055424A4 RG:Z:test_input_1_b
x12 0 ref2 6 30 9M4I13M * 0 0 TTATAAAACAAATAATTAAGTCTACA ?????????????????????????? PG:Z:bull-6ADB4A65 RG:Z:test_input_1_c
x4 0 ref2 10 30 25M * 0 0 CAAATAATTAAGTCTACAGAGCAAC ????????????????????????? PG:Z:bull RG:Z:test_input_1_a
x10 0 ref2 10 30 25M * 0 0 CAAATAATTAAGTCTACAGAGCAAC ????????????????????????? PG:Z:bull-055424A4 RG:Z:test_input_1_b
x13 0 ref2 10 30 25M * 0 0 CAAATAATTAAGTCTACAGAGCAAC ????????????????????????? PG:Z:bull-6ADB4A65 RG:Z:test_input_1_c
x5 0 ref2 12 30 24M * 0 0 AATAATTAAGTCTACAGAGCAACT ???????????????????????? PG:Z:bull RG:Z:test_input_1_a
x11 0 ref2 12 30 24M * 0 0 AATAATTAAGTCTACAGAGCAACT ???????????????????????? PG:Z:bull-055424A4 RG:Z:test_input_1_b
x14 0 ref2 12 30 24M * 0 0 AATAATTAAGTCTACAGAGCAACT ???????????????????????? PG:Z:bull-6ADB4A65 RG:Z:test_input_1_c
x6 0 ref2 14 30 23M * 0 0 TAATTAAGTCTACAGAGCAACTA ??????????????????????? RG:Z:test_input_1_a
x12 0 ref2 14 30 23M * 0 0 TAATTAAGTCTACAGAGCAACTA ??????????????????????? PG:Z:bull-055424A4 RG:Z:test_input_1_b
x15 0 ref2 14 30 23M * 0 0 TAATTAAGTCTACAGAGCAACTA ??????????????????????? PG:Z:bull-6ADB4A65 RG:Z:test_input_1_c
u1 4 * 0 30 23M * 0 0 TAATTAAGTCTACAGAAAAAAAA ??????????????????????? RG:Z:test_input_1_a
u2 4 * 0 30 * * 0 0 TAATTAAGTCTACAGAAAAAAAA ??????????????????????? RG:Z:test_input_1_c
|
