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:FAE@B.EAABCBFCE5D>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>+g...@gde.aede+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>+g...@gde.aede+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:FAE@B.EAABCBFCE5D>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 checkin-list@monetdb.org https://www.monetdb.org/mailman/listinfo/checkin-list