Changeset: 0a3325139576 for MonetDB
URL: http://dev.monetdb.org/hg/MonetDB?cmd=changeset;node=0a3325139576
Modified Files:
sql/backends/monet5/bam/Tests/test.SQL.py
sql/backends/monet5/bam/Tests/test.stable.out
Branch: bamloader
Log Message:
Included tests for sam_export function
diffs (truncated from 508 to 300 lines):
diff --git a/sql/backends/monet5/bam/Tests/test.SQL.py
b/sql/backends/monet5/bam/Tests/test.SQL.py
--- a/sql/backends/monet5/bam/Tests/test.SQL.py
+++ b/sql/backends/monet5/bam/Tests/test.SQL.py
@@ -43,6 +43,8 @@ for match in p.finditer(out):
# Now we will execute all benchmark queries on all BAM files in the bam.files
table
# and output all data contained in the aux table.
+# Furthermore, we transfer every file to the export table, use sam_export to
write
+# the contents to a SAM file and then print the raw contents of this SAM file
c = new_client()
set_var(c, 'rname_1_3', 'chr22', numeric=False)
set_var(c, 'pos_1_3_1', 1000000)
@@ -56,6 +58,7 @@ set_var(c, 'rname_2_10', 'chr22', numeri
set_var(c, 'pos_2_10', 80000000)
set_var(c, 'distance_2_12', 10000000)
+output_files = []
for f in files_to_test:
# benchmark 1
for uc in range(1, 6):
@@ -71,6 +74,15 @@ for f in files_to_test:
#write all aux data
c.stdin.write("SELECT * FROM bam.alignments_extra_%d;"% f[0]);
+ #load into export table
+ c.stdin.write("INSERT INTO bam.export (SELECT qname, flag, rname, pos,
mapq, cigar, rnext, pnext, tlen, seq, qual FROM bam.%salignments_%d);"\
+ % (('unpaired_all_' if f[1] == 1 else ''), f[0]))
+ output_files.append(os.path.join(SRCDIR, ("output_%d.sam"% f[0])))
+ c.stdin.write("CALL sam_export('%s');"% output_files[-1])
+
+
+
+
out, err = c.communicate()
# The output will contain explicit file_ids in the table names, making the
table names
@@ -83,3 +95,11 @@ out = (re.subn(p, replace, out))[0]
sys.stdout.write(out)
sys.stderr.write(err)
+
+
+# All that is left to do is write the contents of the exported SAM files and
delete them
+for f in output_files:
+ sys.stdout.write("\n\nContents of exported SAM file '%s':\n"% f)
+ for ln in open(f):
+ sys.stdout.write(ln)
+ os.remove(f)
\ No newline at end of file
diff --git a/sql/backends/monet5/bam/Tests/test.stable.out
b/sql/backends/monet5/bam/Tests/test.stable.out
--- a/sql/backends/monet5/bam/Tests/test.stable.out
+++ b/sql/backends/monet5/bam/Tests/test.stable.out
@@ -1,9 +1,9 @@
stdout of test 'test` in directory 'sql/backends/monet5/bam` itself:
-# 21:48:24 >
-# 21:48:24 > "mserver5" "--debug=10" "--set" "gdk_nr_threads=0" "--set"
"mapi_open=true" "--set" "mapi_port=30475" "--set"
"mapi_usock=/var/tmp/mtest-32164/.s.monetdb.30475" "--set" "monet_prompt="
"--forcemito" "--set" "mal_listing=2"
"--dbpath=/export/scratch2/robin/Monet/PREFIX/DFT/var/MonetDB/mTests_sql_backends_monet5_bam"
"--set" "mal_listing=0"
-# 21:48:24 >
+# 23:17:55 >
+# 23:17:55 > "mserver5" "--debug=10" "--set" "gdk_nr_threads=0" "--set"
"mapi_open=true" "--set" "mapi_port=35794" "--set"
"mapi_usock=/var/tmp/mtest-20301/.s.monetdb.35794" "--set" "monet_prompt="
"--forcemito" "--set" "mal_listing=2"
"--dbpath=/export/scratch2/robin/Monet/PREFIX/DFT/var/MonetDB/mTests_sql_backends_monet5_bam"
"--set" "mal_listing=0"
+# 23:17:55 >
# MonetDB 5 server v11.18.0
# This is an unreleased version
@@ -13,8 +13,8 @@ stdout of test 'test` in directory 'sql/
# Copyright (c) 1993-July 2008 CWI.
# Copyright (c) August 2008-2014 MonetDB B.V., all rights reserved
# Visit http://www.monetdb.org/ for further information
-# Listening for connection requests on mapi:monetdb://essen.ins.cwi.nl:30475/
-# Listening for UNIX domain connection requests on
mapi:monetdb:///var/tmp/mtest-32164/.s.monetdb.30475
+# Listening for connection requests on mapi:monetdb://essen.ins.cwi.nl:35794/
+# Listening for UNIX domain connection requests on
mapi:monetdb:///var/tmp/mtest-20301/.s.monetdb.35794
# MonetDB/GIS module loaded
# MonetDB/JAQL module loaded
# MonetDB/SQL module loaded
@@ -61,9 +61,9 @@ Ready.
#WARNING To speedup calc.min a bulk operator implementation is needed
#Interpreting lng as wrd.
-# 21:48:24 >
-# 21:48:24 > "/usr/bin/python2" "test.SQL.py" "test"
-# 21:48:24 >
+# 23:17:55 >
+# 23:17:55 > "/usr/bin/python2" "test.SQL.py" "test"
+# 23:17:55 >
#DECLARE rname_1_3 varchar(32);SET rname_1_3='chr22';DECLARE pos_1_3_1
integer;SET pos_1_3_1=1000000;DECLARE pos_1_3_2 integer;SET
pos_1_3_2=20000000;DECLARE qname_1_4 varchar(32);SET
qname_1_4='sim_22_1_1';DECLARE mapq_1_5 integer;SET mapq_1_5=200;DECLARE
mapq_2_1 integer;SET mapq_2_1=100;DECLARE rname_2_9 varchar(32);SET
rname_2_9='chr22';DECLARE pos_2_9 integer;SET pos_2_9=39996433;DECLARE
rname_2_10 varchar(32);SET rname_2_10='chr22';DECLARE pos_2_10 integer;SET
pos_2_10=80000000;DECLARE distance_2_12 integer;SET
distance_2_12=10000000;SELECT qname, flag, rname, pos, mapq, cigar, rnext,
pnext, tlen, seq, qual
#DECLARE rname_1_3 varchar(32);SET rname_1_3='chr22';DECLARE pos_1_3_1
integer;SET pos_1_3_1=1000000;DECLARE pos_1_3_2 integer;SET
pos_1_3_2=20000000;DECLARE qname_1_4 varchar(32);SET
qname_1_4='sim_22_1_1';DECLARE mapq_1_5 integer;SET mapq_1_5=200;DECLARE
mapq_2_1 integer;SET mapq_2_1=100;DECLARE rname_2_9 varchar(32);SET
rname_2_9='chr22';DECLARE pos_2_9 integer;SET pos_2_9=39996433;DECLARE
rname_2_10 varchar(32);SET rname_2_10='chr22';DECLARE pos_2_10 integer;SET
pos_2_10=80000000;DECLARE distance_2_12 integer;SET
distance_2_12=10000000;SELECT qname, flag, rname, pos, mapq, cigar, rnext,
pnext, tlen, seq, qual
@@ -639,7 +639,7 @@ Ready.
% qname, rname, l_flag, l_pos, l_mapq, l_cigar, l_rnext,
l_pnext, l_tlen, l_seq, l_qual, r_flag, r_pos, r_mapq, r_cigar,
r_rnext, r_pnext, r_tlen, r_seq, r_qual, distance # name
% clob, clob, smallint, int, smallint, clob, clob,
int, int, clob, clob, smallint, int, smallint, clob,
clob, int, int, clob, clob, bigint # type
% 0, 0, 1, 1, 1, 0, 0, 1, 1, 0,
0, 1, 1, 1, 0, 0, 1, 1, 0, 0,
1 # length
-#SELECT * FROM bam.alignments_extra_i;SELECT qname, flag, rname, pos, mapq,
cigar, rnext, pnext, tlen, seq, qual
+#SELECT * FROM bam.alignments_extra_i;INSERT INTO bam.export (SELECT qname,
flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual FROM
bam.alignments_i);CALL
sam_export('/export/scratch2/robin/Monet/SOURCE/MonetDB/sql/backends/monet5/bam/Tests/output_1.sam');SELECT
qname, flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual
% bam.alignments_extra_i, bam.alignments_extra_i, bam.alignments_extra_i,
bam.alignments_extra_i # table_name
% tag, virtual_offset, type, value # name
% char, bigint, char, clob # type
@@ -872,7 +872,9 @@ Ready.
[ "AS", 71, "i", "53" ]
[ "YM", 71, "i", "53" ]
[ "YP", 71, "f", "5.01182e-06" ]
-#SELECT * FROM bam.alignments_extra_i;SELECT qname, flag, rname, pos, mapq,
cigar, rnext, pnext, tlen, seq, qual
+#SELECT * FROM bam.alignments_extra_i;INSERT INTO bam.export (SELECT qname,
flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual FROM
bam.alignments_i);CALL
sam_export('/export/scratch2/robin/Monet/SOURCE/MonetDB/sql/backends/monet5/bam/Tests/output_1.sam');SELECT
qname, flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual
+[ 71 ]
+#SELECT * FROM bam.alignments_extra_i;INSERT INTO bam.export (SELECT qname,
flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual FROM
bam.alignments_i);CALL
sam_export('/export/scratch2/robin/Monet/SOURCE/MonetDB/sql/backends/monet5/bam/Tests/output_1.sam');SELECT
qname, flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual
#FROM bam.alignments_i
#WHERE bam_flag(flag, 'seco_alig') = False
#ORDER BY qname;
@@ -1427,7 +1429,7 @@ Ready.
% qname, rname, l_flag, l_pos, l_mapq, l_cigar, l_rnext,
l_pnext, l_tlen, l_seq, l_qual, r_flag, r_pos, r_mapq, r_cigar,
r_rnext, r_pnext, r_tlen, r_seq, r_qual, distance # name
% clob, clob, smallint, int, smallint, clob, clob,
int, int, clob, clob, smallint, int, smallint, clob,
clob, int, int, clob, clob, bigint # type
% 0, 0, 1, 1, 1, 0, 0, 1, 1, 0,
0, 1, 1, 1, 0, 0, 1, 1, 0, 0,
1 # length
-#SELECT * FROM bam.alignments_extra_i;SELECT qname, flag, rname, pos, mapq,
cigar, rnext, pnext, tlen, seq, qual
+#SELECT * FROM bam.alignments_extra_i;INSERT INTO bam.export (SELECT qname,
flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual FROM
bam.alignments_i);CALL
sam_export('/export/scratch2/robin/Monet/SOURCE/MonetDB/sql/backends/monet5/bam/Tests/output_2.sam');SELECT
qname, flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual
% bam.alignments_extra_i, bam.alignments_extra_i, bam.alignments_extra_i,
bam.alignments_extra_i # table_name
% tag, virtual_offset, type, value # name
% char, bigint, char, clob # type
@@ -1660,13 +1662,14 @@ Ready.
[ "AS", 34556849, "i", "53" ]
[ "YM", 34556849, "i", "53" ]
[ "YP", 34556849, "f", "5.01182e-06" ]
-#SELECT * FROM bam.alignments_extra_i;SELECT qname, flag, rname, pos, mapq,
cigar, rnext, pnext, tlen, seq, qual
+#SELECT * FROM bam.alignments_extra_i;INSERT INTO bam.export (SELECT qname,
flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual FROM
bam.alignments_i);CALL
sam_export('/export/scratch2/robin/Monet/SOURCE/MonetDB/sql/backends/monet5/bam/Tests/output_2.sam');SELECT
qname, flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual
+[ 71 ]
+#SELECT * FROM bam.alignments_extra_i;INSERT INTO bam.export (SELECT qname,
flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual FROM
bam.alignments_i);CALL
sam_export('/export/scratch2/robin/Monet/SOURCE/MonetDB/sql/backends/monet5/bam/Tests/output_2.sam');SELECT
qname, flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual
#FROM bam.unpaired_primary_alignments_i
#UNION
#SELECT qname, flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual
#FROM bam.unpaired_alignments_i
#WHERE bam_flag(flag, 'seco_alig') = False
-#ORDER BY qname;
% .L, .L, .L, .L, .L, .L, .L, .L, .L, .L,
.L # table_name
% qname, flag, rname, pos, mapq, cigar, rnext, pnext, tlen,
seq, qual # name
% clob, smallint, clob, int, smallint, clob, clob,
int, int, clob, clob # type
@@ -2220,7 +2223,7 @@ Ready.
% qname, rname, l_flag, l_pos, l_mapq, l_cigar, l_rnext,
l_pnext, l_tlen, l_seq, l_qual, r_flag, r_pos, r_mapq, r_cigar,
r_rnext, r_pnext, r_tlen, r_seq, r_qual, distance # name
% clob, clob, smallint, int, smallint, clob, clob,
int, int, clob, clob, smallint, int, smallint, clob,
clob, int, int, clob, clob, bigint # type
% 0, 0, 1, 1, 1, 0, 0, 1, 1, 0,
0, 1, 1, 1, 0, 0, 1, 1, 0, 0,
1 # length
-#SELECT * FROM bam.alignments_extra_i;SELECT qname, flag, rname, pos, mapq,
cigar, rnext, pnext, tlen, seq, qual
+#SELECT * FROM bam.alignments_extra_i;INSERT INTO bam.export (SELECT qname,
flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual FROM
bam.unpaired_all_alignments_i);CALL
sam_export('/export/scratch2/robin/Monet/SOURCE/MonetDB/sql/backends/monet5/bam/Tests/output_3.sam');SELECT
qname, flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual
% bam.alignments_extra_i, bam.alignments_extra_i, bam.alignments_extra_i,
bam.alignments_extra_i # table_name
% tag, virtual_offset, type, value # name
% char, bigint, char, clob # type
@@ -2453,13 +2456,14 @@ Ready.
[ "YM", 71, "i", "36" ]
[ "YA", 71, "Z",
"91,36,14=1X76=1056D2=1X2=1D4=;92,37,14=1X77=1058D1I1=1X5=;92,37,14=1X76=1058D1=1I1=1X5=;92,37,14=1X75=1058D2=1I1=1X5=;92,37,14=1X74=1058D3=1I1=1X5=;92,37,14=1X73=1058D4=1I1=1X5=;92,37,14=1X72=1058D5=1I1=1X5=;92,37,14=1X71=1058D6=1I1=1X5=;92,37,14=1X70=1058D7=1I1=1X5=;94,69,14=1X77=1057D1X1=1X5="
]
[ "YP", 71, "f", "1.25893e-09" ]
-#SELECT * FROM bam.alignments_extra_i;SELECT qname, flag, rname, pos, mapq,
cigar, rnext, pnext, tlen, seq, qual
+#SELECT * FROM bam.alignments_extra_i;INSERT INTO bam.export (SELECT qname,
flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual FROM
bam.unpaired_all_alignments_i);CALL
sam_export('/export/scratch2/robin/Monet/SOURCE/MonetDB/sql/backends/monet5/bam/Tests/output_3.sam');SELECT
qname, flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual
+[ 64 ]
+#SELECT * FROM bam.alignments_extra_i;INSERT INTO bam.export (SELECT qname,
flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual FROM
bam.unpaired_all_alignments_i);CALL
sam_export('/export/scratch2/robin/Monet/SOURCE/MonetDB/sql/backends/monet5/bam/Tests/output_3.sam');SELECT
qname, flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual
#FROM bam.unpaired_primary_alignments_i
#UNION
#SELECT qname, flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual
#FROM bam.unpaired_alignments_i
#WHERE bam_flag(flag, 'seco_alig') = False
-#ORDER BY qname;
% .L, .L, .L, .L, .L, .L, .L, .L, .L, .L,
.L # table_name
% qname, flag, rname, pos, mapq, cigar, rnext, pnext, tlen,
seq, qual # name
% clob, smallint, clob, int, smallint, clob, clob,
int, int, clob, clob # type
@@ -3041,7 +3045,7 @@ Ready.
% qname, rname, l_flag, l_pos, l_mapq, l_cigar, l_rnext,
l_pnext, l_tlen, l_seq, l_qual, r_flag, r_pos, r_mapq, r_cigar,
r_rnext, r_pnext, r_tlen, r_seq, r_qual, distance # name
% clob, clob, smallint, int, smallint, clob, clob,
int, int, clob, clob, smallint, int, smallint, clob,
clob, int, int, clob, clob, bigint # type
% 0, 0, 1, 1, 1, 0, 0, 1, 1, 0,
0, 1, 1, 1, 0, 0, 1, 1, 0, 0,
1 # length
-#SELECT * FROM bam.alignments_extra_i;
+#SELECT * FROM bam.alignments_extra_i;INSERT INTO bam.export (SELECT qname,
flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual FROM
bam.unpaired_all_alignments_i);CALL
sam_export('/export/scratch2/robin/Monet/SOURCE/MonetDB/sql/backends/monet5/bam/Tests/output_4.sam');
% bam.alignments_extra_i, bam.alignments_extra_i, bam.alignments_extra_i,
bam.alignments_extra_i # table_name
% tag, virtual_offset, type, value # name
% char, bigint, char, clob # type
@@ -3274,8 +3278,334 @@ Ready.
[ "YM", 34753142, "i", "36" ]
[ "YA", 34753142, "Z",
"91,36,14=1X76=1056D2=1X2=1D4=;92,37,14=1X77=1058D1I1=1X5=;92,37,14=1X76=1058D1=1I1=1X5=;92,37,14=1X75=1058D2=1I1=1X5=;92,37,14=1X74=1058D3=1I1=1X5=;92,37,14=1X73=1058D4=1I1=1X5=;92,37,14=1X72=1058D5=1I1=1X5=;92,37,14=1X71=1058D6=1I1=1X5=;92,37,14=1X70=1058D7=1I1=1X5=;94,69,14=1X77=1057D1X1=1X5="
]
[ "YP", 34753142, "f", "1.25893e-09" ]
+#SELECT * FROM bam.alignments_extra_i;INSERT INTO bam.export (SELECT qname,
flag, rname, pos, mapq, cigar, rnext, pnext, tlen, seq, qual FROM
bam.unpaired_all_alignments_i);CALL
sam_export('/export/scratch2/robin/Monet/SOURCE/MonetDB/sql/backends/monet5/bam/Tests/output_4.sam');
+[ 71 ]
-# 21:48:28 >
-# 21:48:28 > "Done."
-# 21:48:28 >
+Contents of exported SAM file
'/export/scratch2/robin/Monet/SOURCE/MonetDB/sql/backends/monet5/bam/Tests/output_1.sam':
+@HD VN:1.0 SO:coordinate
+@SQ SN:chr13 LN:0
+@SQ SN:chr14 LN:0
+@SQ SN:chr18 LN:0
+@SQ SN:chr22 LN:0
+@SQ SN:chr2 LN:0
+@SQ SN:chr3 LN:0
+@SQ SN:chr9 LN:0
+@CO SAM file generated by MonetDB
+sim_22_1_a 385 chr13 17922987 0 3=1X1=1X43=1X16=1X33=
* 0 0
TATACTTATAGAACAAATGAACCCAAAACCACATAAGGTAAACAACAAAGCTACTGGTTCAAAATTAAGCCTAACTTCAACAGTACCAGGCAAAAACCAT
CGF?BB;DAEFEEEGEFGE@5EABGD/>DGDDB@5EA@B=EDEDGC?CE5B3GBAED5D?GEB>D?-EC@##?:DBDA3#2G:@E#F>A5D####A#A##
+sim_22_1_a 419 chr14 18225784 0 5=1X31=1X62= =
18226024 0
TATACTTATAGAACAAATGAACCCAAAACCACATAAGGTAAACAACAAAGCTACTGGTTCAAAATTAAGCCTAACTTCAACAGTACCAGGCAAAAACCAT
CGF?BB;DAEFEEEGEFGE@5EABGD/>DGDDB@5EA@B=EDEDGC?CE5B3GBAED5D?GEB>D?-EC@##?:DBDA3#2G:@E#F>A5D####A#A##
+sim_22_1_a 339 chr14 18226024 0 100= =
18225784 0
GTACAAAAGTGGGAAATTTTAAAGCTGGGCATCCAGGGGAGACATCACGTGTCAGTAGGTTCCGTGATGCCCCCCAAGCCGCAAAACCAGCAAGTTTTTA
####F#EFCGEC5BE?EF#F=@?DF658;DEDG=@FBG@?E=DE=F#FDB?CCGE1GFFEDEFFDGGFE>FFCGEGGG:;B@FFG@DDGGGFGB7FBC?F
+sim_22_1_a 401 chr18 15278551 0 31=1X1=1X16=1X43=1X5=
* 0 0
ATGGTTTTTGCCTGGTACTGTTGAAGTTAGGCTTAATTTTGAACCAGTAGCTTTGTTGTTTACCTTATGTGGTTTTGGGTTCATTTGTTCTATAAGTATA
##A#A####D5A>F#E@:G2#3ADBD:?##@CE-?D>BEG?D5DEABG3B5EC?CGDEDE=B@AE5@BDDGD>/DGBAE5@EGFEGEEEFEAD;BB?FGC
+sim_22_1_a 99 chr22 15047509 255 100= =
15047749 0
TAAAAACTTGCTGGTTTTGCGGCTTGGGGGGCATCACGGAACCTACTGACACGTGATGTCTCCCCTGGATGCCCAGCTTTAAAATTTCCCACTTTTGTAC
F?CBF7BGFGGGDD@GFF@B;:GGGEGCFF>EFGGDFFEDEFFG1EGCC?BDF#F=ED=E?@GBF@=GDED;856FD?@=F#FE?EB5CEGCFE#F####
+sim_22_1_a 147 chr22 15047749 255 100= =
15047509 0
ATGGTTTTTGCCTGGTACTGTTGAAGTTAGGCTTAATTTTGAACCAGTAGCTTTGTTGTTTACCTTATGTGGTTTTGGGTTCATTTGTTCTATAAGTATA
##A#A####D5A>F#E@:G2#3ADBD:?##@CE-?D>BEG?D5DEABG3B5EC?CGDEDE=B@AE5@BDDGD>/DGBAE5@EGFEGEEEFEAD;BB?FGC
+sim_22_1_9 163 chr22 15600021 255 3=1X96= =
15600266 0
CACATTTTCAAAAAACAAAAAAAAGTCTGAGCTCCTACTGTTGATTTAAATTCTTTTATAAATCTCTATCAAACTTTTCATGTTTACAGTTCTTATGCAA
FGDEGFA:GGG@?BE?F;FGCFF5@EF3;5G#FA.G=D?D?AGAE?AFF=>CBDE32?<E?ED#GG?G?E#E#EEE>5#B####EFF##=D####6#:#,
+sim_22_1_9 385 chr22 15600023 0 1=1I1=1I96= *
0 0
CACATTTTCAAAAAACAAAAAAAAGTCTGAGCTCCTACTGTTGATTTAAATTCTTTTATAAATCTCTATCAAACTTTTCATGTTTACAGTTCTTATGCAA
FGDEGFA:GGG@?BE?F;FGCFF5@EF3;5G#FA.G=D?D?AGAE?AFF=>CBDE32?<E?ED#GG?G?E#E#EEE>5#B####EFF##=D####6#:#,
+sim_22_1_9 83 chr22 15600266 255 100= =
15600021 0
TTTAAAATATTAAAAGATGAATTACTATCAATTGTTTTGAATTTTAAACTAAAAATCAGTAGTTACTATAAAATTATTATTAAATGTTCTAATAATTGTA
5<###E?##AB#E=DB#<==;#E#?=AFEE@C@AEA6GDGFFG?DE??::?=@BC#G>GB>GEEDA=G>ADEFDEGABEDDEBCGGEFFF7;=GEGGFGC
+sim_22_1_6 99 chr22 18665735 255 14=1X10=1X74= =
18665943 0
ATTATTCAGCCTCACATTTCAACATCGCAGAACAGGTCAGGGAGTGCAAGGTACCTCCTGCTCTTTCCAGCAGGTGCAAGGTGCACAGGGAGAAGCAGGT
CDDEFG;?DGEDFE+CECEEGFG>C,EG5?DGE87GBDFGCEDD=BGGDD?GDFGDD?DAGGC<G=GF@A@CE@>?>.6#F#FFDEE?@CE==<9?B?>;
+sim_22_1_6 147 chr22 18665943 255 9=1X90= =
18665735 0
CAGGCCGCCCGGGGCGGTTCCACCAAGAAAGTGACATCTGGGAGGGGCATTGAGGCAAAGGGGCCGCGAGAGGGTAGGTGCTCTAGAGGTGGCCCACGGG
?##??3G#9##=5B#>A#.#>#>CFBEC@:BGC>.A5@6-F4:5?;EGGAA>?6@D<?EAEDE6CGGDGGE@9F?GGGGDD5EDDGFD?:>?FED?DEBA
+sim_22_1_16 99 chr22 20264219 255 84=1X4=1X8=1X1= =
20264427 0
TCCAGCTCCATCTGCCAGGTTTTGTTTTGCCTGAGCAGTAGGAGCCTAGAAATATTCCTTGCACAACACAGGCACATAATGACTCGCTGAGCAGGACTTG
FBGGFDGEFA?GGFGEF;AEGCEEGCFE?FE;GA:EAEEGDD>F8E=E=FABFCDCEG@@@E#BDBCBE;=#D5#,55DCAG:B8EGE#AAGBF#?#E#A
+sim_22_1_16 147 chr22 20264427 255 42=1X57= =
20264219 0
CTCTGTCCGTTTCTAATGGGAGCATTAACTCACTGCTTTGGGGGGCTGAGATAGTCAAAAAAGAGCAATTTCAGTGCAGGTTCTAGTCCTGGGCATGCTC
#;>6##<@#C@#8B#>#:GF2G@A2;BFA?B:-<ECG--)GE7F9DGD;/:ACDABDEBD5EG:EB@DDGB#@:AEA?=GG=DE?D5GD<AFEC;DF=GA
+sim_22_1_f 99 chr22 23113643 255 100= =
23113851 0
CTTTTCTATATAAAACAAATAATTTATAATACTGTATTTGTTTTCTATTACTATTGAAATTACTACAAATGTAGTAGCTAAAACAACACATATTTTTCAT
DE>FDG5FDGDGFDAGEFABDGGEGADG?BEDDFGGEDGA?GG=9FEDDGBD58=@AFDFD=CBD9GGD=DCC?EE?F#A6GC#>?F#F@2E?CDF?#B?
+sim_22_1_f 147 chr22 23113851 255 14=1X85= =
23113643 0
TGCCTTTTCCAGCTGATAGAGACAGCCATGTTTCATGGCTTGTAGCCCCTTCCTCCATCTTCAAAGCCAACAACAACTGGTCAAGTTGTCACATTGCATC
####FF-0AD#E>@#?#9??E##A>C=C>D7A,<G?D>=>DDD=F:?G=>CE?:CAFGEDGEFDAGGFEF7CGDCEDEEGA@DECFGE-?DEACD?DC=G
+sim_22_1_f 401 chr22 23113851 0 14=1X77=1056D1=1X2=1D4=
* 0 0
TGCCTTTTCCAGCTGATAGAGACAGCCATGTTTCATGGCTTGTAGCCCCTTCCTCCATCTTCAAAGCCAACAACAACTGGTCAAGTTGTCACATTGCATC
####FF-0AD#E>@#?#9??E##A>C=C>D7A,<G?D>=>DDD=F:?G=>CE?:CAFGEDGEFDAGGFEF7CGDCEDEEGA@DECFGE-?DEACD?DC=G
+sim_22_1_d 321 chr22 23964184 0 33=1X52=1X1=1X11=
* 0 0
CCTGTCAAGCCTACCTTGTCACCTGTTGCCAGGAAGTCCTGCTCCCTCCTAGGGCAGCCCCTGGTTAGTGGAGCTGAGAGAACCCAGGGTGCTGAGCCGT
G5FGBEFEEGFED=EGDABFGGGFGGCE?EF=EEDGG9G?GEAECD9FGBFCEF;GGEC=GC5>F1?CDCCD#FEB=EE5B#CFEB###:ECBA-CE:@A
+sim_22_1_d 99 chr22 24189222 255 86=1X1=1X11= =
24189415 0
CCTGTCAAGCCTACCTTGTCACCTGTTGCCAGGAAGTCCTGCTCCCTCCTAGGGCAGCCCCTGGTTAGTGGAGCTGAGAGAACCCAGGGTGCTGAGCCGT
G5FGBEFEEGFED=EGDABFGGGFGGCE?EF=EEDGG9G?GEAECD9FGBFCEF;GGEC=GC5>F1?CDCCD#FEB=EE5B#CFEB###:ECBA-CE:@A
+sim_22_1_d 147 chr22 24189415 255 100= =
24189222 0
CGGGGGGCGGAGCGCAGGCAGGTGAAACTGGAGGAAATTGGAGATACAGCGGGAGGGAGGGGGAAAAGGTAAAACCTCTCACCACCTCCTCAACGCTTAT
##??BG#C##5AB#D;=G#BF#E#:;@?#<(4BGCG-#EA@<G?FC*GDD@FD:BEBEGFGFDDBDEAD9GGD37DDFDCDFDGEFGFCGEE:?AE>EDG
+sim_22_1_b 99 chr22 24215590 255 23=1X76= =
24215802 0
ACAGTCACCCACAGATTTGCACACCTGTCCAACTTGAGCCAGTGATGAATAAAGACAGATATGGTTGCTGTTCTTTCATGGAGCTTACATCCTAGGGGAG
GGDEGGEAGGD?ED?>GDGDGFG.-GEFG:FE=FFGGFDGDDE==C==:@C7DEE=#GBC,BFFC?<@:A#:6AC>D?C#;EFD/FE#@CF#C=#B?#6=
+sim_22_1_b 147 chr22 24215802 255 20=1X79= =
24215590 0
AGTGGGCAGGGCAGCTTTCTACGAGGAGGTGACTAAATTAAGATCTTGCTGTGAAGATGTTTAGTTGCCTGCTACCCTCCACTTCTCCCTGAATAAAACA
A=#=C?A5.E##?#A##=#C#D?-A#E+E?CEFFCCB1;=-DG:DGFCECF=GA@FFFAECACF-;=DDDFDFB??F?F5?GGF:;FEBGFC;(CGDGG:
+sim_22_1_10 99 chr22 24428205 255 28=1X54=1X16= =
24428413 0
TAATAATCTCAAAAGATAGTTCAATGTAGAACCTAATTAATTCAGTAGGATGACCAGCTCAGTAGCAGTAAGTAGCAGTGCAATATAAGTGGCATTTTAA
EFE:FC@D?5FGGCAEDAFEGG3GEDA5(CFFFEGD#GEFFG?EDCFC?GBAGEFADCFBGCG@FGCF2F0ED>C@E:DAF:A-BF;:/#;#E##A#A?C
+sim_22_1_10 147 chr22 24428413 255 66=1X33= =
24428205 0
ATGCCTTATAAATAAATTTATAAGATGACTACTTTGTAGGCACATCTTGTTATAAGAAAAGTCAATGATATGTGTCCCTTTAATTTTCAAAATGAAAAAT
#GD###A#5ED#B5.#BE==5FG#E3-@#AB7CA)@DEFF=FF#@AEFDFD4>DF2CCCGGA:DGB'FFDG9AEF?3;DGAE?CGCFD2GAG?GB@GFGF
+sim_22_1_3 99 chr22 25791286 255 40=2X31=1X26= =
25791479 0
ATTCTCTGCCCTGAGTAAGGGGAATCAAGAGAGATTTCGATGCAGGAAACAGCACCCGAGAGAATCAATCCAGTTGGGCCTTTTCAATTAGTCTGTTTTT
?GFGFFAEEEGEF-GGEEGGEEFFGDF;??DGDGD/=DBD09EGG5EDG@DEC,FBGD@BF5:FBE6:[email protected]>EF#=#:EE#A>###
+sim_22_1_3 147 chr22 25791479 255
1=1X12=2X12=1X17=1X4=1X48= = 25791286 0
CGTCTCTCCTGCCCCCCACCCTCTGTCTGTGAAAACATTAGAGCTGAAAGACGACTCCCCTGGCGATAAACATCTACGCAAAGCCAATTGTCAGGCGATG
4####B#?C#???#6,D####-BAEF3<)DE>AD(#DD+BFEFFBBEAFBC1GABFAA2GEBGEDF8G7G=AAAD=EGDD:A4ABE?AFGDGD0DD?DB>
+sim_22_1_15 385 chr22 25992044 0 2=1X95=2I *
0 0
GAGAGGCAGCCATAGGTTCGAGTCCTAGCTACTGGCTAGGGGACCTTGGACAAGGCACTGATAGTCCCAAGACTCAGTTTACCCATCTGTGCAGTGATGT
CG.GCEGGGFDGGGGEDCE<BEG5GDGEE=##EEFEFABFDE:#G??>:55;FE:5=A?E;E9BF#E@5=AE;BEGCA:#>EAC#DED#B##A;###B#B
+sim_22_1_15 163 chr22 25992044 255 2=1X95=1X1= =
25992258 0
GAGAGGCAGCCATAGGTTCGAGTCCTAGCTACTGGCTAGGGGACCTTGGACAAGGCACTGATAGTCCCAAGACTCAGTTTACCCATCTGTGCAGTGATGT
CG.GCEGGGFDGGGGEDCE<BEG5GDGEE=##EEFEFABFDE:#G??>:55;FE:5=A?E;E9BF#E@5=AE;BEGCA:#>EAC#DED#B##A;###B#B
+sim_22_1_15 83 chr22 25992258 255 100= =
25992044 0
ATGTGCTCTATGGTCACCTAACCTGTCTGAGACTCAATATACTCATCTGATAAATAGGCAGAATAAAAATAAGAGGTGCCCTGCCTCTTTTCAGAGTCAC
##AEG#=GE>75CADDFEB=DF;<7=CF?<B?CBBC6:GDCDE=EBFFGFDE=BDDEEDEB?DD'AAE?CEEDFGFDDG>EGEGDA>FGD;EEEF:@GF<
+sim_22_1_5 99 chr22 28913667 255 44=1X26=1X28= =
28913875 0
TGGAATGGAGTGGGGTAGAGGACTGCTGATTTAGCAGATGCATTATTTCCCCAAGCCAAAGAGCTGCCCTGGAAATGGGCCAAATATAATTTCTTTTAAG
FGECEGDGGGGG?DGGEBG=FGGFD@EA9:EGGDGCFDDCEC??=DF:DECAEFC:A@EFEGE6>#GBE5C#BF<EDC=D&?G;E#=#B###G=A#@?5#
+sim_22_1_5 147 chr22 28913875 255 100= =
28913667 0
AGGGGCAGTGGGGGTGGTGTGGGGGAGCCGCTGTGAGTCTTAATGCAAAGAGGAATGTAGGCTAGTCAACAGGAAACCCAAAGTGCTTGATACAGCACAG
##DAA#C-#BE#D##F#9:B#@AB#DF#6D6G8F-D8CFF8E?E:#=>B>AG=6EG#CGEDG9B3:AEEG6AC:BE5DDCDFBBG5GCFG4GAEEFBGG?
+sim_22_1_2 99 chr22 29469417 255 100= =
29469608 0
TGATGCTGAGCTTGACGGCGGTCCCCCGGGTACCACTCCATCACCCAGTGAAGGAGGCAGGTGGGACCATCTCAGACATGGGCTGTGTTGGGACCCCCAG
@*G@F:7DEDG?BFGFFGA:DBGG:EAGDCDAGG?E=FGGGGD?EGG5FEFE?#EA@E:DCD?>E-?C?BE<BE?CE:GB=CGC#+?#CC###>C#D@F#
+sim_22_1_2 147 chr22 29469608 255 11=2X13=1X73= =
29469417 0
TGGCCAAGTAATCATGCTAACTAGAGCTCTAGCTTCCTCTCCCTGGTTTCGCTTTTCTAAGTATAAATAAATAAATAAATGATAAATTGCCTTATTTTTT
?#D#A##FBE@-#0EA/?A>#?=EB9%;<-EC>G3EA2ADDF9EGCACFAFC?:GAB5AGE=CGBE5?DEEEFBG:DFDECFDB5AGEGGGGGFD?G?DD
+sim_22_1_8 163 chr22 29509837 255 48=1X51= =
29510064 0
TCCTAGGGCAATGCCTAGTTAGGTAAACCAACTGGTGCTGAAGAGGAGCGCTTGGGTCTGGGAAATGCCCATTGTAATCACTTGGGACCGTGTTACTTGA
F?DE>FGCFDDGGEFAD@GFGEGEGF?DABCF@>@E?<D55ABGEAE#(-#EDE?=D8&@DGEA@#8DDD5DDF#F##E4=@=<BD-#AA?##<A###.#
+sim_22_1_8 83 chr22 29510064 255 66=1X33= =
29509837 0
AGGCCTCTGTTTTTTAAAATACTGCAGATACAAATCCTAGATTTTCATTCTTAGCACTAGAATGGCTTCTCCAGAAACAGATTTATGTATTTGTAGGAAT
>BA#EAB#A?#=DCE#DDEACCDGA?@B?D;AC5EAFA=D)#EG8EG##DGD=?#FFGAABECGFA,G<EGEE:@FDEFE?DGDGGFGFGGGGEGGEFGA
+sim_22_1_13 163 chr22 32901856 255 100= =
32902057 0
TTATTACTCATAACAATCTTATGGGGTACATACTCTCATTATCCCGATTTGACAGATTAGTAAACTGAGGTGCAGACAGGCAAGTGATTCTTCCAGAGAG
F:CADE=EBGDDAEGE>DE@@D=CDDE#EG35CG?E;7CEEBDGGEB@-F=GEGD8CG>=5A=BBAFB#?D:G?E3A##?#F-?#DC#:###E#G8##B#
+sim_22_1_13 337 chr22 32902044 0 1=1I1=1I10=15D86=
* 0 0
AGAGCAAAACATTCTGGCAGATGCACAGCAGGTGGTGGGGCAAAAAGGGACCCGGATGGGGATCAAGAGGCTCACAAGATAGTCTCAGCTCTGCCGCTTA
;5#GFC#@@#AE@FE@E@4B=BE:6;ADA>E#:GF=?#=:?/EA?G>GF3CFCDEB:EFGDG?G@GGFEG=F5FDGE?#GDFGGFEE@AGADGFGEDEEG
+sim_22_1_13 83 chr22 32902057 255 100= =
32901856 0
AGAGCAAAACATTCTGGCAGATGCACAGCAGGTGGTGGGGCAAAAAGGGACCCGGATGGGGATCAAGAGGCTCACAAGATAGTCTCAGCTCTGCCGCTTA
;5#GFC#@@#AE@FE@E@4B=BE:6;ADA>E#:GF=?#=:?/EA?G>GF3CFCDEB:EFGDG?G@GGFEG=F5FDGE?#GDFGGFEE@AGADGFGEDEEG
+sim_22_1_14 99 chr22 37296079 255 100= =
37296332 0
TCGGGCCTAGAACACGGAGCCTGGAGACCCGCGCGCCGTTGACCACAGGGTTTGGCGCCAAACGAGGAGGGCGGGGCGAGGACCAAGGGGCCGGGGCGGG
FDFGECDG?GDF=2BGEGEGEE@GGG=EEG:GFF@>EADDAFEE?::FGEGB>F>E@6EGGAE#DG#=B?(#@BCDDAEEGGFBAFFE-D#:>#:#CD7#
+sim_22_1_14 355 chr22 37296079 0 92=172D6=1X1= =
37296228 0
TCGGGCCTAGAACACGGAGCCTGGAGACCCGCGCGCCGTTGACCACAGGGTTTGGCGCCAAACGAGGAGGGCGGGGCGAGGACCAAGGGGCCGGGGCGGG
FDFGECDG?GDF=2BGEGEGEE@GGG=EEG:GFF@>EADDAFEE?::FGEGB>F>E@6EGGAE#DG#=B?(#@BCDDAEEGGFBAFFE-D#:>#:#CD7#
+sim_22_1_14 403 chr22 37296228 0
1=1X7=104D18=1X9=1X12=1X49= = 37296079 0
GGCGAGGGGGCGGGGCGAGGGCTGAGGGGCGGAAAACTAAGAGGAGCCCCCGCGGGTGGCGAGCTGGGATGACGGAGGGAGGGAGGGAGGGGCGCTGTGG
E##G#A=CD:.E-E.DEDE71#-5EEAGC;5-=#EE5'?DDED?:CFDFG&=?ACFB=DFGG=EBFDGCAGFGG>FG>GFGFFDFGE=AG6DAGEFFEGA
+sim_22_1_14 147 chr22 37296332 255 27=1X9=1X12=1X49=
= 37296079 0
GGCGAGGGGGCGGGGCGAGGGCTGAGGGGCGGAAAACTAAGAGGAGCCCCCGCGGGTGGCGAGCTGGGATGACGGAGGGAGGGAGGGAGGGGCGCTGTGG
E##G#A=CD:.E-E.DEDE71#-5EEAGC;5-=#EE5'?DDED?:CFDFG&=?ACFB=DFGG=EBFDGCAGFGG>FG>GFGFFDFGE=AG6DAGEFFEGA
+sim_22_1_e 163 chr22 39608988 255 92=1X7= =
39609212 0
CGTAGTGAGACCCCATCTCCAAAAAATAAATAAATTAAGGACAGGATGGGCATGGACCACTCTTCAAGATTACTGACCAGCTTCTTGATATTGGACAACT
BC-GDFADDFFFCG@FAF6GCFBD==@=ACG>=G>EC5<>B#E0<CG4,=GBB>DDGAGEDE7BBD6:G#55DD=EEF#B5#;#GEC#BE####A=C?A5
+sim_22_1_e 337 chr22 39609186 0 1=1I1=2X9=27D86=
* 0 0
AATATAAATTAAAATCTCAGGCCTCCATCCCCTATCTACTGAATCAGAATCTGCATTTTAACAGATCTTTTGTGGTTCATATACAGTTGAGTTCTTTTTT
4D@6##D=>G#<E#F#?G#;BC:EE#DB:C#>DE>B>:#ED?CFAFDF>DDC<CG:FEGGF5GEE1EFBFDGDG=GDBFDEDFBEGGGFDEDEDFGC?FF
+sim_22_1_e 83 chr22 39609212 255 4=1X95= =
39608988 0
AATATAAATTAAAATCTCAGGCCTCCATCCCCTATCTACTGAATCAGAATCTGCATTTTAACAGATCTTTTGTGGTTCATATACAGTTGAGTTCTTTTTT
4D@6##D=>G#<E#F#?G#;BC:EE#DB:C#>DE>B>:#ED?CFAFDF>DDC<CG:FEGGF5GEE1EFBFDGDG=GDBFDEDFBEGGGFDEDEDFGC?FF
+sim_22_1_1 419 chr22 39994332 0 5=1X11=1885D83= =
39995061 0
AGTAGATGGGATTACAGGCACCCGCCACCGCGCCCAGCTAATTTTTATATATTTTTTAAGTAGAGATGGGGTTTCACCATGTTGGCCAGGCTGGTCTTGA
EEEEEFAG=GCDAGBF8F:EFGGDEGCFFBAAGCE?DD?5?DGF#DC##BD;D=BD@G?A0F9G6AFC?5#CCCAE6AD-F#CA0#.ADCD;#-#6#A##
+sim_22_1_1 385 chr22 39994467 0 3=1I1=1X14=1751D80=
* 0 0
AGTAGATGGGATTACAGGCACCCGCCACCGCGCCCAGCTAATTTTTATATATTTTTTAAGTAGAGATGGGGTTTCACCATGTTGGCCAGGCTGGTCTTGA
EEEEEFAG=GCDAGBF8F:EFGGDEGCFFBAAGCE?DD?5?DGF#DC##BD;D=BD@G?A0F9G6AFC?5#CCCAE6AD-F#CA0#.ADCD;#-#6#A##
+sim_22_1_1 339 chr22 39995061 0 11=1368D51=1X37=
= 39994332 0
ACTCTGTCGCCCAGGCCGGGGCGTAGCCATGCAAACACGGCTCACTGCAGCCTCGACCTCCCCGGCTCAAGCCATCCTCCCACCTCAACCTCCCAAGTAG
##?#-@=#?D5?EBC##=2#E6A=E8A:F5+F>=CF8#BDFDF?7B8BAGC2EF?FFFG>EE*F=GDBGDGFEFEFG=FBGFFBFFBGEDCF5EFGGGDG
+sim_22_1_1 419 chr22 39995910 0 5=1X12=307D82= =
39996123 0
AGTAGATGGGATTACAGGCACCCGCCACCGCGCCCAGCTAATTTTTATATATTTTTTAAGTAGAGATGGGGTTTCACCATGTTGGCCAGGCTGGTCTTGA
EEEEEFAG=GCDAGBF8F:EFGGDEGCFFBAAGCE?DD?5?DGF#DC##BD;D=BD@G?A0F9G6AFC?5#CCCAE6AD-F#CA0#.ADCD;#-#6#A##
+sim_22_1_1 385 chr22 39996047 0 3=1I1=1X13=171D81=
* 0 0
AGTAGATGGGATTACAGGCACCCGCCACCGCGCCCAGCTAATTTTTATATATTTTTTAAGTAGAGATGGGGTTTCACCATGTTGGCCAGGCTGGTCTTGA
EEEEEFAG=GCDAGBF8F:EFGGDEGCFFBAAGCE?DD?5?DGF#DC##BD;D=BD@G?A0F9G6AFC?5#CCCAE6AD-F#CA0#.ADCD;#-#6#A##
+sim_22_1_1 339 chr22 39996123 0 4=1X1I10=307D46=1X37=
= 39995910 0
ACTCTGTCGCCCAGGCCGGGGCGTAGCCATGCAAACACGGCTCACTGCAGCCTCGACCTCCCCGGCTCAAGCCATCCTCCCACCTCAACCTCCCAAGTAG
##?#-@=#?D5?EBC##=2#E6A=E8A:F5+F>=CF8#BDFDF?7B8BAGC2EF?FFFG>EE*F=GDBGDGFEFEFG=FBGFFBFFBGEDCF5EFGGGDG
+sim_22_1_1 163 chr22 39996217 255 100= =
39996429 0
AGTAGATGGGATTACAGGCACCCGCCACCGCGCCCAGCTAATTTTTATATATTTTTTAAGTAGAGATGGGGTTTCACCATGTTGGCCAGGCTGGTCTTGA
EEEEEFAG=GCDAGBF8F:EFGGDEGCFFBAAGCE?DD?5?DGF#DC##BD;D=BD@G?A0F9G6AFC?5#CCCAE6AD-F#CA0#.ADCD;#-#6#A##
+sim_22_1_1 83 chr22 39996429 255 62=1X37= =
39996217 0
ACTCTGTCGCCCAGGCCGGGGCGTAGCCATGCAAACACGGCTCACTGCAGCCTCGACCTCCCCGGCTCAAGCCATCCTCCCACCTCAACCTCCCAAGTAG
##?#-@=#?D5?EBC##=2#E6A=E8A:F5+F>=CF8#BDFDF?7B8BAGC2EF?FFFG>EE*F=GDBGDGFEFEFG=FBGFFBFFBGEDCF5EFGGGDG
+sim_22_1_11 163 chr22 41987203 255 68=1X31= =
41987433 0
GCCTGGATGGAGCCCTGGGGCTGGAGGCTGGCACTTCCAATCGAGACTGCCCGTTTGGTCCTGCCCTGTGTGACAGTATACAGGCTTATCTTGGCCACCC
GGGF:GFGG?EGFEG?DAFBEF=FA+>EGG=E:EDBEB7DCCE=DG@?:?CAE4E=DCBGFG@FE8EB>FDAE#CE=ECA@#D-D?:CEA#BD#E#B#EA
+sim_22_1_11 337 chr22 41987433 0 19=1X50=1X26=1I1=1I
* 0 0
TGACTATTCTAGGAGACCACGAGTTGGGAAAGGTATGGCAAGGCACCCACAGGAGAGCTGATGACAACTGTGAGGATCAGTGACTGTGGATGATGGAACA
=?GF3A:EA#GF6#D5=EF##AG>*EE==A#GB=GA.DE'B@AACDECGB3G>DGFDCGFBFE?EFFD;E7GBFGDCAEGBGDAEFB;5GD=GEECF0FE
+sim_22_1_11 83 chr22 41987433 255 19=1X50=1X26=1X2=
= 41987203 0
TGACTATTCTAGGAGACCACGAGTTGGGAAAGGTATGGCAAGGCACCCACAGGAGAGCTGATGACAACTGTGAGGATCAGTGACTGTGGATGATGGAACA
=?GF3A:EA#GF6#D5=EF##AG>*EE==A#GB=GA.DE'B@AACDECGB3G>DGFDCGFBFE?EFFD;E7GBFGDCAEGBGDAEFB;5GD=GEECF0FE
+sim_22_1_17 337 chr22 42395458 0 2=1D9=617D9=1X79=
* 0 0
TAATTTTACCTCTCAGAAAGCTCTGGCTGCTGCTGTTCATCTTTTTGTAAATGTTCACTTAACAGCTGATCCTTTGGAACAAGAACTGGAGAGACAGTGG
C=#EDEC#5####E#D-E3E#ECGA?EE=GCDD:#@E:B?F?E@GDBDFFG5GC:FE=BEBF:GA?EDG@B<>A?G:FFEFCFECFDCDC<DGEDDF;EG
+sim_22_1_17 385 chr22 42395633 0 5=1I3=232D71=1X12=1X6=
* 0 0
TCAGTGTTATCATTTCTAAATTTAAAAGGTAAATTCAATAATTCCTAAGATCTCTGCTTGCTCTAAAATTCCCTGACTCTCATTACTGTCTTTATAGTGT
G>AEE=EF-GB?EGG@GGEEAB-GCE>E9;DEACG?:CCGDA(E7DFF.GFGEBEEED;D#FG4G:DDGGC=GG.;D@#A#B????#C#B#CB##AFE##
+sim_22_1_17 163 chr22 42395864 255 8=1X71=1X12=1X6=
= 42396076 0
TCAGTGTTATCATTTCTAAATTTAAAAGGTAAATTCAATAATTCCTAAGATCTCTGCTTGCTCTAAAATTCCCTGACTCTCATTACTGTCTTTATAGTGT
G>AEE=EF-GB?EGG@GGEEAB-GCE>E9;DEACG?:CCGDA(E7DFF.GFGEBEEED;D#FG4G:DDGGC=GG.;D@#A#B????#C#B#CB##AFE##
+sim_22_1_17 83 chr22 42396076 255 20=1X79= =
42395864 0
TAATTTTACCTCTCAGAAAGCTCTGGCTGCTGCTGTTCATCTTTTTGTAAATGTTCACTTAACAGCTGATCCTTTGGAACAAGAACTGGAGAGACAGTGG
C=#EDEC#5####E#D-E3E#ECGA?EE=GCDD:#@E:B?F?E@GDBDFFG5GC:FE=BEBF:GA?EDG@B<>A?G:FFEFCFECFDCDC<DGEDDF;EG
+sim_22_1_7 337 chr22 45059241 0 1=1X8=1134D31=1X56=1I1=
* 0 0
ACCAGCCTGGCCAACATGGTGAAACGCTGTCTCTACTAAAAGTACAAAAAAATTAGCTAGGCGTGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGGAG
5#>?CG=E,A=#AG#?#>GAEA@=DB<>ED@?CFC@A?=EG65GGD:C@EE(BFABFG5B?G;EFE?@5EG5GAGD5@;DFGG?F?7GAGGAFEFA>E-E
+sim_22_1_7 163 chr22 45060162 255 89=1X5=1X4= =
45060375 0
TAGGAATTTTTTTTTTTTTTGAGACAGAGTTTCGCTCTTGTTGCCCAGGCTGGAGTGCAATGGTGCGATTGGCTGGAATTTTTTTTCTTGGCTGTCATAA
FGFCEEEGDE?CB::GE=FADFCEABEGBF>D;FGGGFFCDFGDGABCCADCDDDC;EF=:ACGGCEC#C=?=>EG@9#A@?=E;=:#B##A##@%C#D#
+sim_22_1_7 419 chr22 45060235 0 1=1I12=72I3=1X5=1X4=
= 45060375 0
TAGGAATTTTTTTTTTTTTTGAGACAGAGTTTCGCTCTTGTTGCCCAGGCTGGAGTGCAATGGTGCGATTGGCTGGAATTTTTTTTCTTGGCTGTCATAA
FGFCEEEGDE?CB::GE=FADFCEABEGBF>D;FGGGFFCDFGDGABCCADCDDDC;EF=:ACGGCEC#C=?=>EG@9#A@?=E;=:#B##A##@%C#D#
+sim_22_1_7 339 chr22 45060375 0 41=1X56=1I1= =
45060235 0
ACCAGCCTGGCCAACATGGTGAAACGCTGTCTCTACTAAAAGTACAAAAAAATTAGCTAGGCGTGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGGAG
5#>?CG=E,A=#AG#?#>GAEA@=DB<>ED@?CFC@A?=EG65GGD:C@EE(BFABFG5B?G;EFE?@5EG5GAGD5@;DFGG?F?7GAGGAFEFA>E-E
+sim_22_1_7 83 chr22 45060375 255 41=1X56=1X1= =
45060162 0
ACCAGCCTGGCCAACATGGTGAAACGCTGTCTCTACTAAAAGTACAAAAAAATTAGCTAGGCGTGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGGAG
5#>?CG=E,A=#AG#?#>GAEA@=DB<>ED@?CFC@A?=EG65GGD:C@EE(BFABFG5B?G;EFE?@5EG5GAGD5@;DFGG?F?7GAGGAFEFA>E-E
+sim_22_1_4 99 chr22 46558150 255 100= =
46558385 0
CCCTGGCTGGAGTACAAGTTACTGGGGCTCAGTTCTTAGGATTCCTACAAAGCCCAATTTGGGTCACCTGACTTGTTAAGTCAGGTGAAGGTGACTTACC
G6EG5GGGGDGGADFGEB?GFDGGAEGGCFFBECGGB*BEA9BFG6E=?=)C5FGAFE6G2FB=CFECF#DAEAG<AGD7F*44D:#B>+?F##;E;#DD
+sim_22_1_4 147 chr22 46558385 255 1X49=1X9=1X39= =
46558150 0
TGACTCCAGGGAAGTTAAATAAATGCCAAGAGGATTCACGTGGCAGGGCCATGGAGGTGGGGAGGAAGGAGCACAGGAGACACCAGGTGCTCACTGTCCA
#E###BF#?;;#>@=1#@=E=B#E-EDF?DBC#EDEF?;6FE>[email protected]+DAF:4E>FF=EDA;DAEGEGAEFD?GD=DG=G=C5FEGDGGFAG
+sim_22_1_4 401 chr22 46558386 0 1I49=1X9=1X39= *
0 0
TGACTCCAGGGAAGTTAAATAAATGCCAAGAGGATTCACGTGGCAGGGCCATGGAGGTGGGGAGGAAGGAGCACAGGAGACACCAGGTGCTCACTGTCCA
#E###BF#?;;#>@=1#@=E=B#E-EDF?DBC#EDEF?;6FE>[email protected]+DAF:4E>FF=EDA;DAEGEGAEFD?GD=DG=G=C5FEGDGGFAG
+sim_22_1_a 401 chr2 132593958 0 33=1X16=1X3=1X39=1X5=
* 0 0
ATGGTTTTTGCCTGGTACTGTTGAAGTTAGGCTTAATTTTGAACCAGTAGCTTTGTTGTTTACCTTATGTGGTTTTGGGTTCATTTGTTCTATAAGTATA
##A#A####D5A>F#E@:G2#3ADBD:?##@CE-?D>BEG?D5DEABG3B5EC?CGDEDE=B@AE5@BDDGD>/DGBAE5@EGFEGEEEFEAD;BB?FGC
+sim_22_1_9 321 chr3 76011277 0 41=1X34=1X2=1X20=
* 0 0
TACAATTATTAGAACATTTAATAATAATTTTATAGTAACTACTGATTTTTAGTTTAAAATTCAAAACAATTGATAGTAATTCATCTTTTAATATTTTAAA
CGFGGEG=;7FFFEGGCBEDDEBAGEDFEDA>G=ADEEG>BG>G#CB@=?::??ED?GFFGDG6AEA@C@EEFA=?#E#;==<#BD=E#BA##?E###<5
+sim_22_1_a 401 chr9 66420362 0 50=1X43=1X5= *
0 0
ATGGTTTTTGCCTGGTACTGTTGAAGTTAGGCTTAATTTTGAACCAGTAGCTTTGTTGTTTACCTTATGTGGTTTTGGGTTCATTTGTTCTATAAGTATA
##A#A####D5A>F#E@:G2#3ADBD:?##@CE-?D>BEG?D5DEABG3B5EC?CGDEDE=B@AE5@BDDGD>/DGBAE5@EGFEGEEEFEAD;BB?FGC
+sim_22_1_a 401 chr9 69153222 0 50=1X43=1X5= *
0 0
ATGGTTTTTGCCTGGTACTGTTGAAGTTAGGCTTAATTTTGAACCAGTAGCTTTGTTGTTTACCTTATGTGGTTTTGGGTTCATTTGTTCTATAAGTATA
##A#A####D5A>F#E@:G2#3ADBD:?##@CE-?D>BEG?D5DEABG3B5EC?CGDEDE=B@AE5@BDDGD>/DGBAE5@EGFEGEEEFEAD;BB?FGC
+
+
+Contents of exported SAM file
'/export/scratch2/robin/Monet/SOURCE/MonetDB/sql/backends/monet5/bam/Tests/output_2.sam':
+@HD VN:1.0 SO:coordinate
+@SQ SN:chr13 LN:0
+@SQ SN:chr14 LN:0
+@SQ SN:chr18 LN:0
+@SQ SN:chr22 LN:0
+@SQ SN:chr2 LN:0
+@SQ SN:chr3 LN:0
+@SQ SN:chr9 LN:0
+@CO SAM file generated by MonetDB
+sim_22_1_a 385 chr13 17922987 0 3=1X1=1X43=1X16=1X33=
* 0 0
TATACTTATAGAACAAATGAACCCAAAACCACATAAGGTAAACAACAAAGCTACTGGTTCAAAATTAAGCCTAACTTCAACAGTACCAGGCAAAAACCAT
CGF?BB;DAEFEEEGEFGE@5EABGD/>DGDDB@5EA@B=EDEDGC?CE5B3GBAED5D?GEB>D?-EC@##?:DBDA3#2G:@E#F>A5D####A#A##
+sim_22_1_a 419 chr14 18225784 0 5=1X31=1X62= =
18226024 0
TATACTTATAGAACAAATGAACCCAAAACCACATAAGGTAAACAACAAAGCTACTGGTTCAAAATTAAGCCTAACTTCAACAGTACCAGGCAAAAACCAT
CGF?BB;DAEFEEEGEFGE@5EABGD/>DGDDB@5EA@B=EDEDGC?CE5B3GBAED5D?GEB>D?-EC@##?:DBDA3#2G:@E#F>A5D####A#A##
+sim_22_1_a 339 chr14 18226024 0 100= =
18225784 0
GTACAAAAGTGGGAAATTTTAAAGCTGGGCATCCAGGGGAGACATCACGTGTCAGTAGGTTCCGTGATGCCCCCCAAGCCGCAAAACCAGCAAGTTTTTA
####F#EFCGEC5BE?EF#F=@?DF658;DEDG=@FBG@?E=DE=F#FDB?CCGE1GFFEDEFFDGGFE>FFCGEGGG:;B@FFG@DDGGGFGB7FBC?F
+sim_22_1_a 401 chr18 15278551 0 31=1X1=1X16=1X43=1X5=
* 0 0
ATGGTTTTTGCCTGGTACTGTTGAAGTTAGGCTTAATTTTGAACCAGTAGCTTTGTTGTTTACCTTATGTGGTTTTGGGTTCATTTGTTCTATAAGTATA
##A#A####D5A>F#E@:G2#3ADBD:?##@CE-?D>BEG?D5DEABG3B5EC?CGDEDE=B@AE5@BDDGD>/DGBAE5@EGFEGEEEFEAD;BB?FGC
+sim_22_1_a 99 chr22 15047509 255 100= =
15047749 0
TAAAAACTTGCTGGTTTTGCGGCTTGGGGGGCATCACGGAACCTACTGACACGTGATGTCTCCCCTGGATGCCCAGCTTTAAAATTTCCCACTTTTGTAC
F?CBF7BGFGGGDD@GFF@B;:GGGEGCFF>EFGGDFFEDEFFG1EGCC?BDF#F=ED=E?@GBF@=GDED;856FD?@=F#FE?EB5CEGCFE#F####
+sim_22_1_a 147 chr22 15047749 255 100= =
15047509 0
ATGGTTTTTGCCTGGTACTGTTGAAGTTAGGCTTAATTTTGAACCAGTAGCTTTGTTGTTTACCTTATGTGGTTTTGGGTTCATTTGTTCTATAAGTATA
##A#A####D5A>F#E@:G2#3ADBD:?##@CE-?D>BEG?D5DEABG3B5EC?CGDEDE=B@AE5@BDDGD>/DGBAE5@EGFEGEEEFEAD;BB?FGC
+sim_22_1_9 163 chr22 15600021 255 3=1X96= =
15600266 0
CACATTTTCAAAAAACAAAAAAAAGTCTGAGCTCCTACTGTTGATTTAAATTCTTTTATAAATCTCTATCAAACTTTTCATGTTTACAGTTCTTATGCAA
FGDEGFA:GGG@?BE?F;FGCFF5@EF3;5G#FA.G=D?D?AGAE?AFF=>CBDE32?<E?ED#GG?G?E#E#EEE>5#B####EFF##=D####6#:#,
+sim_22_1_9 385 chr22 15600023 0 1=1I1=1I96= *
0 0
CACATTTTCAAAAAACAAAAAAAAGTCTGAGCTCCTACTGTTGATTTAAATTCTTTTATAAATCTCTATCAAACTTTTCATGTTTACAGTTCTTATGCAA
FGDEGFA:GGG@?BE?F;FGCFF5@EF3;5G#FA.G=D?D?AGAE?AFF=>CBDE32?<E?ED#GG?G?E#E#EEE>5#B####EFF##=D####6#:#,
+sim_22_1_9 83 chr22 15600266 255 100= =
15600021 0
TTTAAAATATTAAAAGATGAATTACTATCAATTGTTTTGAATTTTAAACTAAAAATCAGTAGTTACTATAAAATTATTATTAAATGTTCTAATAATTGTA
5<###E?##AB#E=DB#<==;#E#?=AFEE@C@AEA6GDGFFG?DE??::?=@BC#G>GB>GEEDA=G>ADEFDEGABEDDEBCGGEFFF7;=GEGGFGC
+sim_22_1_6 99 chr22 18665735 255 14=1X10=1X74= =
18665943 0
ATTATTCAGCCTCACATTTCAACATCGCAGAACAGGTCAGGGAGTGCAAGGTACCTCCTGCTCTTTCCAGCAGGTGCAAGGTGCACAGGGAGAAGCAGGT
CDDEFG;?DGEDFE+CECEEGFG>C,EG5?DGE87GBDFGCEDD=BGGDD?GDFGDD?DAGGC<G=GF@A@CE@>?>.6#F#FFDEE?@CE==<9?B?>;
+sim_22_1_6 147 chr22 18665943 255 9=1X90= =
18665735 0
CAGGCCGCCCGGGGCGGTTCCACCAAGAAAGTGACATCTGGGAGGGGCATTGAGGCAAAGGGGCCGCGAGAGGGTAGGTGCTCTAGAGGTGGCCCACGGG
?##??3G#9##=5B#>A#.#>#>CFBEC@:BGC>.A5@6-F4:5?;EGGAA>?6@D<?EAEDE6CGGDGGE@9F?GGGGDD5EDDGFD?:>?FED?DEBA
+sim_22_1_16 99 chr22 20264219 255 84=1X4=1X8=1X1= =
20264427 0
TCCAGCTCCATCTGCCAGGTTTTGTTTTGCCTGAGCAGTAGGAGCCTAGAAATATTCCTTGCACAACACAGGCACATAATGACTCGCTGAGCAGGACTTG
FBGGFDGEFA?GGFGEF;AEGCEEGCFE?FE;GA:EAEEGDD>F8E=E=FABFCDCEG@@@E#BDBCBE;=#D5#,55DCAG:B8EGE#AAGBF#?#E#A
+sim_22_1_16 147 chr22 20264427 255 42=1X57= =
20264219 0
CTCTGTCCGTTTCTAATGGGAGCATTAACTCACTGCTTTGGGGGGCTGAGATAGTCAAAAAAGAGCAATTTCAGTGCAGGTTCTAGTCCTGGGCATGCTC
#;>6##<@#C@#8B#>#:GF2G@A2;BFA?B:-<ECG--)GE7F9DGD;/:ACDABDEBD5EG:EB@DDGB#@:AEA?=GG=DE?D5GD<AFEC;DF=GA
+sim_22_1_f 99 chr22 23113643 255 100= =
23113851 0
CTTTTCTATATAAAACAAATAATTTATAATACTGTATTTGTTTTCTATTACTATTGAAATTACTACAAATGTAGTAGCTAAAACAACACATATTTTTCAT
DE>FDG5FDGDGFDAGEFABDGGEGADG?BEDDFGGEDGA?GG=9FEDDGBD58=@AFDFD=CBD9GGD=DCC?EE?F#A6GC#>?F#F@2E?CDF?#B?
+sim_22_1_f 147 chr22 23113851 255 14=1X85= =
23113643 0
TGCCTTTTCCAGCTGATAGAGACAGCCATGTTTCATGGCTTGTAGCCCCTTCCTCCATCTTCAAAGCCAACAACAACTGGTCAAGTTGTCACATTGCATC
####FF-0AD#E>@#?#9??E##A>C=C>D7A,<G?D>=>DDD=F:?G=>CE?:CAFGEDGEFDAGGFEF7CGDCEDEEGA@DECFGE-?DEACD?DC=G
+sim_22_1_f 401 chr22 23113851 0 14=1X77=1056D1=1X2=1D4=
* 0 0
TGCCTTTTCCAGCTGATAGAGACAGCCATGTTTCATGGCTTGTAGCCCCTTCCTCCATCTTCAAAGCCAACAACAACTGGTCAAGTTGTCACATTGCATC
####FF-0AD#E>@#?#9??E##A>C=C>D7A,<G?D>=>DDD=F:?G=>CE?:CAFGEDGEFDAGGFEF7CGDCEDEEGA@DECFGE-?DEACD?DC=G
+sim_22_1_d 321 chr22 23964184 0 33=1X52=1X1=1X11=
* 0 0
CCTGTCAAGCCTACCTTGTCACCTGTTGCCAGGAAGTCCTGCTCCCTCCTAGGGCAGCCCCTGGTTAGTGGAGCTGAGAGAACCCAGGGTGCTGAGCCGT
G5FGBEFEEGFED=EGDABFGGGFGGCE?EF=EEDGG9G?GEAECD9FGBFCEF;GGEC=GC5>F1?CDCCD#FEB=EE5B#CFEB###:ECBA-CE:@A
+sim_22_1_d 99 chr22 24189222 255 86=1X1=1X11= =
24189415 0
CCTGTCAAGCCTACCTTGTCACCTGTTGCCAGGAAGTCCTGCTCCCTCCTAGGGCAGCCCCTGGTTAGTGGAGCTGAGAGAACCCAGGGTGCTGAGCCGT
G5FGBEFEEGFED=EGDABFGGGFGGCE?EF=EEDGG9G?GEAECD9FGBFCEF;GGEC=GC5>F1?CDCCD#FEB=EE5B#CFEB###:ECBA-CE:@A
+sim_22_1_d 147 chr22 24189415 255 100= =
24189222 0
CGGGGGGCGGAGCGCAGGCAGGTGAAACTGGAGGAAATTGGAGATACAGCGGGAGGGAGGGGGAAAAGGTAAAACCTCTCACCACCTCCTCAACGCTTAT
##??BG#C##5AB#D;=G#BF#E#:;@?#<(4BGCG-#EA@<G?FC*GDD@FD:BEBEGFGFDDBDEAD9GGD37DDFDCDFDGEFGFCGEE:?AE>EDG
+sim_22_1_b 99 chr22 24215590 255 23=1X76= =
24215802 0
ACAGTCACCCACAGATTTGCACACCTGTCCAACTTGAGCCAGTGATGAATAAAGACAGATATGGTTGCTGTTCTTTCATGGAGCTTACATCCTAGGGGAG
GGDEGGEAGGD?ED?>GDGDGFG.-GEFG:FE=FFGGFDGDDE==C==:@C7DEE=#GBC,BFFC?<@:A#:6AC>D?C#;EFD/FE#@CF#C=#B?#6=
+sim_22_1_b 147 chr22 24215802 255 20=1X79= =
24215590 0
AGTGGGCAGGGCAGCTTTCTACGAGGAGGTGACTAAATTAAGATCTTGCTGTGAAGATGTTTAGTTGCCTGCTACCCTCCACTTCTCCCTGAATAAAACA
A=#=C?A5.E##?#A##=#C#D?-A#E+E?CEFFCCB1;=-DG:DGFCECF=GA@FFFAECACF-;=DDDFDFB??F?F5?GGF:;FEBGFC;(CGDGG:
+sim_22_1_10 99 chr22 24428205 255 28=1X54=1X16= =
24428413 0
TAATAATCTCAAAAGATAGTTCAATGTAGAACCTAATTAATTCAGTAGGATGACCAGCTCAGTAGCAGTAAGTAGCAGTGCAATATAAGTGGCATTTTAA
EFE:FC@D?5FGGCAEDAFEGG3GEDA5(CFFFEGD#GEFFG?EDCFC?GBAGEFADCFBGCG@FGCF2F0ED>C@E:DAF:A-BF;:/#;#E##A#A?C
+sim_22_1_10 147 chr22 24428413 255 66=1X33= =
24428205 0
ATGCCTTATAAATAAATTTATAAGATGACTACTTTGTAGGCACATCTTGTTATAAGAAAAGTCAATGATATGTGTCCCTTTAATTTTCAAAATGAAAAAT
#GD###A#5ED#B5.#BE==5FG#E3-@#AB7CA)@DEFF=FF#@AEFDFD4>DF2CCCGGA:DGB'FFDG9AEF?3;DGAE?CGCFD2GAG?GB@GFGF
+sim_22_1_3 99 chr22 25791286 255 40=2X31=1X26= =
25791479 0
ATTCTCTGCCCTGAGTAAGGGGAATCAAGAGAGATTTCGATGCAGGAAACAGCACCCGAGAGAATCAATCCAGTTGGGCCTTTTCAATTAGTCTGTTTTT
?GFGFFAEEEGEF-GGEEGGEEFFGDF;??DGDGD/=DBD09EGG5EDG@DEC,FBGD@BF5:FBE6:[email protected]>EF#=#:EE#A>###
+sim_22_1_3 147 chr22 25791479 255
1=1X12=2X12=1X17=1X4=1X48= = 25791286 0
CGTCTCTCCTGCCCCCCACCCTCTGTCTGTGAAAACATTAGAGCTGAAAGACGACTCCCCTGGCGATAAACATCTACGCAAAGCCAATTGTCAGGCGATG
4####B#?C#???#6,D####-BAEF3<)DE>AD(#DD+BFEFFBBEAFBC1GABFAA2GEBGEDF8G7G=AAAD=EGDD:A4ABE?AFGDGD0DD?DB>
+sim_22_1_15 385 chr22 25992044 0 2=1X95=2I *
0 0
GAGAGGCAGCCATAGGTTCGAGTCCTAGCTACTGGCTAGGGGACCTTGGACAAGGCACTGATAGTCCCAAGACTCAGTTTACCCATCTGTGCAGTGATGT
CG.GCEGGGFDGGGGEDCE<BEG5GDGEE=##EEFEFABFDE:#G??>:55;FE:5=A?E;E9BF#E@5=AE;BEGCA:#>EAC#DED#B##A;###B#B
_______________________________________________
checkin-list mailing list
[email protected]
https://www.monetdb.org/mailman/listinfo/checkin-list
