[galaxy-dev] [hg] galaxy 2698: Updated BWA wrapper tool: replaced threads opt...

details: http://www.bx.psu.edu/hg/galaxy/rev/14d54c9dbd59 changeset: 2698:14d54c9dbd59 user: Kelly Vincent <kpvincent@bx.psu.edu> date: Wed Sep 16 11:14:39 2009 -0400 description: Updated BWA wrapper tool: replaced threads option to hard-coded value, changed output and test file output datatype to sam, changed test file names to standard, better cleanup of temp files, added dbkey parameter and exec_before_job method to set it. Also changed name of a FASTQ Conversions tool test file 17 file(s) affected in this change: test-data/bwa_wrapper_in1.fastq test-data/bwa_wrapper_in2.fastq test-data/bwa_wrapper_in3.fastq test-data/bwa_wrapper_in4.fastq test-data/bwa_wrapper_in5.fastq test-data/bwa_wrapper_in6.fastq test-data/bwa_wrapper_out0.sam test-data/bwa_wrapper_out0b.sam test-data/bwa_wrapper_out1.sam test-data/bwa_wrapper_out2.sam test-data/bwa_wrapper_out3.sam test-data/bwa_wrapper_out4.sam test-data/fastq_conv_in1.fastq tools/next_gen_conversion/fastq_conversions.xml tools/sr_mapping/bwa_wrapper.py tools/sr_mapping/bwa_wrapper.xml tools/sr_mapping/bwa_wrapper_code.py diffs (1880 lines): diff -r 1a56a541f9f9 -r 14d54c9dbd59 test-data/bwa_wrapper_in1.fastq --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/bwa_wrapper_in1.fastq Wed Sep 16 11:14:39 2009 -0400 @@ -0,0 +1,4 @@ +@081017-and-081020:1:1:1715:1759 +GGACTCAGATAGTAATCCACGCTCCTTTAAAATATC ++ +II#IIIIIII$5+.(9IIIIIII$%*$G$A31I&&B diff -r 1a56a541f9f9 -r 14d54c9dbd59 test-data/bwa_wrapper_in2.fastq --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/bwa_wrapper_in2.fastq Wed Sep 16 11:14:39 2009 -0400 @@ -0,0 +1,576 @@ +@1831_573_1004/1 +AATACTTTCGGCGCCCTAAACCAGCTCACTGGGG ++ +><C&&9952+C>5<.?<79,=42<292:<(9/-7 +@1831_573_1050/1 +TTTATGGGTATGGCCGCTCACAGGCCAGCGGCCT ++ +;@@17?@=>7??@A8?==@4A?A4)&+.'&+'1, +@1831_573_1067/1 +AAGTATGAGTCATTTACCGGGATGCGAAGAAAAG ++ +++#%(',%/$,+&1#&),)&+'*'*%#$&#%('$ +@1831_573_1219/1 +CGCCCTATAATAATACAGCGCGTTAGACCGCTCA ++ +@@=4/+)5)408?'665>*/5?<61';<3,:,5- +@1831_573_1242/1 +GCTGAATCGCACCTAGCAACTATAGCCGGGCCGC ++ +>>::>17=A5?@@=;7A=;2.60>82<8=74+;; +@1831_573_1333/1 +AGAATCGTTACCACACACTGCGTCTGGGTATCCG ++ +='@#%3=.>)/34*117,,/6-4+.9742456<) +@1831_573_1362/1 +CGATCTCAACCAGGTCCGCGCCCACCCCTGCCTC ++ +B@@?@@?@C@BA?@>@<@8A@?@'.8'?17:,+? +@1831_573_1448/1 +TCACGCCGGTCCTTGACTGGCGTTCTCTTCGAGG ++ +7=A)?,.@A@@:@@<=/7)@<-#8662%9613&+ +@1831_573_1490/1 +CTCGTCATGTTACGCAAAGGCACGTCACAGCACC ++ +@@@?,@??@$><=>+%/*>*1,)?26&12'9%7. +@1831_573_1523/1 +ATGGAACGGAACGGGTGAGGAGGGGAACGTAGGG ++ +,//*,#8''.6''6'2&27/80)%�'%#*,2' +@1831_573_1578/1 +CGAGTAGCAAACAAGACGCCAATCCAGGCGACCC ++ +?=@<@C@?@@??=>@;7<<8<694)3511(+1<. +@1831_573_1647/1 +AGGGGTTTACACTATTCGACTGGGTGAGAGGCGT ++ +?@@<@@???;A==@@;9?@?3=;6.4/;6:=/;4 +@1831_573_1684/1 +TTCAACTGCGTCGACGTAGCGCACAGGCGTCCGT ++ +>1:8<@39@@??9=;@@<>?@@<;<?7>76?9>? +@1831_573_1769/1 +TGGACGTATAGTGGCGATGAGCATGTAGGTTCTC ++ +54749&&767%/7;3$-7;3#,3//#,45/#&06 +@1831_573_1853/1 +CAAAACGCCCGGGGCCTCACATGCGCGGCAGTTC ++ +4./B@@72B+/4?@?7?+@9/+99.')2<2&)2& +@1831_573_1943/1 +ATAACGTATGGCAGTGAACGGGCGGAACCTGCCC ++ +>/*1,),@-)'0*>5'$/?6@/.&?8/(')A43, +@1831_573_1977/1 +GGCGTAGGGCTCATTGTGCAAGTATCCGACCTCC ++ +@?6@8@@9A@=4633A7?9?<3:+7>'@%@?%32 +@1831_574_109/1 +TCGGTTGCGTTACTTCATGGGAGGGCTTTACATT ++ +7651#4817>@65<2:84716=788<1995.7?6 +@1831_574_148/1 +CGAACCTCGTATAACGGAGTAGTCGGAAACAGTC ++ +;->3;;3>//>0;;7=,=4-4/+>029<9<624; +@1831_574_185/1 +CCGTTTTGCCTAGTAATGCTCGGCGCAGCGTCGC ++ +<4;%%@7>#*>97%6;+8$#<07%1372$%<54# +@1831_574_243/1 +AGGCACCGTAACTCAGGACATTCTCGATTAGTTA ++ +@@@95@=9?*><45$<@<4%@9?7#=?>7&>?97 +@1831_574_257/1 +ATACCTTCCAAAGCAATAGAATAAAAAACAGTAC ++ +>>B>)@?>B@B?@1%*%,#+00.'(+&5&%#$3& +@1831_574_293/1 +TGCTGCAAATAAACATACAGCCTTCTAATGACTA ++ +')7'.*%#&$'/.$##$.$&*)+*$#8%&%#&&% +@1831_574_389/1 +CATGGCTATGCACCGGTTTGTAGCGTACTCGAGA ++ +)$>/>8):2@:213;;.1)@3%64%+)<7)+$92 +@1831_574_575/1 +TTCTTGGCAAGCGCAGATTATGCGTTCCGCCTAG ++ +?@>=?@>8@@@=@6*%&<='5@<<',0<=4*5/< +@1831_574_592/1 +TCATTTACCACGTCAGGGTCGGAGTCATTCATTA ++ +3-5)19,,<@>4067<2.-864346;9<1/;212 +@1831_574_617/1 +AAGCATCGGCGGGAGCGCAAGCTGGGAAGGTGCC ++ +>@?;>?=@==@>@=5?@+@:@40@><:0)')64- +@1831_574_725/1 +GACAAGATGGCTATTATTTACAATCCGATCTGCA ++ +'>?52+5=;4$6;<8.1<47*9+?7(+5;-).)4 +@1831_574_734/1 +CCTGTACGAAAGAACGTAGGCATGGGCTGGGGGG ++ +-&.$2&<*(-%):/)%5:)/%7+,9034*A.?70 +@1831_574_824/1 +AGCGCAAATTATGCGTTCCGCCTAGCGGAGAACT ++ +@@@;;),+@<-<7;<&89??&29?<$2:A<);8? +@1831_574_959/1 +CGCGCTAGGACTCGGCCCCAAGAGACGTTCCGCC ++ +A=;+?@9;/<@;<8#&&/*%(%.0$0$2)+1$-1 +@1831_574_1062/1 +ACCGGTAATATAAGGCAACATGATTACGGCCACG ++ +@821<//:136/8<1/5.427.;>288/0<1$)/ +@1831_574_1092/1 +GACTGGCGAAATCATCGCGGAAAAACCCCTATCA ++ +8?.,98)03%#,,,(,++$'69'#9'7#';6.%2 +@1831_574_1103/1 +ATCTCCTGATTAGACATATCTCCGTAGCTCACGC ++ +??>?@=>@@A<>?=@==<:?<;58'&&):+35'1 +@1831_574_1116/1 +CACCTCACGTGAGTATAGCAGCCCGAGCGTCACC ++ +@++,)/)?./)%<)2>.==:8?&;44&/&,)/*4 +@1831_574_1194/1 +TTATCACATTTGGGGATCGGAAGGGACTACTTCG ++ +5>;,??9@9@?@1@@=@3=@>@4?@@94?75/,6 +@1831_574_1204/1 +CTTACTGGTCTGCTGGACATATAGTGGCGATGAA ++ +A5A5B@5:;@@,.9?A1?#.?;3),>82%)='7* +@1831_574_1306/1 +ATTGCTTAGATCCAGTGGCGCTCAATACAACGGA ++ +2.7&%''.$&#'%##$#'#%')')$'$%$,)'.1 +@1831_574_1387/1 +GTACTTCTCAATGCTGCACTACTATGTACGCCCC ++ +?A??:@?@@9?@@?9C;:>&).,,);&'<&7(/' +@1831_574_1431/1 +GACCAGTTTCAGGGCTAACCGTCCCTACTCGACC ++ +7?>9BA897;+8*=6#6.+5&-#26$0.,5&'$1 +@1831_574_1560/1 +GGCGTCTTAGGATTGAAGAGGGCCTCCCACCCCC ++ +))&65<4=>99@:597(1&;+&(%&4#,#&&&*$ +@1831_574_1591/1 +TGAGCACTTATGGCTAGGCGTAGGGGACCGTGAG ++ +=<A><?@78@<=;@&@6>8:?,><&&<7>8(8+9 +@1831_574_1624/1 +ACGGGAAGGGCTGGAATCTACCCAGTAGGCATTG ++ +<6.&(,7-+))7..'/;43.)927.(<85&13-) +@1831_574_1826/1 +TACGTCGCGACCGAGCGTTATATAGTCTGACCCC ++ +@<<,))459<?.+57'-'&6?1)(;1;0&%4,#$ +@1831_574_1903/1 +AGTGCAACATCTGCTTTGCTTATCAGCACACGGC ++ +956?@>;7@@@8:@>8><=9>9>@2=71?1:6;/ +@1831_574_1961/1 +GTTTCACTTCGGTTATTAAGAAACCCAAATGGAA ++ +>-07-@>(2'@<#064@%5%@@6.&<<(9)=:47 +@1831_575_54/1 +TTTCTTATGGGTAGAACAGCTGCCACTGACTGAA ++ +=@@68<2;<8;?@<>+9>7,,6972&57:4&9<8 +@1831_575_80/1 +TCTTTGGGTTTGGGGCAATTTGGTATGTTCGTCT ++ +?5><@?/=5%(&.5)+#/0'/'>#$/'37#()6& +@1831_575_192/1 +AACTACGCCCCTTAATTACACAGCGCGTTAGACC ++ +:8<,9(.=/)#%8$,4'<,:)3/.'.5<,>-;7' +@1831_575_197/1 +TTCGCCTACACTTAGATACCTCTTAAACTCAATG ++ +-:?463@A:88?6?437<=.2&1.;@67.4@%2, +@1831_575_223/1 +ACGCACAAAGGAGCTCGGCGCATAGCAACGCAGA ++ +?>?8?@:?;>+$'%&,+#$.--#%$1'&#-7$&' +@1831_575_420/1 +CCCACATGGAAAACACTCACCGCCGAACAGAGCG ++ +2?.>?(-''2+..@),$%1@;')1)?+)7.').. +@1831_575_434/1 +ATCGCTGCGAGGTCACCCTGGTTACGCCCCTTCC ++ +@@><@B@>=@@?>><@?;)@?;6.?=4,'';/3) +@1831_575_444/1 +TAGGCGACCGTGAGGACAAGAGCTGTTGCCTTGA ++ +2@@;<.=?@<+=:A<<@?<</>@><(<8>4&5;> +@1831_575_459/1 +CTTATCAGCACACGGTTTACCAGTCCGACTCCAA ++ +@@C@?A@@A@??A@@'6@@B'=;C@'4<=9,@1. +@1831_575_506/1 +CGCAGAACCCGCACGCTTGTGCTCACCACTGTAC ++ +?)2*5@3+.7=?A:8@,/&):@37/<89'+4*)* +@1831_575_569/1 +AGCAGACTGCTGTAACACGGTGTGGTGTAAGGAT ++ +'9-,@29+479%&>9#&08=>7#9>,69/+'*3) +@1831_575_622/1 +CCAAATCCGGCCCAGTAAGTGTCCTGTCGCACCC ++ +A8B=@@?<@?<>>=?96<497)<4/<':4+;74+ +@1831_575_644/1 +CTTGAATGGCGATCTCGTCGAGGAAATACCACTA ++ +<8>.>&,+)8.10':5(/-62&*.=/%(.<&&18 +@1831_575_663/1 +GGCAACTTATCCGCATTGGTCCCACTGGAGCGCA ++ +@@@@A</;A>@6>0792<5966:?;25+:><1,+ +@1831_575_681/1 +TCTCCTGATTAGACATATCACTGTTGAACACAGC ++ +72?.)25;96-14'*.03'#4#++$)/8+%/((& +@1831_575_711/1 +TATGTTCGTCCACGTCAGACGCGCAAAGTTGCGC ++ +%+:2(#'1+.)&*45-#,<A-%4;>+))<;##*& +@1831_575_730/1 +CACACAGGAACCATAGCGTATGTTATTCACCCCC ++ +%8<%-')5,*1&=#)(1*1)->7$.,4'4'&'*& +@1831_575_904/1 +ACCCGCTTAAAGACGTGAATTTTGCCTCCGCGCC ++ +<><.689@@7@<54:4=@>69=:+6-=42;8<:# +@1831_575_938/1 +TCATCAGGGAAGGCTAGGGGTTTACACTATTCGA ++ +:@@:?@;=8;>=;;:=<93<4.;<6:2<9;4)87 +@1831_575_970/1 +TGACTCCTACAGTCTTTATAGTACCGAGGGAGGC ++ +@=@=;@@><@@:9<:=?;>@=?9>>?@:2===<5 +@1831_575_991/1 +TTCGGCGATCCCCCCCACGGCGCGATGCCGCGCA ++ +>>;7???/4<:::,#&*%'$--5-#(7)&&,-)% +@1831_575_1138/1 +TTGAAAGACCTGAACGCGATTTCATGTTTACTGC ++ +=0;%1936+?96?,?5<:2>@6@<:;4>.7;*:1 +@1831_575_1157/1 +TCGCTGTTTAGATTTCGGGAGGGTACTTCAGTGC ++ +B@B>?A?A@?9'?;?&8?@@6;<@=+<:99+965 +@1831_575_1180/1 +GAATTCACGGCAGTGTTATCACCGTTTCCTTCCA ++ +<B>5?':/78?;1;55:997@@>94=7<186->7 +@1831_575_1283/1 +GGTGGAATACTAAGGACTAATGTGCTGTCTCTTT ++ +C@A54><?A78@<?2>:#>=?@2>80:'?69.'7 +@1831_575_1302/1 +CGACTATTCGTTTCGTCTAGAACGTGACACTAGC ++ +12>/#2'&-2&%'7+$%*#$/1(&&12/$,2,,* +@1831_575_1310/1 +CTTGCTCTCGAGCTATGCCTCAGGACACGCCCTT ++ +<<>663?99@(@7=11@:@@/7??<@<;>'3-7; +@1831_575_1321/1 +CAACACAGTTGAACGGCGGAGGAGTAAAGATGCG ++ +;B@?@>=A@<A?>@>?6>@5>'4@9?)7@7?/89 +@1831_575_1373/1 +CGCTACCGGTTCCAACGGCTGCCTGACTCGCGGA ++ +@@@>@;=9<@??7@>@:@;@;17?630/7=4,2< +@1831_575_1419/1 +TGGGGAATATAACAGCGTAGCGTTGAACACTAGA ++ +//3->&-.->%11/@&/&7:<,-1=,+2<6'%&< +@1831_575_1436/1 +TGCATCTAGCGCGTATTGCGAGCCTTACGGATTC ++ +2@=;39*3<*4<:>52>>9'7:<='0=>1'136< +@1831_575_1442/1 +TCTGACAAGGTTCCTGCACCTGCGTCTGAGAGGG ++ +>58<?7>+1B6-/9<<:5?>A:1@?539.@4/1> +@1831_575_1454/1 +CCTCCTAACCACGAGCCGAGGGGTCTCTGCCCCT ++ +4=?8=6589=0:A=:?>6<;@/,4+-6+/'.,(, +@1831_575_1500/1 +CACAAAAGGTCCCTACCTGTCTACCCTACATAGC ++ +?@@?A@?8<6<:985=99999==B3<5<<,7?>9 +@1831_575_1535/1 +CTCGACGATATGACCGCCAGCCACTTAACTCCGC ++ +B=@<@@B?:A@@=>@?68@=?@5?7;:6<<;>5< +@1831_575_1724/1 +TCGTAAGTGTTAAGGAGCTGTGTACAAAACAACA ++ +<B.+@),2*.%)))4%2@;7#%(%+$8))85%&5 +@1831_575_1829/1 +CATTTGCTGACCCTGCGTAGTTTAGTCTCACAGC ++ +2(*'%.3$+7)@&%$'3*+*#/#/*+0.=&#)+0 +@1831_575_1898/1 +CTTACCATATCATCTCAACCCATAACAGCACTTA ++ +B>@@&1/))'40)%#8/.%#8$((#;4'$'63,, +@1831_575_1964/1 +GACAGACCATGAGGCTGAAGACTACTAAGTGCGT ++ +>39)@2<2/@+9?2=&)>>@*62=5&2<42.'?+ +@1831_576_32/1 +TACGCAACGATTTATGGCCTTATAATTGAGGCCA ++ +<>?>?=7<2)522;><<@40@>704<>5=23@+& +@1831_576_74/1 +ACATTCTGCAGTGGGACAGAGCGGTACGCCGCAA ++ +=@@?@>(??<B5?@@@9<@><+><;@';>6961? +@1831_576_86/1 +ATGAAAACGCATTAGGACAACCATAATGGCCTCA ++ +811)2:*.++5</:3+43924*))/:,6&29)2/ +@1831_576_89/1 +GCTGTTTGATTTGAGAAGAGGAATTAAGCGCCGA ++ +)69+,'.4=-,>/>(*$#)3030*'(,%)2##$, +@1831_576_266/1 +ATGGGGTCACTCGACCTAATCCCGCGGCTTTGGT ++ +@@?=;<?7=@>9,>@1$&&89$/:>7'3178%&6 +@1831_576_327/1 +GCCGTTCTACTCTAGCTGCAACTTGCGATTGCTA ++ +@?@6@@;<@?><@>==3@:==<3@@>53<9><6= +@1831_576_331/1 +GACGCTTTACTCCGGTAGTACCGTGCCGTTTATA ++ +-2)+(*.*1/;5%.-9&#/1'+($*$##()%/$. +@1831_576_387/1 +ACACGCCATGATCCGATAAGAAGGGAACGTAAGG ++ +8).$5#1#*%.$##.*#$%##-%,+,1#&%.%)) +@1831_576_406/1 +AGGTCTTACAGCACGGGGCTGAGCGCATCTGACC ++ +4<>%%14:*4656)&<251&2+3#&19,6&4>5( +@1831_576_449/1 +CTCGAACCGCGGGGTCCAAAGACTGCTGCAAGGA ++ ++'&'.,,$/+.)$$8&%#+?&,#)-&###7,+#* +@1831_576_519/1 +TACCTGCCTACTACTTGCTCTCGAGCTATGCCTC ++ +=46<97@>2/6?;2<4A881>9121+<1/4.9+7 +@1831_576_603/1 +CAATATGTCTTAGTCGTGACTCGGCAACTTATCC ++ +@46=@C.??<A79@@;-<@@>29B?>55<B7598 +@1831_576_655/1 +GAACAGTCTATAGTGGCGGGAATCTCGTCGTCAG ++ +@AA9@@<6*>@@5/<@>9'=;>7+@?9>/9;+,% +@1831_576_677/1 +TTTACTCAGTTGATACATCACTGTAGCACATAGG ++ +290&/*0#&'&,.2'#&*$&('#-%($*#%$)#% +@1831_576_718/1 +CGTGCCTTTCAGGGTCTTTTCTGGTCTGGTCGCT ++ +4###$0###%(#########,####%####$### +@1831_576_722/1 +CGTATGATGGCGAGTCTTTATAATCCAATCTGAA ++ +*,-,##.)*&(*1%*(%(-2#+)-#.&-#%%$') +@1831_576_754/1 +AGGCGTCCTGCATCGACCGTTCTATCCGCTTAGA ++ +4//#(.$)'',>(
$<,##%((,#5?#0*%1* +@1831_576_815/1 +TAGGCCTGATATGACACGATCACAGTGCAACATC ++ +3>@?9>@?B=>=;>A???=>:25=4.25?6<57. +@1831_576_882/1 +TGTAAGACGGTGAGGTGTAAGGATCAGACCACGG ++ +26?8?@:4>@>96??<<=5'1<>9846=<9<1>8 +@1831_576_898/1 +AGTACTGTCGCGCATTGGGGTCCTGGTCGTTGCT ++ +@@C2?C>?<>@886B?;?.??87=B<8<15??=. +@1831_576_923/1 +CTGGACATGAGAGACTGCACGCGGTACACGTCGG ++ +6.51=;.699<96>;;49<;;11;<@59:9=647 +@1831_576_930/1 +CTGGCATGTACGTCCATGTCAGACGAGCAGAACT ++ +#*0.2&.:((#'14'##-)#%$$2%#$/1&#%/# +@1831_576_1019/1 +GATGCGCGCTGTCATGGCATCGAACCATCGGTCG ++ +<$>$,1,&++&@,.)'+/+#9'69/6'2(+-'9- +@1831_576_1068/1 +AAGAGTGACTCACTTACCGGGATGCGATGGAGCC ++ +>><=4?;@7=??9?;9>5@9?:8@:=5';7;'#4 +@1831_576_1131/1 +AGTTCGGGAAGGGCTGGAATCTACCCAGTAGGCA ++ +04@#@/@,4>9?2+?1571@',>=;(759;*92< +@1831_576_1168/1 +CACTTAACTCCGCTGTCGGAAGCCTTACAAGACA ++ +@CC6@@2?9=>7?;76<;467@;9,0%26'',4$ +@1831_576_1207/1 +CAACCTGACTAAACGGGGATACGCTGGCGCTACA ++ +$>2,(&?4?(,@:<&,@>?$&:8A%%=0.%,597 +@1831_576_1289/1 +TAGCGCAAGTCCAGAATGTTCATAGACTCGCGAT ++ +9<>1B@@;7@@663==28,5':8<<,.=2>>.50 +@1831_576_1329/1 +CCAATAGCAGAGACCTAATAGGAAACGAAAGCAA ++ +>&0,870/A@@/5.;=;:'&@3'&$$%,+#($&2 +@1831_576_1367/1 +GGTCTCATCCGTTCCAATCGGGACTTTGACCAGT ++ +?@@<>@@=@<@<=@>@??9:?<=>=<8;59@787 +@1831_576_1416/1 +TAGCGTTCAACGTCGATCTCATCTTGCCGATGGC ++ +';?..1.<@'.=+6.5?7<0-?7;(%=>56.98@ +@1831_576_1461/1 +GAGGGGCGGCCCGGTTCAAGCAGGTAAGCAACAA ++ +1&74#(#(&##$#,'###%#%%#*#&%%##%%%' +@1831_576_1605/1 +AGTGCAACATCTGCTTTGCTTATCAGCACACGGC ++ +<@C@@>1;@:;<<@@9@???9<3?5-21=4877, +@1831_576_1664/1 +CGCGCACAACTCGCCATGATACGACAAGACCCGA ++ +@?>?@=A:@@>@='=<>=).>=9*8)(7#/++;) +@1831_576_1671/1 +GTCTTTGTAAGCCTGGCCTGGTCAGGTCTGGTCT ++ +@?@5>@@8&8?25#&)&5&&)15&&,&4%&,&6% +@1831_576_1729/1 +CGTTTCGTCTACAACGTGACACTAGCACAGTATA ++ +?@?<?<@>A?/==?>?7>?@8<?@>99;><+0=> +@1831_576_1880/1 +TATGCGCTGTTGAGCTTACAACTCACTAGCGAAT ++ +515/792,:,7/%/05,%$):+#8%2(1754))3 +@1831_576_1982/1 +AATGTCGTCAGACGACTTTGGCGCGAAAACCATA ++ +@@?@,A><@9@=?@9B=8<6@@6@2<A?>.7<+@ +@1831_576_1987/1 +AAGGTCTTCTGTCGGCGCTGGACTGAACACCTGA ++ +##()(#&##)%)/-#%$11#%3>'##&$,#$$'+ +@1831_576_2014/1 +CCGTGACACACAATGCCGGCCCCAGCCTAGCAAT ++ +@@9,B29-5>'?,+?79+/A';'2@'5&/9,6&1 +@1831_576_2028/1 +ACTCGCCGCATCCCCGAGTGACGCTCGAGACCAA ++ +<>9>5'2(&707.8#&&39,0%7/#(#,*%&5*+ +@1831_577_40/1 +CCCCGCGTTACGAACGAGAGAAATCTCTTATAAT ++ +2/4('.')(<*?#$)%&<,/39<(.2,+<=@611 +@1831_577_119/1 +TCCCACAAGCCATTGACATGCTCGCTCTAAACAG ++ +>?C11)-1/)#;#/,850*+.+$$5550+%-.40 +@1831_577_133/1 +TGCTTGTACGGTCTAACGGGGTATGGGTTTCTGG ++ +4@>>?=1?:+>@07@@><>@*<+1@15)96'2$) +@1831_577_255/1 +ATTGAGGCCAAGATAATTGAGGAGAGAGAAGGTG ++ +<=87@>69%**#&#-+$.#&&#+$-+%&%.,%5# +@1831_577_281/1 +TATGTACGTCGCGTACACTCCGGGGCCCGCAAAA ++ +-:B;8@=;>7@>4?=?+659?;5<7?;9@8(>:? +@1831_577_288/1 +CATCCGAGGCTATCAAGGCGTAAGCACTGACCTA ++ +;/79??&=B::298*6.7/+4&21,7,6?.7#'6 +@1831_577_322/1 +GAATGCTGGAGTACATTATGACAAACTATCGGAG ++ +=;>4506;255464-<#7+194&2<?65968)7/ +@1831_577_362/1 +CGATTAGTTACCACTCGTATTCGCACGCCCAGGA ++ +2B?:@?@<?958=,3:,90:&'-99,6<5.($+6 +@1831_577_382/1 +GTCGCGTATTCCCCGATGCTATGTAGACTTGCAA ++ +2957/:1))=76(*24;,3+:<.&.&-=1=2/5* +@1831_577_464/1 +TAGAGGCACCCTAACTCAGGGCTTTCTCGATTAG ++ +=<7/<826)>#.'&4204+5#/041.7*91&756 +@1831_577_488/1 +TGAATAGTTATGGCCAGAATGTCTGCACCGATAC ++ +48.,*>6566<?8=<=<2>6;94>;=9>@8924@ +@1831_577_511/1 +GGTGCTTATCAGTTCTTTCTCGGGACTTGTATTT ++ +7/+&$:<7%6,$$%'%/+)#$7((&*3>16'0/+ +@1831_577_545/1 +ACCGCTCTTTGGGTATGGGGCAATCTGGCATGTT ++ +><<9?<>?A<481@<@8==@76/61<95.5988- +@1831_577_559/1 +GTGCCACTATGTTCGACAGACCCTAAGGCGGAAG ++ +<@;??@>@?/2<.@1=>12=61/;=?.&2+92)' +@1831_577_562/1 +GTTCCACTACGTTCCACGCAAAGGAATCCGAATC ++ +((&*&0%)1%)#($2-,***%/-,,))&,-.1'* +@1831_577_637/1 +GCCTTCGCGGGAGCATATCAGTCGATACATCCCA ++ +-,4#>:-6+:8,&(5;3=0>7=68&1/9&'?;4, +@1831_577_641/1 +TATCTACCACCGCGGTGGCGCGAGAATGCTCCCT ++ +47;/':A;;5?:72,(=),#*?+.#&7$8#%7/' +@1831_577_692/1 +CCGGTGAGAATTACATCGCGAGTACGCCCAAGGA ++ +B,@?@B(@A?@+@@>@+?=>@'=<@<=<9=?75< diff -r 1a56a541f9f9 -r 14d54c9dbd59 test-data/bwa_wrapper_in3.fastq --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/bwa_wrapper_in3.fastq Wed Sep 16 11:14:39 2009 -0400 @@ -0,0 +1,192 @@ +@1378_11_329/2 +GTTCGTGGCCGGTGGGTGTTTGGG ++ +###$$#$#$&#####$'$#$###$ +@1378_17_1788/2 +TGGGTGGATGTGTTGTCGTTCATG ++ +#$#$###$#$#######$#$#### +@1378_25_2035/2 +GTGCGTCGGGGAGGGTGCTGTCGG ++ +######%#$%#$$###($###&&% +@1378_28_770/2 +GGTTGGGAGTTGGCGCGCGTGGTT ++ +###$##%%$####%#%$#%%#($$ +@1378_33_1945/2 +GGGGGGGGGGGGGGGGGGGGGGGG ++ +.6;?;55697>8779457986264 +@1378_34_789/2 +GGGCTTGCGTTAGTGAGAGGTTGT ++ +###%$%$%%###$####$###$#& +@1378_35_263/2 +GGCGGGGGGGGGGGGGGGGGGGGG ++ +$%%&,)3/)%$-6'$&5&#$''5. +@1378_43_186/2 +GCTAGGGTTTGGGTTTGCGGTGGG ++ +$%#$########%##%#$###'#' +@1378_51_1671/2 +GGTGTTGTTCGGATGTAGCGTGGG ++ +##%$$####$########$##0## +@1378_56_324/2 +GCTGGGTGATCTAAGGAATGTAAT ++ +#$#####$$$#$$$$##%###*(# +@1378_56_773/2 +TCTGGTCGGTTTCGGGGAGTGGAA ++ +##%%#&$###$#$##%$####%%$ +@1378_62_2027/2 +GTTGGCCTGGCCTGCCGTGCTGCG ++ +*##),/%##$)#%##1$#'%.#&# +@1378_62_2029/2 +GGCGGTGTGTGGTGCGGCTGTGCG ++ +/$$$=(####%####)$$%$-&%# +@1378_67_1795/2 +CGTGCCAACGACCGTCCTACGGAG ++ +@@.?7@@@.@:@A<<>>=>70<@7 +@1378_68_466/2 +GTGTCATCTGAGGTAAAGCATTGT ++ +/##$09#$#.=$#$76+$%1'### +@1378_68_1692/2 +TGTTCGATTTAGGTCTTGCCTCGC ++ +###%%$%%##%$#$##$###$$#$ +@1378_80_664/2 +TGTCTGCGTTGTATCTCTGGTGTA ++ +%##%,%$$#&$$###$#$%##'%# +@1378_85_1786/2 +CCCTAGGAGCGTATACCGGACGAG ++ +,'&/%/@,&1,&'/)&,6&&1)(( +@1378_86_1011/2 +AGGCGATGGGATATTATTTTACTT ++ +:$###)%##$9$###1$$#$2### +@1378_86_1789/2 +GGCCAGGTACTAATAAATTCAACA ++ +1.51/&%#'$,*?%#6&0$$5,20 +@1378_91_1596/2 +GCTTTTTCATTCGGTGCCTTTGGA ++ +'>%/3%=()8'#.%?50$&5>%)% +@1378_94_1595/2 +GCGGGTTAACCCAGTACCTCTGGC ++ +,8@6'>8<76>@@(->9)7'*941 +@1378_95_1039/2 +GTTCTGTGCCAGGTGAGGTACGGA ++ +&##,./#$&)6##+,'#$$0(##$ +@1378_95_1767/2 +CTAATTGACCGGGCAAGCTATTAA ++ +',,2'(10/>:=/7).3&%6&)&6 +@1378_96_1037/2 +CTGCTGGGCCATTTGACTTACTCA ++ +'$#+#(##-%5##+*&###-.$$$ +@1378_96_1764/2 +TCTTTGTAACCCACTTAGTATTTC ++ +4='28>.@4881B807:822>%)1 +@1378_98_1574/2 +AGGCGAGTGTGGGGGTTGTTTGAG ++ ++%%$#)##%##$####%###$%$# +@1378_107_1647/2 +GGTCTGGTTCTATGTTGGTCGACT ++ +###'$$#$$$(#%###(#$##$%# +@1378_111_829/2 +GGTCTTCTATGGTGCTAAATTAGT ++ +$##'&###$##)#%($#$%#$'%# +@1378_111_1900/2 +GCACGCCTTTGGGCTAAGCCGTAA ++ +)$)'#%$########$'#&%$#(# +@1378_112_1483/2 +TGGAGTGGTGTGTTTGCTGAGCCA ++ +#$#)#############$#%#%'% +@1378_125_1287/2 +TGACGTGGGTTGTCCCGTGAGATT ++ +##$%%#$###$##$$#&%##$(%% +@1378_126_468/2 +AAAAAACGGGCCGGGCGGACTAGG ++ +@@B@@@>@2>1+?:92)>@8&&<+ +@1378_127_664/2 +TCGCTTTGCCTATGTTTGTTCGGA ++ +#%$%#&##$%#%%###$$###)-' +@1378_129_463/2 +AAAAAAAAAAAAAAAAAAAAAAAA ++ +,)&20%8'(&35-***1-2&+510 +@1378_129_875/2 +GACCTTTACGTATTGGGGGTTGGC ++ +###)###+###$##$#&%##$,#$ +@1378_140_1251/2 +TTTCCTTCGTGTGCGTGCGGAGTG ++ +#%#%$##$$$######.$$$%#%( +@1378_141_809/2 +TCTCGTGGTTTCTTTTTTATGTGT ++ +##%)##$$#####%$#$#%%#'## +@1378_144_983/2 +GTTCGTTCGTGGTGTACGAGGGTG ++ +#(#%#####($#%##$$#%##%#) +@1378_153_270/2 +GGCCGTGTGCGGGTGTAGATTGGA ++ +%$##($######&##$&$$$$%## +@1378_155_1689/2 +CGCGTTCGGACAAAGCTAGCACCT ++ +@4@@=6=@@<4:@@681@@@.24. +@1378_157_1580/2 +GGGATTGAAGGGATGTATGCTAGG ++ +#%$&%#$$'%$%#$##*#%$$$$# +@1378_161_317/2 +AATCCATACCCACAAAAGCAGGCC ++ +.&%','(@''?7//+&)+2.+)0) +@1378_177_735/2 +CGAGCCCTAAACCATGAGATCGGA ++ +@@B>@AC9@A<B8@=9>7@5)>,0 +@1378_181_1684/2 +TTTCTGTTGTGGTTTTGTTGGGGT ++ +$##'$%'##%##$%$#$$####$* +@1378_187_1407/2 +TTGGGTGAAATCTTGTCGAGTGGA ++ +####&##$$###$#####%##%%) +@1378_203_721/2 +CCTTTACGATCATAAACCATGACC ++ +9<&./1&@;12')?<1',?/)&/. +@1378_206_2039/2 +CCCCCCCCCCCCCCCCCCCCCCCC ++ +?=>>?=:;<>>4::<=56199(05 diff -r 1a56a541f9f9 -r 14d54c9dbd59 test-data/bwa_wrapper_in4.fastq --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/bwa_wrapper_in4.fastq Wed Sep 16 11:14:39 2009 -0400 @@ -0,0 +1,192 @@ +@1378_11_329/1 +AGACCGGGCGGGGTGGCGTTCGGT ++ +%##+'#######%###$#$##$(# +@1378_17_1788/1 +TGCCGTGTCTTGCTAACGCCGATT ++ +#'#$$#$###%%##$$$$###### +@1378_25_2035/1 +CTGCGTGTTGGTGTCTACTGGGGT ++ +#%#'##$#$##&%#%$$$%#%#'# +@1378_28_770/1 +CTGCTATGCCGGCCGCCTGCCCCT ++ +%##(#$$#%$#$%###$%%$##$# +@1378_33_1945/1 +GGGGGGGGGGGGGGGGGGGGGGGG ++ ++9;64287788752796/576352 +@1378_34_789/1 +ATGGTGGCTGACGCGTTTGACTGT ++ +#$##%#$##$&$#%##$##$###$ +@1378_35_263/1 +GGGCGGGGGGCGGCGGGGGGGGGG ++ +),,#%'$)'-(%&%'%',%%&&(& +@1378_43_186/1 +ATACTAGTTGGGACGCGTTGTGCT ++ +#$(4%$########$#$###$$$# +@1378_51_1671/1 +TGGGTAGGTGCGGCTCTTGCGGTC ++ +##$#$%###%%#%#$##$$##### +@1378_56_324/1 +AGGCTTGGATTGTTGCGCTGACCT ++ +/+#**###%%##%#####$$$#$# +@1378_56_773/1 +TGTCGTGAGGTCACTTATCCCCAT ++ +&%#%##%%#####&#$%##$%##$ +@1378_62_2027/1 +CTTCCACGATCTGCTCGCTGTGGT ++ +(#&&$##$$#$%#%$$$#$###'# +@1378_62_2029/1 +TCTGGGCTGTCTTCGGGTCGGTGT ++ +$%$$####$##$$#)##%%#$### +@1378_67_1795/1 +TAGTGGGCTAGTCCGGGATCTACC ++ +B@C2@@?;-@@1>'A@@91@.@1' +@1378_68_466/1 +GTGATCGTCGGTGCCAGTCCCTGT ++ +#(%)+##$#$#%#+$%##$##### +@1378_68_1692/1 +GGTCTTCTGGGCACGGGCCAAGTT ++ +###)%'&%##$###%#######$# +@1378_80_664/1 +CTGCTTTGATCCCCGGTGGAGCAC ++ +7#%###$$6#######$##$$$## +@1378_85_1786/1 +ATACTATGTCGATCTGTAAAAAAA ++ +)&.)#3%@$&%-,2#&+.-%0&./ +@1378_86_1011/1 +CTACGTTATTGCTCTGTTTGTCCT ++ +######$%##$$$%###%#$#### +@1378_86_1789/1 +CGAATATCCCCTGGTTTAGAAGTA ++ +61?1;@/>@<=/(.7'(<5@=%$= +@1378_91_1596/1 +TTAGCGGTTGACTATCTGCTGACA ++ +*&+'#9'(%*'#//,&<),/)'*# +@1378_94_1595/1 +CGTGCGACAGCCCATGTTTTCAGA ++ +-=..5,3826&*+.+#+#%%6;%# +@1378_95_1039/1 +CGGCGTCCATCTTCGCCTTGAGAT ++ +$##.#$##$$#%$#$%%$###$)$ +@1378_95_1767/1 +ACACGCACCATCTGGCGGCTAACC ++ +86<65<4<C27/$5+#;-$+#2%$ +@1378_96_1037/1 +ATCCCCCAAGATGCCTGTTGATTG ++ +$#$'##$$$#%$$#%###+##$#$ +@1378_96_1764/1 +TACACCGGAACCTTGCACTAAAGC ++ +6-9'<#+81?+*<8-=69&6;*3) +@1378_98_1574/1 +GTTCTGCCGGTGTCTGTGGCGGGC ++ +$$#+&$$####%$$$###$%#%%# +@1378_107_1647/1 +AGGCCTACTACGCGTCATTGATAG ++ +&#$$#$(.#%#$$####&$%##($ +@1378_111_829/1 +TGCGGCACTTGCTTCTTCGTATTT ++ +%#%##%#$%#$#%###$$##&#$$ +@1378_111_1900/1 +TCCCCTCGCTCGGCTCTGTGCTGT ++ +$&%*$#(#)##$#'##%(##$#$% +@1378_112_1483/1 +TGTCCAGCTATGCGGCTTCCTCCT ++ +%#$+#%#&#$#####%####%$## +@1378_125_1287/1 +TGTCTCTGGGGGGCCTGGTTAGGT ++ +$##13$'%#$###$$###$$$#&# +@1378_126_468/1 +TCTCCTGATTAGACATATCTCCGT ++ +7>@318?=,=8@49579?..7,7) +@1378_127_664/1 +AGAGGTTGGTGTCTTGTCGCAGCT ++ +##'#$######$$%######$$$# +@1378_129_463/1 +CCTGTGGGCCAAGCCCAATGAAAG ++ +8),,<6(--)<()5&.&/8+;.2+ +@1378_129_875/1 +TTTCTATGGCTTACGCTGTCTGCC ++ +#$($##%####%$#$#####$### +@1378_140_1251/1 +ATCCTAGCGCGGTGTCTTGGGGAC ++ +#$%1#$$$##$##$#$#$##$%$$ +@1378_141_809/1 +TGTCCTCCAGTGTCTGTTGGGTGT ++ +%&,-##$$#(%###$#$$'###'# +@1378_144_983/1 +AGCGCCCGGTTGGTGCGGCTCGTC ++ +-$(&%*$#*#))#$$$#%%$#$## +@1378_153_270/1 +AGTCCTTGTCCCCTGGGTTTTCCC ++ ++''$#&%$%#$##&$$($#&#$$# +@1378_155_1689/1 +TTGGGAGGGAGAGAGACTAGACCG ++ +?4@@4<<??<?-@4@=4<5<?:/> +@1378_157_1580/1 +TGGGCCTCGGTGCCCTTGGTCTGT ++ +#%)$##'#$$$&#####%#$#$## +@1378_161_317/1 +TTGGCCGGCAACCCCGGTACCTAA ++ +7<,<'@)@>.)2@/')'&(?/-<( +@1378_177_735/1 +AGTACGCCATGTATTTGCGACCAG ++ +=?@?C<88@=>:7>@55/.,416, +@1378_181_1684/1 +CGACTCCCGCATTCACGGTCAAGT ++ +&*#,##$#&$*$$#$#$$$#%$## +@1378_187_1407/1 +TGGCGTCCACTCGTGGGTCTATCG ++ +$#$'%#$%$%&$%#####$#$#%# +@1378_203_721/1 +CCTAAATAACCCAGGGTAAAAGAT ++ +9<1+2+1@=259@+;2.71&@+2& +@1378_206_2039/1 +CCCCCCCCCCCCCCCCCCCCCCCC ++ +9@?8>;>>><<994;=895895-8 diff -r 1a56a541f9f9 -r 14d54c9dbd59 test-data/bwa_wrapper_in5.fastq --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/bwa_wrapper_in5.fastq Wed Sep 16 11:14:39 2009 -0400 @@ -0,0 +1,4 @@ +@081017-and-081020:1:1:1715:1759/2 +ACGCTCCTTTAAAATATC ++/2 +IIIII$%*$G$A31I&&B diff -r 1a56a541f9f9 -r 14d54c9dbd59 test-data/bwa_wrapper_in6.fastq --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/bwa_wrapper_in6.fastq Wed Sep 16 11:14:39 2009 -0400 @@ -0,0 +1,4 @@ +@081017-and-081020:1:1:1715:1759/1 +GGACTCAGATAGTAATCC ++/1 +II#IIIIIII$5+.(9II diff -r 1a56a541f9f9 -r 14d54c9dbd59 test-data/bwa_wrapper_out0.sam --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/bwa_wrapper_out0.sam Wed Sep 16 11:14:39 2009 -0400 @@ -0,0 +1,2 @@ +@SQ SN:PHIX174 LN:5386 +081017-and-081020:1:1:1715:1759 16 PHIX174 322 25 36M * 0 0 GATATTTTAAAGGAGCGTGGATTACTATCTGAGTCC B&&I13A$G$*%$IIIIIII9(.+5$IIIIIII#II XT:A:U NM:i:2 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:2C8A24 diff -r 1a56a541f9f9 -r 14d54c9dbd59 test-data/bwa_wrapper_out0b.sam --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/bwa_wrapper_out0b.sam Wed Sep 16 11:14:39 2009 -0400 @@ -0,0 +1,145 @@ +@SQ SN:phiX LN:5386 +1831_573_1004 0 phiX 278 25 33M1S * 0 0 TTAGATATGAGTCACATTTTGTTCATGGTAGAGG ]]>!4[WPFW]\ZSV]\YNRZOWWTTU]MJQEM XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_573_1050 4 * 0 0 * * 0 0 TTTATGGGTATGGCCGCTCACAGGCCAGCGGCCT ;@@17?@=>7??@A8?==@4A?A4)&+.'&+'1, +1831_573_1067 4 * 0 0 * * 0 0 AAGTATGAGTCATTTACCGGGATGCGAAGAAAAG ++#%(',%/$,+&1#&),)&+'*'*%#$&#%('$ +1831_573_1219 16 phiX 2202 37 1S33M * 0 0 TACGTCACCTTATCAGTCCAATTTAAATTACACT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_573_1242 0 phiX 648 37 33M1S * 0 0 CATCCCGTCAACATTCAAACGGCCTGTCTCATCC ]]]]XQ]]]]]]]][]]]VIMOW]SW]]]T*1] XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_573_1333 4 * 0 0 * * 0 0 AGAATCGTTACCACACACTGCGTCTGGGTATCCG ='@#%3=.>)/34*117,,/6-4+.9742456<) +1831_573_1362 4 * 0 0 * * 0 0 CGATCTCAACCAGGTCCGCGCCCACCCCTGCCTC B@@?@@?@C@BA?@>@<@8A@?@'.8'?17:,+? +1831_573_1448 16 phiX 2563 25 1S33M * 0 0 CGAAGTAACGTAAGTCTACCTATGTAGACAGTTG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_573_1490 0 phiX 1767 25 33M1S * 0 0 CGTCGTTAGGCCAGTTTTCTGGTCGTGTTCAACC ]]]TU]]]=;]]]R9=BQQDF>QZQ1,LBIGEN XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_573_1523 4 * 0 0 * * 0 0 ATGGAACGGAACGGGTGAGGAGGGGAACGTAGGG ,//*,#8''.6''6'2&27/80)%�'%#*,2' +1831_573_1578 4 * 0 0 * * 0 0 CGAGTAGCAAACAAGACGCCAATCCAGGCGACCC ?=@<@C@?@@??=>@;7<<8<694)3511(+1<. +1831_573_1647 16 phiX 1118 37 1S33M * 0 0 ACTGAGGAAGCTCTACCTGCGGCAACCGCGAGAG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_573_1684 16 phiX 3696 37 1S33M * 0 0 ATCACGACTCCAACTGAATCAAGTAGTCGTTTGC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_573_1769 4 * 0 0 * * 0 0 TGGACGTATAGTGGCGATGAGCATGTAGGTTCTC 54749&&767%/7;3$-7;3#,3//#,45/#&06 +1831_573_1853 4 * 0 0 * * 0 0 CAAAACGCCCGGGGCCTCACATGCGCGGCAGTTC 4./B@@72B+/4?@?7?+@9/+99.')2<2&)2& +1831_573_1943 4 * 0 0 * * 0 0 ATAACGTATGGCAGTGAACGGGCGGAACCTGCCC >/*1,),@-)'0*>5'$/?6@/.&?8/(')A43, +1831_573_1977 4 * 0 0 * * 0 0 GGCGTAGGGCTCATTGTGCAAGTATCCGACCTCC @?6@8@@9A@=4633A7?9?<3:+7>'@%@?%32 +1831_574_109 0 phiX 3729 37 33M1S * 0 0 CAGAATCAGCGGTATGGCTCTTCTCATATTGGCT VTO/2URQ]]]TZWU[UTQP\]XY]VS[WLN]] XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_574_148 16 phiX 2369 37 1S33M * 0 0 GATCCAAAAGACGAATCCTCAAATTAGTACAAAG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_574_185 4 * 0 0 * * 0 0 CCGTTTTGCCTAGTAATGCTCGGCGCAGCGTCGC <4;%%@7>#*>97%6;+8$#<07%1372$%<54# +1831_574_243 0 phiX 2129 37 33M1S * 0 0 TCTGGTGATTTGCAAGAACGCGTACTTATTCGCA ]]]W]]]]RQ]YRBI]]Y0<]]]CI]]]FM]]Y XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_574_257 4 * 0 0 * * 0 0 ATACCTTCCAAAGCAATAGAATAAAAAACAGTAC >>B>)@?>B@B?@1%*%,#+00.'(+&5&%#$3& +1831_574_293 4 * 0 0 * * 0 0 TGCTGCAAATAAACATACAGCCTTCTAATGACTA ')7'.*%#&$'/.$##$.$&*)+*$#8%&%#&&% +1831_574_389 4 * 0 0 * * 0 0 CATGGCTATGCACCGGTTTGTAGCGTACTCGAGA )$>/>8):2@:213;;.1)@3%64%+)<7)+$92 +1831_574_575 16 phiX 4006 37 1S33M * 0 0 CTTACAGTGCGACTAATAAAACTGAAACTCGCAT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_574_592 16 phiX 2950 37 1S33M * 0 0 TCGCCATAACGAAGACGAGAACGACCACCGCGGT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_574_617 4 * 0 0 * * 0 0 AAGCATCGGCGGGAGCGCAAGCTGGGAAGGTGCC >@?;>?=@==@>@=5?@+@:@40@><:0)')64- +1831_574_725 4 * 0 0 * * 0 0 GACAAGATGGCTATTATTTACAATCCGATCTGCA '>?52+5=;4$6;<8.1<47*9+?7(+5;-).)4 +1831_574_734 16 phiX 5263 25 1S33M * 0 0 CCTCTCTACTCTAACTCCGACCCTTTTCAATGAC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_574_824 16 phiX 3997 37 1S33M * 0 0 ATGGGAAGACTTACAGTGCGACTAATAAAACTGA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_574_959 4 * 0 0 * * 0 0 CGCGCTAGGACTCGGCCCCAAGAGACGTTCCGCC A=;+?@9;/<@;<8#&&/*%(%.0$0$2)+1$-1 +1831_574_1062 4 * 0 0 * * 0 0 ACCGGTAATATAAGGCAACATGATTACGGCCACG @821<//:136/8<1/5.427.;>288/0<1$)/ +1831_574_1092 4 * 0 0 * * 0 0 GACTGGCGAAATCATCGCGGAAAAACCCCTATCA 8?.,98)03%#,,,(,++$'69'#9'7#';6.%2 +1831_574_1103 4 * 0 0 * * 0 0 ATCTCCTGATTAGACATATCTCCGTAGCTCACGC ??>?@=>@@A<>?=@==<:?<;58'&&):+35'1 +1831_574_1116 4 * 0 0 * * 0 0 CACCTCACGTGAGTATAGCAGCCCGAGCGTCACC @++,)/)?./)%<)2>.==:8?&;44&/&,)/*4 +1831_574_1194 0 phiX 1136 37 33M1S * 0 0 GCCGTTGGCGCTCTCCGTCTTTCTCCATTGCGTG \]PT]]]]]]]ZZ]]]\Y]]]]\]]]V\]UMDK XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_574_1204 16 phiX 1445 25 1S33M * 0 0 TGGATTACTCGAATTAGTTCTACTACGAGCAATA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_574_1306 4 * 0 0 * * 0 0 ATTGCTTAGATCCAGTGGCGCTCAATACAACGGA 2.7&%''.$&#'%##$#'#%')')$'$%$,)'.1 +1831_574_1387 4 * 0 0 * * 0 0 GTACTTCTCAATGCTGCACTACTATGTACGCCCC ?A??:@?@@9?@@?9C;:>&).,,);&'<&7(/' +1831_574_1431 4 * 0 0 * * 0 0 GACCAGTTTCAGGGCTAACCGTCCCTACTCGACC 7?>9BA897;+8*=6#6.+5&-#26$0.,5&'$1 +1831_574_1560 4 * 0 0 * * 0 0 GGCGTCTTAGGATTGAAGAGGGCCTCCCACCCCC ))&65<4=>99@:597(1&;+&(%&4#,#&&&*$ +1831_574_1591 16 phiX 5118 37 1S33M * 0 0 CTTCGACAAGTCTTAGTCTTACTCGGCGTTGAAG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_574_1624 16 phiX 4058 37 1S33M * 0 0 CATAACTCCGAACACCGTAAAGATGAGAAAGAGT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_574_1826 4 * 0 0 * * 0 0 TACGTCGCGACCGAGCGTTATATAGTCTGACCCC @<<,))459<?.+57'-'&6?1)(;1;0&%4,#$ +1831_574_1903 0 phiX 5209 37 33M1S * 0 0 AGCTGGGTTACGACGCGACGCCGTTCAACCAGAC WT]]]][]]]][]]]]]]]]]]][X]QYYTYZ- XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_574_1961 0 phiX 3555 25 33M1S * 0 0 CGCGTTGCGTCTATTATGGAAAACACCAATCTTA TFPMV]OCBP]H<OS]<1CN]]M=K]52KO]WT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_575_54 16 phiX 3587 37 1S33M * 0 0 TAAGGTTCGTTGTCGTCCAAAGGCTCTAATACGC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_80 4 * 0 0 * * 0 0 TCTTTGGGTTTGGGGCAATTTGGTATGTTCGTCT ?5><@?/=5%(&.5)+#/0'/'>#$/'37#()6& +1831_575_192 4 * 0 0 * * 0 0 AACTACGCCCCTTAATTACACAGCGCGTTAGACC :8<,9(.=/)#%8$,4'<,:)3/.'.5<,>-;7' +1831_575_197 16 phiX 4861 37 1S33M * 0 0 CCGAACGTTTTATGCACCGGAATACCAATGTCAT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_223 4 * 0 0 * * 0 0 ACGCACAAAGGAGCTCGGCGCATAGCAACGCAGA ?>?8?@:?;>+$'%&,+#$.--#%$1'&#-7$&' +1831_575_420 4 * 0 0 * * 0 0 CCCACATGGAAAACACTCACCGCCGAACAGAGCG 2?.>?(-''2+..@),$%1@;')1)?+)7.').. +1831_575_434 4 * 0 0 * * 0 0 ATCGCTGCGAGGTCACCCTGGTTACGCCCCTTCC @@><@B@>=@@?>><@?;)@?;6.?=4,'';/3) +1831_575_444 16 phiX 5098 37 1S33M * 0 0 TTCGACAGCGATGAAGGGTTCTTCGACAAGTCTT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_459 4 * 0 0 * * 0 0 CTTATCAGCACACGGTTTACCAGTCCGACTCCAA @@C@?A@@A@??A@@'6@@B'=;C@'4<=9,@1. +1831_575_506 4 * 0 0 * * 0 0 CGCAGAACCCGCACGCTTGTGCTCACCACTGTAC ?)2*5@3+.7=?A:8@,/&):@37/<89'+4*)* +1831_575_569 0 phiX 5043 25 33M1S * 0 0 GACCTTGCTGCTAAAGGTCTAGGAGCTAAAGAAT IOBU[TMHTYG4M]7$?Q]]]57]SKXQC;:FE XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_575_622 0 phiX 2485 37 33M1S * 0 0 GTTTTACAGACACCTAAAGCTACATCGTCAACGC ]]]]]]]]]]]]]]]X[YVYINYLTLJWHO[$! XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_644 4 * 0 0 * * 0 0 CTTGAATGGCGATCTCGTCGAGGAAATACCACTA <8>.>&,+)8.10':5(/-62&*.=/%(.<&&18 +1831_575_663 16 phiX 1554 25 1S33M * 0 0 TGCTGAAGATGGTGTAGATAACTGTAATACCCAG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_575_681 4 * 0 0 * * 0 0 TCTCCTGATTAGACATATCACTGTTGAACACAGC 72?.)25;96-14'*.03'#4#++$)/8+%/((& +1831_575_711 4 * 0 0 * * 0 0 TATGTTCGTCCACGTCAGACGCGCAAAGTTGCGC %+:2(#'1+.)&*45-#,<A-%4;>+))<;##*& +1831_575_730 4 * 0 0 * * 0 0 CACACAGGAACCATAGCGTATGTTATTCACCCCC %8<%-')5,*1&=#)(1*1)->7$.,4'4'&'*& +1831_575_904 4 * 0 0 * * 0 0 ACCCGCTTAAAGACGTGAATTTTGCCTCCGCGCC <><.689@@7@<54:4=@>69=:+6-=42;8<:# +1831_575_938 16 phiX 1133 37 1S33M * 0 0 TCTGCGGCAACCGCGAGAGGCAGAAAGAGGTAAC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_970 16 phiX 232 37 1S33M * 0 0 GCTCCTCTTCACCGAATTATACGAACCGTGCAAG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_991 4 * 0 0 * * 0 0 TTCGGCGATCCCCCCCACGGCGCGATGCCGCGCA >>;7???/4<:::,#&*%'$--5-#(7)&&,-)% +1831_575_1138 0 phiX 1913 37 33M1S * 0 0 GCTAAAGGTGCTTTGACTTATACCGATATTGCTC VT7-SURJS]X]TT]Z]UY]]]]]]X[UN[NMT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1157 16 phiX 3736 25 1S33M * 0 0 GTCGCCATACCGAGAAGAGTATAACCGCGATGAC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1180 0 phiX 4169 37 33M1S * 0 0 GGGCGTTGAGTTCGATAATGGTGATATGTATGTA ]]\]OJROX]]UUYSX\[Y]]]]VZ]\VRWLT] XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1283 4 * 0 0 * * 0 0 GGTGGAATACTAAGGACTAATGTGCTGTCTCTTT C@A54><?A78@<?2>:#>=?@2>80:'?69.'7 +1831_575_1302 16 phiX 5202 37 1S33M * 0 0 GGAATGGTTCGACCCAATGCTGCGCTGCGGCAAG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1310 16 phiX 1612 37 1S33M * 0 0 AATGACTGGTTCTTGCACTAATGAAGTACGTCGC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1321 0 phiX 2664 37 33M1S * 0 0 TTTGGTTCGCTTTGAGTCTTCTTCGGTTCCGACG ]]]]]]]]]]]]]]]]]]]\ND]]]7/]]]WPZ XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1373 16 phiX 4221 37 1S33M * 0 0 TAGACTGCAAGCACTACTCAAACATAGACAATGA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1419 4 * 0 0 * * 0 0 TGGGGAATATAACAGCGTAGCGTTGAACACTAGA //3->&-.->%11/@&/&7:<,-1=,+2<6'%&< +1831_575_1436 0 phiX 3460 37 33M1S * 0 0 CTGGCATTCAGTCGGCGACTTCACGCCAGAATAC []]WULFXOGY]]\PY]]IGZ]]M@V]XAAMR[ XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1442 0 phiX 3317 37 33M1S * 0 0 TGCTTGGGAGCGTGCTGGTGCTGATGCTTCCTCG \V]]]]RE\]LEQ]]]X]]]]TZ]]QUPW]LIX XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1454 4 * 0 0 * * 0 0 CCTCCTAACCACGAGCCGAGGGGTCTCTGCCCCT 4=?8=6589=0:A=:?>6<;@/,4+-6+/'.,(, +1831_575_1500 16 phiX 4269 37 1S33M * 0 0 GCTTAACCGTGTTACGATGTTACACGAGGGGGTT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1535 16 phiX 1225 37 1S33M * 0 0 GAGTGCAAATACCACTTGTCACCTAATTCAAGTA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1724 4 * 0 0 * * 0 0 TCGTAAGTGTTAAGGAGCTGTGTACAAAACAACA <B.+@),2*.%)))4%2@;7#%(%+$8))85%&5 +1831_575_1829 4 * 0 0 * * 0 0 CATTTGCTGACCCTGCGTAGTTTAGTCTCACAGC 2(*'%.3$+7)@&%$'3*+*#/#/*+0.=&#)+0 +1831_575_1898 4 * 0 0 * * 0 0 CTTACCATATCATCTCAACCCATAACAGCACTTA B>@@&1/))'40)%#8/.%#8$((#;4'$'63,, +1831_575_1964 4 * 0 0 * * 0 0 GACAGACCATGAGGCTGAAGACTACTAAGTGCGT >39)@2<2/@+9?2=&)>>@*62=5&2<42.'?+ +1831_576_32 16 phiX 988 25 1S33M * 0 0 TAACTCCTATTTAATACAGATTATAAGTTTGACC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_576_74 4 * 0 0 * * 0 0 ACATTCTGCAGTGGGACAGAGCGGTACGCCGCAA =@@?@>(??<B5?@@@9<@><+><;@';>6961? +1831_576_86 16 phiX 2981 37 1S33M * 0 0 TGTACAGATTTAACAAACCTCCGCCAGTTTTTCG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_89 4 * 0 0 * * 0 0 GCTGTTTGATTTGAGAAGAGGAATTAAGCGCCGA )69+,'.4=-,>/>(*$#)3030*'(,%)2##$, +1831_576_266 16 phiX 280 37 1S33M * 0 0 ATCTATACTCAGTGTAAAACAAGTACCATCTCTA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_576_327 16 phiX 4517 25 1S33M * 0 0 TATGATAAGTCGCAAACTACTTACGTTACGCTGT XT:A:U CM:i:1 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:1G31 +1831_576_331 4 * 0 0 * * 0 0 GACGCTTTACTCCGGTAGTACCGTGCCGTTTATA -2)+(*.*1/;5%.-9&#/1'+($*$##()%/$. +1831_576_387 4 * 0 0 * * 0 0 ACACGCCATGATCCGATAAGAAGGGAACGTAAGG 8).$5#1#*%.$##.*#$%##-%,+,1#&%.%)) +1831_576_406 4 * 0 0 * * 0 0 AGGTCTTACAGCACGGGGCTGAGCGCATCTGACC 4<>%%14:*4656)&<251&2+3#&19,6&4>5( +1831_576_449 4 * 0 0 * * 0 0 CTCGAACCGCGGGGTCCAAAGACTGCTGCAAGGA +'&'.,,$/+.)$$8&%#+?&,#)-&###7,+#* +1831_576_519 16 phiX 1625 25 1S33M * 0 0 GGCACTAATGAAGTACGTCGCAATGGCACTACAA XT:A:U CM:i:1 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:25A7 +1831_576_603 16 phiX 1576 37 1S33M * 0 0 GGTAATACCCAGACGTTCGACGAATACGATTAAA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_655 16 phiX 3292 37 1S33M * 0 0 CAACGACGACGTAAAGGACTCGAATTACGAACCC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_576_677 4 * 0 0 * * 0 0 TTTACTCAGTTGATACATCACTGTAGCACATAGG 290&/*0#&'&,.2'#&*$&('#-%($*#%$)#% +1831_576_718 4 * 0 0 * * 0 0 CGTGCCTTTCAGGGTCTTTTCTGGTCTGGTCGCT 4###$0###%(#########,####%####$### +1831_576_722 4 * 0 0 * * 0 0 CGTATGATGGCGAGTCTTTATAATCCAATCTGAA *,-,##.)*&(*1%*(%(-2#+)-#.&-#%%$') +1831_576_754 4 * 0 0 * * 0 0 AGGCGTCCTGCATCGACCGTTCTATCCGCTTAGA 4//#(.$)'',>(
$<,##%((,#5?#0*%1* +1831_576_815 0 phiX 5186 37 33M1S * 0 0 GGAGTGCTTAATCCAACTTACCAAGCTGGGTTAC Z]]]]]]]]]]]]]]]]]]]UP[ZKIP]][ZUN XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_882 0 phiX 5053 37 33M1S * 0 0 CTAAAGGTCTAGGAGCTAAAGAATGGAACAACTG Q]]]]]W[]]]X]]]]][EAV]]ZUS\]]]VX] XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_898 0 phiX 518 37 33M1S * 0 0 GATTGCTACTGACCGCTCTCGTGCTCGTCGCTGT ]]]Z]]]]]]]YW]]]]VV]]X]]]]]VO]]]T XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_923 0 phiX 1050 37 33M1S * 0 0 ATCTTGGCTTCCTTGCTGGTCAGATTGGTCGTCG MLOW]RMX[]]X]]]XV]]]UKU]]]W\\]\ST XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_930 4 * 0 0 * * 0 0 CTGGCATGTACGTCCATGTCAGACGAGCAGAACT #*0.2&.:((#'14'##-)#%$$2%#$/1&#%/# +1831_576_1019 4 * 0 0 * * 0 0 GATGCGCGCTGTCATGGCATCGAACCATCGGTCG <$>$,1,&++&@,.)'+/+#9'69/6'2(+-'9- +1831_576_1068 4 * 0 0 * * 0 0 AAGAGTGACTCACTTACCGGGATGCGATGGAGCC >><=4?;@7=??9?;9>5@9?:8@:=5';7;'#4 +1831_576_1131 16 phiX 4061 37 1S33M * 0 0 TACTCCGAACACCGTAAAGATGAGAAAGAGTTAG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_576_1168 4 * 0 0 * * 0 0 CACTTAACTCCGCTGTCGGAAGCCTTACAAGACA @CC6@@2?9=>7?;76<;467@;9,0%26'',4$ +1831_576_1207 4 * 0 0 * * 0 0 CAACCTGACTAAACGGGGATACGCTGGCGCTACA $>2,(&?4?(,@:<&,@>?$&:8A%%=0.%,597 +1831_576_1289 16 phiX 1899 37 1S33M * 0 0 AAAGTCATGGAATTGCGATTTCCACGAAACTGAA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_1329 4 * 0 0 * * 0 0 CCAATAGCAGAGACCTAATAGGAAACGAAAGCAA >&0,870/A@@/5.;=;:'&@3'&$$%,+#($&2 +1831_576_1367 0 phiX 1839 37 33M1S * 0 0 AGCATGGCACTATGTTTACTCTTGCGCTTGTTCT ]]]]]]]]]]]]]]]]]]\]]]]]]]\YW]]XX XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_1416 0 phiX 4522 37 33M1S * 0 0 TTCAGCGTTTGATGAATGCAATGCGACAGGCTCC K]VEHHS]P>TQJML]]\UFU][L6K]\TMPZ] XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_1461 4 * 0 0 * * 0 0 GAGGGGCGGCCCGGTTCAAGCAGGTAAGCAACAA 1&74#(#(&##$#,'###%#%%#*#&%%##%%%' +1831_576_1605 0 phiX 5209 37 33M1S * 0 0 AGCTGGGTTACGACGCGACGCCGTTCAACCAGAC ]]]]]XU]]]]]]]]]]]]]]X[]KHLWZUXW, XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_576_1664 4 * 0 0 * * 0 0 CGCGCACAACTCGCCATGATACGACAAGACCCGA @?>?@=A:@@>@='=<>=).>=9*8)(7#/++;) +1831_576_1671 4 * 0 0 * * 0 0 GTCTTTGTAAGCCTGGCCTGGTCAGGTCTGGTCT @?@5>@@8&8?25#&)&5&&)15&&,&4%&,&6% +1831_576_1729 16 phiX 5194 37 1S33M * 0 0 TATTAGGTTGAATGGTTCGACCCAATGCTGCGCT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_1880 4 * 0 0 * * 0 0 TATGCGCTGTTGAGCTTACAACTCACTAGCGAAT 515/792,:,7/%/05,%$):+#8%2(1754))3 +1831_576_1982 16 phiX 5024 37 1S33M * 0 0 TATTGTTTTTCAGTCTATACCTGGAACGACGATT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_1987 4 * 0 0 * * 0 0 AAGGTCTTCTGTCGGCGCTGGACTGAACACCTGA ##()(#&##)%)/-#%$11#%3>'##&$,#$$'+ +1831_576_2014 4 * 0 0 * * 0 0 CCGTGACACACAATGCCGGCCCCAGCCTAGCAAT @@9,B29-5>'?,+?79+/A';'2@'5&/9,6&1 +1831_576_2028 4 * 0 0 * * 0 0 ACTCGCCGCATCCCCGAGTGACGCTCGAGACCAA <>9>5'2(&707.8#&&39,0%7/#(#,*%&5*+ +1831_577_40 0 phiX 2772 25 33M1S * 0 0 TGTGACTATTGACGTCCTTCCCCGTACGCCGGGT JLE8>>9:MOR=!&74KQDKU]M?IG@P]]]PK XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_577_119 4 * 0 0 * * 0 0 TCCCACAAGCCATTGACATGCTCGCTCTAAACAG >?C11)-1/)#;#/,850*+.+$$5550+%-.40 +1831_577_133 16 phiX 1714 37 1S33M * 0 0 CGTACGCGAGATTAGAGACCCGTAGACCGATACT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_577_255 4 * 0 0 * * 0 0 ATTGAGGCCAAGATAATTGAGGAGAGAGAAGGTG <=87@>69%**#&#-+$.#&&#+$-+%&%.,%5# +1831_577_281 16 phiX 4577 37 1S33M * 0 0 TAAAACTGTGAGAGTGCAACCGACTGCTGGCTAA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_577_288 0 phiX 3201 25 33M1S * 0 0 GGCACTTCTGCCGTTTCTGATAAGTTGCTTGATA SOY]]NL]]]UTZKIMNOCH/-LFLLK]VN5%F XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_577_322 16 phiX 458 37 1S33M * 0 0 CAAGACGGCAAAACCTAAATTGGCTTCTACTAAA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_577_362 0 phiX 2153 37 33M1S * 0 0 CTTATTCGCCACCATGATTATGACCAGTGTTTCA ]]]]]]]]]WV]RHVONRSI6=O[NK[ZL?%(J XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_577_382 16 phiX 4686 37 1S33M * 0 0 TCCAGCGTTCCGATTACTAAGTGTGCGGCTGACG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_577_464 0 phiX 2126 25 33M1S * 0 0 CCTTCTGGTGATTTGCAAGAACGCGTACTTATTG ]\OT]SQHP<,>6COKMHI3-HMNHNJLS@FUT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_577_488 16 phiX 4948 37 1S33M * 0 0 GCCGGACAACTACGATTTCCACTCGGCGAATTTC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_577_511 4 * 0 0 * * 0 0 GGTGCTTATCAGTTCTTTCTCGGGACTTGTATTT 7/+&$:<7%6,$$%'%/+)#$7((&*3>16'0/+ +1831_577_545 0 phiX 1711 37 33M1S * 0 0 TGTCATGCGCTCTAATCTCTGGGCATCTGGCTAT ]]]]]]]]]YURZ]]]]]]]VNNPV]WLLWZYN XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_577_559 0 phiX 3840 25 33M1S * 0 0 GCTGTTGCCGATACTTGGAACAATTTCTGGAAAG ]]]]]]]]WJWSWZW]XLX\P#-]]V=A(/TD9 XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_577_562 4 * 0 0 * * 0 0 GTTCCACTACGTTCCACGCAAAGGAATCCGAATC ((&*&0%)1%)#($2-,***%/-,,))&,-.1'* +1831_577_637 0 phiX 4247 37 33M1S * 0 0 TGTTACTGAGAAGTTAATGGATGAATTGGCACAA BI2<]PLJN[M;7FYWYVW]]\WG@IQH6O]XI XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_577_641 4 * 0 0 * * 0 0 TATCTACCACCGCGGTGGCGCGAGAATGCTCCCT 47;/':A;;5?:72,(=),#*?+.#&7$8#%7/' +1831_577_692 0 phiX 3905 37 33M1S * 0 0 GTCTAGGAAATAACCGTCAGGATTGACACCCTCA WU]]]SQ]]]TT]]]TS]]]PM]]]]]]]]]UZ XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 diff -r 1a56a541f9f9 -r 14d54c9dbd59 test-data/bwa_wrapper_out1.sam --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/bwa_wrapper_out1.sam Wed Sep 16 11:14:39 2009 -0400 @@ -0,0 +1,145 @@ +@SQ SN:PHIX174 LN:5386 +1831_573_1004 0 PHIX174 278 25 33M1S * 0 0 TTAGATATGAGTCACATTTTGTTCATGGTAGAGG ]]>!4[WPFW]\ZSV]\YNRZOWWTTU]MJQEM XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_573_1050 4 * 0 0 * * 0 0 TTTATGGGTATGGCCGCTCACAGGCCAGCGGCCT ;@@17?@=>7??@A8?==@4A?A4)&+.'&+'1, +1831_573_1067 4 * 0 0 * * 0 0 AAGTATGAGTCATTTACCGGGATGCGAAGAAAAG ++#%(',%/$,+&1#&),)&+'*'*%#$&#%('$ +1831_573_1219 16 PHIX174 2202 37 1S33M * 0 0 TACGTCACCTTATCAGTCCAATTTAAATTACACT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_573_1242 0 PHIX174 648 37 33M1S * 0 0 CATCCCGTCAACATTCAAACGGCCTGTCTCATCC ]]]]XQ]]]]]]]][]]]VIMOW]SW]]]T*1] XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_573_1333 4 * 0 0 * * 0 0 AGAATCGTTACCACACACTGCGTCTGGGTATCCG ='@#%3=.>)/34*117,,/6-4+.9742456<) +1831_573_1362 4 * 0 0 * * 0 0 CGATCTCAACCAGGTCCGCGCCCACCCCTGCCTC B@@?@@?@C@BA?@>@<@8A@?@'.8'?17:,+? +1831_573_1448 16 PHIX174 2563 25 1S33M * 0 0 CGAAGTAACGTAAGTCTACCTATGTAGACAGTTG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_573_1490 0 PHIX174 1767 25 33M1S * 0 0 CGTCGTTAGGCCAGTTTTCTGGTCGTGTTCAACC ]]]TU]]]=;]]]R9=BQQDF>QZQ1,LBIGEN XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_573_1523 4 * 0 0 * * 0 0 ATGGAACGGAACGGGTGAGGAGGGGAACGTAGGG ,//*,#8''.6''6'2&27/80)%�'%#*,2' +1831_573_1578 4 * 0 0 * * 0 0 CGAGTAGCAAACAAGACGCCAATCCAGGCGACCC ?=@<@C@?@@??=>@;7<<8<694)3511(+1<. +1831_573_1647 16 PHIX174 1118 37 1S33M * 0 0 ACTGAGGAAGCTCTACCTGCGGCAACCGCGAGAG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_573_1684 16 PHIX174 3696 37 1S33M * 0 0 ATCACGACTCCAACTGAATCAAGTAGTCGTTTGC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_573_1769 4 * 0 0 * * 0 0 TGGACGTATAGTGGCGATGAGCATGTAGGTTCTC 54749&&767%/7;3$-7;3#,3//#,45/#&06 +1831_573_1853 4 * 0 0 * * 0 0 CAAAACGCCCGGGGCCTCACATGCGCGGCAGTTC 4./B@@72B+/4?@?7?+@9/+99.')2<2&)2& +1831_573_1943 4 * 0 0 * * 0 0 ATAACGTATGGCAGTGAACGGGCGGAACCTGCCC >/*1,),@-)'0*>5'$/?6@/.&?8/(')A43, +1831_573_1977 4 * 0 0 * * 0 0 GGCGTAGGGCTCATTGTGCAAGTATCCGACCTCC @?6@8@@9A@=4633A7?9?<3:+7>'@%@?%32 +1831_574_109 0 PHIX174 3729 37 33M1S * 0 0 CAGAATCAGCGGTATGGCTCTTCTCATATTGGCT VTO/2URQ]]]TZWU[UTQP\]XY]VS[WLN]] XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_574_148 16 PHIX174 2369 37 1S33M * 0 0 GATCCAAAAGACGAATCCTCAAATTAGTACAAAG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_574_185 4 * 0 0 * * 0 0 CCGTTTTGCCTAGTAATGCTCGGCGCAGCGTCGC <4;%%@7>#*>97%6;+8$#<07%1372$%<54# +1831_574_243 0 PHIX174 2129 37 33M1S * 0 0 TCTGGTGATTTGCAAGAACGCGTACTTATTCGCA ]]]W]]]]RQ]YRBI]]Y0<]]]CI]]]FM]]Y XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_574_257 4 * 0 0 * * 0 0 ATACCTTCCAAAGCAATAGAATAAAAAACAGTAC >>B>)@?>B@B?@1%*%,#+00.'(+&5&%#$3& +1831_574_293 4 * 0 0 * * 0 0 TGCTGCAAATAAACATACAGCCTTCTAATGACTA ')7'.*%#&$'/.$##$.$&*)+*$#8%&%#&&% +1831_574_389 4 * 0 0 * * 0 0 CATGGCTATGCACCGGTTTGTAGCGTACTCGAGA )$>/>8):2@:213;;.1)@3%64%+)<7)+$92 +1831_574_575 16 PHIX174 4006 37 1S33M * 0 0 CTTACAGTGCGACTAATAAAACTGAAACTCGCAT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_574_592 16 PHIX174 2950 37 1S33M * 0 0 TCGCCATAACGAAGACGAGAACGACCACCGCGGT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_574_617 4 * 0 0 * * 0 0 AAGCATCGGCGGGAGCGCAAGCTGGGAAGGTGCC >@?;>?=@==@>@=5?@+@:@40@><:0)')64- +1831_574_725 4 * 0 0 * * 0 0 GACAAGATGGCTATTATTTACAATCCGATCTGCA '>?52+5=;4$6;<8.1<47*9+?7(+5;-).)4 +1831_574_734 16 PHIX174 5263 25 1S33M * 0 0 CCTCTCTACTCTAACTCCGACCCTTTTCAATGAC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_574_824 16 PHIX174 3997 37 1S33M * 0 0 ATGGGAAGACTTACAGTGCGACTAATAAAACTGA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_574_959 4 * 0 0 * * 0 0 CGCGCTAGGACTCGGCCCCAAGAGACGTTCCGCC A=;+?@9;/<@;<8#&&/*%(%.0$0$2)+1$-1 +1831_574_1062 4 * 0 0 * * 0 0 ACCGGTAATATAAGGCAACATGATTACGGCCACG @821<//:136/8<1/5.427.;>288/0<1$)/ +1831_574_1092 4 * 0 0 * * 0 0 GACTGGCGAAATCATCGCGGAAAAACCCCTATCA 8?.,98)03%#,,,(,++$'69'#9'7#';6.%2 +1831_574_1103 4 * 0 0 * * 0 0 ATCTCCTGATTAGACATATCTCCGTAGCTCACGC ??>?@=>@@A<>?=@==<:?<;58'&&):+35'1 +1831_574_1116 4 * 0 0 * * 0 0 CACCTCACGTGAGTATAGCAGCCCGAGCGTCACC @++,)/)?./)%<)2>.==:8?&;44&/&,)/*4 +1831_574_1194 0 PHIX174 1136 37 33M1S * 0 0 GCCGTTGGCGCTCTCCGTCTTTCTCCATTGCGTG \]PT]]]]]]]ZZ]]]\Y]]]]\]]]V\]UMDK XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_574_1204 16 PHIX174 1445 25 1S33M * 0 0 TGGATTACTCGAATTAGTTCTACTACGAGCAATA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_574_1306 4 * 0 0 * * 0 0 ATTGCTTAGATCCAGTGGCGCTCAATACAACGGA 2.7&%''.$&#'%##$#'#%')')$'$%$,)'.1 +1831_574_1387 4 * 0 0 * * 0 0 GTACTTCTCAATGCTGCACTACTATGTACGCCCC ?A??:@?@@9?@@?9C;:>&).,,);&'<&7(/' +1831_574_1431 4 * 0 0 * * 0 0 GACCAGTTTCAGGGCTAACCGTCCCTACTCGACC 7?>9BA897;+8*=6#6.+5&-#26$0.,5&'$1 +1831_574_1560 4 * 0 0 * * 0 0 GGCGTCTTAGGATTGAAGAGGGCCTCCCACCCCC ))&65<4=>99@:597(1&;+&(%&4#,#&&&*$ +1831_574_1591 16 PHIX174 5118 37 1S33M * 0 0 CTTCGACAAGTCTTAGTCTTACTCGGCGTTGAAG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_574_1624 16 PHIX174 4058 37 1S33M * 0 0 CATAACTCCGAACACCGTAAAGATGAGAAAGAGT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_574_1826 4 * 0 0 * * 0 0 TACGTCGCGACCGAGCGTTATATAGTCTGACCCC @<<,))459<?.+57'-'&6?1)(;1;0&%4,#$ +1831_574_1903 0 PHIX174 5209 37 33M1S * 0 0 AGCTGGGTTACGACGCGACGCCGTTCAACCAGAC WT]]]][]]]][]]]]]]]]]]][X]QYYTYZ- XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_574_1961 0 PHIX174 3555 25 33M1S * 0 0 CGCGTTGCGTCTATTATGGAAAACACCAATCTTA TFPMV]OCBP]H<OS]<1CN]]M=K]52KO]WT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_575_54 16 PHIX174 3587 37 1S33M * 0 0 TAAGGTTCGTTGTCGTCCAAAGGCTCTAATACGC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_80 4 * 0 0 * * 0 0 TCTTTGGGTTTGGGGCAATTTGGTATGTTCGTCT ?5><@?/=5%(&.5)+#/0'/'>#$/'37#()6& +1831_575_192 4 * 0 0 * * 0 0 AACTACGCCCCTTAATTACACAGCGCGTTAGACC :8<,9(.=/)#%8$,4'<,:)3/.'.5<,>-;7' +1831_575_197 16 PHIX174 4861 37 1S33M * 0 0 CCGAACGTTTTATGCACCGGAATACCAATGTCAT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_223 4 * 0 0 * * 0 0 ACGCACAAAGGAGCTCGGCGCATAGCAACGCAGA ?>?8?@:?;>+$'%&,+#$.--#%$1'&#-7$&' +1831_575_420 4 * 0 0 * * 0 0 CCCACATGGAAAACACTCACCGCCGAACAGAGCG 2?.>?(-''2+..@),$%1@;')1)?+)7.').. +1831_575_434 4 * 0 0 * * 0 0 ATCGCTGCGAGGTCACCCTGGTTACGCCCCTTCC @@><@B@>=@@?>><@?;)@?;6.?=4,'';/3) +1831_575_444 16 PHIX174 5098 37 1S33M * 0 0 TTCGACAGCGATGAAGGGTTCTTCGACAAGTCTT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_459 4 * 0 0 * * 0 0 CTTATCAGCACACGGTTTACCAGTCCGACTCCAA @@C@?A@@A@??A@@'6@@B'=;C@'4<=9,@1. +1831_575_506 4 * 0 0 * * 0 0 CGCAGAACCCGCACGCTTGTGCTCACCACTGTAC ?)2*5@3+.7=?A:8@,/&):@37/<89'+4*)* +1831_575_569 0 PHIX174 5043 25 33M1S * 0 0 GACCTTGCTGCTAAAGGTCTAGGAGCTAAAGAAT IOBU[TMHTYG4M]7$?Q]]]57]SKXQC;:FE XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_575_622 0 PHIX174 2485 37 33M1S * 0 0 GTTTTACAGACACCTAAAGCTACATCGTCAACGC ]]]]]]]]]]]]]]]X[YVYINYLTLJWHO[$! XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_644 4 * 0 0 * * 0 0 CTTGAATGGCGATCTCGTCGAGGAAATACCACTA <8>.>&,+)8.10':5(/-62&*.=/%(.<&&18 +1831_575_663 16 PHIX174 1554 25 1S33M * 0 0 TGCTGAAGATGGTGTAGATAACTGTAATACCCAG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_575_681 4 * 0 0 * * 0 0 TCTCCTGATTAGACATATCACTGTTGAACACAGC 72?.)25;96-14'*.03'#4#++$)/8+%/((& +1831_575_711 4 * 0 0 * * 0 0 TATGTTCGTCCACGTCAGACGCGCAAAGTTGCGC %+:2(#'1+.)&*45-#,<A-%4;>+))<;##*& +1831_575_730 4 * 0 0 * * 0 0 CACACAGGAACCATAGCGTATGTTATTCACCCCC %8<%-')5,*1&=#)(1*1)->7$.,4'4'&'*& +1831_575_904 4 * 0 0 * * 0 0 ACCCGCTTAAAGACGTGAATTTTGCCTCCGCGCC <><.689@@7@<54:4=@>69=:+6-=42;8<:# +1831_575_938 16 PHIX174 1133 37 1S33M * 0 0 TCTGCGGCAACCGCGAGAGGCAGAAAGAGGTAAC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_970 16 PHIX174 232 37 1S33M * 0 0 GCTCCTCTTCACCGAATTATACGAACCGTGCAAG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_991 4 * 0 0 * * 0 0 TTCGGCGATCCCCCCCACGGCGCGATGCCGCGCA >>;7???/4<:::,#&*%'$--5-#(7)&&,-)% +1831_575_1138 0 PHIX174 1913 37 33M1S * 0 0 GCTAAAGGTGCTTTGACTTATACCGATATTGCTC VT7-SURJS]X]TT]Z]UY]]]]]]X[UN[NMT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1157 16 PHIX174 3736 25 1S33M * 0 0 GTCGCCATACCGAGAAGAGTATAACCGCGATGAC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1180 0 PHIX174 4169 37 33M1S * 0 0 GGGCGTTGAGTTCGATAATGGTGATATGTATGTA ]]\]OJROX]]UUYSX\[Y]]]]VZ]\VRWLT] XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1283 4 * 0 0 * * 0 0 GGTGGAATACTAAGGACTAATGTGCTGTCTCTTT C@A54><?A78@<?2>:#>=?@2>80:'?69.'7 +1831_575_1302 16 PHIX174 5202 37 1S33M * 0 0 GGAATGGTTCGACCCAATGCTGCGCTGCGGCAAG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1310 16 PHIX174 1612 37 1S33M * 0 0 AATGACTGGTTCTTGCACTAATGAAGTACGTCGC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1321 0 PHIX174 2664 37 33M1S * 0 0 TTTGGTTCGCTTTGAGTCTTCTTCGGTTCCGACG ]]]]]]]]]]]]]]]]]]]\ND]]]7/]]]WPZ XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1373 16 PHIX174 4221 37 1S33M * 0 0 TAGACTGCAAGCACTACTCAAACATAGACAATGA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1419 4 * 0 0 * * 0 0 TGGGGAATATAACAGCGTAGCGTTGAACACTAGA //3->&-.->%11/@&/&7:<,-1=,+2<6'%&< +1831_575_1436 0 PHIX174 3460 37 33M1S * 0 0 CTGGCATTCAGTCGGCGACTTCACGCCAGAATAC []]WULFXOGY]]\PY]]IGZ]]M@V]XAAMR[ XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1442 0 PHIX174 3317 37 33M1S * 0 0 TGCTTGGGAGCGTGCTGGTGCTGATGCTTCCTCG \V]]]]RE\]LEQ]]]X]]]]TZ]]QUPW]LIX XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1454 4 * 0 0 * * 0 0 CCTCCTAACCACGAGCCGAGGGGTCTCTGCCCCT 4=?8=6589=0:A=:?>6<;@/,4+-6+/'.,(, +1831_575_1500 16 PHIX174 4269 37 1S33M * 0 0 GCTTAACCGTGTTACGATGTTACACGAGGGGGTT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1535 16 PHIX174 1225 37 1S33M * 0 0 GAGTGCAAATACCACTTGTCACCTAATTCAAGTA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_575_1724 4 * 0 0 * * 0 0 TCGTAAGTGTTAAGGAGCTGTGTACAAAACAACA <B.+@),2*.%)))4%2@;7#%(%+$8))85%&5 +1831_575_1829 4 * 0 0 * * 0 0 CATTTGCTGACCCTGCGTAGTTTAGTCTCACAGC 2(*'%.3$+7)@&%$'3*+*#/#/*+0.=&#)+0 +1831_575_1898 4 * 0 0 * * 0 0 CTTACCATATCATCTCAACCCATAACAGCACTTA B>@@&1/))'40)%#8/.%#8$((#;4'$'63,, +1831_575_1964 4 * 0 0 * * 0 0 GACAGACCATGAGGCTGAAGACTACTAAGTGCGT >39)@2<2/@+9?2=&)>>@*62=5&2<42.'?+ +1831_576_32 16 PHIX174 988 25 1S33M * 0 0 TAACTCCTATTTAATACAGATTATAAGTTTGACC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_576_74 4 * 0 0 * * 0 0 ACATTCTGCAGTGGGACAGAGCGGTACGCCGCAA =@@?@>(??<B5?@@@9<@><+><;@';>6961? +1831_576_86 16 PHIX174 2981 37 1S33M * 0 0 TGTACAGATTTAACAAACCTCCGCCAGTTTTTCG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_89 4 * 0 0 * * 0 0 GCTGTTTGATTTGAGAAGAGGAATTAAGCGCCGA )69+,'.4=-,>/>(*$#)3030*'(,%)2##$, +1831_576_266 16 PHIX174 280 37 1S33M * 0 0 ATCTATACTCAGTGTAAAACAAGTACCATCTCTA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_576_327 16 PHIX174 4517 37 1S33M * 0 0 TATGATAAGTCGCAAACTACTTACGTTACGCTGT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_331 4 * 0 0 * * 0 0 GACGCTTTACTCCGGTAGTACCGTGCCGTTTATA -2)+(*.*1/;5%.-9&#/1'+($*$##()%/$. +1831_576_387 4 * 0 0 * * 0 0 ACACGCCATGATCCGATAAGAAGGGAACGTAAGG 8).$5#1#*%.$##.*#$%##-%,+,1#&%.%)) +1831_576_406 4 * 0 0 * * 0 0 AGGTCTTACAGCACGGGGCTGAGCGCATCTGACC 4<>%%14:*4656)&<251&2+3#&19,6&4>5( +1831_576_449 4 * 0 0 * * 0 0 CTCGAACCGCGGGGTCCAAAGACTGCTGCAAGGA +'&'.,,$/+.)$$8&%#+?&,#)-&###7,+#* +1831_576_519 16 PHIX174 1625 37 1S33M * 0 0 GGCACTAATGAAGTACGTCGCAATGGCACTACAA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_603 16 PHIX174 1576 37 1S33M * 0 0 GGTAATACCCAGACGTTCGACGAATACGATTAAA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_655 16 PHIX174 3292 37 1S33M * 0 0 CAACGACGACGTAAAGGACTCGAATTACGAACCC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_576_677 4 * 0 0 * * 0 0 TTTACTCAGTTGATACATCACTGTAGCACATAGG 290&/*0#&'&,.2'#&*$&('#-%($*#%$)#% +1831_576_718 4 * 0 0 * * 0 0 CGTGCCTTTCAGGGTCTTTTCTGGTCTGGTCGCT 4###$0###%(#########,####%####$### +1831_576_722 4 * 0 0 * * 0 0 CGTATGATGGCGAGTCTTTATAATCCAATCTGAA *,-,##.)*&(*1%*(%(-2#+)-#.&-#%%$') +1831_576_754 4 * 0 0 * * 0 0 AGGCGTCCTGCATCGACCGTTCTATCCGCTTAGA 4//#(.$)'',>(
$<,##%((,#5?#0*%1* +1831_576_815 0 PHIX174 5186 37 33M1S * 0 0 GGAGTGCTTAATCCAACTTACCAAGCTGGGTTAC Z]]]]]]]]]]]]]]]]]]]UP[ZKIP]][ZUN XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_882 0 PHIX174 5053 37 33M1S * 0 0 CTAAAGGTCTAGGAGCTAAAGAATGGAACAACTG Q]]]]]W[]]]X]]]]][EAV]]ZUS\]]]VX] XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_898 0 PHIX174 518 37 33M1S * 0 0 GATTGCTACTGACCGCTCTCGTGCTCGTCGCTGT ]]]Z]]]]]]]YW]]]]VV]]X]]]]]VO]]]T XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_923 0 PHIX174 1050 37 33M1S * 0 0 ATCTTGGCTTCCTTGCTGGTCAGATTGGTCGTCG MLOW]RMX[]]X]]]XV]]]UKU]]]W\\]\ST XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_930 4 * 0 0 * * 0 0 CTGGCATGTACGTCCATGTCAGACGAGCAGAACT #*0.2&.:((#'14'##-)#%$$2%#$/1&#%/# +1831_576_1019 4 * 0 0 * * 0 0 GATGCGCGCTGTCATGGCATCGAACCATCGGTCG <$>$,1,&++&@,.)'+/+#9'69/6'2(+-'9- +1831_576_1068 4 * 0 0 * * 0 0 AAGAGTGACTCACTTACCGGGATGCGATGGAGCC >><=4?;@7=??9?;9>5@9?:8@:=5';7;'#4 +1831_576_1131 16 PHIX174 4061 37 1S33M * 0 0 TACTCCGAACACCGTAAAGATGAGAAAGAGTTAG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_576_1168 4 * 0 0 * * 0 0 CACTTAACTCCGCTGTCGGAAGCCTTACAAGACA @CC6@@2?9=>7?;76<;467@;9,0%26'',4$ +1831_576_1207 4 * 0 0 * * 0 0 CAACCTGACTAAACGGGGATACGCTGGCGCTACA $>2,(&?4?(,@:<&,@>?$&:8A%%=0.%,597 +1831_576_1289 16 PHIX174 1899 37 1S33M * 0 0 AAAGTCATGGAATTGCGATTTCCACGAAACTGAA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_1329 4 * 0 0 * * 0 0 CCAATAGCAGAGACCTAATAGGAAACGAAAGCAA >&0,870/A@@/5.;=;:'&@3'&$$%,+#($&2 +1831_576_1367 0 PHIX174 1839 37 33M1S * 0 0 AGCATGGCACTATGTTTACTCTTGCGCTTGTTCT ]]]]]]]]]]]]]]]]]]\]]]]]]]\YW]]XX XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_1416 0 PHIX174 4522 37 33M1S * 0 0 TTCAGCGTTTGATGAATGCAATGCGACAGGCTCC K]VEHHS]P>TQJML]]\UFU][L6K]\TMPZ] XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_1461 4 * 0 0 * * 0 0 GAGGGGCGGCCCGGTTCAAGCAGGTAAGCAACAA 1&74#(#(&##$#,'###%#%%#*#&%%##%%%' +1831_576_1605 0 PHIX174 5209 37 33M1S * 0 0 AGCTGGGTTACGACGCGACGCCGTTCAACCAGAC ]]]]]XU]]]]]]]]]]]]]]X[]KHLWZUXW, XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_576_1664 4 * 0 0 * * 0 0 CGCGCACAACTCGCCATGATACGACAAGACCCGA @?>?@=A:@@>@='=<>=).>=9*8)(7#/++;) +1831_576_1671 4 * 0 0 * * 0 0 GTCTTTGTAAGCCTGGCCTGGTCAGGTCTGGTCT @?@5>@@8&8?25#&)&5&&)15&&,&4%&,&6% +1831_576_1729 16 PHIX174 5194 37 1S33M * 0 0 TATTAGGTTGAATGGTTCGACCCAATGCTGCGCT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_1880 4 * 0 0 * * 0 0 TATGCGCTGTTGAGCTTACAACTCACTAGCGAAT 515/792,:,7/%/05,%$):+#8%2(1754))3 +1831_576_1982 16 PHIX174 5024 37 1S33M * 0 0 TATTGTTTTTCAGTCTATACCTGGAACGACGATT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_576_1987 4 * 0 0 * * 0 0 AAGGTCTTCTGTCGGCGCTGGACTGAACACCTGA ##()(#&##)%)/-#%$11#%3>'##&$,#$$'+ +1831_576_2014 4 * 0 0 * * 0 0 CCGTGACACACAATGCCGGCCCCAGCCTAGCAAT @@9,B29-5>'?,+?79+/A';'2@'5&/9,6&1 +1831_576_2028 4 * 0 0 * * 0 0 ACTCGCCGCATCCCCGAGTGACGCTCGAGACCAA <>9>5'2(&707.8#&&39,0%7/#(#,*%&5*+ +1831_577_40 0 PHIX174 2772 25 33M1S * 0 0 TGTGACTATTGACGTCCTTCCCCGTACGCCGGGT JLE8>>9:MOR=!&74KQDKU]M?IG@P]]]PK XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_577_119 4 * 0 0 * * 0 0 TCCCACAAGCCATTGACATGCTCGCTCTAAACAG >?C11)-1/)#;#/,850*+.+$$5550+%-.40 +1831_577_133 16 PHIX174 1714 37 1S33M * 0 0 CGTACGCGAGATTAGAGACCCGTAGACCGATACT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_577_255 4 * 0 0 * * 0 0 ATTGAGGCCAAGATAATTGAGGAGAGAGAAGGTG <=87@>69%**#&#-+$.#&&#+$-+%&%.,%5# +1831_577_281 16 PHIX174 4577 37 1S33M * 0 0 TAAAACTGTGAGAGTGCAACCGACTGCTGGCTAA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_577_288 0 PHIX174 3201 25 33M1S * 0 0 GGCACTTCTGCCGTTTCTGATAAGTTGCTTGATA SOY]]NL]]]UTZKIMNOCH/-LFLLK]VN5%F XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_577_322 16 PHIX174 458 37 1S33M * 0 0 CAAGACGGCAAAACCTAAATTGGCTTCTACTAAA XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_577_362 0 PHIX174 2153 37 33M1S * 0 0 CTTATTCGCCACCATGATTATGACCAGTGTTTCA ]]]]]]]]]WV]RHVONRSI6=O[NK[ZL?%(J XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_577_382 16 PHIX174 4686 37 1S33M * 0 0 TCCAGCGTTCCGATTACTAAGTGTGCGGCTGACG XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_577_464 0 PHIX174 2126 25 33M1S * 0 0 CCTTCTGGTGATTTGCAAGAACGCGTACTTATTG ]\OT]SQHP<,>6COKMHI3-HMNHNJLS@FUT XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_577_488 16 PHIX174 4948 37 1S33M * 0 0 GCCGGACAACTACGATTTCCACTCGGCGAATTTC XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_577_511 4 * 0 0 * * 0 0 GGTGCTTATCAGTTCTTTCTCGGGACTTGTATTT 7/+&$:<7%6,$$%'%/+)#$7((&*3>16'0/+ +1831_577_545 0 PHIX174 1711 37 33M1S * 0 0 TGTCATGCGCTCTAATCTCTGGGCATCTGGCTAT ]]]]]]]]]YURZ]]]]]]]VNNPV]WLLWZYN XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 +1831_577_559 0 PHIX174 3840 25 33M1S * 0 0 GCTGTTGCCGATACTTGGAACAATTTCTGGAAAG ]]]]]]]]WJWSWZW]XLX\P#-]]V=A(/TD9 XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:33 +1831_577_562 4 * 0 0 * * 0 0 GTTCCACTACGTTCCACGCAAAGGAATCCGAATC ((&*&0%)1%)#($2-,***%/-,,))&,-.1'* +1831_577_637 0 PHIX174 4247 37 33M1S * 0 0 TGTTACTGAGAAGTTAATGGATGAATTGGCACAA BI2<]PLJN[M;7FYWYVW]]\WG@IQH6O]XI XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0 MD:Z:33 +1831_577_641 4 * 0 0 * * 0 0 TATCTACCACCGCGGTGGCGCGAGAATGCTCCCT 47;/':A;;5?:72,(=),#*?+.#&7$8#%7/' +1831_577_692 0 PHIX174 3905 37 33M1S * 0 0 GTCTAGGAAATAACCGTCAGGATTGACACCCTCA WU]]]SQ]]]TT]]]TS]]]PM]]]]]]]]]UZ XT:A:U CM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:33 diff -r 1a56a541f9f9 -r 14d54c9dbd59 test-data/bwa_wrapper_out2.sam --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/bwa_wrapper_out2.sam Wed Sep 16 11:14:39 2009 -0400 @@ -0,0 +1,97 @@ +@SQ SN:PHIX174 LN:5386 +1378_11_329 77 * 0 0 * * 0 0 GTTCGTGGCCGGTGGGTGTTTGGG ###$$#$#$&#####$'$#$###$ +1378_11_329 141 * 0 0 * * 0 0 AGACCGGGCGGGGTGGCGTTCGGT %##+'#######%###$#$##$(# +1378_17_1788 77 * 0 0 * * 0 0 TGGGTGGATGTGTTGTCGTTCATG #$#$###$#$#######$#$#### +1378_17_1788 141 * 0 0 * * 0 0 TGCCGTGTCTTGCTAACGCCGATT #'#$$#$###%%##$$$$###### +1378_25_2035 77 * 0 0 * * 0 0 GTGCGTCGGGGAGGGTGCTGTCGG ######%#$%#$$###($###&&% +1378_25_2035 141 * 0 0 * * 0 0 CTGCGTGTTGGTGTCTACTGGGGT #%#'##$#$##&%#%$$$%#%#'# +1378_28_770 77 * 0 0 * * 0 0 GGTTGGGAGTTGGCGCGCGTGGTT ###$##%%$####%#%$#%%#($$ +1378_28_770 141 * 0 0 * * 0 0 CTGCTATGCCGGCCGCCTGCCCCT %##(#$$#%$#$%###$%%$##$# +1378_33_1945 77 * 0 0 * * 0 0 GGGGGGGGGGGGGGGGGGGGGGGG .6;?;55697>8779457986264 +1378_33_1945 141 * 0 0 * * 0 0 GGGGGGGGGGGGGGGGGGGGGGGG +9;64287788752796/576352 +1378_34_789 77 * 0 0 * * 0 0 GGGCTTGCGTTAGTGAGAGGTTGT ###%$%$%%###$####$###$#& +1378_34_789 141 * 0 0 * * 0 0 ATGGTGGCTGACGCGTTTGACTGT #$##%#$##$&$#%##$##$###$ +1378_35_263 77 * 0 0 * * 0 0 GGCGGGGGGGGGGGGGGGGGGGGG $%%&,)3/)%$-6'$&5&#$''5. +1378_35_263 141 * 0 0 * * 0 0 GGGCGGGGGGCGGCGGGGGGGGGG ),,#%'$)'-(%&%'%',%%&&(& +1378_43_186 77 * 0 0 * * 0 0 GCTAGGGTTTGGGTTTGCGGTGGG $%#$########%##%#$###'#' +1378_43_186 141 * 0 0 * * 0 0 ATACTAGTTGGGACGCGTTGTGCT #$(4%$########$#$###$$$# +1378_51_1671 77 * 0 0 * * 0 0 GGTGTTGTTCGGATGTAGCGTGGG ##%$$####$########$##0## +1378_51_1671 141 * 0 0 * * 0 0 TGGGTAGGTGCGGCTCTTGCGGTC ##$#$%###%%#%#$##$$##### +1378_56_324 77 * 0 0 * * 0 0 GCTGGGTGATCTAAGGAATGTAAT #$#####$$$#$$$$##%###*(# +1378_56_324 141 * 0 0 * * 0 0 AGGCTTGGATTGTTGCGCTGACCT /+#**###%%##%#####$$$#$# +1378_56_773 77 * 0 0 * * 0 0 TCTGGTCGGTTTCGGGGAGTGGAA ##%%#&$###$#$##%$####%%$ +1378_56_773 141 * 0 0 * * 0 0 TGTCGTGAGGTCACTTATCCCCAT &%#%##%%#####&#$%##$%##$ +1378_62_2027 77 * 0 0 * * 0 0 GTTGGCCTGGCCTGCCGTGCTGCG *##),/%##$)#%##1$#'%.#&# +1378_62_2027 141 * 0 0 * * 0 0 CTTCCACGATCTGCTCGCTGTGGT (#&&$##$$#$%#%$$$#$###'# +1378_62_2029 77 * 0 0 * * 0 0 GGCGGTGTGTGGTGCGGCTGTGCG /$$$=(####%####)$$%$-&%# +1378_62_2029 141 * 0 0 * * 0 0 TCTGGGCTGTCTTCGGGTCGGTGT $%$$####$##$$#)##%%#$### +1378_67_1795 77 * 0 0 * * 0 0 CGTGCCAACGACCGTCCTACGGAG @@.?7@@@.@:@A<<>>=>70<@7 +1378_67_1795 141 * 0 0 * * 0 0 TAGTGGGCTAGTCCGGGATCTACC B@C2@@?;-@@1>'A@@91@.@1' +1378_68_466 77 * 0 0 * * 0 0 GTGTCATCTGAGGTAAAGCATTGT /##$09#$#.=$#$76+$%1'### +1378_68_466 141 * 0 0 * * 0 0 GTGATCGTCGGTGCCAGTCCCTGT #(%)+##$#$#%#+$%##$##### +1378_68_1692 77 * 0 0 * * 0 0 TGTTCGATTTAGGTCTTGCCTCGC ###%%$%%##%$#$##$###$$#$ +1378_68_1692 141 * 0 0 * * 0 0 GGTCTTCTGGGCACGGGCCAAGTT ###)%'&%##$###%#######$# +1378_80_664 77 * 0 0 * * 0 0 TGTCTGCGTTGTATCTCTGGTGTA %##%,%$$#&$$###$#$%##'%# +1378_80_664 141 * 0 0 * * 0 0 CTGCTTTGATCCCCGGTGGAGCAC 7#%###$$6#######$##$$$## +1378_85_1786 77 * 0 0 * * 0 0 CCCTAGGAGCGTATACCGGACGAG ,'&/%/@,&1,&'/)&,6&&1)(( +1378_85_1786 141 * 0 0 * * 0 0 ATACTATGTCGATCTGTAAAAAAA )&.)#3%@$&%-,2#&+.-%0&./ +1378_86_1011 77 * 0 0 * * 0 0 AGGCGATGGGATATTATTTTACTT :$###)%##$9$###1$$#$2### +1378_86_1011 141 * 0 0 * * 0 0 CTACGTTATTGCTCTGTTTGTCCT ######$%##$$$%###%#$#### +1378_86_1789 77 * 0 0 * * 0 0 GGCCAGGTACTAATAAATTCAACA 1.51/&%#'$,*?%#6&0$$5,20 +1378_86_1789 141 * 0 0 * * 0 0 CGAATATCCCCTGGTTTAGAAGTA 61?1;@/>@<=/(.7'(<5@=%$= +1378_91_1596 77 * 0 0 * * 0 0 GCTTTTTCATTCGGTGCCTTTGGA '>%/3%=()8'#.%?50$&5>%)% +1378_91_1596 141 * 0 0 * * 0 0 TTAGCGGTTGACTATCTGCTGACA *&+'#9'(%*'#//,&<),/)'*# +1378_94_1595 77 * 0 0 * * 0 0 GCGGGTTAACCCAGTACCTCTGGC ,8@6'>8<76>@@(->9)7'*941 +1378_94_1595 141 * 0 0 * * 0 0 CGTGCGACAGCCCATGTTTTCAGA -=..5,3826&*+.+#+#%%6;%# +1378_95_1039 77 * 0 0 * * 0 0 GTTCTGTGCCAGGTGAGGTACGGA &##,./#$&)6##+,'#$$0(##$ +1378_95_1039 141 * 0 0 * * 0 0 CGGCGTCCATCTTCGCCTTGAGAT $##.#$##$$#%$#$%%$###$)$ +1378_95_1767 77 * 0 0 * * 0 0 CTAATTGACCGGGCAAGCTATTAA ',,2'(10/>:=/7).3&%6&)&6 +1378_95_1767 141 * 0 0 * * 0 0 ACACGCACCATCTGGCGGCTAACC 86<65<4<C27/$5+#;-$+#2%$ +1378_96_1037 77 * 0 0 * * 0 0 CTGCTGGGCCATTTGACTTACTCA '$#+#(##-%5##+*&###-.$$$ +1378_96_1037 141 * 0 0 * * 0 0 ATCCCCCAAGATGCCTGTTGATTG $#$'##$$$#%$$#%###+##$#$ +1378_96_1764 77 * 0 0 * * 0 0 TCTTTGTAACCCACTTAGTATTTC 4='28>.@4881B807:822>%)1 +1378_96_1764 141 * 0 0 * * 0 0 TACACCGGAACCTTGCACTAAAGC 6-9'<#+81?+*<8-=69&6;*3) +1378_98_1574 77 * 0 0 * * 0 0 AGGCGAGTGTGGGGGTTGTTTGAG +%%$#)##%##$####%###$%$# +1378_98_1574 141 * 0 0 * * 0 0 GTTCTGCCGGTGTCTGTGGCGGGC $$#+&$$####%$$$###$%#%%# +1378_107_1647 77 * 0 0 * * 0 0 GGTCTGGTTCTATGTTGGTCGACT ###'$$#$$$(#%###(#$##$%# +1378_107_1647 141 * 0 0 * * 0 0 AGGCCTACTACGCGTCATTGATAG &#$$#$(.#%#$$####&$%##($ +1378_111_829 77 * 0 0 * * 0 0 GGTCTTCTATGGTGCTAAATTAGT $##'&###$##)#%($#$%#$'%# +1378_111_829 141 * 0 0 * * 0 0 TGCGGCACTTGCTTCTTCGTATTT %#%##%#$%#$#%###$$##&#$$ +1378_111_1900 77 * 0 0 * * 0 0 GCACGCCTTTGGGCTAAGCCGTAA )$)'#%$########$'#&%$#(# +1378_111_1900 141 * 0 0 * * 0 0 TCCCCTCGCTCGGCTCTGTGCTGT $&%*$#(#)##$#'##%(##$#$% +1378_112_1483 77 * 0 0 * * 0 0 TGGAGTGGTGTGTTTGCTGAGCCA #$#)#############$#%#%'% +1378_112_1483 141 * 0 0 * * 0 0 TGTCCAGCTATGCGGCTTCCTCCT %#$+#%#&#$#####%####%$## +1378_125_1287 77 * 0 0 * * 0 0 TGACGTGGGTTGTCCCGTGAGATT ##$%%#$###$##$$#&%##$(%% +1378_125_1287 141 * 0 0 * * 0 0 TGTCTCTGGGGGGCCTGGTTAGGT $##13$'%#$###$$###$$$#&# +1378_126_468 77 * 0 0 * * 0 0 AAAAAACGGGCCGGGCGGACTAGG @@B@@@>@2>1+?:92)>@8&&<+ +1378_126_468 141 * 0 0 * * 0 0 TCTCCTGATTAGACATATCTCCGT 7>@318?=,=8@49579?..7,7) +1378_127_664 77 * 0 0 * * 0 0 TCGCTTTGCCTATGTTTGTTCGGA #%$%#&##$%#%%###$$###)-' +1378_127_664 141 * 0 0 * * 0 0 AGAGGTTGGTGTCTTGTCGCAGCT ##'#$######$$%######$$$# +1378_129_463 77 * 0 0 * * 0 0 AAAAAAAAAAAAAAAAAAAAAAAA ,)&20%8'(&35-***1-2&+510 +1378_129_463 141 * 0 0 * * 0 0 CCTGTGGGCCAAGCCCAATGAAAG 8),,<6(--)<()5&.&/8+;.2+ +1378_129_875 77 * 0 0 * * 0 0 GACCTTTACGTATTGGGGGTTGGC ###)###+###$##$#&%##$,#$ +1378_129_875 141 * 0 0 * * 0 0 TTTCTATGGCTTACGCTGTCTGCC #$($##%####%$#$#####$### +1378_140_1251 77 * 0 0 * * 0 0 TTTCCTTCGTGTGCGTGCGGAGTG #%#%$##$$$######.$$$%#%( +1378_140_1251 141 * 0 0 * * 0 0 ATCCTAGCGCGGTGTCTTGGGGAC #$%1#$$$##$##$#$#$##$%$$ +1378_141_809 77 * 0 0 * * 0 0 TCTCGTGGTTTCTTTTTTATGTGT ##%)##$$#####%$#$#%%#'## +1378_141_809 141 * 0 0 * * 0 0 TGTCCTCCAGTGTCTGTTGGGTGT %&,-##$$#(%###$#$$'###'# +1378_144_983 77 * 0 0 * * 0 0 GTTCGTTCGTGGTGTACGAGGGTG #(#%#####($#%##$$#%##%#) +1378_144_983 141 * 0 0 * * 0 0 AGCGCCCGGTTGGTGCGGCTCGTC -$(&%*$#*#))#$$$#%%$#$## +1378_153_270 77 * 0 0 * * 0 0 GGCCGTGTGCGGGTGTAGATTGGA %$##($######&##$&$$$$%## +1378_153_270 141 * 0 0 * * 0 0 AGTCCTTGTCCCCTGGGTTTTCCC +''$#&%$%#$##&$$($#&#$$# +1378_155_1689 77 * 0 0 * * 0 0 CGCGTTCGGACAAAGCTAGCACCT @4@@=6=@@<4:@@681@@@.24. +1378_155_1689 141 * 0 0 * * 0 0 TTGGGAGGGAGAGAGACTAGACCG ?4@@4<<??<?-@4@=4<5<?:/> +1378_157_1580 77 * 0 0 * * 0 0 GGGATTGAAGGGATGTATGCTAGG #%$&%#$$'%$%#$##*#%$$$$# +1378_157_1580 141 * 0 0 * * 0 0 TGGGCCTCGGTGCCCTTGGTCTGT #%)$##'#$$$&#####%#$#$## +1378_161_317 77 * 0 0 * * 0 0 AATCCATACCCACAAAAGCAGGCC .&%','(@''?7//+&)+2.+)0) +1378_161_317 141 * 0 0 * * 0 0 TTGGCCGGCAACCCCGGTACCTAA 7<,<'@)@>.)2@/')'&(?/-<( +1378_177_735 77 * 0 0 * * 0 0 CGAGCCCTAAACCATGAGATCGGA @@B>@AC9@A<B8@=9>7@5)>,0 +1378_177_735 141 * 0 0 * * 0 0 AGTACGCCATGTATTTGCGACCAG =?@?C<88@=>:7>@55/.,416, +1378_181_1684 77 * 0 0 * * 0 0 TTTCTGTTGTGGTTTTGTTGGGGT $##'$%'##%##$%$#$$####$* +1378_181_1684 141 * 0 0 * * 0 0 CGACTCCCGCATTCACGGTCAAGT &*#,##$#&$*$$#$#$$$#%$## +1378_187_1407 77 * 0 0 * * 0 0 TTGGGTGAAATCTTGTCGAGTGGA ####&##$$###$#####%##%%) +1378_187_1407 141 * 0 0 * * 0 0 TGGCGTCCACTCGTGGGTCTATCG $#$'%#$%$%&$%#####$#$#%# +1378_203_721 77 * 0 0 * * 0 0 CCTTTACGATCATAAACCATGACC 9<&./1&@;12')?<1',?/)&/. +1378_203_721 141 * 0 0 * * 0 0 CCTAAATAACCCAGGGTAAAAGAT 9<1+2+1@=259@+;2.71&@+2& +1378_206_2039 77 * 0 0 * * 0 0 CCCCCCCCCCCCCCCCCCCCCCCC ?=>>?=:;<>>4::<=56199(05 +1378_206_2039 141 * 0 0 * * 0 0 CCCCCCCCCCCCCCCCCCCCCCCC 9@?8>;>>><<994;=895895-8 diff -r 1a56a541f9f9 -r 14d54c9dbd59 test-data/bwa_wrapper_out3.sam --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/bwa_wrapper_out3.sam Wed Sep 16 11:14:39 2009 -0400 @@ -0,0 +1,2 @@ +@SQ SN:PHIX174 LN:5386 +081017-and-081020:1:1:1715:1759 16 PHIX174 322 25 36M * 0 0 GATATTTTAAAGGAGCGTGGATTACTATCTGAGTCC B&&I13A$G$*%$IIIIIII9(.+5$IIIIIII#II XT:A:U NM:i:2 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:2C8A24 diff -r 1a56a541f9f9 -r 14d54c9dbd59 test-data/bwa_wrapper_out4.sam --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/bwa_wrapper_out4.sam Wed Sep 16 11:14:39 2009 -0400 @@ -0,0 +1,3 @@ +@SQ SN:PHIX174 LN:5386 +081017-and-081020:1:1:1715:1759 113 PHIX174 322 25 18M = 340 18 GATATTTTAAAGGAGCGT B&&I13A$G$*%$IIIII XT:A:U NM:i:2 SM:i:25 AM:i:25 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0 MD:Z:2C8A6 +081017-and-081020:1:1:1715:1759 177 PHIX174 340 37 18M = 322 -18 GGATTACTATCTGAGTCC II9(.+5$IIIIIII#II XT:A:U NM:i:0 SM:i:37 AM:i:25 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0 MD:Z:18 diff -r 1a56a541f9f9 -r 14d54c9dbd59 test-data/fastq_conv_in1.fastq --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/fastq_conv_in1.fastq Wed Sep 16 11:14:39 2009 -0400 @@ -0,0 +1,4 @@ +@081017-and-081020:1:1:1715:1759 +GGACTCAGATAGTAATCCACGCTCCTTTAAAATATC ++ +II#IIIIIII$5+.(9IIIIIII$%*$G$A31I&&B diff -r 1a56a541f9f9 -r 14d54c9dbd59 tools/next_gen_conversion/fastq_conversions.xml --- a/tools/next_gen_conversion/fastq_conversions.xml Tue Sep 15 14:26:20 2009 -0400 +++ b/tools/next_gen_conversion/fastq_conversions.xml Wed Sep 16 11:14:39 2009 -0400 @@ -52,7 +52,7 @@ <tests> <test> <param name="type" value="sol2std" /> - <param name="input" value="bwa_phiX_sanger.fastq" ftype="fastqsolexa" /> + <param name="input" value="fastq_conv_in1.fastq" ftype="fastqsolexa" /> <output name="outputFastqsanger" file="fastq_conv_out1.fastqsanger" /> </test> <test> diff -r 1a56a541f9f9 -r 14d54c9dbd59 tools/sr_mapping/bwa_wrapper.py --- a/tools/sr_mapping/bwa_wrapper.py Tue Sep 15 14:26:20 2009 -0400 +++ b/tools/sr_mapping/bwa_wrapper.py Wed Sep 16 11:14:39 2009 -0400 @@ -3,8 +3,6 @@ """ Runs BWA on single-end or paired-end data. Produces a SAM file containing the mappings. - -usage: python bwa_wrapper.py reference_sequence indexing_algorithm(is_or_bwtsw) forward_fastq_file reverse_fastq_file(or_None) output alignment_type(single_or_paired) parameters(pre_set_or_full) file_type(solexa_or_solid) file_source(indexed_or_history) maxEditDist fracMissingAligns maxGapOpens maxGapExtens disallowLongDel disallowIndel seed maxEditDistSeed numThreads mismatchPenalty gapOpenPenalty gapExtensPenalty colorSpaceRev suboptAlign noIterSearch outputTopN maxInsertSize maxOccurPairing\nThe last eighteen need to all be specified, or all be None """ import optparse, os, sys, tempfile @@ -16,6 +14,7 @@ def __main__(): #Parse Command Line parser = optparse.OptionParser() + parser.add_option('', '--threads', dest='threads', help='The number of threads to use') parser.add_option('', '--ref', dest='ref', help='The reference genome to use or index') parser.add_option('', '--fastq', dest='fastq', help='The (forward) fastq file to use for the mapping') parser.add_option('', '--rfastq', dest='rfastq', help='The reverse fastq file to use for mapping if paired-end data') @@ -32,7 +31,6 @@ parser.add_option('-i', '--disallowIndel', dest='disallowIndel', help='Disallow indel within specified bps') parser.add_option('-l', '--seed', dest='seed', help='Take the first specified subsequences') parser.add_option('-k', '--maxEditDistSeed', dest='maxEditDistSeed', help='Maximum edit distance to the seed') - parser.add_option('-t', '--numThreads', dest='numThreads', help='Number of threads') parser.add_option('-M', '--mismatchPenalty', dest='mismatchPenalty', help='Mismatch penalty') parser.add_option('-O', '--gapOpenPenalty', dest='gapOpenPenalty', help='Gap open penalty') parser.add_option('-E', '--gapExtensPenalty', dest='gapExtensPenalty', help='Gap extension penalty') @@ -42,8 +40,9 @@ parser.add_option('', '--outputTopN', dest='outputTopN', help='Output top specified hits') parser.add_option('', '--maxInsertSize', dest='maxInsertSize', help='Maximum insert size for a read pair to be considered mapped good') parser.add_option('', '--maxOccurPairing', dest='maxOccurPairing', help='Maximum occurrences of a read for pairings') + parser.add_option('', '--dbkey', dest='dbkey', help='') (options, args) = parser.parse_args() - + # index if necessary if options.fileSource == 'history': # make temp directory for placement of indices and copy reference file there @@ -65,46 +64,42 @@ else: indexing_cmds = '-a %s' % indexingAlg options.ref = os.path.join(tmp_dir,os.path.split(options.ref)[1]) - cmd1 = 'bwa index %s %s 2> /dev/null' % (indexing_cmds, options.ref) + cmd1 = 'cd %s; bwa index %s %s 2> /dev/null' % (tmp_dir, indexing_cmds, options.ref) try: os.system(cmd1) except Exception, erf: stop_err('Error indexing reference sequence\n' + str(erf)) - # set up aligning and generate aligning command options if options.params == 'pre_set': if options.fileType == 'solid': - aligning_cmds = '-c' + aligning_cmds = '-c -t %s' % options.threads else: - aligning_cmds = '' + aligning_cmds = '-t %s' % options.threads gen_alignment_cmds = '' else: aligning_cmds = '-n %s -o %s -e %s -d %s -i %s %s -k %s -t %s -M %s -O %s -E %s %s %s %s' % \ ((options.fracMissingAligns, options.maxEditDist)[options.maxEditDist != '0'], options.maxGapOpens, options.maxGapExtens, options.disallowLongDel, - options.disallowIndel, ('',' -l %s'%options.seed)[options.seed!=-1], - options.maxEditDistSeed, options.numThreads, options.mismatchPenalty, + options.disallowIndel, ('',' -l %s'%options.seed)[options.seed!='-1'], + options.maxEditDistSeed, options.threads, options.mismatchPenalty, options.gapOpenPenalty, options.gapExtensPenalty, ('',' -c')[options.colorSpaceRev=='true'], ('',' -R')[options.suboptAlign=='true'], ('',' -N')[options.noIterSearch=='true']) if options.genAlignType == 'single': gen_alignment_cmds = '-n %s' % options.outputTopN elif options.genAlignType == 'paired': gen_alignment_cmds = '-a %s -o %s' % (options.maxInsertSize, options.maxOccurPairing) - - # set up output file - file(options.output,'w').write('QNAME\tFLAG\tRNAME\tPOS\tMAPQ\tCIGAR\tMRNM\tMPOS\tISIZE\tSEQ\tQUAL\tOPT\n') + # set up output files tmp_align_out = tempfile.NamedTemporaryFile() + tmp_align_out2 = tempfile.NamedTemporaryFile() # prepare actual aligning and generate aligning commands cmd2 = 'bwa aln %s %s %s > %s 2> /dev/null' % (aligning_cmds, options.ref, options.fastq, tmp_align_out.name) cmd2b = '' if options.genAlignType == 'paired': - tmp_align_out2 = tempfile.NamedTemporaryFile() cmd2b = 'bwa aln %s %s %s > %s 2> /dev/null' % (aligning_cmds, options.ref, options.rfastq, tmp_align_out2.name) cmd3 = 'bwa sampe %s %s %s %s %s %s >> %s 2> /dev/null' % (gen_alignment_cmds, options.ref, tmp_align_out.name, tmp_align_out2.name, options.fastq, options.rfastq, options.output) else: cmd3 = 'bwa samse %s %s %s %s >> %s 2> /dev/null' % (gen_alignment_cmds, options.ref, tmp_align_out.name, options.fastq, options.output) - # align try: os.system(cmd2) @@ -116,11 +111,13 @@ os.system(cmd2b) except Exception, erf: stop_err("Error aligning second sequence\n" + str(erf)) - # generate align try: os.system(cmd3) except Exception, erf: stop_err("Error sequence aligning sequence\n" + str(erf)) - + # clean up temp files + tmp_align_out.close() + tmp_align_out2.close() + if __name__=="__main__": __main__() diff -r 1a56a541f9f9 -r 14d54c9dbd59 tools/sr_mapping/bwa_wrapper.xml --- a/tools/sr_mapping/bwa_wrapper.xml Tue Sep 15 14:26:20 2009 -0400 +++ b/tools/sr_mapping/bwa_wrapper.xml Wed Sep 16 11:14:39 2009 -0400 @@ -1,7 +1,8 @@ -<tool id="bwa_wrapper" name="BWA" version="1.0.0"> +<tool id="bwa_wrapper" name="BWA" version="1.0.1"> <description> fast mapping of reads against reference sequence</description> <command interpreter="python"> bwa_wrapper.py + --threads="8" #if $solidOrSolexa.solidRefGenomeSource.refGenomeSource == "history": --ref=$solidOrSolexa.solidRefGenomeSource.ownFile #else: @@ -27,7 +28,6 @@ --disallowIndel="None" --seed="None" --maxEditDistSeed="None" - --numThreads="None" --mismatchPenalty="None" --gapOpenPenalty="None" --gapExtensPenalty="None" @@ -46,7 +46,6 @@ --disallowIndel=$params.disallowIndel --seed=$params.seed --maxEditDistSeed=$params.maxEditDistSeed - --numThreads=$params.numThreads --mismatchPenalty=$params.mismatchPenalty --gapOpenPenalty=$params.gapOpenPenalty --gapExtensPenalty=$params.gapExtensPenalty @@ -57,9 +56,14 @@ --maxInsertSize=$params.maxInsertSize --maxOccurPairing=$params.maxOccurPairing #end if + #if $solidOrSolexa.solidRefGenomeSource.refGenomeSource == "history": + --dbkey=$dbkey + #else: + --dbkey="None" + #end if </command> <inputs> - <conditional name="solidOrSolexa"> + <conditional name="solidOrSolexa"> <param name="solidSolexa" type="select" label="Select SOLiD or Solexa format for the original dataset"> <option value="solexa">Solexa</option> <option value="solid">SOLiD</option> @@ -70,18 +74,18 @@ <option value="indexed">Use a built-in index</option> <option value="history">Use one from the history</option> </param> - <when value="history"> - <param name="ownFile" type="data" format="fasta" label="Select a reference genome" /> - </when> <when value="indexed"> <param name="indices" type="select" label="Select a reference genome"> <options from_file="sequence_index_color.loc"> <column name="value" index="1" /> <column name="name" index="0" /> - <filter type="sort_by" column="0" /> + <filter type="sort_by" column="0" /> </options> </param> </when> + <when value="history"> + <param name="ownFile" type="data" format="fasta" metadata_name="dbkey" label="Select a reference genome" /> + </when> </conditional> </when> <when value="solexa"> @@ -90,39 +94,39 @@ <option value="indexed">Use a built-in index</option> <option value="history">Use one from the history</option> </param> - <when value="history"> - <param name="ownFile" type="data" format="fasta" label="Select a reference genome" /> - </when> <when value="indexed"> <param name="indices" type="select" label="Select a reference genome"> <options from_file="sequence_index_base.loc"> <column name="value" index="1" /> <column name="name" index="0" /> - <filter type="sort_by" column="0" /> + <filter type="sort_by" column="0" /> </options> </param> + </when> + <when value="history"> + <param name="ownFile" type="data" format="fasta" metadata_name="dbkey" label="Select a reference genome" /> </when> </conditional> </when> </conditional> - <conditional name="paired"> - <param name="sPaired" type="select" label="Is this library mate-paired?"> - <option value="single">Single-end</option> - <option value="paired">Paired-end</option> - </param> - <when value="single"> - <param name="input1" type="data" format="fastqsanger" label="FASTQ file" /> - </when> - <when value="paired"> - <param name="input1" type="data" format="fastqsanger" label="Forward FASTQ file" /> - <param name="input2" type="data" format="fastqsanger" label="Reverse FASTQ file" /> - </when> - </conditional> - <conditional name="params"> + <conditional name="paired"> + <param name="sPaired" type="select" label="Is this library mate-paired?"> + <option value="single">Single-end</option> + <option value="paired">Paired-end</option> + </param> + <when value="single"> + <param name="input1" type="data" format="fastqsanger" label="FASTQ file" /> + </when> + <when value="paired"> + <param name="input1" type="data" format="fastqsanger" label="Forward FASTQ file" /> + <param name="input2" type="data" format="fastqsanger" label="Reverse FASTQ file" /> + </when> + </conditional> + <conditional name="params"> <param name="source_select" type="select" label="BWA settings to use" help="For most mapping needs use Commonly used settings. If you want full control use Full List"> - <option value="pre_set">Commonly used</option> - <option value="full">Full Parameter List</option> - </param> + <option value="pre_set">Commonly used</option> + <option value="full">Full Parameter List</option> + </param> <when value="pre_set" /> <when value="full"> <param name="maxEditDist" type="integer" value="0" label="Maximum edit distance" help="Enter this value OR a fraction of missing alignments, not both" /> @@ -133,9 +137,8 @@ <param name="disallowIndel" type="integer" value="5" label="Disallow insertion/deletion within [value] bp towards the end" /> <param name="seed" type="integer" value="-1" label="Number of first subsequences to take as seed" help="Enter -1 for infinity" /> <param name="maxEditDistSeed" type="integer" value="2" label="Maximum edit distance in the seed" /> - <param name="numThreads" type="integer" value="1" label="Number of threads, in multi-threading mode" /> <param name="mismatchPenalty" type="integer" value="3" label="Mismatch penalty" help="BWA will not search for suboptimal hits with a score lower than [value]" /> - <param name="gapOpenPenalty" type="integer" value="1" label="Gap open penalty" /> + <param name="gapOpenPenalty" type="integer" value="11" label="Gap open penalty" /> <param name="gapExtensPenalty" type="integer" value="4" label="Gap extension penalty" /> <param name="colorSpaceRev" type="select" label="Reverse query but don't compement it" help="Reverse query for all alignment in color space"> <option value="false">Don't reverse query</option> @@ -150,34 +153,33 @@ </conditional> </inputs> <outputs> - <data format="tabular" name="output" /> + <data format="sam" name="output" /> </outputs> -<!-- Tests all fail because of problem with nested conditionals in test framework <tests> <test> <param name="solidSolexa" value="solexa" /> <param name="refGenomeSource" value="indexed" /> <param name="indices" value="phiX" /> <param name="sPaired" value="single" /> - <param name="input1" value="bwa_phiX_sanger.fastq" /> + <param name="input1" value="bwa_wrapper_in1.fastq" /> <param name="source_select" value="pre_set" /> - <output name="output" file="bwa_wrapper_out0.tabular" /> + <output name="output" file="bwa_wrapper_out0.sam" ftype="sam" /> </test> <test> <param name="solidSolexa" value="solid" /> <param name="refGenomeSource" value="history" /> <param name="ownFile" value="phiX.fa" /> <param name="sPaired" value="single" /> - <param name="input1" value="bwa_phiX_sanger.fastq" /> + <param name="input1" value="bwa_wrapper_in0.fastq" /> <param name="source_select" value="pre_set" /> - <output name="output" file="bwa_wrapper_out0b.tabular" /> + <output name="output" file="bwa_wrapper_out0b.sam" ftype="sam" /> </test> <test> <param name="solidSolexa" value="solid" /> <param name="refGenomeSource" value="indexed" /> <param name="indices" value="phiX" /> <param name="sPaired" value="single" /> - <param name="input1" value="bwa_solid.fastq" /> + <param name="input1" value="bwa_wrapper_in2.fastq" /> <param name="source_select" value="full" /> <param name="maxEditDist" value="0" /> <param name="fracMissingAligns" value="0.04" /> @@ -187,9 +189,8 @@ <param name="disallowIndel" value="5" /> <param name="seed" value="-1" /> <param name="maxEditDistSeed" value="2" /> - <param name="numThreads" value="1" /> <param name="mismatchPenalty" value="3" /> - <param name="gapOpenPenalty" value="1" /> + <param name="gapOpenPenalty" value="11" /> <param name="gapExtensPenalty" value="4" /> <param name="colorSpaceRev" value="true" /> <param name="suboptAlign" value="true" /> @@ -197,15 +198,15 @@ <param name="outputTopN" value="-1" /> <param name="maxInsertSize" value="500" /> <param name="maxOccurPairing" value="100000" /> - <output name="output" file="bwa_wrapper_out1.tabular" /> + <output name="output" file="bwa_wrapper_out1.sam" ftype="sam" /> </test> <test> <param name="solidSolexa" value="solid" /> <param name="refGenomeSource" value="indexed" /> <param name="indices" value="phiX" /> <param name="sPaired" value="paired" /> - <param name="input1" value="bwa_solid_f.fastq" /> - <param name="input2" value="bwa_solid_r.fastq" /> + <param name="input1" value="bwa_wrapper_in3.fastq" /> + <param name="input2" value="bwa_wrapper_in4.fastq" /> <param name="source_select" value="full" /> <param name="maxEditDist" value="0" /> <param name="fracMissingAligns" value="0.04" /> @@ -215,9 +216,8 @@ <param name="disallowIndel" value="5" /> <param name="seed" value="-1" /> <param name="maxEditDistSeed" value="2" /> - <param name="numThreads" value="1" /> <param name="mismatchPenalty" value="3" /> - <param name="gapOpenPenalty" value="1" /> + <param name="gapOpenPenalty" value="11" /> <param name="gapExtensPenalty" value="4" /> <param name="colorSpaceRev" value="true" /> <param name="suboptAlign" value="true" /> @@ -225,14 +225,14 @@ <param name="outputTopN" value="-1" /> <param name="maxInsertSize" value="500" /> <param name="maxOccurPairing" value="100000" /> - <output name="output" file="bwa_wrapper_out2.tabular" /> + <output name="output" file="bwa_wrapper_out2.sam" ftype="sam" /> </test> <test> <param name="solidSolexa" value="solexa" /> <param name="refGenomeSource" value="indexed" /> <param name="indices" value="phiX" /> <param name="sPaired" value="single" /> - <param name="input1" value="bwa_phiX_sanger.fastq" /> + <param name="input1" value="bwa_wrapper_in1.fastq" /> <param name="source_select" value="full" /> <param name="maxEditDist" value="0" /> <param name="fracMissingAligns" value="0.04" /> @@ -242,9 +242,8 @@ <param name="disallowIndel" value="5" /> <param name="seed" value="-1" /> <param name="maxEditDistSeed" value="2" /> - <param name="numThreads" value="1" /> <param name="mismatchPenalty" value="3" /> - <param name="gapOpenPenalty" value="1" /> + <param name="gapOpenPenalty" value="11" /> <param name="gapExtensPenalty" value="4" /> <param name="colorSpaceRev" value="false" /> <param name="suboptAlign" value="true" /> @@ -252,15 +251,15 @@ <param name="outputTopN" value="-1" /> <param name="maxInsertSize" value="500" /> <param name="maxOccurPairing" value="100000" /> - <output name="output" file="bwa_wrapper_out3.tabular" /> + <output name="output" file="bwa_wrapper_out3.sam" ftype="sam" /> </test> <test> <param name="solidSolexa" value="solexa" /> <param name="refGenomeSource" value="indexed" /> <param name="indices" value="phiX" /> <param name="sPaired" value="paired" /> - <param name="input1" value="bwa_phiX_sanger_f.fastq" /> - <param name="input2" value="bwa_phiX_sanger_r.fastq" /> + <param name="input1" value="bwa_wrapper_in5.fastq" /> + <param name="input2" value="bwa_wrapper_in6.fastq" /> <param name="source_select" value="full" /> <param name="maxEditDist" value="0" /> <param name="fracMissingAligns" value="0.04" /> @@ -270,9 +269,8 @@ <param name="disallowIndel" value="5" /> <param name="seed" value="-1" /> <param name="maxEditDistSeed" value="2" /> - <param name="numThreads" value="1" /> <param name="mismatchPenalty" value="3" /> - <param name="gapOpenPenalty" value="1" /> + <param name="gapOpenPenalty" value="11" /> <param name="gapExtensPenalty" value="4" /> <param name="colorSpaceRev" value="false" /> <param name="suboptAlign" value="true" /> @@ -280,11 +278,10 @@ <param name="outputTopN" value="-1" /> <param name="maxInsertSize" value="500" /> <param name="maxOccurPairing" value="100000" /> - <output name="output" file="bwa_wrapper_out4.tabular" /> - </test> + <output name="output" file="bwa_wrapper_out4.sam" ftype="sam" /> + </test> </tests> ---> - <help> + <help> **What it does** @@ -294,7 +291,7 @@ **Input formats** -BWA accepts files in FASTQ format. +BWA accepts files in Sanger FASTQ format. ------ @@ -336,7 +333,7 @@ ------- -**BWA Settings** +**BWA settings** All of the options have a default value. You can change any of them. All of the options in BWA have been implemented here. @@ -390,9 +387,9 @@ -o INT Maximum occurrences of a read for pairing. A read with more occurrences will be treated as a single-end read. Reducing this parameter helps faster pairing. [100000] - - </help> + </help> + <code file="bwa_wrapper_code.py" /> </tool> diff -r 1a56a541f9f9 -r 14d54c9dbd59 tools/sr_mapping/bwa_wrapper_code.py --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tools/sr_mapping/bwa_wrapper_code.py Wed Sep 16 11:14:39 2009 -0400 @@ -0,0 +1,8 @@ +import os + +def exec_before_job(app, inp_data, out_data, param_dict, tool): + try: + refFile = param_dict['solidOrSolexa']['solidRefGenomeSource']['indices'].value + out_data['output'].set_dbkey(os.path.split(refFile)[1].split('.')[0]) + except Exception, eq: + out_data['output'].set_dbkey(param_dict['dbkey'])