Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,855	0	A	C	44.6%	138.2 / 23.1	83	G247G (GGT→GGG)	gpmA	phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base A (41/5);  new base C (0/37);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 1.83e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTC  >  NC_000913/786708‑787002                                                                                                                                                      |                                                                                                                                                     tCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTAcc                                                                                                                                                    <  2:216335/150‑1 (MQ=255)  cAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCt                                                                                                                                                  <  2:384902/151‑1 (MQ=255)   agagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCt                                                                                                                                                  <  2:647682/150‑1 (MQ=255)    gagTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑gg                                                                                                                                               <  1:837570/151‑1 (MQ=255)        ggCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTg                                                                                                                                           >  2:1308604/1‑151 (MQ=255)             aTCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAAcc                                                                                                                                      <  2:976477/151‑1 (MQ=255)             aTCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAAcc                                                                                                                                      >  1:1076385/1‑151 (MQ=255)               ccGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAAccgc                                                                                                                                    >  2:18520/1‑151 (MQ=255)                         accacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTc                                                                                                                          >  1:1251939/1‑151 (MQ=255)                         accacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTc                                                                                                                          >  2:1170358/1‑151 (MQ=255)                          ccacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCg                                                                                                                         >  2:662857/1‑151 (MQ=255)                             caGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTg                                                                                                                      >  2:463582/1‑151 (MQ=255)                              aGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGc                                                                                                                     >  2:1056009/1‑151 (MQ=255)                              aGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGc                                                                                                                     >  2:565032/1‑151 (MQ=255)                               gcccTGATAAGTCGCGCGGCGTCGCATCGGTCAATGTGCTCCTTTGTT‑AACAAAAAAAAACCCGACTCATTTGCAGTCGGTTTTCTCTTTTTAACC‑AAATGGCTTTTACTTCGTTTTCCCTT‑GTTTTGCAACCGCCGCTGCTTTCGCTg                                                                                                                      <  1:1011010/147‑1 (MQ=255)                                   tGATAGGTCGCGCAGCGTCGCATCAGGCAATGTGCCCCATTGTT‑AGCAACAAAAAAGCCGACTCCCTTGCAGTCGGTTTTCTCATTTTAAAC‑AAATGACGTTTCCTTCGTTTTCCCCT‑GGTTTGCACCCGCCGCTGCTTTCGCTGCGAt                                                                                                                  <  1:828330/149‑1 (MQ=255)                                   tGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGAtct                                                                                                                >  2:231425/1‑151 (MQ=255)                                    gATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                               >  1:1085758/1‑151 (MQ=255)                                     aaaggTCGCGCGGCGTCGCATCAGGCAATGTGCTCCATTTTA‑ACCAAAAAAAAACCCGACTCCCTTCCGGTCGGTTTTCTCATTTAAAAC‑GAAGGACTTTTACTTCGCTTTCCCCG‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                               <  1:67364/146‑1 (MQ=255)                                     aaaaGCCGCGCAGCGTCGCATCAGGCAGTGTCCTCCATTGTT‑AGCAACAAAAAAGCCGACTCCCTTGCAGTCGCTTTTCTCATTTTAAAC‑GAATGGCTTTTCCTTCTCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCg                                                                                                              <  1:992907/149‑1 (MQ=255)                                       aaGTCGCGCAGCGTCGCATCGGGCAATGTGCTCAATTTTT‑AGCAAAAAAAAAGCCGCCTCCCTTGCAGTCGTTTTTCTCATTTTAACC‑GAATGCCGTTTACTTCGTTTTCCCCT‑GTTTTCCAACCGCCGCTGCTTTCGCt                                                                                                                       <  1:73514/140‑1 (MQ=255)                                       aaGTCGCGCAGCGTCGCATCGGGAAATGTCCTCCTTTGTT‑AGAAAAAAAAAAGCCGACTCTCTTGCAGTGGGTTTTCTCATTTTAAAC‑GAATGACTTTTCCTTCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCg                                                                                                              <  1:404324/149‑1 (MQ=255)                                       aaGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc                                                                                                            >  2:540587/1‑151 (MQ=255)                                       aaGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTA‑AGCAAAAAAAAACCCGCCTCACTTGCAGTCGGCTTTCTCTTTTAAAAC‑GAATGACGTTTACTTCGCTTTCCCCT‑GGTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCg                                                                                                              <  2:759827/149‑1 (MQ=255)                                        aGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc                                                                                                            >  2:431576/1‑150 (MQ=255)                                          tCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGc                                                                                                         >  1:425772/1‑151 (MQ=255)                                           cgcgCAGCGTCGCATCAGGAAATGTTCCCCTTTTTT‑AAAAAAAAAAAGCCCCACTCCCTTGCAGTCGGTTTTCTCATTTAAAAC‑AAATGCTTTTTTCTTCGTTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCa                                                                                                           <  1:793079/148‑1 (MQ=255)                                               cAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAc                                                                                                    >  1:150415/1‑151 (MQ=255)                                                 gCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                  >  2:1302674/1‑151 (MQ=255)                                                  cGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                  >  1:668398/1‑150 (MQ=255)                                                   gtcacttggggCAATGTGCT‑CTTTGTTAAGCAAAAAAAAGGCCGACCCCTGTGCAGTCGGCTTTCTCTTTTAAAAC‑GAATGCCTTTTCCTCGGTTTTCCCCT‑GTTTGCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAg                                                                                                <  1:737839/142‑1 (MQ=255)                                                   gTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAGAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                  >  1:459530/1‑149 (MQ=255)                                                   gTCGCATCAGGAAGTGTGCTCCATTGTA‑AACAAAAAAAAAGCCGACTCACTTGCAGTGGGCTTTCTCATTTTAAAC‑GAATGACTTTTACTTCGCTTTCCCCG‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAg                                                                                                <  1:1205285/151‑1 (MQ=255)                                                          cAGGTAATGTCCTCCATTTTA‑ACCAAAAAAAAACCCCACTCCCTCGCATTCCTTTTTCTCTTTTAAACC‑AAATGACGTTTCCTTCCCTTCCCCCT‑TTTTTGAACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                           <  1:29215/149‑1 (MQ=255)                                                          cAGGCATTGTCCCCCATTTTT‑GCCAACAAAAAACCCGACTCATTTGCGTTGGGCTTTCTC‑TTTTAAACGGAGGGACGTTGCCTTCGCTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGAGATCTCGTCAGCATTACCCAGata                                                                                             <  2:1213233/147‑1 (MQ=255)                                                          cAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                           >  2:1012715/1‑149 (MQ=255)                                                            gacAATTTGCTCCTTTGTT‑AAAAAAAAAAAAGCCGACTCATTTGCATTCGCCTTCCTAATTTTAAAC‑AAATGCCTTTTATTTCGCTTCCCCCT‑GGTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGat                                                                                              <  2:95339/142‑1 (MQ=255)                                                                 tgtgCTCCTTTGTT‑AAAAAAAAAAAGGCCGCCTCCCTTGCGGTCGGCTTTCTCTTTTAAACC‑GATTGCCTTTTCCCTCGCTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTc                                                                                   <  1:132955/150‑1 (MQ=255)                                                                    gCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCg                                                                               >  2:963057/1‑151 (MQ=255)                                                                     cTCCATTTTT‑AGCAAAAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCa                                                                                  <  2:69808/147‑1 (MQ=255)                                                                      tCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGt                                                                             >  1:726862/1‑151 (MQ=255)                                                                       cctttGTT‑ACAAACAAAAAAGCCGACTCACTTGCAGTCGGTTTTCTCATTTAAAAC‑GAAGGACGTTTACTCCGCTTTCCCCT‑GGTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCa                                                                                  <  2:548407/142‑1 (MQ=255)                                                                       ccATTGTA‑AGCAACAAAAAAGCCGCCTCACTTGCATTCGGCTTTCTCTTTTAAAAC‑AAATGACTTTTACTTCGTTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGc                                                                                <  2:178689/147‑1 (MQ=255)                                                                          ttttt‑AGAAAAAAAAAAGCCGACTCCCTTGCAGTCGGTTTTCTAATTTAAACC‑GAATGACTTTTACTTCCCTTTCCCTT‑GTTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGc                                                                                <  1:18520/141‑1 (MQ=255)                                                                           ttttaaaaaaaaaaaaaCCCGACTCACTTGCAGTCGTTTTTCTCTTTTAAAAC‑GAATGACGTTTCCTTCTCTTTCCCCT‑GTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTg                                                                          <  1:723945/139‑1 (MQ=255)                                                                           tGTT‑AGAAAAAAAAAACCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGCCTTTTACTTCGTTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                           <  1:39459/148‑1 (MQ=255)                                                                            tttaaccaataaAAAAGCCGACTCCCTTGCATTGGTCTTACTCATTTAAAAG‑GAAAGATTTTTCCTTCGTTTTCCCTT‑GGTGTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGa                                                                         <  2:1190929/144‑1 (MQ=255)                                                                             tt‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAAt                                                                      >  2:77658/1‑151 (MQ=255)                                                                              taaacaacaaAAAACCCGCCCCATTTGCATTCGGTTTTCTATTTTAAAAC‑GAATATCTTTTACTCCG‑TTTCCCCTGGGTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaaa                                                                       <  1:464378/146‑1 (MQ=255)                                                                                  caacaaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtctc                                                                  >  2:806600/1‑151 (MQ=255)                                                                                  caacaaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtct                                                                   >  2:665573/1‑150 (MQ=255)                                                                                   aacaaAAAAGCCGACTCCCTTCCCGTCGCTTTTCTTTTTTAAAAC‑GAATGACTTTTACTTCTCTTTCCCCT‑TGTTTACAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtctc                                                                  <  1:348145/150‑1 (MQ=255)                                                                                   aacaaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcg                                                                 >  2:101363/1‑151 (MQ=255)                                                                                   aacaaAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcg                                                                 >  2:971278/1‑151 (MQ=255)                                                                                    aaaaaaacGCCGCCTCATTTGCGGTCGGCTTTCTCATTTAAAAC‑AAAGGAGGTTTACTTCGCTTTCCCTG‑GGTTAACACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                                <  1:150318/149‑1 (MQ=255)                                                                                      aaaaagGCCGCCTCTCTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgtc                                                               <  2:1210429/150‑1 (MQ=255)                                                                                          acgccgtcccccTTGCGGTCGTTTTTCTCTTTTTAAAC‑AAATGACTTTTCCTTCTCTTCCCCCG‑TTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAAc                                                           <  1:862877/139‑1 (MQ=255)                                                                                               gACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATac                                                     >  1:164005/1‑151 (MQ=255)                                                                                                   cACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                  >  1:181154/1‑150 (MQ=255)                                                                                                     cTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCa                                                 <  2:317775/149‑1 (MQ=255)                                                                                                          aGTCGGTTTTCTCATTTAAAAC‑GAATGGCGTTTACTTCGCTTTCCCCG‑TGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCa                                           <  2:612318/150‑1 (MQ=255)                                                                                                              ggCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGcc                                       >  2:622836/1‑150 (MQ=255)                                                                                                                cTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                                    >  1:1154999/1‑151 (MQ=255)                                                                                                                cTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                                    >  1:260182/1‑151 (MQ=255)                                                                                                                  ttCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                                  >  2:633513/1‑151 (MQ=255)                                                                                                                  ttCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                                  >  1:6770/1‑151 (MQ=255)                                                                                                                    ctcATTTAAAAC‑GAATGACTTTTACTTCGCTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACCCATACACCAGCGGCACGCCAGTCg                                  <  2:884970/149‑1 (MQ=255)                                                                                                                      cATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGat                              >  2:137515/1‑151 (MQ=255)                                                                                                                        tttaAAAC‑AAATGACTTTTACTTCGTTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGata                             <  1:985335/150‑1 (MQ=255)                                                                                                                             aaaggaagagcGTTTTCTCCGCTTCCCCTT‑GTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAg                        <  1:223077/141‑1 (MQ=255)                                                                                                                                cgaagtacttttACTTCGTTTTCCCCT‑GGTTTGCATCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCt                      <  2:701779/141‑1 (MQ=255)                                                                                                                                 gAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTc                     >  1:244041/1‑150 (MQ=255)                                                                                                                                  aaTGACTTTTCCTTCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCaa                   <  2:618274/151‑1 (MQ=255)                                                                                                                                    tGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCa                    >  2:158160/1‑148 (MQ=255)                                                                                                                                    tGACG‑TTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGa                 >  2:798950/1‑150 (MQ=255)                                                                                                                                      aCGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAt               >  1:1045/1‑151 (MQ=255)                                                                                                                                        gTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtct             >  2:1131745/1‑151 (MQ=255)                                                                                                                                         tttCCTTCGCTTTCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGa                 <  2:138471/146‑1 (MQ=255)                                                                                                                                          ttccTTCGCTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtctc            <  1:254861/147‑1 (MQ=255)                                                                                                                                            aCTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttc         >  1:464997/1‑151 (MQ=255)                                                                                                                                              ttCGCTTTCCCCT‑GTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtctt          <  2:1076445/148‑1 (MQ=255)                                                                                                                                                ccttttCCCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCttcttc      <  2:1192820/147‑1 (MQ=255)                                                                                                                                                   tttACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTc  >  1:1248115/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTT‑AGCAACAAAAAAGCCGACTCACTTGCAGTCGGCTTTCTCATTTTAAAC‑GAATGACGTTTACTTCGCTTTACCCT‑GGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTC  >  NC_000913/786708‑787002

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