Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,947,481	1	.	A	29.7%	125.2 / 99.0	64	intergenic (+18/‑16)	metZ/metW	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base . (6/39);  new base A (0/19);  total (6/58)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-01

GGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGG  >  NC_000913/2947331‑2947613                                                                                                                                                           |                                                                                                                                     ggcgTAAAAAAACGGGTGCATGAAAACTCCTGGGTAGTATTGTCGCCA‑TCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑a                                                                                                                                     <  5:78560/148‑1 (MQ=11)  gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑a                                                                                                                                     <  3:915369/151‑1 (MQ=255)  gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑a                                                                                                                                     <  6:1080575/151‑1 (MQ=255)  gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑a                                                                                                                                     <  1:715204/151‑1 (MQ=255)  gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTT‑AAGAAG‑T‑a                                                                                                                                     >  2:1024140/1‑150 (MQ=255)  gggTAAAAAAAAGGGTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑a                                                                                                                                     <  4:1365016/151‑1 (MQ=255)  ggctAAAAACAAGCGTCCATGAAAACTACATCGTAGTATATTCCCCA‑TCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aa                                                                                                                                    <  4:875356/148‑1 (MQ=1)   ggTAAAACAAACGGTGCAGGAAAACCCGTGCGGAGAATAGTCCCCC‑TCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aaa                                                                                                                                   <  5:1151998/151‑1 (MQ=11)   ggTAAAAAAACGGCTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aaa                                                                                                                                   <  7:1538602/151‑1 (MQ=255)   gcaaaactccacggggCAGGAAATCGCGAGTGCTGTAGACTGGGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAGCCGAAGGTC‑GGCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aaa                                                                                                                                   <  3:526338/136‑1 (MQ=11)    gTAACAACAAGGCTGAATGAAAACAAGATCGGAGTATACTGCGAA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aaaa                                                                                                                                  <  3:1456643/151‑2 (MQ=14)    gTAAACAAAAGGCTGCATCAAATCACGAGCGGAGTATACTACGAA‑TCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aaaa                                                                                                                                  <  5:82883/151‑2 (MQ=14)    gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑aaa                                                                                                                                   <  8:513332/151‑1 (MQ=255)    gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑aaa                                                                                                                                   <  8:1664086/151‑1 (MQ=255)    gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑aaa                                                                                                                                   <  5:1166939/151‑1 (MQ=255)     tAAAAATACGGTGGCATGAAAAAGCGAGCGGAGTATAGTGCGCC‑TCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aaaa                                                                                                                                  <  3:598397/150‑2 (MQ=17)     tAAAAAACAGGGTGAATGAAAACCATATGGGAGTATACTCCGCA‑TCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aaaa                                                                                                                                  <  4:1595515/150‑2 (MQ=11)     tAAAAAAACCGTTGCATGAAAACTCTAGCGGAGTATTGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aaaa                                                                                                                                  <  1:466772/150‑2 (MQ=17)      aaaataaaGGGTGCATGAAAAACCGTTCGGAGTATACTGCGCC‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aaaag                                                                                                                                 <  4:764116/150‑3 (MQ=14)      aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAg                                                                                                                                  <  4:542221/150‑1 (MQ=255)      aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGc                                                                                                                                 <  7:1030850/151‑1 (MQ=255)      aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGc                                                                                                                                 <  4:393722/151‑1 (MQ=255)      aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGG‑CCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGc                                                                                                                                 <  8:1418390/150‑1 (MQ=255)       aaatcaaGGTGGCATGAAAAATCCAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aaaagca                                                                                                                               <  6:1018834/146‑5 (MQ=14)       aaaaaacGGCTGAATGAAAACTCTAGCGGAGTATATTCGGCT‑TCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑T‑aaaagca                                                                                                                               <  8:1570572/151‑5 (MQ=11)       aaaaaaaGGTTGCATGAAACCGCGAGCGGAGTATAGTGCGCA‑TCCACGGCCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCGCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑ACAGCa                                                                                                                                <  8:960535/151‑1 (MQ=255)       aaaaaaaGGGTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCa                                                                                                                                <  4:50046/151‑1 (MQ=255)         atacaccttcaaTGAAAACGAGATGGCAGTATACTGCGCC‑TCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGt                                                                                                                              <  3:761932/140‑1 (MQ=11)         aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGt                                                                                                                              <  6:60150/151‑1 (MQ=255)          ataaGGTGGCATGAATACTCTTTGGGAGTATTTTTCGCT‑TCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTa                                                                                                                             <  8:1635761/149‑1 (MQ=11)          acaagcctgaaTGAAAACTACAGCGCAGTATATTGCGAA‑TACCCTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTa                                                                                                                             <  6:87837/141‑1 (MQ=11)           aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTAc                                                                                                                            <  6:169819/151‑1 (MQ=255)           aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTAc                                                                                                                            <  5:9444/151‑1 (MQ=255)           aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTAc                                                                                                                            <  3:1617109/151‑1 (MQ=255)           aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTAc                                                                                                                            <  8:1480358/151‑1 (MQ=255)           aaaGGTTGAATGAAAACTCGATCGGAGTATAGTCCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTAc                                                                                                                            <  7:625152/151‑1 (MQ=14)            aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACg                                                                                                                           <  2:942551/151‑1 (MQ=255)            aaGGGTGCATGAATAAGCGAGCGGAGTATAGTGCGCA‑TCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACg                                                                                                                           <  8:1039296/151‑1 (MQ=14)             cggTTGCATGAAAACGCTTGTGTTGTATATTCGGAT‑TCCCGGGACGCGGGGTGGAGCAGCCTGGTAGGTCGGCGGGATCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATT‑AAATTTGATGAAG‑TAAAAGCAGTACgg                                                                                                                          <  1:580622/150‑1 (MQ=11)             aGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATT‑AAATTTTATGAAGTT‑AAAGCAGTACgg                                                                                                                          >  3:1529161/1‑151 (MQ=255)               gTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTg                                                                                                                        <  3:387958/151‑1 (MQ=255)                gtGCATGAAAACTACAGCGTAGTATAGAG‑GCATTACCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGTGa                                                                                                                       <  6:60798/150‑1 (MQ=14)                gtGCATCAAAAAGAGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGTGa                                                                                                                       <  4:1658782/150‑1 (MQ=14)                 tGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGa                                                                                                                       <  1:1359167/150‑1 (MQ=255)                 tGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGa                                                                                                                       <  4:1422114/150‑1 (MQ=255)                 ggCATGAAAACGAGTGGGTAGTATTTTGCGAT‑TCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GGCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGTGAc                                                                                                                      <  2:1701989/150‑1 (MQ=11)                 gccATGAAAACTCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGTGAc                                                                                                                      <  3:1025148/149‑1 (MQ=14)                 gcaatgaattctctaGCGGTGTATACTGCGCA‑TCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGTGa                                                                                                                       <  3:525295/136‑1 (MQ=11)                  gaaTGAAAACGCAAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGTGAcg                                                                                                                     <  7:110221/149‑1 (MQ=14)                  gCATGAAAACTCTAGCGGAGTATAGTTCGCT‑TCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGTGAcg                                                                                                                     <  2:1527216/151‑1 (MQ=11)                  gCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTG‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGAc                                                                                                                      <  6:495684/150‑1 (MQ=255)                  gCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGAcg                                                                                                                     <  4:1222571/151‑1 (MQ=255)                  gCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGAcg                                                                                                                     <  7:870164/151‑1 (MQ=255)                  ccagcaaaaAGCTAGAGGAGTATAGTGCGCA‑TCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGTGAcg                                                                                                                     <  4:1263301/146‑1 (MQ=11)                  caaggAAAACTCCAGCGTAGTATAGTCCGCA‑TCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGTGAc                                                                                                                      <  1:1411224/146‑1 (MQ=11)                   aatgaaatctagtgCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGTGAcgc                                                                                                                    <  5:640425/138‑1 (MQ=11)                    aTGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAGATCCGGACCCCGCACCCAATTAATAATTGATGAAT‑T‑AAAGCAGTACGCTGCcgcg                                                                                                                   >  7:1387942/1‑151 (MQ=18)                      gAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCgg                                                                                                                  <  4:247685/150‑1 (MQ=255)                       tcctatatatcGGAGTATAGTGCGCC‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGTGAcgc                                                                                                                    <  4:936632/137‑1 (MQ=11)                       aaaaCGCGAGCTGAGTATAGTGCGCA‑TCCACTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCgggg                                                                                                                <  4:1091750/151‑1 (MQ=255)                       aaaaCGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCgggg                                                                                                                <  2:138498/151‑1 (MQ=255)                       aaaaCCCGTGCGGAGTATTCTGGCCT‑TCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGTGACGCgggg                                                                                                                <  2:413983/151‑1 (MQ=11)                          aCGCGAGCGGATTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTgg                                                                                                             <  1:393146/151‑1 (MQ=255)                           cgggAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGGAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGa                                                                                                            <  2:659824/148‑1 (MQ=255)                           cgcgAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGa                                                                                                            <  3:842731/151‑1 (MQ=255)                           cgcgAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGa                                                                                                            <  4:495595/151‑1 (MQ=255)                             cgAGCGTAGTATAGTGAGCA‑TCGACGGACGAGGTGTGGAGCAGCCTGGTAGCTCGTCGGTCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCGAGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGagc                                                                                                          <  6:1324786/151‑1 (MQ=17)                               aGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGagcag                                                                                                        <  4:796614/151‑1 (MQ=255)                               aGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTAATAAAT‑A‑AAAGCAGTACGGTGAcgc                                                                                                                    >  6:935324/1‑139 (MQ=17)                                 gcgAGTATAGTGCGAT‑TGGTCTGACGCGGGGGGGAGTAGCCTGGTAGTTCGGTGGGATCATAACCCGAAGGGG‑GTAGGTTCTAATCCGGCCGCCGCAACGAATTAAAATTTGTTGTAG‑T‑AATGGAGTACGGTGACGCGGGGTGGAGCAGcc                                                                                                      <  8:1007289/149‑1 (MQ=9)                                 cGGACTATAGTGCTCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGATCATAACCCGAAGGTC‑GGCGGTTCAAATCCTGCCCCCGCAACCAATTAAAATTTGAAGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGcc                                                                                                      <  1:1276032/151‑1 (MQ=17)                                  ggAGTATAGTTATCA‑TCTCCGGACTCTGTGTGGAGCAGTCTTGTATCTCGTTTGTCTCTTAACCCGAATGTG‑GTCTGTTCTTATAATTCAGCCTCAACGAATTAAAATTTGATGAAG‑T‑CAAGGAGTACTGTGACGCGGGGTGGAGCAGCCt                                                                                                     <  6:171843/151‑1 (MQ=9)                                  ggAGTATAGTGCTCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCt                                                                                                     <  6:642141/151‑1 (MQ=255)                                  ggAGTATAGTGCGCA‑TCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑TAAAAGCAGTACGGt                                                                                                                         <  8:622343/131‑1 (MQ=11)                                         agggcccA‑TCCACGGACCCGGTCTGGATCAGCGTGGTCGCTCGTATGTCTCATAACCCTAATGGG‑GTCTGTTGAAATCAGTCACGCCCAAAGAATTAAAATTTGTTGAAT‑T‑TTTGCAGTACCGTGACGCGGGGTGGAGCAGCCTGGTAGc                                                                                               <  4:71670/144‑1 (MQ=9)                                              cgattcgcctgACGCGGTGTGGATCAGCCTGGTATCTCGGCGGCCTAATAACCCGAATGGC‑GGCGGTGCTAATCATTCCCCCGCAACCAATTACATTTTGATGAAG‑T‑AACGCACTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgtc                                                                                          <  4:321431/141‑1 (MQ=2)                                                   cctccgaatcttggTGGATCAGCCTGTTCTCTCTTCGGCATCATAAGCCCAAGGTCTGATGGTCCAAAAAC‑TCCCCCGCAACCAAATAAAATTTGCTGTAT‑T‑CACGGACTACCGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCt                                                                                     <  5:1304463/139‑1 (MQ=2)                                                                      cagcCTGGTAGCTCGTCGCTCTCATAACCCGAAGGTC‑GTCGGTGCAAATCCCGCCCCCGCAACCAAATAAAAATTGATGAAG‑T‑ATCGCACTACCGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGtcgtcg                                                                  <  3:1225276/151‑1 (MQ=11)                                                                                                    gAAGCTC‑GTCTGTTCAAATACGGCCCCCTCAACCAATTAAAAATGGATGAAG‑T‑AAAGCACTCCCCTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaa                                    <  7:1108262/151‑1 (MQ=14)                                                                                                          gcggcgcgtccaaTCCCCCCCCCGCAACCAATTATACTTTGATGTAG‑A‑AAACCCCTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtga                               <  3:951730/140‑1 (MQ=14)                                                                                                                                ggccagcAATTATAATTTGATGAAG‑T‑TTCGGCATACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGa         <  2:66793/145‑1 (MQ=17)                                                                                                                                   aaCCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATTTAGATGAAGGAAAAGccgaggggtgcgg       >  2:95197/1‑138 (MQ=14)                                                                                                                                   aaCCAATTAAAATTTCATTAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAttattgaac                                <  2:920835/125‑9 (MQ=11)                                                                                                                                       aaTTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGATGATGATCTAAAAGCAGTACGGTGACGCgggg  >  3:80754/1‑151 (MQ=17)                                                                                                                                       aaTTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAATTTTCAGGAAGTAAAGGCTGTACTGTGACGCgggg  >  3:89573/1‑151 (MQ=17)                                                                                                                                                           |                                                                                                                                     GGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCA‑TCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTC‑GTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑T‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGG  >  NC_000913/2947331‑2947613

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑