.
+
+=head1 SEE ALSO
+
+constant(3), perl(1).
+
+=cut
diff -r 000000000000 -r acc8d8bfeb9a example/germline.bam
Binary file example/germline.bam has changed
diff -r 000000000000 -r acc8d8bfeb9a example/germline.bam.bai
Binary file example/germline.bam.bai has changed
diff -r 000000000000 -r acc8d8bfeb9a example/germline.cover
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/example/germline.cover Wed Feb 08 16:59:24 2012 -0500
@@ -0,0 +1,2 @@
+4 125892484 + 1 17
+4 125894262 + 2 11
diff -r 000000000000 -r acc8d8bfeb9a example/germline.sclip.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/example/germline.sclip.txt Wed Feb 08 16:59:24 2012 -0500
@@ -0,0 +1,3 @@
+4 125892484 + HWUSI-EAS729_614T6:6:61:16691:2909#0 A C
+4 125894262 + HWUSI-EAS1543_614TM:1:62:2828:15248#0 AA CC
+4 125894262 + HWUSI-EAS729_614T6:3:68:9882:5425#0 AAAGAGTT CCCCCCCC
diff -r 000000000000 -r acc8d8bfeb9a example/tumor.bam
Binary file example/tumor.bam has changed
diff -r 000000000000 -r acc8d8bfeb9a example/tumor.bam.bai
Binary file example/tumor.bam.bai has changed
diff -r 000000000000 -r acc8d8bfeb9a example/tumor.cover
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/example/tumor.cover Wed Feb 08 16:59:24 2012 -0500
@@ -0,0 +1,11 @@
+4 125892327 + 1 30
+4 125892458 + 1 27
+4 125893225 + 1 28
+4 125893227 + 5 29
+4 125893365 - 1 26
+4 125893979 - 1 16
+10 66301086 - 2 33
+10 66301858 + 4 15
+10 66301865 - 9 22
+10 66301871 - 1 24
+10 66302136 + 1 51
diff -r 000000000000 -r acc8d8bfeb9a example/tumor.cover.somatic.cover
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/example/tumor.cover.somatic.cover Wed Feb 08 16:59:24 2012 -0500
@@ -0,0 +1,11 @@
+4 125892327 + 1 30
+4 125892458 + 1 27
+4 125893225 + 1 28
+4 125893227 + 5 29
+4 125893365 - 1 26
+4 125893979 - 1 16
+10 66301086 - 2 33
+10 66301858 + 4 15
+10 66301865 - 9 22
+10 66301871 - 1 24
+10 66302136 + 1 51
diff -r 000000000000 -r acc8d8bfeb9a example/tumor.predSV.html
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/example/tumor.predSV.html Wed Feb 08 16:59:24 2012 -0500
@@ -0,0 +1,379 @@
+4 125893227 + 5 10 66301858 - 4 CTX 29 14 83 71 0.895173453996983 0.230769230769231 0.735384615384615 0.5 1 4 125893135 176 10 66301773 TTATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTTCAGTAGCTGTCACCTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCTCCCGAGTAGCTGGGATTACAGGTGCCCACCACCATGCCTGGCTAATTTT
+
+ 125893140 125893160 125893180 125893200 125893220 125893240 125893260 125893280 125893300 125893320
+ . | . | . | . | . | . | . | . | . | . |
+4:125893129-125893326 CTAAAATTATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCATCAAAATTTTAAAATATCCCTT
+
+HWI-EAS90_614M9:6:67:12631:17754#0 + CTAAAATTATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACT
+HWUSI-EAS1591_6113C:4:21:6779:18838#0 + CTAAAATTATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACT
+HWI-EAS90_614M9:1:32:10389:16117#0 + attttGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*Cagtagctgtcactgtaggcaaaaat
+HWUSI-EAS1591_6113C:2:27:1355:21017#0 + ATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCaacaat
+HWUSI-EAS1591_6113C:1:93:16865:6051#0 + AAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCAT
+HWI-EAS90_614M9:3:78:1961:11755#0 + ATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAaAAtGCCCAtTATTC*TGC
+HWI-EAS231_614TA:5:6:17800:6620#0 + TTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATAT
+HWI-EAS90_614M9:3:95:10317:10581#0 + ATATATCATATTCTAAATTATGAAAAGAGAATATGA*Tctcttttccag*agccggccccggaagccaaaaagccccttattcttcctaattatctagggga
+HWI-EAS292_102081428:2:115:1681:2531#0 + aaattatgaaaagagaatatgattctctttt*taggagctgttactgtagggaaaaatgcccattattctgctaaatatctaagtgaggggtgcttctatt
+HWI-EAS90_614M9:1:67:1331:13704#0 + AATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATA
+HWI-EAS231_614TA:2:73:4389:2239#0 + TGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAG
+HWI-EAS231_614TA:7:88:14296:6002#0 + TGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAG
+HWI-EAS231_614TA:1:26:7960:13504#0 + TTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGT
+HWI-EAS90_614M9:7:18:13917:3050#0 + CTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCA
+HWI-EAS292_102081428:1:66:16549:20224#0+ TCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TgCTAAATATCTAAGTGATGCTTGCTACTATTAACTGaAtGTATaTACTAGCGTCAA
+HWI-EAS90_614M9:5:77:4642:6832#0 + CTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTaACTGAATGTATATACTAGCgtcaaa
+HWI-EAS90_614M9:8:34:7077:4638#0 + CTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAA
+HWI-EAS231_614TA:2:33:16646:21114#0 + CTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAAAttttaaaatatccctt
+HWUSI-EAS1591_6113C:4:9:4513:12290#0 - TTATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACCTCCTGG
+HWI-EAS231_614TA:7:49:16816:4213#0 - TATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCA
+HWI-EAS231_614TA:1:118:16290:13237#0 - ATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAAT
+HWI-EAS231_614TA:6:113:8101:10568#0 - TTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATG
+HWI-EAS90_614M9:6:48:8352:11907#0 - ATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATT
+HWUSI-EAS1591_6113C:4:39:14328:11130#0 - ATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACCTCCTGGGTTCAAGTGATTCTCCTG
+HWI-EAS90_614M9:1:85:18838:8742#0 - TTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACCTCCtgggttcaggTGATTCTCCTGCCTCTACCT
+HWI-EAS90_614M9:4:32:9763:4986#0 - TTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCaCCTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCT
+HWI-EAS231_614TA:4:103:18855:9257#0 - ATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGT
+HWI-EAS292_102081428:2:23:3944:1922#0 - aTAtaGATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGT
+HWUSI-EAS1591_6113C:2:117:3095:2051#0 - CTCTTTT*CAGtagctgtCaCCTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCTCCCGAGTAGCTGGGATTACAGGTGCCCACCACCATGCCTGGCTAAT
+HWI-EAS90_614M9:5:28:11006:5910#0 - TCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAA
+HWI-EAS90_614M9:5:1:12903:5535#0 - TTTT*CAGTAgCTGTCACCTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCTCCCGAGTAGCTGGGATTACAGGTGCCCACCACCATGCCTGGCTAATTTT
+HWI-EAS90_614M9:2:56:9827:2778#0 - tcactgtaGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAAATTTTAAAATATCC
+HWI-EAS90_614M9:4:47:5786:19724#0 - tcactGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAAATTTTAAAATATCC
+HWI-EAS90_614M9:6:67:11670:7706#0 - TCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAAATTTTAAAATATCC
+HWUSI-EAS1591_6113C:1:117:14883:9756#0 - TCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAAATTTTAAAATATCC
+HWI-EAS90_614M9:6:37:7790:16256#0 - CACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAAATTTTAAAATATCCC
+HWUSI-EAS1591_6113C:4:38:15097:15951#0 - ACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAAATTTTAAAATATCCCT
+_________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
+HWUSI-EAS729_614T6:8:86:3208:16885#0 + TTATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGC
+HWUSI-EAS729_614T6:3:97:10823:18653#0 + ATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAA
+HWUSI-EAS729_614T6:4:113:3054:5730#0 + TGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGgcaag
+HWUSI-EAS729_614T6:2:63:5049:3137#0 + AATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAA
+HWUSI-EAS729_614T6:5:95:4777:9034#0 + AATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAA
+HWUSI-EAS729_614T6:8:79:9632:17904#0 + AATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAA
+HWUSI-EAS1543_614TM:8:50:14524:11101#0 + ATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAAT
+HWUSI-EAS729_614T6:3:79:10477:10666#0 + TTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCC
+HWUSI-EAS729_614T6:4:80:3846:17138#0 + TTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCC
+HWUSI-EAS729_614T6:4:83:3398:15440#0 + TTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCC
+HWUSI-EAS729_614T6:5:37:19001:5491#0 + TGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCC
+HWUSI-EAS729_614T6:8:113:1681:2586#0 + TCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCT
+HWUSI-EAS729_614T6:8:3:18130:10917#0 + TATTTTGAAATATATATCATATtCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAA
+HWUSI-EAS729_614T6:7:32:3784:11386#0 + TGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCT
+HWUSI-EAS729_614T6:7:117:8731:4955#0 + ATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGA
+HWUSI-EAS729_614T6:8:51:17611:11416#0 + TATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGAT
+HWUSI-EAS1543_614TM:1:102:4785:3367#0 + TATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCT
+HWUSI-EAS1543_614TM:6:60:2952:14857#0 + TATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCT
+HWUSI-EAS729_614T6:5:55:11097:17856#0 + ATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTG
+HWUSI-EAS729_614T6:2:8:16102:20401#0 + AAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAAT
+HWUSI-EAS1543_614TM:1:75:3691:5075#0 + GATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGC
+HWUSI-EAS729_614T6:8:52:15200:10784#0 + CACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCATCAAAattntAAAATATCCC
+HWUSI-EAS1543_614TM:7:17:14845:19479#0 - ATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATAATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAAT
+HWUSI-EAS729_614T6:3:13:15991:15051#0 - TATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TG
+HWUSI-EAS729_614T6:7:65:3141:2489#0 - TATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TG
+HWUSI-EAS729_614T6:8:54:5081:15134#0 - TATCATATTTTcAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TG
+HWUSI-EAS1543_614TM:4:38:16709:14763#0 - CATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTA
+HWUSI-EAS1543_614TM:8:50:14524:11101#0 - TATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCT
+HWUSI-EAS729_614T6:1:56:9586:15810#0 - TAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATaTAAGTGATGCTTGCTACTAT
+HWUSI-EAS729_614T6:1:37:11212:13638#0 - TAAATTATGAAtAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTAT
+HWUSI-EAS729_614T6:1:29:1825:18616#0 - AATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTA
+HWUSI-EAS729_614T6:5:48:14044:12231#0 - AATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTA
+HWUSI-EAS729_614T6:4:115:6035:15619#0 - ATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAA
+HWUSI-EAS729_614T6:5:72:1269:20227#0 - aagagaatatgattctctTtT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGT
+HWUSI-EAS729_614T6:7:32:3784:11386#0 - AAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGT
+HWUSI-EAS1543_614TM:1:102:4785:3367#0 - AGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAAATTTTA
+HWUSI-EAS1543_614TM:6:60:2952:14857#0 - agtagCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGctactATTAACTGAATGTATATACTAGCGTCAAAATTTTA
+_________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
+
+
+ 66301780 66301800 66301820 66301840 66301860 66301880 66301900 66301920 66301940
+ . | . | . | . | . | . | . | . | . | .
+10:66301762-66301954 GTAAAAACACAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGA*AGAGGATATGTA
+
+HWI-EAS231_614TA:7:82:16847:8457#0 + CACAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGgagaatcacttgaacccagggggtggaggttggag
+HWI-EAS231_614TA:3:12:6209:13521#0 + CAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGagaatcacttgaacccaggagggggagggtggagtg
+HWI-EAS90_614M9:8:23:6627:13184#0 + GAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAaa
+HWI-EAS90_614M9:6:39:4439:17404#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAagtggga*agaggat
+HWI-EAS292_102081428:4:40:3376:13334#0 - AAAAAtTAgcCAGGCATGGtGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGA
+HWI-EAS90_614M9:4:10:8095:5678#0 - AAAAAttagCcaGgcatGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGA
+HWI-EAS90_614M9:6:117:3466:13678#0 - gccaggggtggtggtgggcgcctgtaatcccagctactcgGGAGGTAGAGGCaGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGAT
+HWI-EAS231_614TA:6:38:5961:6062#0 - aggcgtggtggtgggcacctgtaatcccagctactcGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCAT
+HWI-EAS231_614TA:1:37:3390:19595#0 - acacaaaaattagccaggcaTGGTGGTGGGCaCCTGTAATCCCAGCTaCTcGGGAGGtAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGACAGCTACT
+HWI-EAS90_614M9:1:98:9476:6999#0 - ggTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCAgTTGAACCCAGGAGGTGACAGCTACTGAAAAGAGAATCATATTCTCTTTT
+HWI-EAS231_614TA:8:4:10123:13264#0 - gTAatCCCAGCTaCtCGGGAGGTAGaggcaGGAGAATCaCTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGG
+HWUSI-EAS1591_6113C:2:96:6055:2122#0 - actcgGGAGGTAGAGGCAGGAGAATCaCTTGAACCCAGGAGGTGGAaGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAA
+HWI-EAS90_614M9:7:110:15112:12998#0 - GTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGACAGCTACTGAAAAGAGAATCATATTCTCTTTTCATAATTTAGAATATGATATATATTTCAAAAT
+HWI-EAS90_614M9:5:12:2797:13435#0 - GCAGGAGAATCACTTGAACCCAGGAGGTGACAGCTACTGAAAAGAGAATCATATTCTCTTTTCATAATTTAGAATATGATATATATTTCAAAATATGATA
+HWI-EAS231_614TA:7:82:16847:8457#0 - ggtggaggttggagtgagccaagatcatgccacttcactcTAGCCTGGGCAACAGAGGAAGACTCCACCTc*AAAAAAAAAAAGTGGGA*AGAGGATATGTA
+____________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
+HWUSI-EAS1543_614TM:3:9:10897:10640#0 + GTAAAAACACAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAaccccggaggtggag
+HWUSI-EAS1543_614TM:5:85:1510:15876#0 + GTAAAAACACAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGtAGAGGCAGGggaatcacttgaacccaggagggggcg
+HWUSI-EAS1543_614TM:6:36:6112:14330#0 + ACAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGagaatcacttgaacccaggagggggggggtggggg
+HWUSI-EAS1543_614TM:5:9:12006:6856#0 + AAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGaatcacttgaaccccggagggggagggtgggggga
+HWUSI-EAS729_614T6:4:52:7979:7186#0 + GCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGggggatcactttaaaccca*gagggggggggtgggggggggcaaggt
+HWUSI-EAS729_614T6:7:102:2873:18725#0 + GCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGgagaatcacttgaacccaggagggggagggtggggggagcgaagaa
+HWUSI-EAS729_614T6:8:57:13075:5535#0 + GCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGagaatcacttgaaacccggagggggggggtggagggagccaagaa
+HWUSI-EAS1543_614TM:5:116:13350:19305#0+ CCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGgtggagggtggggtgagccaagatc
+HWUSI-EAS729_614T6:2:14:10861:9808#0 + AGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAAtCACTTGAACCCaggagggggagggtggagtgagccaagatcat
+HWUSI-EAS1543_614TM:6:116:11405:10100#0+ CATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAAtcacttgaacccaggaggtggagggtggggggagccaagatcaagcc
+HWUSI-EAS1543_614TM:8:103:11766:8746#0 + GGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGgTGGaGTGAGCCAAGATCATGCCACTGCA
+HWUSI-EAS1543_614TM:2:81:13732:11884#0 + GTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACtTGAACCCaggaggtggaggttggagtgagccaagatcatgccaccgcac
+HWUSI-EAS729_614T6:6:27:13980:13345#0 + CTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCaggagggggaggttggagtgagccaagatcatgccactgccctctagccc
+HWUSI-EAS729_614T6:7:80:16372:13691#0 + CTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAAtCACTTGAAcCCAGGaggtggagggtggagtgagccaagatcatgccactgcactctagcct
+HWUSI-EAS1543_614TM:4:24:17238:9026#0 + TCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGtGGAGGtTGGAGtGaGCCAAGATCATGCCACTGCACTCTAGCCtGGGCaACAGAGGaAGA
+HWUSI-EAS729_614T6:4:85:2284:10568#0 + TGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTccacctc*acaacaaaa
+HWUSI-EAS729_614T6:8:88:7531:4372#0 + TGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAaaaaagtgggaa
+HWUSI-EAS1543_614TM:6:91:4071:17532#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACtCCACCtC*AaAAAAaaaaagcggga*agaggat
+HWUSI-EAS1543_614TM:7:25:16711:4592#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAagaGGGA*AGaggat
+HWUSI-EAS729_614T6:3:34:6821:21373#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGaGGAAGACTCCACCTC*AAAAAAagaaagtggga*agaggat
+HWUSI-EAS729_614T6:6:117:5396:15508#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAaaaaaagtggga*agaggat
+HWUSI-EAS1543_614TM:8:35:4484:12519#0 - CACAAAaAttAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAG
+HWUSI-EAS1543_614TM:4:40:15892:1713#0 - gggcgtgtgggtgggcacctgtaatcccagctactcGGGAGGTAGAGGCAGGAGAATCACTTGAaCCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCAT
+HWUSI-EAS1543_614TM:5:90:12703:4722#0 - aTGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCA
+HWUSI-EAS1543_614TM:3:102:1564:19915#0 - gggaggtagaggcgggaGAATCACTtgaaCcCaGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACT
+HWUSI-EAS1543_614TM:5:5:18169:18490#0 - GGTAGAGGCaGGAGaAtCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAtGAAGACTCCAC
+HWUSI-EAS729_614T6:4:44:11568:10646#0 - GGTAGAGGCAGGAGAAtcACTtgaaCCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCAC
+HWUSI-EAS729_614T6:8:38:2365:13463#0 - ggtagaggcaggagaatcacttgaaCCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCAC
+HWUSI-EAS1543_614TM:8:59:4470:6454#0 - GTAGAGGcAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACC
+HWUSI-EAS729_614T6:1:58:5165:11774#0 - agaggCAgGAGAATCaCTtgAaCCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC
+HWUSI-EAS729_614T6:6:64:8262:4344#0 - AGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC
+HWUSI-EAS1543_614TM:3:9:10897:10640#0 - AGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAA
+HWUSI-EAS1543_614TM:5:85:1510:15876#0 - aggagaatcacttgaacccgggaggtggagGTtGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAA
+HWUSI-EAS1543_614TM:7:109:3095:2034#0 - aggagaatcacttgagcccgggaggtggaggTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAA
+HWUSI-EAS1543_614TM:4:24:17238:9026#0 - tcactTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGT
+HWUSI-EAS1543_614TM:6:36:6112:14330#0 - acttgaacccaggaggtggaGGTTGGAGTGAGCCAAGATCATgCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGG
+HWUSI-EAS1543_614TM:4:67:1427:8600#0 - ttgaccccagggggtgGAGGTtGGAGTGAGCCaaGatCatgccacTGCcCTCTAGCCTGGGcAACAGAGGaAGaCTCCACCTC*AAAAAAAAAAAGTGGGA
+HWUSI-EAS1543_614TM:6:77:17645:6273#0 - ttgaacccaggaggtggaggTTGGAGtGAGCcAAgATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGA
+HWUSI-EAS1543_614TM:7:100:16004:10107#0- TTGAACCCaGGaGGTGGaGGTTGGAgTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGaCTCCACCTC*AAAAAAAAAAAGTGGGA
+HWUSI-EAS729_614T6:2:93:11989:6538#0 - tgagcccaggaggtggaGGTTGGAGTGAGCcAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGAA
+HWUSI-EAS1543_614TM:2:81:13732:11884#0 - caggaggggggggttggggtgagccaagatcatgccactgcactctagcCtGGGCAAcaGAGGAAgactccacctC*AAAAAAAAAAAGTGGGA*AGAGGAT
+HWUSI-EAS729_614T6:3:11:5262:1432#0 - caggagggggggggtggagtgagccaagatcatgcccctgcactctagcCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGA*AGAGGAT
+____________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
+
+
5 7052198 - 0 10 66301865 + 8 CTX 0 22 0 81 0.761379310344828 0.482758620689655 0 0 1 5 7052278 164 10 66301947 AGCCATGGACCTTGTGGTGGGTTCTTAACAATGGTGAGTCCGGAGTTCTTAACGATGGTGAGTCCGTAGTTTGTTCCTTCAGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTCAAAAAAAAAAAGTGGGAAGAGG
+
+
+
+5:7052198-7052198 T
+
+____________________
+____________________
+
+
+ 66301780 66301800 66301820 66301840 66301860 66301880 66301900 66301920 66301940 66301960
+ . | . | . | . | . | . | . | . | . | . |
+10:66301769-66301962 CACAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGA*AGAGGATATGTATGAGCAG*A
+
+HWI-EAS231_614TA:7:82:16847:8457#0 + CACAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGgagaatcacttgaacccagggggtggaggttggag
+HWI-EAS231_614TA:3:12:6209:13521#0 + CAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGagaatcacttgaacccaggagggggagggtggagtg
+HWI-EAS90_614M9:8:23:6627:13184#0 + GAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAaa
+HWI-EAS90_614M9:6:39:4439:17404#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAagtggga*agaggat
+HWI-EAS231_614TA:2:4:4943:8854#0 + AGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGA*AGAGGAtATGtATGAGC
+HWI-EAS231_614TA:4:78:6983:13721#0 + AGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGA*AGAGGATATGTATGAGC
+HWUSI-EAS1591_6113C:2:82:8140:13005#0 + TTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGA*AGAGGATATGTATGAGCAGA
+HWI-EAS231_614TA:4:36:3528:11249#0 + TTTGTTCCTTCAGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTCAAAAAAAAAAGtGggAAGAGGAtatgt
+HWI-EAS231_614TA:5:52:9612:8954#0 + CGTAGTTTGTTCCTTCAGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTCAAAAAAAAAAAGTGGGAAGAGG
+HWI-EAS90_614M9:1:118:15320:11068#0 + AGCCATGGACCTTGTGGTGGGTTCTTAACAATGGTGAGTCCGGAGTTCTTAACGATGGTGAGTCCGTAGTTTGTTCCTTCAGGAGTGAGCCAAGATCATG
+HWI-EAS90_614M9:3:104:3717:11676#0 + GAGTCCGTAGTTTGTTCCTTCAGGAGTGAGCCAAGATCAtGCCACtgcactccagcctgggcaacagaggaagactccacctc*caaaaaaaaaaaaggga
+HWI-EAS90_614M9:4:115:17151:16115#0 + AGTCCGGAGTTCTCAACGATGGTGAGTCCGTAGTTTGTTCCTTCAGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTC
+HWI-EAS90_614M9:5:27:7334:16234#0 + CTTGTGGTGGGTTCTTAACAATGGTGAGTCCGGAGTTCTTAACGATGGTGAGTCCGTAGTTTGTTCCTTCAGGAGTGAGCCAAGATCATGCCACTGCACT
+HWI-EAS90_614M9:7:21:19162:4940#0 + TTCAGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAgtGGga*AGaggataTGtatgagca
+HWUSI-EAS1591_6113C:1:24:6705:7145#0 + AGCCATGGACCTTGTGGTGGGTTCTTAACAATGGTGAGTCCGGAGTTCTTAACGATGGCGAGTCCGTAGTTTGTTCCTTCAGGAGTGAGCcAAGATCATG
+HWUSI-EAS1591_6113C:2:55:7685:7454#0 + AGTCCGTAGTTTGTTCCTTCAGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTCAAAAAAAAAAAGtGGGaA
+HWI-EAS292_102081428:4:40:3376:13334#0 - AAAAAtTAgcCAGGCATGGtGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGA
+HWI-EAS90_614M9:4:10:8095:5678#0 - AAAAAttagCcaGgcatGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGA
+HWI-EAS90_614M9:6:117:3466:13678#0 - gccaggggtggtggtgggcgcctgtaatcccagctactcgGGAGGTAGAGGCaGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGAT
+HWI-EAS231_614TA:6:38:5961:6062#0 - aggcgtggtggtgggcacctgtaatcccagctactcGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCAT
+HWI-EAS231_614TA:8:4:10123:13264#0 - gTAatCCCAGCTaCtCGGGAGGTAGaggcaGGAGAATCaCTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGG
+HWUSI-EAS1591_6113C:2:96:6055:2122#0 - actcgGGAGGTAGAGGCAGGAGAATCaCTTGAACCCAGGAGGTGGAaGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAA
+HWI-EAS231_614TA:7:82:16847:8457#0 - ggtggaggttggagtgagccaagatcatgccacttcactcTAGCCTGGGCAACAGAGGAAGACTCCACCTc*AAAAAAAAAAAGTGGGA*AGAGGATATGTA
+HWI-EAS231_614TA:3:12:6209:13521#0 - gaggttggagtgagcCAAGAtCaTGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGA*AGAGGATATGTATGAG
+______________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
+HWUSI-EAS1543_614TM:6:36:6112:14330#0 + ACAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGagaatcacttgaacccaggagggggggggtggggg
+HWUSI-EAS1543_614TM:5:9:12006:6856#0 + AAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGaatcacttgaaccccggagggggagggtgggggga
+HWUSI-EAS729_614T6:4:52:7979:7186#0 + GCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGggggatcactttaaaccca*gagggggggggtgggggggggcaaggt
+HWUSI-EAS729_614T6:7:102:2873:18725#0 + GCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGgagaatcacttgaacccaggagggggagggtggggggagcgaagaa
+HWUSI-EAS729_614T6:8:57:13075:5535#0 + GCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGagaatcacttgaaacccggagggggggggtggagggagccaagaa
+HWUSI-EAS1543_614TM:5:116:13350:19305#0+ CCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGgtggagggtggggtgagccaagatc
+HWUSI-EAS729_614T6:2:14:10861:9808#0 + AGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAAtCACTTGAACCCaggagggggagggtggagtgagccaagatcat
+HWUSI-EAS1543_614TM:6:116:11405:10100#0+ CATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAAtcacttgaacccaggaggtggagggtggggggagccaagatcaagcc
+HWUSI-EAS1543_614TM:8:103:11766:8746#0 + GGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGgTGGaGTGAGCCAAGATCATGCCACTGCA
+HWUSI-EAS1543_614TM:2:81:13732:11884#0 + GTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACtTGAACCCaggaggtggaggttggagtgagccaagatcatgccaccgcac
+HWUSI-EAS729_614T6:6:27:13980:13345#0 + CTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCaggagggggaggttggagtgagccaagatcatgccactgccctctagccc
+HWUSI-EAS729_614T6:7:80:16372:13691#0 + CTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAAtCACTTGAAcCCAGGaggtggagggtggagtgagccaagatcatgccactgcactctagcct
+HWUSI-EAS1543_614TM:4:24:17238:9026#0 + TCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGtGGAGGtTGGAGtGaGCCAAGATCATGCCACTGCACTCTAGCCtGGGCaACAGAGGaAGA
+HWUSI-EAS729_614T6:4:85:2284:10568#0 + TGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTccacctc*acaacaaaa
+HWUSI-EAS729_614T6:8:88:7531:4372#0 + TGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAaaaaagtgggaa
+HWUSI-EAS1543_614TM:6:91:4071:17532#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACtCCACCtC*AaAAAAaaaaagcggga*agaggat
+HWUSI-EAS1543_614TM:7:25:16711:4592#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAagaGGGA*AGaggat
+HWUSI-EAS729_614T6:3:34:6821:21373#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGaGGAAGACTCCACCTC*AAAAAAagaaagtggga*agaggat
+HWUSI-EAS729_614T6:6:117:5396:15508#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAaaaaaagtggga*agaggat
+HWUSI-EAS1543_614TM:4:21:14965:3057#0 + GTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTCAAAAAAAAAAAagaggga*agaggatatgaatgggca
+HWUSI-EAS729_614T6:4:1:12928:15062#0 + GTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAgTGGGA*AGAGGATATGTATGAGCAG
+HWUSI-EAS1543_614TM:8:35:4484:12519#0 - CACAAAaAttAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAG
+HWUSI-EAS1543_614TM:4:40:15892:1713#0 - gggcgtgtgggtgggcacctgtaatcccagctactcGGGAGGTAGAGGCAGGAGAATCACTTGAaCCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCAT
+HWUSI-EAS1543_614TM:5:90:12703:4722#0 - aTGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCA
+HWUSI-EAS1543_614TM:3:102:1564:19915#0 - gggaggtagaggcgggaGAATCACTtgaaCcCaGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACT
+HWUSI-EAS1543_614TM:5:5:18169:18490#0 - GGTAGAGGCaGGAGaAtCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAtGAAGACTCCAC
+HWUSI-EAS729_614T6:4:44:11568:10646#0 - GGTAGAGGCAGGAGAAtcACTtgaaCCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCAC
+HWUSI-EAS729_614T6:8:38:2365:13463#0 - ggtagaggcaggagaatcacttgaaCCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCAC
+HWUSI-EAS1543_614TM:8:59:4470:6454#0 - GTAGAGGcAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACC
+HWUSI-EAS729_614T6:1:58:5165:11774#0 - agaggCAgGAGAATCaCTtgAaCCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC
+HWUSI-EAS729_614T6:6:64:8262:4344#0 - AGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC
+HWUSI-EAS1543_614TM:3:9:10897:10640#0 - AGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAA
+HWUSI-EAS1543_614TM:5:85:1510:15876#0 - aggagaatcacttgaacccgggaggtggagGTtGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAA
+HWUSI-EAS1543_614TM:7:109:3095:2034#0 - aggagaatcacttgagcccgggaggtggaggTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAA
+HWUSI-EAS1543_614TM:4:24:17238:9026#0 - tcactTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGT
+HWUSI-EAS1543_614TM:6:36:6112:14330#0 - acttgaacccaggaggtggaGGTTGGAGTGAGCCAAGATCATgCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGG
+HWUSI-EAS1543_614TM:4:67:1427:8600#0 - ttgaccccagggggtgGAGGTtGGAGTGAGCCaaGatCatgccacTGCcCTCTAGCCTGGGcAACAGAGGaAGaCTCCACCTC*AAAAAAAAAAAGTGGGA
+HWUSI-EAS1543_614TM:6:77:17645:6273#0 - ttgaacccaggaggtggaggTTGGAGtGAGCcAAgATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGA
+HWUSI-EAS1543_614TM:7:100:16004:10107#0- TTGAACCCaGGaGGTGGaGGTTGGAgTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGaCTCCACCTC*AAAAAAAAAAAGTGGGA
+HWUSI-EAS729_614T6:2:93:11989:6538#0 - tgagcccaggaggtggaGGTTGGAGTGAGCcAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGAA
+HWUSI-EAS1543_614TM:2:81:13732:11884#0 - caggaggggggggttggggtgagccaagatcatgccactgcactctagcCtGGGCAAcaGAGGAAgactccacctC*AAAAAAAAAAAGTGGGA*AGAGGAT
+HWUSI-EAS729_614T6:3:11:5262:1432#0 - caggagggggggggtggagtgagccaagatcatgcccctgcactctagcCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGA*AGAGGAT
+______________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
+
+
10 66301858 + 4 4 125893225 - 1 CTX 15 28 71 81 0.735384615384615 0.5 0.889507154213037 0.243243243243243 1 10 66301777 153 4 125893154 TTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGACAGCTACTGAAAAGAGAATCATATTCTCTTTTCATAATTTAGAATATGATATATATTTCAAAATATGATA
+
+ 66301780 66301800 66301820 66301840 66301860 66301880 66301900 66301920 66301940
+ . | . | . | . | . | . | . | . | . | .
+10:66301762-66301954 GTAAAAACACAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGA*AGAGGATATGTA
+
+HWI-EAS231_614TA:7:82:16847:8457#0 + CACAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGgagaatcacttgaacccagggggtggaggttggag
+HWI-EAS231_614TA:3:12:6209:13521#0 + CAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGagaatcacttgaacccaggagggggagggtggagtg
+HWI-EAS90_614M9:8:23:6627:13184#0 + GAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAaa
+HWI-EAS90_614M9:6:39:4439:17404#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAagtggga*agaggat
+HWI-EAS292_102081428:4:40:3376:13334#0 - AAAAAtTAgcCAGGCATGGtGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGA
+HWI-EAS90_614M9:4:10:8095:5678#0 - AAAAAttagCcaGgcatGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGA
+HWI-EAS90_614M9:6:117:3466:13678#0 - gccaggggtggtggtgggcgcctgtaatcccagctactcgGGAGGTAGAGGCaGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGAT
+HWI-EAS231_614TA:6:38:5961:6062#0 - aggcgtggtggtgggcacctgtaatcccagctactcGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCAT
+HWI-EAS231_614TA:1:37:3390:19595#0 - acacaaaaattagccaggcaTGGTGGTGGGCaCCTGTAATCCCAGCTaCTcGGGAGGtAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGACAGCTACT
+HWI-EAS90_614M9:1:98:9476:6999#0 - ggTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCAgTTGAACCCAGGAGGTGACAGCTACTGAAAAGAGAATCATATTCTCTTTT
+HWI-EAS231_614TA:8:4:10123:13264#0 - gTAatCCCAGCTaCtCGGGAGGTAGaggcaGGAGAATCaCTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGG
+HWUSI-EAS1591_6113C:2:96:6055:2122#0 - actcgGGAGGTAGAGGCAGGAGAATCaCTTGAACCCAGGAGGTGGAaGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAA
+HWI-EAS90_614M9:7:110:15112:12998#0 - GTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGACAGCTACTGAAAAGAGAATCATATTCTCTTTTCATAATTTAGAATATGATATATATTTCAAAAT
+HWI-EAS90_614M9:5:12:2797:13435#0 - GCAGGAGAATCACTTGAACCCAGGAGGTGACAGCTACTGAAAAGAGAATCATATTCTCTTTTCATAATTTAGAATATGATATATATTTCAAAATATGATA
+HWI-EAS231_614TA:7:82:16847:8457#0 - ggtggaggttggagtgagccaagatcatgccacttcactcTAGCCTGGGCAACAGAGGAAGACTCCACCTc*AAAAAAAAAAAGTGGGA*AGAGGATATGTA
+____________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
+HWUSI-EAS1543_614TM:3:9:10897:10640#0 + GTAAAAACACAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAaccccggaggtggag
+HWUSI-EAS1543_614TM:5:85:1510:15876#0 + GTAAAAACACAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGtAGAGGCAGGggaatcacttgaacccaggagggggcg
+HWUSI-EAS1543_614TM:6:36:6112:14330#0 + ACAAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGagaatcacttgaacccaggagggggggggtggggg
+HWUSI-EAS1543_614TM:5:9:12006:6856#0 + AAAAATTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGaatcacttgaaccccggagggggagggtgggggga
+HWUSI-EAS729_614T6:4:52:7979:7186#0 + GCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGggggatcactttaaaccca*gagggggggggtgggggggggcaaggt
+HWUSI-EAS729_614T6:7:102:2873:18725#0 + GCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGgagaatcacttgaacccaggagggggagggtggggggagcgaagaa
+HWUSI-EAS729_614T6:8:57:13075:5535#0 + GCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGagaatcacttgaaacccggagggggggggtggagggagccaagaa
+HWUSI-EAS1543_614TM:5:116:13350:19305#0+ CCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGgtggagggtggggtgagccaagatc
+HWUSI-EAS729_614T6:2:14:10861:9808#0 + AGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAAtCACTTGAACCCaggagggggagggtggagtgagccaagatcat
+HWUSI-EAS1543_614TM:6:116:11405:10100#0+ CATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAAtcacttgaacccaggaggtggagggtggggggagccaagatcaagcc
+HWUSI-EAS1543_614TM:8:103:11766:8746#0 + GGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGgTGGaGTGAGCCAAGATCATGCCACTGCA
+HWUSI-EAS1543_614TM:2:81:13732:11884#0 + GTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACtTGAACCCaggaggtggaggttggagtgagccaagatcatgccaccgcac
+HWUSI-EAS729_614T6:6:27:13980:13345#0 + CTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCaggagggggaggttggagtgagccaagatcatgccactgccctctagccc
+HWUSI-EAS729_614T6:7:80:16372:13691#0 + CTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAAtCACTTGAAcCCAGGaggtggagggtggagtgagccaagatcatgccactgcactctagcct
+HWUSI-EAS1543_614TM:4:24:17238:9026#0 + TCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGtGGAGGtTGGAGtGaGCCAAGATCATGCCACTGCACTCTAGCCtGGGCaACAGAGGaAGA
+HWUSI-EAS729_614T6:4:85:2284:10568#0 + TGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTccacctc*acaacaaaa
+HWUSI-EAS729_614T6:8:88:7531:4372#0 + TGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAaaaaagtgggaa
+HWUSI-EAS1543_614TM:6:91:4071:17532#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACtCCACCtC*AaAAAAaaaaagcggga*agaggat
+HWUSI-EAS1543_614TM:7:25:16711:4592#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAagaGGGA*AGaggat
+HWUSI-EAS729_614T6:3:34:6821:21373#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGaGGAAGACTCCACCTC*AAAAAAagaaagtggga*agaggat
+HWUSI-EAS729_614T6:6:117:5396:15508#0 + CAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAaaaaaagtggga*agaggat
+HWUSI-EAS1543_614TM:8:35:4484:12519#0 - CACAAAaAttAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAG
+HWUSI-EAS1543_614TM:4:40:15892:1713#0 - gggcgtgtgggtgggcacctgtaatcccagctactcGGGAGGTAGAGGCAGGAGAATCACTTGAaCCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCAT
+HWUSI-EAS1543_614TM:5:90:12703:4722#0 - aTGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCA
+HWUSI-EAS1543_614TM:3:102:1564:19915#0 - gggaggtagaggcgggaGAATCACTtgaaCcCaGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACT
+HWUSI-EAS1543_614TM:5:5:18169:18490#0 - GGTAGAGGCaGGAGaAtCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAtGAAGACTCCAC
+HWUSI-EAS729_614T6:4:44:11568:10646#0 - GGTAGAGGCAGGAGAAtcACTtgaaCCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCAC
+HWUSI-EAS729_614T6:8:38:2365:13463#0 - ggtagaggcaggagaatcacttgaaCCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCAC
+HWUSI-EAS1543_614TM:8:59:4470:6454#0 - GTAGAGGcAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACC
+HWUSI-EAS729_614T6:1:58:5165:11774#0 - agaggCAgGAGAATCaCTtgAaCCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC
+HWUSI-EAS729_614T6:6:64:8262:4344#0 - AGAGGCAGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC
+HWUSI-EAS1543_614TM:3:9:10897:10640#0 - AGGAGAATCACTTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAA
+HWUSI-EAS1543_614TM:5:85:1510:15876#0 - aggagaatcacttgaacccgggaggtggagGTtGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAA
+HWUSI-EAS1543_614TM:7:109:3095:2034#0 - aggagaatcacttgagcccgggaggtggaggTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAA
+HWUSI-EAS1543_614TM:4:24:17238:9026#0 - tcactTGAACCCAGGAGGTGGAGGTTGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGT
+HWUSI-EAS1543_614TM:6:36:6112:14330#0 - acttgaacccaggaggtggaGGTTGGAGTGAGCCAAGATCATgCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGG
+HWUSI-EAS1543_614TM:4:67:1427:8600#0 - ttgaccccagggggtgGAGGTtGGAGTGAGCCaaGatCatgccacTGCcCTCTAGCCTGGGcAACAGAGGaAGaCTCCACCTC*AAAAAAAAAAAGTGGGA
+HWUSI-EAS1543_614TM:6:77:17645:6273#0 - ttgaacccaggaggtggaggTTGGAGtGAGCcAAgATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGA
+HWUSI-EAS1543_614TM:7:100:16004:10107#0- TTGAACCCaGGaGGTGGaGGTTGGAgTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGaCTCCACCTC*AAAAAAAAAAAGTGGGA
+HWUSI-EAS729_614T6:2:93:11989:6538#0 - tgagcccaggaggtggaGGTTGGAGTGAGCcAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGAA
+HWUSI-EAS1543_614TM:2:81:13732:11884#0 - caggaggggggggttggggtgagccaagatcatgccactgcactctagcCtGGGCAAcaGAGGAAgactccacctC*AAAAAAAAAAAGTGGGA*AGAGGAT
+HWUSI-EAS729_614T6:3:11:5262:1432#0 - caggagggggggggtggagtgagccaagatcatgcccctgcactctagcCTGGGCAACAGAGGAAGACTCCACCTC*AAAAAAAAAAAGTGGGA*AGAGGAT
+____________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
+
+
+ 125893140 125893160 125893180 125893200 125893220 125893240 125893260 125893280 125893300 125893320
+ . | . | . | . | . | . | . | . | . | . |
+4:125893129-125893324 CTAAAATTATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCATCAAAATTTTAAAATATCCC
+
+HWI-EAS90_614M9:6:67:12631:17754#0 + CTAAAATTATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACT
+HWUSI-EAS1591_6113C:4:21:6779:18838#0 + CTAAAATTATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACT
+HWI-EAS90_614M9:1:32:10389:16117#0 + attttGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*Cagtagctgtcactgtaggcaaaaat
+HWUSI-EAS1591_6113C:2:27:1355:21017#0 + ATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCaacaat
+HWUSI-EAS1591_6113C:1:93:16865:6051#0 + AAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCAT
+HWI-EAS90_614M9:3:78:1961:11755#0 + ATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAaAAtGCCCAtTATTC*TGC
+HWI-EAS231_614TA:5:6:17800:6620#0 + TTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATAT
+HWI-EAS90_614M9:3:95:10317:10581#0 + ATATATCATATTCTAAATTATGAAAAGAGAATATGA*Tctcttttccag*agccggccccggaagccaaaaagccccttattcttcctaattatctagggga
+HWI-EAS292_102081428:2:115:1681:2531#0 + aaattatgaaaagagaatatgattctctttt*taggagctgttactgtagggaaaaatgcccattattctgctaaatatctaagtgaggggtgcttctatt
+HWI-EAS90_614M9:1:67:1331:13704#0 + AATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATA
+HWI-EAS231_614TA:2:73:4389:2239#0 + TGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAG
+HWI-EAS231_614TA:7:88:14296:6002#0 + TGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAG
+HWI-EAS231_614TA:1:26:7960:13504#0 + TTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGT
+HWI-EAS90_614M9:7:18:13917:3050#0 + CTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCA
+HWI-EAS292_102081428:1:66:16549:20224#0+ TCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TgCTAAATATCTAAGTGATGCTTGCTACTATTAACTGaAtGTATaTACTAGCGTCAA
+HWI-EAS90_614M9:5:77:4642:6832#0 + CTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTaACTGAATGTATATACTAGCgtcaaa
+HWI-EAS90_614M9:8:34:7077:4638#0 + CTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAA
+HWUSI-EAS1591_6113C:4:9:4513:12290#0 - TTATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACCTCCTGG
+HWI-EAS231_614TA:7:49:16816:4213#0 - TATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCA
+HWI-EAS231_614TA:1:118:16290:13237#0 - ATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAAT
+HWI-EAS231_614TA:6:113:8101:10568#0 - TTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATG
+HWI-EAS90_614M9:6:48:8352:11907#0 - ATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATT
+HWUSI-EAS1591_6113C:4:39:14328:11130#0 - ATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACCTCCTGGGTTCAAGTGATTCTCCTG
+HWI-EAS90_614M9:1:85:18838:8742#0 - TTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACCTCCtgggttcaggTGATTCTCCTGCCTCTACCT
+HWI-EAS90_614M9:4:32:9763:4986#0 - TTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCaCCTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCT
+HWI-EAS231_614TA:4:103:18855:9257#0 - ATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGT
+HWI-EAS292_102081428:2:23:3944:1922#0 - aTAtaGATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGT
+HWI-EAS90_614M9:5:18:17924:10181#0 - tctctttt*cagtagctgtcccctcctGGGTTCAAGtGaTtCTCCTGCCtCTACCTCCCGAGTAGCTGGGATTACAGGTGCCCACCACCATGCCTGGCTAA
+HWUSI-EAS1591_6113C:2:117:3095:2051#0 - CTCTTTT*CAGtagctgtCaCCTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCTCCCGAGTAGCTGGGATTACAGGTGCCCACCACCATGCCTGGCTAAT
+HWI-EAS90_614M9:5:28:11006:5910#0 - TCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAA
+HWI-EAS90_614M9:5:1:12903:5535#0 - TTTT*CAGTAgCTGTCACCTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCTCCCGAGTAGCTGGGATTACAGGTGCCCACCACCATGCCTGGCTAATTTT
+HWI-EAS90_614M9:2:56:9827:2778#0 - tcactgtaGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAAATTTTAAAATATCC
+HWI-EAS90_614M9:4:47:5786:19724#0 - tcactGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAAATTTTAAAATATCC
+HWI-EAS90_614M9:6:67:11670:7706#0 - TCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAAATTTTAAAATATCC
+HWUSI-EAS1591_6113C:1:117:14883:9756#0 - TCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAAATTTTAAAATATCC
+HWI-EAS90_614M9:6:37:7790:16256#0 - CACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAAATTTTAAAATATCCC
+_______________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
+HWUSI-EAS729_614T6:8:86:3208:16885#0 + TTATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGC
+HWUSI-EAS729_614T6:3:97:10823:18653#0 + ATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAA
+HWUSI-EAS729_614T6:4:113:3054:5730#0 + TGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGgcaag
+HWUSI-EAS729_614T6:2:63:5049:3137#0 + AATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAA
+HWUSI-EAS729_614T6:5:95:4777:9034#0 + AATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAA
+HWUSI-EAS729_614T6:8:79:9632:17904#0 + AATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAA
+HWUSI-EAS1543_614TM:8:50:14524:11101#0 + ATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAAT
+HWUSI-EAS729_614T6:3:79:10477:10666#0 + TTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCC
+HWUSI-EAS729_614T6:4:80:3846:17138#0 + TTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCC
+HWUSI-EAS729_614T6:4:83:3398:15440#0 + TTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCC
+HWUSI-EAS729_614T6:5:37:19001:5491#0 + TGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCC
+HWUSI-EAS729_614T6:8:113:1681:2586#0 + TCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCT
+HWUSI-EAS729_614T6:8:3:18130:10917#0 + TATTTTGAAATATATATCATATtCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAA
+HWUSI-EAS729_614T6:7:32:3784:11386#0 + TGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCT
+HWUSI-EAS729_614T6:7:117:8731:4955#0 + ATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGA
+HWUSI-EAS729_614T6:8:51:17611:11416#0 + TATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGAT
+HWUSI-EAS1543_614TM:1:102:4785:3367#0 + TATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCT
+HWUSI-EAS1543_614TM:6:60:2952:14857#0 + TATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCT
+HWUSI-EAS729_614T6:5:55:11097:17856#0 + ATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTG
+HWUSI-EAS729_614T6:2:8:16102:20401#0 + AAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAAT
+HWUSI-EAS1543_614TM:1:75:3691:5075#0 + GATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGC
+HWUSI-EAS729_614T6:8:52:15200:10784#0 + CACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCATCAAAattntAAAATATCCC
+HWUSI-EAS1543_614TM:7:17:14845:19479#0 - ATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATAATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAAT
+HWUSI-EAS729_614T6:3:13:15991:15051#0 - TATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TG
+HWUSI-EAS729_614T6:7:65:3141:2489#0 - TATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TG
+HWUSI-EAS729_614T6:8:54:5081:15134#0 - TATCATATTTTcAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TG
+HWUSI-EAS1543_614TM:4:38:16709:14763#0 - CATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTA
+HWUSI-EAS1543_614TM:8:50:14524:11101#0 - TATTCTAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCT
+HWUSI-EAS729_614T6:1:56:9586:15810#0 - TAAATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATaTAAGTGATGCTTGCTACTAT
+HWUSI-EAS729_614T6:1:37:11212:13638#0 - TAAATTATGAAtAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTAT
+HWUSI-EAS729_614T6:1:29:1825:18616#0 - AATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTA
+HWUSI-EAS729_614T6:5:48:14044:12231#0 - AATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTA
+HWUSI-EAS729_614T6:4:115:6035:15619#0 - ATTATGAAAAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAA
+HWUSI-EAS729_614T6:5:72:1269:20227#0 - aagagaatatgattctctTtT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGT
+HWUSI-EAS729_614T6:7:32:3784:11386#0 - AAGAGAATATGATTCTCTTTT*CAGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGT
+HWUSI-EAS1543_614TM:1:102:4785:3367#0 - AGTAGCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGCTACTATTAACTGAATGTATATACTAGCGTCAAAATTTTA
+HWUSI-EAS1543_614TM:6:60:2952:14857#0 - agtagCTGTCACTGTAGGCAAAAATGCCCATTATTC*TGCTAAATATCTAAGTGATGCTTGctactATTAACTGAATGTATATACTAGCGTCAAAATTTTA
+_______________________________________________________________________________________________________________________________________________________________________________________________________________________________________________
+
+
\ No newline at end of file
diff -r 000000000000 -r acc8d8bfeb9a example/tumor.predSV.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/example/tumor.predSV.txt Wed Feb 08 16:59:24 2012 -0500
@@ -0,0 +1,3 @@
+4 125893227 + 5 10 66301858 - 4 CTX 29 14 83 71 0.895173453996983 0.230769230769231 0.735384615384615 0.5 1 4 125893135 176 10 66301773 TTATGAATTTTGAAATATATATCATATTTTGAAATATATATCATATTCTAAATTATGAAAAGAGAATATGATTCTCTTTTCAGTAGCTGTCACCTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCTCCCGAGTAGCTGGGATTACAGGTGCCCACCACCATGCCTGGCTAATTTT
+5 7052198 - 0 10 66301865 + 8 CTX 0 22 0 81 0.761379310344828 0.482758620689655 0 0 1 5 7052278 164 10 66301947 AGCCATGGACCTTGTGGTGGGTTCTTAACAATGGTGAGTCCGGAGTTCTTAACGATGGTGAGTCCGTAGTTTGTTCCTTCAGGAGTGAGCCAAGATCATGCCACTGCACTCTAGCCTGGGCAACAGAGGAAGACTCCACCTCAAAAAAAAAAAGTGGGAAGAGG
+10 66301858 + 4 4 125893225 - 1 CTX 15 28 71 81 0.735384615384615 0.5 0.889507154213037 0.243243243243243 1 10 66301777 153 4 125893154 TTAGCCAGGCATGGTGGTGGGCACCTGTAATCCCAGCTACTCGGGAGGTAGAGGCAGGAGAATCACTTGAACCCAGGAGGTGACAGCTACTGAAAAGAGAATCATATTCTCTTTTCATAATTTAGAATATGATATATATTTCAAAATATGATA
diff -r 000000000000 -r acc8d8bfeb9a example/tumor.sclip.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/example/tumor.sclip.txt Wed Feb 08 16:59:24 2012 -0500
@@ -0,0 +1,27 @@
+4 125892327 + HWUSI-EAS1591_6113C:3:17:12332:19420#0 CC CC
+4 125892458 + HWUSI-EAS1591_6113C:4:91:6281:9961#0 GACTAACCACCACGGTACATGTTTTCCTATGTAAAAAACCTGCACATTCTACACATGTATCCCAGAACTTAAAGTAAAACAC B@C@?:CC>CCBCCCCACBCDCCCCCC;<:<9CCCCC@CCCCCBCCCCCCCCCCCCCCCACCCCCCCCCCCCCCCCCCCCAC
+4 125893225 + HWI-EAS90_614M9:5:18:17924:10181#0 CCCTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCTCCCGAGTAGCTGGGATTACAGGTGCCCACCACCATGCCTGGCTAA #######@@7@:8@><16+6(B>AABCAA3AB@CC6CCCCCCCDCCCCCCBCCDCCCCCCCCCCCCDCCCCCCCCCCCCCC
+4 125893227 + HWI-EAS90_614M9:4:32:9763:4986#0 CTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCT CCCC@CCCCCCCCCCCCCCCCCCCCCCCCCCCCC
+4 125893227 + HWI-EAS90_614M9:5:1:12903:5535#0 CTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCTCCCGAGTAGCTGGGATTACAGGTGCCCACCACCATGCCTGGCTAATTTT B;@@>DABBB>DBAAAB@BABDDDDCCBCDACCBCCCCCCCBCCCCCCCCCCCCCCCCCCCCCCBCCCCCCCCCCCCCCCCCC
+4 125893227 + HWUSI-EAS1591_6113C:2:117:3095:2051#0 CTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCTCCCGAGTAGCTGGGATTACAGGTGCCCACCACCATGCCTGGCTAAT 65977A::??BB>?BA8???=BDCCCCBCCCCCBCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCDCCCCCCCCCCCCCCC
+4 125893227 + HWUSI-EAS1591_6113C:4:39:14328:11130#0 CTCCTGGGTTCAAGTGATTCTCCTG CCCCCCCCCCCCCCCCCCCCCCCCC
+4 125893227 + HWUSI-EAS1591_6113C:4:9:4513:12290#0 CTCCTGG CCCCCCC
+4 125893365 - HWI-EAS292_102081428:1:66:17910:13190#0 TAACTACTT GFGFDGGGG
+4 125893979 - HWUSI-EAS1591_6113C:1:81:16173:11673#0 ATAAAATATTTATATTTTATATATTTGTATTCAAT CCCCCCCCCCCCCCCCCCCCCCCCCCCDCCCCCCC
+10 66301086 - HWI-EAS231_614TA:6:108:5679:13132#0 AAACATGGTTCAAATAAATTATTTATAAAGAAG CCCCCCCCCCCCCCCCCBBCCCCCCCCCCCCCC
+10 66301086 - HWI-EAS292_102081428:1:76:17776:13382#0 AAACATGGTTCAAATAAATTATTTATAAAGAAG GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG
+10 66301858 + HWI-EAS231_614TA:1:37:3390:19595#0 ACAGCTACT CCCCCCCCC
+10 66301858 + HWI-EAS90_614M9:1:98:9476:6999#0 ACAGCTACTGAAAAGAGAATCATATTCTCTTTT CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
+10 66301858 + HWI-EAS90_614M9:5:12:2797:13435#0 ACAGCTACTGAAAAGAGAATCATATTCTCTTTTCATAATTTAGAATATGATATATATTTCAAAATATGATA CCCCCCDCCCCCCCCCCCCACCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
+10 66301858 + HWI-EAS90_614M9:7:110:15112:12998#0 ACAGCTACTGAAAAGAGAATCATATTCTCTTTTCATAATTTAGAATATGATATATATTTCAAAAT CCCCCCBCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
+10 66301865 - HWI-EAS231_614TA:4:36:3528:11249#0 TTTGTTCCTTCA CCCCCCCCCCCC
+10 66301865 - HWI-EAS231_614TA:5:52:9612:8954#0 CGTAGTTTGTTCCTTCA CCCCCCCCCCCCCCCCC
+10 66301865 - HWI-EAS90_614M9:1:118:15320:11068#0 AGCCATGGACCTTGTGGTGGGTTCTTAACAATGGTGAGTCCGGAGTTCTTAACGATGGTGAGTCCGTAGTTTGTTCCTTCA CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCACCCDCCCCCCCCCCCCCCCC@CCBCCCDDCCDBCDBCCDBCBCDDD
+10 66301865 - HWI-EAS90_614M9:3:104:3717:11676#0 GAGTCCGTAGTTTGTTCCTTCA BBBBBB@9B@878;69?B9AB>
+10 66301865 - HWI-EAS90_614M9:4:115:17151:16115#0 AGTCCGGAGTTCTCAACGATGGTGAGTCCGTAGTTTGTTCCTTCA CCCCCCCCCCCCCACCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
+10 66301865 - HWI-EAS90_614M9:5:27:7334:16234#0 CTTGTGGTGGGTTCTTAACAATGGTGAGTCCGGAGTTCTTAACGATGGTGAGTCCGTAGTTTGTTCCTTCA CCCCCCCDCCCCCCCCCCCCCBCCBCCCDCCCCCCBCCCCCCCCC@CC>CAC@CCBABCADDDBDBBBBBB
+10 66301865 - HWI-EAS90_614M9:7:21:19162:4940#0 TTCA CCCC
+10 66301865 - HWUSI-EAS1591_6113C:1:24:6705:7145#0 AGCCATGGACCTTGTGGTGGGTTCTTAACAATGGTGAGTCCGGAGTTCTTAACGATGGCGAGTCCGTAGTTTGTTCCTTCA CCCCCCCCCCCCCCCCCDCCCCCCCCCCCCCCCCDCCCDCCCCCCACCCCCCDCCDC98A@C 33;
+
+my $rmdup = 0;
+my $print_read = 1;
+$use_scratch = 0;
+$min_percent_id = 90;
+my ($work_dir, $out_dir);
+my $ref_genome;
+# input/output
+my ($out_prefix, $range, $input_bam );
+my $out_suffix = "sclip.txt";
+my $paired = 1;
+my ( $help, $man, $version, $usage );
+my $optionOK = GetOptions(
+ 'i|in|input=s' => \$input_bam,
+ 'o|out_dir=s' => \$out_dir,
+ 'ref_genome=s' => \$ref_genome,
+ 'p=s' => \$out_prefix,
+ 'scratch!' => \$use_scratch,
+ 'paired!' => \$paired,
+ 'rmdup!' => \$rmdup,
+ 'lq_cutoff=i' => \$lowqual_cutoff,
+ 'min_pct_id=i' => \$min_percent_id,
+ 'min_pct_hq=i' => \$min_percent_hq,
+ 'print_read!' => \$print_read,
+ 'r|range=s' => \$range,
+ 'h|help|?' => \$help,
+ 'man' => \$man,
+ 'usage' => \$usage,
+ 'v|version' => \$version,
+);
+
+pod2usage(-verbose=>2) if($man or $usage);
+pod2usage(1) if($help or $version );
+
+my $start_dir = getcwd;
+if($input_bam) {
+ croak "The bam file you specified does not exist!\n" unless(-e $input_bam);
+ $input_bam = File::Spec->rel2abs($input_bam);
+}
+else{
+ croak "You must specify the input bam file or sample name";
+}
+croak "You must provide the reference genome in fasta format!" if(!$ref_genome);
+croak "The reference genome file you speicified does not exist!\n"
+ unless(-e $ref_genome);
+
+my $input_base = fileparse($input_bam);
+
+#setup output dir and workind directory
+$out_dir = getcwd if(!$out_dir);
+mkdir $out_dir if(!-e $out_dir || ! -d $out_dir);
+$work_dir = get_work_dir() if($use_scratch);
+$work_dir = $out_dir if(!$work_dir);
+$use_scratch = undef if($work_dir eq $out_dir);
+chdir($work_dir);
+
+# figure out output prefix
+$out_prefix = $input_base if(!$out_prefix);
+
+my $sam = Bio::DB::Sam->new( -bam => $input_bam, -fasta => $ref_genome);
+
+my $output_file;
+my $validator = SCValidator->new();
+if(!$paired) {
+ $validator->remove_validator("strand_validator");
+}
+if($range) {
+ my ($chr, $start, $end) = parse_range($range);
+ my $tmp = $chr;
+ $tmp = $tmp . ".$start" if($start);
+ $tmp = $tmp . ".$end" if($end);
+ $output_file = join('.', $out_prefix, $tmp, $out_suffix);
+ $output_file = File::Spec->catfile($out_dir, $output_file);
+
+ my($pcover, $ncover) = extract_range_sclip(
+ -SAM => $sam,
+ -RANGE => $range,
+ -WORK_DIR => $work_dir,
+ -OUTPUT => $output_file,
+ -VALIDATOR => $validator);
+ $output_file = join('.', $out_prefix, $tmp, "cover");
+ $output_file = File::Spec->catfile($out_dir, $output_file);
+ warn "$output_file file exists and it will be replaced!" if(-e $output_file);
+ $rmdup = 0;
+ #$start = 0 unless $start;
+ #my ($tid) = $sam->header->parse_region($chr);
+ #my $chr_len = $sam->header->target_len->[$tid];
+ #$end = $chr_len unless $end;
+
+ #my ($coverage) = $sam->features(-type => 'coverage', -seq_id => $chr,
+ # -start => $start, -end => $end);
+ #my @cc = $coverage->coverage;
+
+ $rmdup = 0;
+ my($fh, $fname) = tempfile( DIR => $work_dir);
+ foreach my $p (keys(%{$pcover})) {
+ my $c = ($rmdup ? scalar(@{$pcover->{$p}}) : $pcover->{$p});
+ print $fh join("\t", $chr, $p, "+", $c, count_coverage($sam, $chr, $p)), "\n";
+ }
+ foreach my $p (keys(%{$ncover})) {
+ my $c = ($rmdup ? scalar(@{$ncover->{$p}}) : $ncover->{$p});
+ print $fh join("\t", $chr, $p, "-", $c , count_coverage($sam, $chr, $p)), "\n";
+ }
+ system("sort -k 2 -n $fname -o $output_file");
+ system("rm $fname");
+}
+else{
+ $output_file = join('.', $out_prefix, $out_suffix);
+ $output_file = File::Spec->catfile($out_dir, $output_file);
+ my $header = $sam->header;
+ my $target_names = $header->target_name;
+ $output_file = join('.', $out_prefix, "cover");
+ my $read_file = join('.', $out_prefix, "sclip.txt");
+ $output_file = File::Spec->catfile($out_dir, $output_file);
+ $read_file = File::Spec->catfile($out_dir, $read_file);
+
+ if(-e $output_file) {
+ warn "$output_file file exists and it will be replaced!";
+ system("rm $output_file");
+ }
+ my @t_names = uniq @{$target_names};
+ foreach my $chr (@t_names){
+ my($fh, $fname) = tempfile( DIR => $work_dir);
+ my($pcover, $ncover) = extract_range_sclip(
+ -SAM => $sam,
+ -RANGE =>$chr,
+ -WORK_DIR => $work_dir,
+ -OUTPUT => $read_file,
+ -VALIDATOR => $validator);
+ foreach my $p (keys(%{$pcover})) {
+ my $c = ($rmdup ? scalar(@{$pcover->{$p}}) : $pcover->{$p});
+ print $fh join("\t", $chr, $p, "+", $c, count_coverage($sam, $chr, $p) ), "\n";
+ }
+ foreach my $p (keys(%{$ncover})) {
+ my $c = ($rmdup ? scalar(@{$ncover->{$p}}) : $ncover->{$p});
+ print $fh join("\t", $chr, $p, "-", $c, count_coverage($sam, $chr, $p)), "\n";
+ }
+ system("sort -k 2 -n $fname -o $fname.sorted");
+ system("cat $fname.sorted >> $output_file");
+ system("rm $fname");
+ system("rm $fname.sorted");
+ }
+}
+chdir $start_dir;
+exit(0);
+
+sub count_coverage {
+ my ($sam, $chr, $pos, $clip) = @_;
+ if($rmdup) {
+ my @pairs;
+ my $seg = $sam->segment(-seq_id => $chr, -start => $pos, -end => $pos);
+ my $n = 0;
+ my $itr = $seg->features(-iterator => 1);
+ while( my $a = $itr->next_seq) {
+ my $sclip_len = 0;
+ if($clip) {
+ my @cigar_array = @{$a->cigar_array};
+ $sclip_len = $cigar_array[0]->[1] if($cigar_array[0]->[0] eq 'S' && $clip == RIGHT_CLIP);
+ $sclip_len = $cigar_array[$#cigar_array]->[1] if($cigar_array[$#cigar_array]->[0] eq 'S' && $clip == RIGHT_CLIP);
+ }
+ next if(@pairs > 0 && is_PCR_dup($a, \@pairs, $sclip_len));
+ $n++;
+ #return $n if( $n > $max_repetitive_cover);
+ if($a->mpos) {
+ push @pairs, [$a->start, $a->end, $a->mate_start, $a->mate_end, $sclip_len];
+ }
+ else {
+ push @pairs, [$a->start, $a->end, 0, 0, $sclip_len];
+ }
+
+ }
+ return $n;
+ }
+ else{
+ my ($c) = $sam->features(-type => 'coverage', -seq_id=> $chr,
+ -start => $pos, -end => $pos);
+ my @c_d = $c->coverage;
+ return $c_d[0];
+ }
+}
+
+sub extract_range_sclip {
+ my %arg = @_;
+ my $sam = $arg{-SAM} || croak "missing -SAM";
+ my $range = $arg{-RANGE} || croak "missing -RANGE";
+ my $output_file = $arg{-OUTPUT} || croak "missing -OUTPUT";
+ my $validator = $arg{-VALIDATOR} || croak "missing -VALIDATOR";
+
+ my($fh, $fname) = tempfile( DIR => $work_dir);
+ my (%plus_cover, %neg_cover);
+ $sam->fetch($range,
+ sub {
+ my $a = shift;
+ my $cigar_str = $a->cigar_str;
+# print STDERR $a->qname, "\t", $cigar_str, "\n";
+ my @cigar_array = @{$a->cigar_array};
+ return if($a->cigar_str !~ m/S/);
+ if($paired && !$a->proper_pair) { #paired but mate is not mapped
+ $validator->remove_validator("strand_validator");
+ }
+ #return if(!$a->proper_pair && $paired);
+ #return if($paired && !$a->mpos);
+ my ($sclip_len, $ort, $pos, $seq, $qual_str, $qual);
+ $qual_str = join( "", (map { chr $_ + FQ_BASE_NUMBER } $a->qscore));
+ if($cigar_array[0]->[0] eq 'S' && $validator->validate($a, LEFT_CLIP) ) {
+ $sclip_len = $cigar_array[0]->[1]; $ort = "-"; $pos = $a->start;
+ $seq = substr($a->query->dna, 0, $sclip_len );
+ $qual = substr($qual_str, 0, $sclip_len);
+
+ my $print = 1;
+ if($rmdup) {
+ if(exists $neg_cover{$pos}) {
+ $print = 0 if(is_PCR_dup($a, $neg_cover{$pos}, $sclip_len));
+ }
+ else {
+ $neg_cover{$pos} = [];
+ }
+ if($print == 1) {
+ if($a->mpos) {
+ push @{$neg_cover{$pos}}, [$a->start, $a->end, $a->mate_start, $a->mate_end, $sclip_len];
+ }
+ else {
+ push @{$neg_cover{$pos}}, [$a->start, $a->end, 0, 0, $sclip_len];
+ }
+ }
+ }
+ else {
+ if(exists $neg_cover{$pos}) {
+ $neg_cover{$pos}++;
+ }
+ else{
+ $neg_cover{$pos} = 1;
+ }
+ }
+ print $fh join("\t", $a->seq_id, $pos, $ort, $a->qname, $seq, $qual), "\n"
+ if($print_read && $print == 1);
+ }
+
+ if($cigar_array[$#cigar_array]->[0] eq 'S' && $validator->validate($a, RIGHT_CLIP)) {
+ $sclip_len = $cigar_array[$#cigar_array]->[1]; $ort = '+'; $pos = $a->end;
+ my $l = length($a->query->dna);
+ $seq = substr($a->query->dna, $l - $sclip_len);
+ $qual = substr($qual_str, $l - $sclip_len);
+
+ my $print = 1;
+ if($rmdup) {
+ if(exists $plus_cover{$pos}) {
+ $print = 0 if(is_PCR_dup($a, $plus_cover{$pos}, $sclip_len));
+ }
+ else {
+ $plus_cover{$pos} = [];
+ }
+ if($print ==1 ) {
+ if($a->mpos) {
+ push @{$plus_cover{$pos}}, [$a->start, $a->end, $a->mate_start, $a->mate_end, $sclip_len];
+ }
+ else {
+ push @{$plus_cover{$pos}}, [$a->start, 0, 0, $a->mate_end, $sclip_len];
+ }
+ }
+ }
+ else {
+ if(exists $plus_cover{$pos}) {
+ $plus_cover{$pos}++;
+ }
+ else{
+ $plus_cover{$pos} = 1;
+ }
+ }
+ $print = is_PCR_dup($a, $plus_cover{$pos}, $sclip_len) if($rmdup);
+ print $fh join("\t", $a->seq_id, $pos, $ort, $a->qname, $seq, $qual), "\n"
+ if($print_read && $print);
+ }
+ if($paired && !$a->proper_pair) { #paired but mate is not mapped, add back
+ $validator->add_validator("strand_validator");
+ }
+ }
+ );
+ if($print_read) {
+ system("sort -k 2 -n $fname -o $fname.sorted");
+ system("cat $fname.sorted >> $output_file");
+ system("rm $fname");
+ system("rm $fname.sorted");
+ }
+ return(\%plus_cover, \%neg_cover);
+}
+
+
+=head1 NAME
+
+extractSClip.pl - extract positions with soft clipped read in bam file.
+
+
+=head1 VERSION
+
+This documentation refers to extractSClip.pl version 0.0.1.
+
+
+=head1 USAGE
+
+ # extract all positions with soft clipped reads in whole genome:
+ ./extractSClip.pl -i sample.bam -g hg18.fa
+ # extract chr1 positions with soft clipped reads
+ ./extractSClip.pl -i sample.bam -g hg18.fa -r chr1
+
+
+=head1 REQUIRED ARGUMENTS
+
+ -i: Input bam file.
+ --ref_genome: The genome file in fa file, must be the same used to map reads.
+
+
+=head1 OPTIONS
+
+ -r: The range of positions need to be extracted. Format: chr1 or chr1:500-5000.
+ -o: The output directory, default is current directory.
+ --scratch: use scracth space, default off.
+ --rmdup: remove PCR dumplicate reads, default on, use --normdup to turn it off.
+ --lq_cutoff: low quality cutoff value, default 20.
+ --min_pct_id: minimum percent identify for the aligned high qual part,default 90.
+ --min_pct_hq: minimum percent high quality for soft clipped part, default 80.
+ --print_read: individual soft-clipped read will be printed, default off.
+ -h, --help: The help page.
+ --man: Print the man page.
+ --usage: Print usage information.
+ -v, --version: print version information.
+
+
+=head1 DESCRIPTION
+
+This is a program to extract all soft-clipped positions such that for each position
+a list of requirements need to be satisfied. More specifically, the orientaion of
+pair-end read should be satisfied, the minimum percent identify of aligned part
+need to be satisfied, and the minimum percent of hiqh quality soft-clipped part
+should be satisfied.
+
+
+=head1 DEPENDENCIES
+
+The program depend on several packages:
+1. Bioperl perl module.
+2. Bio::DB::Sam, version 1.5 or later, it requires samtools lib installed.
+
+=head1 BUGS AND LIMITATIONS
+
+There are no known bugs in this module, but the method is limitted to bam file
+that has soft-clipping cigar string generated.Please report problems to
+Jianmin Wang (Jianmin.Wang@stjude.org)
+Patches are welcome.
+
+=head1 AUTHOR
+
+Jianmin Wang (Jianmin.Wang@stjude.org)
+
+
+=head1 LICENCE AND COPYRIGHT
+
+Copyright (c) 2010 by St. Jude Children's Research Hospital.
+
+This program is free software: you can redistribute it and/or modify
+it under the terms of the GNU General Public License as published by
+the Free Software Foundation, either version 2 of the License, or
+(at your option) any later version.
+
+This program is distributed in the hope that it will be useful,
+but WITHOUT ANY WARRANTY; without even the implied warranty of
+MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+GNU General Public License for more details.
+
+You should have received a copy of the GNU General Public License
+along with this program. If not, see .
diff -r 000000000000 -r acc8d8bfeb9a germline.bam
Binary file germline.bam has changed
diff -r 000000000000 -r acc8d8bfeb9a germline.bam.bai
Binary file germline.bam.bai has changed
diff -r 000000000000 -r acc8d8bfeb9a getSomaticPindel.pl
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/getSomaticPindel.pl Wed Feb 08 16:59:24 2012 -0500
@@ -0,0 +1,142 @@
+#!/usr/bin/perl -w
+
+use strict;
+use Carp;
+use Getopt::Long;
+use English;
+use Pod::Usage;
+use Data::Dumper;
+
+my ( $help, $man, $version, $usage );
+my ($infile, $outfile);
+my $optionOK = GetOptions(
+ 'h|help|?' => \$help,
+ 'man' => \$man,
+ 'usage' => \$usage,
+ 'v|version' => \$version,
+ 'i|input=s' => \$infile,
+ 'o|output=s' => \$outfile,
+);
+
+pod2usage(2) if($man);
+pod2usage( -verbose => 99, -sections => "USAGE|REQUIRED ARGUMENTS|OPTIONS" )
+ if($help or $usage or !$infile or !$outfile);
+pod2usage( -verbose => 99, -sections => "VERSION") if($version);
+
+open my $IN "<$infile" or croak "can't open $infile:$OS_ERROR";
+while(my $line = <$IN>) {
+
+}
+=head1 NAME
+
+ -
+
+
+=head1 VERSION
+
+The initial template usually just has:
+
+This documentation refers to version 0.0.1.
+
+
+=head1 USAGE
+
+ # Brief working invocation example(s) here showing the most common usage(s)
+
+ # This section will be as far as many users ever read,
+ # so make it as educational and exemplary as possible.
+=head1 REQUIRED ARGUMENTS
+
+A complete list of every argument that must appear on the command line.
+when the application is invoked, explaining what each of them does, any
+restrictions on where each one may appear (i.e., flags that must appear
+before or after filenames), and how the various arguments and options
+may interact (e.g., mutual exclusions, required combinations, etc.)
+
+If all of the application's arguments are optional, this section
+may be omitted entirely.
+
+=head1 OPTIONS
+
+A complete list of every available option with which the application
+can be invoked, explaining what each does, and listing any restrictions,
+or interactions.
+
+If the application has no options, this section may be omitted entirely.
+
+
+=head1 DESCRIPTION
+
+A full description of the application and its features.
+May include numerous subsections (i.e., =head2, =head3, etc.).
+
+
+=head1 DIAGNOSTICS
+
+A list of every error and warning message that the application can generate
+(even the ones that will "never happen"), with a full explanation of each
+problem, one or more likely causes, and any suggested remedies. If the
+application generates exit status codes (e.g., under Unix), then list the exit
+status associated with each error.
+
+=head1 CONFIGURATION AND ENVIRONMENT
+
+A full explanation of any configuration system(s) used by the application,
+including the names and locations of any configuration files, and the
+meaning of any environment variables or properties that can be set. These
+descriptions must also include details of any configuration language used.
+
+
+=head1 DEPENDENCIES
+
+A list of all the other modules that this module relies upon, including any
+restrictions on versions, and an indication of whether these required modules are
+part of the standard Perl distribution, part of the module's distribution,
+or must be installed separately.
+
+
+=head1 INCOMPATIBILITIES
+
+A list of any modules that this module cannot be used in conjunction with.
+This may be due to name conflicts in the interface, or competition for
+system or program resources, or due to internal limitations of Perl
+(for example, many modules that use source code filters are mutually
+incompatible).
+
+
+=head1 BUGS AND LIMITATIONS
+
+A list of known problems with the module, together with some indication of
+whether they are likely to be fixed in an upcoming release.
+
+Also a list of restrictions on the features the module does provide:
+data types that cannot be handled, performance issues and the circumstances
+in which they may arise, practical limitations on the size of data sets,
+special cases that are not (yet) handled, etc.
+
+The initial template usually just has:
+
+There are no known bugs in this module.
+Please report problems to ()
+Patches are welcome.
+
+=head1 AUTHOR
+
+ ()
+
+
+
+=head1 LICENCE AND COPYRIGHT
+
+Copyright (c) (). All rights reserved.
+
+followed by whatever licence you wish to release it under.
+For Perl code that is often just:
+
+This module is free software; you can redistribute it and/or
+modify it under the same terms as Perl itself. See L.
+
+This program is distributed in the hope that it will be useful,
+but WITHOUT ANY WARRANTY; without even the implied warranty of
+MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+
diff -r 000000000000 -r acc8d8bfeb9a getUniqSV.pl
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/getUniqSV.pl Wed Feb 08 16:59:24 2012 -0500
@@ -0,0 +1,222 @@
+#!/nfs_exports/apps/64-bit/gnu-apps/perl5.8.9/bin/perl -w
+
+use strict;
+use Carp;
+use Getopt::Long;
+use English;
+use Pod::Usage;
+use Data::Dumper;
+use DBI;
+use DBD::mysql;
+
+my $db = 'hg18';
+my $port = 3306;
+my $host = '10.4.19.125';
+my $USER = "sjXuser";
+
+my ( $help, $man, $version, $usage );
+my ($in, $normal, $out);
+my $optionOK = GetOptions(
+ 'd|i|in|input=s' => \$in,
+ 'g|normal=s' => \$normal,
+ 'o|out|output=s' => \$out,
+ 'h|help|?' => \$help,
+ 'man' => \$man,
+ 'usage' => \$usage,
+ 'v|version' => \$version,
+);
+
+pod2usage(2) if($man);
+pod2usage( -verbose => 99, -sections => "USAGE|REQUIRED ARGUMENTS|OPTIONS" )
+ if($help or $usage);
+pod2usage( -verbose => 99, -sections => "VERSION") if($version);
+
+croak "Missing input file" if(!$in);
+open STDOUT, ">$out" if($out);
+open my $IN, "<$in" or croak "can't open $in: $OS_ERROR";
+open my $NORMAL, "<$normal" or croak "can't open $normal: $OS_ERROR" if($normal);
+
+my $dbh = DBI->connect("dbi:mysql:$db:$host:$port", $USER) or
+ die "Can't connect to MySQL database: $DBI::errstr\n";
+
+my @SV;
+my @normal_SV;
+if($normal) {
+ while(my $line = <$NORMAL>) {
+ chomp $line;
+ my ($chrA, $posA, $ortA, $countA, $chrB, $posB, $ortB, $countB, $type) = split /\t/, $line;
+ next if( search_sv(\@normal_SV, $chrA, $posA, $chrB, $posB, $type));
+ push @normal_SV, [$chrA, $posA, $chrB, $posB, $type];
+ }
+}
+while( my $line = <$IN> ) {
+ chomp $line;
+ my ($chrA, $posA, $ortA, $countA, $chrB, $posB, $ortB, $countB, $type) = split /\t/, $line;
+ next if( search_sv(\@normal_SV, $chrA, $posA, $chrB, $posB, $type)); # the SV exists in normal sample
+ push @normal_SV, [$chrA, $posA, $chrB, $posB, $type];
+ print STDOUT $line;
+ if($type eq 'DEL') {
+ my $tmp = get_gene_info($chrA, $posA, $posB);
+ print STDOUT "\t", $tmp if($tmp);
+ }
+ else {
+ print STDOUT "\t", get_gene_info($chrA, $posA, $posA);
+ print STDOUT "\t", get_gene_info($chrB, $posB, $posB);
+ }
+ print STDOUT "\n";
+}
+close $IN;
+close STDOUT if($out);
+$dbh->disconnect;
+exit(0);
+
+sub search_sv {
+ my ($r_SV, $chrA, $posA, $chrB, $posB, $type) = @_;
+
+ foreach my $sv (@{$r_SV}) {
+ return 1 if( $sv->[0] eq $chrA && abs($sv->[1] - $posA) < 20 &&
+ $sv->[2] eq $chrB && abs($sv->[3] - $posB) < 20 && $sv->[4] eq $type);
+ return 1 if( $sv->[0] eq $chrB && abs($sv->[1] - $posB) < 20 &&
+ $sv->[2] eq $chrA && abs($sv->[3] - $posA) < 20 && $sv->[4] eq $type);
+ return 1 if( $sv->[0] eq $chrA && abs($sv->[1] - $posA) < 100 &&
+ $sv->[2] eq $chrB && abs($sv->[3] - $posB) < 100 && $sv->[4] eq $type && $type eq 'CTX');
+ return 1 if( $sv->[0] eq $chrB && abs($sv->[1]- $posB) < 100 &&
+ $sv->[2] eq $chrA && abs($sv->[3] - $posA) < 100 && $sv->[4] eq $type && $type eq 'CTX');
+ }
+ return;
+}
+
+sub get_gene_info {
+ my ($chr, $st, $ed) = @_;
+ $chr = 'chr' . $chr if($chr !~ m/chr/);
+ my $rtn="";
+ my $sth = $dbh->prepare("select distinct geneName
+ from refFlat
+ where (chrom = '$chr' AND $ed >= txStart AND $st <= txEnd)
+ order by txStart"
+ );
+ $sth->execute();
+ my $tbl_ref = $sth->fetchall_arrayref();
+ foreach my $r (@{$tbl_ref}) {
+ if($rtn eq '') {
+ $rtn = $r->[0];
+ }
+ else {
+ $rtn = $rtn . ", " . $r->[0];
+ }
+ }
+ return $rtn;
+}
+=head1 NAME
+
+ -
+
+
+=head1 VERSION
+
+The initial template usually just has:
+
+This documentation refers to version 0.0.1.
+
+
+=head1 USAGE
+
+ # Brief working invocation example(s) here showing the most common usage(s)
+
+ # This section will be as far as many users ever read,
+ # so make it as educational and exemplary as possible.
+=head1 REQUIRED ARGUMENTS
+
+A complete list of every argument that must appear on the command line.
+when the application is invoked, explaining what each of them does, any
+restrictions on where each one may appear (i.e., flags that must appear
+before or after filenames), and how the various arguments and options
+may interact (e.g., mutual exclusions, required combinations, etc.)
+
+If all of the application's arguments are optional, this section
+may be omitted entirely.
+
+=head1 OPTIONS
+
+A complete list of every available option with which the application
+can be invoked, explaining what each does, and listing any restrictions,
+or interactions.
+
+If the application has no options, this section may be omitted entirely.
+
+
+=head1 DESCRIPTION
+
+A full description of the application and its features.
+May include numerous subsections (i.e., =head2, =head3, etc.).
+
+
+=head1 DIAGNOSTICS
+
+A list of every error and warning message that the application can generate
+(even the ones that will "never happen"), with a full explanation of each
+problem, one or more likely causes, and any suggested remedies. If the
+application generates exit status codes (e.g., under Unix), then list the exit
+status associated with each error.
+
+=head1 CONFIGURATION AND ENVIRONMENT
+
+A full explanation of any configuration system(s) used by the application,
+including the names and locations of any configuration files, and the
+meaning of any environment variables or properties that can be set. These
+descriptions must also include details of any configuration language used.
+
+
+=head1 DEPENDENCIES
+
+A list of all the other modules that this module relies upon, including any
+restrictions on versions, and an indication of whether these required modules are
+part of the standard Perl distribution, part of the module's distribution,
+or must be installed separately.
+
+
+=head1 INCOMPATIBILITIES
+
+A list of any modules that this module cannot be used in conjunction with.
+This may be due to name conflicts in the interface, or competition for
+system or program resources, or due to internal limitations of Perl
+(for example, many modules that use source code filters are mutually
+incompatible).
+
+
+=head1 BUGS AND LIMITATIONS
+
+A list of known problems with the module, together with some indication of
+whether they are likely to be fixed in an upcoming release.
+
+Also a list of restrictions on the features the module does provide:
+data types that cannot be handled, performance issues and the circumstances
+in which they may arise, practical limitations on the size of data sets,
+special cases that are not (yet) handled, etc.
+
+The initial template usually just has:
+
+There are no known bugs in this module.
+Please report problems to ()
+Patches are welcome.
+
+=head1 AUTHOR
+
+ ()
+
+
+
+=head1 LICENCE AND COPYRIGHT
+
+Copyright (c) (). All rights reserved.
+
+followed by whatever licence you wish to release it under.
+For Perl code that is often just:
+
+This module is free software; you can redistribute it and/or
+modify it under the same terms as Perl itself. See L.
+
+This program is distributed in the hope that it will be useful,
+but WITHOUT ANY WARRANTY; without even the implied warranty of
+MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+
diff -r 000000000000 -r acc8d8bfeb9a ptrfinder
Binary file ptrfinder has changed
diff -r 000000000000 -r acc8d8bfeb9a removeDupSV.pl
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/removeDupSV.pl Wed Feb 08 16:59:24 2012 -0500
@@ -0,0 +1,188 @@
+#!/nfs_exports/apps/64-bit/gnu-apps/perl5.8.9/bin/perl -w
+
+use strict;
+use Carp;
+use Getopt::Long;
+use English;
+use Pod::Usage;
+use Data::Dumper;
+use DBI;
+use DBD::mysql;
+
+my ( $help, $man, $version, $usage );
+my ($in, $out);
+my $min_inv = 10000;
+my $min_total_sclip = 6;
+my $add_inv_back;
+my $width = 20;
+my $optionOK = GetOptions(
+ 'i|in|input=s' => \$in,
+ 'o|out|output=s' => \$out,
+ 'min_inv_size=i' => \$min_inv,
+ 'min_total_sclip=i' => \$min_total_sclip,
+ 'add_inv_back!' => \$add_inv_back,
+ 'width=i' => \$width,
+ 'h|help|?' => \$help,
+ 'man' => \$man,
+ 'usage' => \$usage,
+ 'v|version' => \$version,
+);
+
+pod2usage(2) if($man);
+pod2usage( -verbose => 99, -sections => "USAGE|REQUIRED ARGUMENTS|OPTIONS" )
+ if($help or $usage);
+pod2usage( -verbose => 99, -sections => "VERSION") if($version);
+
+croak "Missing input file" if(!$in);
+open STDOUT, ">$out" if($out);
+open my $IN, "<$in" or croak "can't open $in: $OS_ERROR";
+
+my @SV;
+my @discard_SV;
+while(my $line = <$IN>) {
+ chomp $line;
+ my ($chrA, $posA, $ortA, $countA, $chrB, $posB, $ortB, $countB, $type) = split /\t/, $line;
+ next if( search_sv(\@SV, $chrA, $posA, $chrB, $posB, $type));
+ if(($type eq 'ITX' || $type eq 'INV') && (abs($posB - $posA) < $min_inv && ($countA + $countB) < $min_total_sclip)) {
+ print STDERR "$line is removed due to possible circulization\n";
+ push @discard_SV, $line;
+ next;
+ }
+ push @SV, [$chrA, $posA, $chrB, $posB, $type, $line];
+}
+
+close $IN;
+
+foreach my $sv (@SV) {
+ print STDOUT $sv->[5];
+# if($sv->[4] eq 'DEL') {
+# my $tmp = get_gene_info($sv->[0], $sv->[1], $sv->[3]);
+# print STDOUT "\t", $tmp if($tmp);
+# }
+# else {
+# print STDOUT "\t", get_gene_info($sv->[0], $sv->[1], $sv->[1]);
+# print STDOUT "\t", get_gene_info($sv->[2], $sv->[3], $sv->[3]);
+# }
+ print STDOUT "\n";
+}
+
+if($add_inv_back) {
+ foreach my $sv (@discard_SV) {
+ my ($chrA, $posA, $ortA, $countA, $chrB, $posB, $ortB, $countB, $type) = split /\t/, $sv;
+ next if(search_nearby_sv(\@SV, $chrA, $posA, $posB));
+ print STDOUT $sv;
+ print "\n";
+ }
+}
+close STDOUT if($out);
+exit(0);
+
+sub search_nearby_sv {
+ my($r_SV, $chrA, $posA, $posB) = @_;
+ foreach my $sv (@{$r_SV}) {
+ return 1 if($sv->[0] eq $chrA && abs($sv->[1] - $posA) < 10000);
+ return 1 if($sv->[2] eq $chrA && abs($sv->[3] - $posA) < 10000);
+ return 1 if($sv->[0] eq $chrA && abs($sv->[1] - $posB) < 10000);
+ return 1 if($sv->[2] eq $chrA && abs($sv->[3] - $posB) < 10000);
+ }
+ return;
+}
+
+sub search_sv {
+ my ($r_SV, $chrA, $posA, $chrB, $posB, $type) = @_;
+
+ for(my $i = 0; $i < scalar @{$r_SV}; $i++) {
+ my $sv = $r_SV->[$i];
+ if($type eq "INV" && $r_SV->[$i][4] eq "ITX" ) {
+ if(( $sv->[0] eq $chrA && abs($sv->[1] - $posA) < 20 &&
+ $sv->[2] eq $chrB && abs($sv->[3] - $posB) < 20 )||
+ ( $sv->[0] eq $chrB && abs($sv->[1] - $posB) < 20 &&
+ $sv->[2] eq $chrA && abs($sv->[3] - $posA) < 20)) {
+ delete $r_SV->[$i];
+ return;
+ }
+ }
+ else {
+ return 1 if( $sv->[0] eq $chrA && abs($sv->[1] - $posA) < $width &&
+ $sv->[2] eq $chrB && abs($sv->[3] - $posB) < $width && $sv->[4] eq $type);
+ return 1 if( $sv->[0] eq $chrB && abs($sv->[1] - $posB) < $width &&
+ $sv->[2] eq $chrA && abs($sv->[3] - $posA) < $width && $sv->[4] eq $type);
+ return 1 if( $sv->[0] eq $chrA && abs($sv->[1] - $posA) < $width &&
+ $sv->[2] eq $chrB && abs($sv->[3] - $posB) < $width && $sv->[4] eq $type && $type eq 'CTX');
+ return 1 if( $sv->[0] eq $chrB && abs($sv->[1]- $posB) < $width &&
+ $sv->[2] eq $chrA && abs($sv->[3] - $posA) < $width && $sv->[4] eq $type && $type eq 'CTX');
+ }
+ }
+ return;
+}
+
+=head1 NAME
+
+removeDupSV.pl - a program to remove duplicated SVs for CREST.pl
+
+
+=head1 VERSION
+
+This documentation refers to removeDupSV.pl version 0.0.1.
+
+
+=head1 USAGE
+ The program need an input file generated by CREST.pl
+ removeDupSV.pl -i tumor.predSV.txt -o tumor.predSV.txt.nodup
+
+=head1 REQUIRED ARGUMENTS
+
+ To run the program, several parameter must specified.
+ -i The input file generated by CREST.pl
+ -o The output file
+
+=head1 OPTIONS
+
+ The options that can be used for the program.
+ --min_inv_size For ITX/INV, the minimum event size, default 10000
+ --min_total_sclip For ITX/INV with size smaller than min_inv_size, the min
+ number of soft-clipped reads needed to keep it, default 6
+ --(no)add_inv_back Add ITX/INV back if there is another event nearby,
+ default OFF.
+ --width The max breakpoint distance such that it's considered as
+ same breakpoint, default 20
+ -h, --help, -? Help information
+ --man Man page.
+ --usage Usage information.
+ --version Software version.
+
+
+=head1 DESCRIPTION
+
+This is a simple program to remove duplicated SVs, also it remove small INV/ITX
+events since most of those are due to library preparation.
+
+
+=head1 BUGS AND LIMITATIONS
+
+There are no known bugs in this module, but the method is limitted to bam file
+that has soft-clipping cigar string generated.Please report problems to
+Jianmin Wang (Jianmin.Wang@stjude.org)
+Patches are welcome.
+
+=head1 AUTHOR
+
+Jianmin Wang (Jianmin.Wang@stjude.org)
+
+
+=head1 LICENCE AND COPYRIGHT
+
+Copyright (c) 2010 by St. Jude Children's Research Hospital.
+
+This program is free software: you can redistribute it and/or modify
+it under the terms of the GNU General Public License as published by
+the Free Software Foundation, either version 2 of the License, or
+(at your option) any later version.
+
+This program is distributed in the hope that it will be useful,
+but WITHOUT ANY WARRANTY; without even the implied warranty of
+MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+GNU General Public License for more details.
+
+You should have received a copy of the GNU General Public License
+along with this program. If not, see .
diff -r 000000000000 -r acc8d8bfeb9a tumor.bam
Binary file tumor.bam has changed
diff -r 000000000000 -r acc8d8bfeb9a tumor.bam.bai
Binary file tumor.bam.bai has changed
diff -r 000000000000 -r acc8d8bfeb9a tumor.bam.cover
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/tumor.bam.cover Wed Feb 08 16:59:24 2012 -0500
@@ -0,0 +1,11 @@
+4 125892327 + 1 30
+4 125892458 + 1 27
+4 125893225 + 1 28
+4 125893227 + 5 29
+4 125893365 - 1 26
+4 125893979 - 1 16
+10 66301086 - 2 33
+10 66301858 + 4 15
+10 66301865 - 9 22
+10 66301871 - 1 24
+10 66302136 + 1 51
diff -r 000000000000 -r acc8d8bfeb9a tumor.bam.sclip.txt
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/tumor.bam.sclip.txt Wed Feb 08 16:59:24 2012 -0500
@@ -0,0 +1,27 @@
+4 125892327 + HWUSI-EAS1591_6113C:3:17:12332:19420#0 CC CC
+4 125892458 + HWUSI-EAS1591_6113C:4:91:6281:9961#0 GACTAACCACCACGGTACATGTTTTCCTATGTAAAAAACCTGCACATTCTACACATGTATCCCAGAACTTAAAGTAAAACAC B@C@?:CC>CCBCCCCACBCDCCCCCC;<:<9CCCCC@CCCCCBCCCCCCCCCCCCCCCACCCCCCCCCCCCCCCCCCCCAC
+4 125893225 + HWI-EAS90_614M9:5:18:17924:10181#0 CCCTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCTCCCGAGTAGCTGGGATTACAGGTGCCCACCACCATGCCTGGCTAA #######@@7@:8@><16+6(B>AABCAA3AB@CC6CCCCCCCDCCCCCCBCCDCCCCCCCCCCCCDCCCCCCCCCCCCCC
+4 125893227 + HWI-EAS90_614M9:4:32:9763:4986#0 CTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCT CCCC@CCCCCCCCCCCCCCCCCCCCCCCCCCCCC
+4 125893227 + HWI-EAS90_614M9:5:1:12903:5535#0 CTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCTCCCGAGTAGCTGGGATTACAGGTGCCCACCACCATGCCTGGCTAATTTT B;@@>DABBB>DBAAAB@BABDDDDCCBCDACCBCCCCCCCBCCCCCCCCCCCCCCCCCCCCCCBCCCCCCCCCCCCCCCCCC
+4 125893227 + HWUSI-EAS1591_6113C:2:117:3095:2051#0 CTCCTGGGTTCAAGTGATTCTCCTGCCTCTACCTCCCGAGTAGCTGGGATTACAGGTGCCCACCACCATGCCTGGCTAAT 65977A::??BB>?BA8???=BDCCCCBCCCCCBCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCDCCCCCCCCCCCCCCC
+4 125893227 + HWUSI-EAS1591_6113C:4:39:14328:11130#0 CTCCTGGGTTCAAGTGATTCTCCTG CCCCCCCCCCCCCCCCCCCCCCCCC
+4 125893227 + HWUSI-EAS1591_6113C:4:9:4513:12290#0 CTCCTGG CCCCCCC
+4 125893365 - HWI-EAS292_102081428:1:66:17910:13190#0 TAACTACTT GFGFDGGGG
+4 125893979 - HWUSI-EAS1591_6113C:1:81:16173:11673#0 ATAAAATATTTATATTTTATATATTTGTATTCAAT CCCCCCCCCCCCCCCCCCCCCCCCCCCDCCCCCCC
+10 66301086 - HWI-EAS231_614TA:6:108:5679:13132#0 AAACATGGTTCAAATAAATTATTTATAAAGAAG CCCCCCCCCCCCCCCCCBBCCCCCCCCCCCCCC
+10 66301086 - HWI-EAS292_102081428:1:76:17776:13382#0 AAACATGGTTCAAATAAATTATTTATAAAGAAG GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGG
+10 66301858 + HWI-EAS231_614TA:1:37:3390:19595#0 ACAGCTACT CCCCCCCCC
+10 66301858 + HWI-EAS90_614M9:1:98:9476:6999#0 ACAGCTACTGAAAAGAGAATCATATTCTCTTTT CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
+10 66301858 + HWI-EAS90_614M9:5:12:2797:13435#0 ACAGCTACTGAAAAGAGAATCATATTCTCTTTTCATAATTTAGAATATGATATATATTTCAAAATATGATA CCCCCCDCCCCCCCCCCCCACCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
+10 66301858 + HWI-EAS90_614M9:7:110:15112:12998#0 ACAGCTACTGAAAAGAGAATCATATTCTCTTTTCATAATTTAGAATATGATATATATTTCAAAAT CCCCCCBCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
+10 66301865 - HWI-EAS231_614TA:4:36:3528:11249#0 TTTGTTCCTTCA CCCCCCCCCCCC
+10 66301865 - HWI-EAS231_614TA:5:52:9612:8954#0 CGTAGTTTGTTCCTTCA CCCCCCCCCCCCCCCCC
+10 66301865 - HWI-EAS90_614M9:1:118:15320:11068#0 AGCCATGGACCTTGTGGTGGGTTCTTAACAATGGTGAGTCCGGAGTTCTTAACGATGGTGAGTCCGTAGTTTGTTCCTTCA CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCACCCDCCCCCCCCCCCCCCCC@CCBCCCDDCCDBCDBCCDBCBCDDD
+10 66301865 - HWI-EAS90_614M9:3:104:3717:11676#0 GAGTCCGTAGTTTGTTCCTTCA BBBBBB@9B@878;69?B9AB>
+10 66301865 - HWI-EAS90_614M9:4:115:17151:16115#0 AGTCCGGAGTTCTCAACGATGGTGAGTCCGTAGTTTGTTCCTTCA CCCCCCCCCCCCCACCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
+10 66301865 - HWI-EAS90_614M9:5:27:7334:16234#0 CTTGTGGTGGGTTCTTAACAATGGTGAGTCCGGAGTTCTTAACGATGGTGAGTCCGTAGTTTGTTCCTTCA CCCCCCCDCCCCCCCCCCCCCBCCBCCCDCCCCCCBCCCCCCCCC@CC>CAC@CCBABCADDDBDBBBBBB
+10 66301865 - HWI-EAS90_614M9:7:21:19162:4940#0 TTCA CCCC
+10 66301865 - HWUSI-EAS1591_6113C:1:24:6705:7145#0 AGCCATGGACCTTGTGGTGGGTTCTTAACAATGGTGAGTCCGGAGTTCTTAACGATGGCGAGTCCGTAGTTTGTTCCTTCA CCCCCCCCCCCCCCCCCDCCCCCCCCCCCCCCCCDCCCDCCCCCCACCCCCCDCCDC98A@C