Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,621,118	0	A	C	46.0%	66.8 / 13.9	50	V7G (GTG→GGG)	eamA	cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base A (21/6);  new base C (0/23);  total (21/29)
Fisher's exact test for biased strand distribution p-value = 5.11e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCG  >  NC_000913/1620968‑1621262                                                                                                                                                        |                                                                                                                                                 ggTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAaca                                                                                                                                                  <  2:409464/151‑1 (MQ=255)    aCTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAccc                                                                                                                                               >  1:676920/1‑151 (MQ=255)       ttCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCAt                                                                                                                                            >  2:152898/1‑151 (MQ=255)          ggTCGTGCGATAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACCAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCAtttt                                                                                                                                         >  1:771771/1‑151 (MQ=255)             cGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTc                                                                                                                                      <  1:315694/151‑1 (MQ=255)                gCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcgc                                                                                                                                    <  2:776197/150‑1 (MQ=255)                    cAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACAt                                                                                                                               >  1:1018323/1‑151 (MQ=255)                       aaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCt                                                                                                                            >  1:869072/1‑151 (MQ=255)                       aaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCt                                                                                                                            >  2:413882/1‑151 (MQ=255)                              aGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCct                                                                                                                     >  2:309319/1‑151 (MQ=255)                                   gAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGa                                                                                                                 >  1:742894/1‑150 (MQ=255)                                      aaGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAAttt                                                                                                             >  1:817989/1‑151 (MQ=255)                                       aGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTc                                                                                                            >  1:148146/1‑151 (MQ=255)                                             cAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCttt                                                                                                      >  2:903262/1‑151 (MQ=255)                                                 aTACAGCGCAAACCGGCCAGCATCAGAGGTGGCATGTTATGAACACCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTCGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCCTCTTTTGt                                                                                                   >  1:772715/1‑150 (MQ=255)                                                       cgcAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCAAAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCaa                                                                                             <  2:613697/150‑1 (MQ=255)                                                         cAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATc                                                                                           >  1:455987/1‑150 (MQ=255)                                                          aacccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCAAAAAATTTACCCCCCATACGCCCCCTCCCAGTAGCCCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCaa                                                                                             <  1:55713/144‑1 (MQ=255)                                                             ccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGC‑CCAAAAAATTTAGCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                      <  1:920132/151‑1 (MQ=255)                                                             ccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACt                                                                                       <  2:1022377/150‑1 (MQ=255)                                                               ggCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGC‑CCACAAAATTTACCCCCCATCCGCCCACTCCCAGTAGCGCCAACACCCCATTTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAg                                                                                    <  1:113400/151‑1 (MQ=255)                                                               ggCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                      >  1:893317/1‑149 (MQ=255)                                                               ggCCAGCATCAGCGGGGGCATTTTATGACCCCCCTCTTTGATGC‑CCACAAAATTTACCCCCCATACCACCCCTCCCAGTAGCGCCACCCCCCCATTTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTa                                                                                     <  1:1031057/150‑1 (MQ=255)                                                                 ccagattCAGCGGTGCCATTTTAT‑AACCCCCCTTTTTGTGACCCACAAATTTTACCCCCCATCCGCCCACTCCCATTACCCCCACCCCCCTATTTTTCCGCAACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCt                                                                                  <  1:676876/145‑1 (MQ=255)                                                                          gCGGTGGCATTTTATAACGCCCCACTTGACCCC‑CCACAAAATTTAGCCCCATTCCGACCACTACCAGTAGCCCCAACCCCCCCTTTTTCCGCGACATTCTCCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAttt                                                                            <  1:413124/148‑1 (MQ=255)                                                                          gCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCACACGCCCCCTCCCAGTAGCGCCAACCCCCCACCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACt                                                                         <  2:676941/151‑1 (MQ=255)                                                                            ggTGGCATGTTATGAACCCCCACTTTGATGA‑CCAAAAAATTTAGCCCCCATACGCCCACTACCAGTACCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtt                                                                        <  2:983360/150‑1 (MQ=255)                                                                             gTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtttt                                                                      >  1:624688/1‑151 (MQ=255)                                                                               ggCATTTTATGAACCCCCACTTTGATGA‑CCACAAAATTTACCCCCCATACGACCACTCCCAGTAGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACttt                                                                       <  1:171161/148‑1 (MQ=255)                                                                                   tGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTCCCAGAAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCtt                                                                   <  2:384786/148‑1 (MQ=255)                                                                                   tGTTATGAACCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTa                                                                 <  2:512793/150‑1 (MQ=255)                                                                                          aacccccccTTTGATGC‑CCACAAAATTTACCCCCCATACGACCACTACCAGTAGCGCCACCCCCCCATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                          <  1:181251/147‑1 (MQ=255)                                                                                          aaGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                          >  2:526442/1‑150 (MQ=255)                                                                                             ccccACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTc                                                      >  1:616351/1‑151 (MQ=255)                                                                                             ccccACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTc                                                      >  2:207349/1‑151 (MQ=255)                                                                                                       atga‑CCACAAAATTTAGCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCCATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGcc                                            <  2:906148/151‑1 (MQ=255)                                                                                                       atga‑CCACAAAATTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGc                                             <  2:854581/150‑1 (MQ=255)                                                                                                       atga‑CCACAAAATTTACCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCg                                              <  2:879629/149‑1 (MQ=255)                                                                                                                aaaaTTTGCCCCCCATACGCCCACTACCAGTATCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTac                                      <  1:60411/149‑1 (MQ=255)                                                                                                                  aaTTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgacg                                  >  2:628394/1‑151 (MQ=255)                                                                                                                     ttACCCCCCATCCGCCCACTCCCAGTAGCGCCACCCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGc                                <  1:781162/150‑1 (MQ=255)                                                                                                                            ccATACGCCCTCACCCAGTACCCCCACCCCCCCTCCTTTCCCCAAATTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgc                        <  1:332950/151‑1 (MQ=255)                                                                                                                                aCGCCCCCTCCCAGTAGCCCCAACCCCCCATCTTTTCGCGAATTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCaa                    <  1:208068/151‑1 (MQ=255)                                                                                                                                  gACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAAt                  >  1:996670/1‑151 (MQ=255)                                                                                                                                    ccccTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAAt                  <  1:833280/146‑1 (MQ=255)                                                                                                                                     cACTCCCAGTAGCGCAAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAg               <  2:50874/151‑1 (MQ=255)                                                                                                                                      aCTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGt              >  2:107876/1‑151 (MQ=255)                                                                                                                                          ccAGTAGCGCAACCCCCCCATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTc             <  2:470446/148‑1 (MQ=255)                                                                                                                                          ccAGTACCGCCAACCCCCCATCTTTTCCCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCt          <  2:201290/151‑1 (MQ=255)                                                                                                                                                 cgcAAACACCCCTTCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTc   <  2:945461/151‑1 (MQ=255)                                                                                                                                                    cAACACCCCATCTTTCCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCg  <  2:972098/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                 GGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGA‑CCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCG  >  NC_000913/1620968‑1621262

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