Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,349,1450AC45.6% 90.2 / 15.0 57F110V (TTC→GTC) inaAacid‑inducible Kdo/WaaP family putative kinase
Reads supporting (aligned to +/- strand):  ref base A (28/3);  new base C (0/26);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 3.30e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGA  >  NC_000913/2348995‑2349291
                                                                                                                                                       |                                                                                                                                                  
aGATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGaa                                                                                                                                                    <  2:834098/151‑1 (MQ=255)
 gATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAAc                                                                                                                                                   >  1:849274/1‑151 (MQ=255)
   tatGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAAccc                                                                                                                                                 >  2:377399/1‑151 (MQ=255)
    atGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAAcccc                                                                                                                                                <  1:84394/151‑1 (MQ=255)
            cATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATc                                                                                                                                        >  2:350882/1‑151 (MQ=255)
             aTAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCt                                                                                                                                       >  2:998438/1‑151 (MQ=255)
                  aacaGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGt                                                                                                                                  >  1:915863/1‑151 (MQ=255)
                  aacaGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGt                                                                                                                                  >  1:885001/1‑151 (MQ=255)
                   acaGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAg                                                                                                                                   >  2:3104/1‑149 (MQ=255)
                       ccATGCTGACAATTAATGATATCAATTTTCTTAAACGCCGGCCCCACTGCTTCCAA‑CATGGTTTGCCTTATTTCGTAAAAAAAAGGCGATCCTGAATGCTGGGCATACCAGTAAGCAATGCTGATAACCCCCGCCATATCTTCAGTCa                                                                                                                                <  2:440817/148‑1 (MQ=255)
                            ctgaataTTAATGCTATGCTTTTTCTAAAACGCCAGCCCCACTGTTTTAAA‑CAGGGCTTCCCGTACTTCGTCAGAAAAGGGCAATACGGCATGCGGGGCATACCAGCAAGAAAGGCTAATAAACCCCGCCATATCTTCAGTCACCAGCa                                                                                                                          <  2:594526/143‑1 (MQ=255)
                            cTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAAc                                                                                                                        >  1:638479/1‑151 (MQ=255)
                                 gATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCg                                                                                                                   >  2:206883/1‑151 (MQ=255)
                                 gATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCg                                                                                                                   >  1:422134/1‑151 (MQ=255)
                                   ttattccTATGTATTTTCTAAACCGCCCGCCCCCGTGTTTACAC‑CTTGGCTTGCCGTCCTTCTTCAAAATAGGGCGTTACTGTTTGCTGGGCAAACCAGTCAGAAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCa                                                                                                                          <  2:640969/136‑1 (MQ=255)
                                          tATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCAcc                                                                                                          >  2:1004346/1‑151 (MQ=255)
                                           aTGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCAccc                                                                                                         >  1:287742/1‑151 (MQ=255)
                                            tattttttCTAAAACCCCAGCCCCCTTTTTTTAAA‑CATGGTTTGCCGTCCTTCGAAAAAAAAGGGCGATCCGGCGTGCGGGCCATCCCAGTCAGAAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCAtt                                                                                                              <  2:25364/141‑1 (MQ=255)
                                                ttttCTAAACCGCCCGCCCCACTGCTTTCAA‑CGGGGCTTGCCGTACTTCGACAAAAAAAGGCGATACTGCATGCTGGGCATCCCAGCCACCAAGGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCaa                                                                                                     <  2:566849/150‑1 (MQ=255)
                                                       aaccGCCAGCGCCCCTGTTTTCAA‑CATGGCTTGCCGTACTTCGTCAAAATAGGGCGATCCTGCATGCGGGGCATCCCAGTCAGCAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTc                                                                                                <  1:357265/145‑1 (MQ=255)
                                                        aaCGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCg                                                                                            >  1:793070/1‑151 (MQ=255)
                                                                 gcCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCa                                                                                    >  1:1014085/1‑150 (MQ=255)
                                                                   cACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaaa                                                                                  >  2:232702/1‑150 (MQ=255)
                                                                      ttttttCAC‑CAGGGTTTCCCTTTCTTCGAAAAAAAAAGGCGATCCTGCTTGCGGGACATACCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaaaaa                                                                                <  1:76538/146‑1 (MQ=255)
                                                                         tttCAA‑CATGGTTGCCCTTACTTCGTAAAAAAAGGGCGATCCTGCATGCGGGGCATCCCAGTCAGAAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGAt                                                                            <  1:978064/150‑1 (MQ=255)
                                                                         tttCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATc                                                                           >  1:333680/1‑151 (MQ=255)
                                                                          ttCAACCATGGCTTCCCGTCCTTCGCAAAAAAAGGGCGATCCTGCAGGCGGGCCACCCCAGCCAACAAGGCTGAGAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACg                                                                              <  2:56870/148‑1 (MQ=255)
                                                                          tcaaa‑CATGTCTGGCCGTACTTCGTAAAAATAAGGCGATACTGCATGCTGGACACACCAGTCATCAATGCTGAAAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGAt                                                                            <  2:652063/146‑1 (MQ=255)
                                                                                  tGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCAc                                                                   >  2:973663/1‑151 (MQ=255)
                                                                                       tcccGTCCTTCGTAAAAAAAAGGCAATTCTGCATGCTGGCCATCCCGGCCAAAACTGTTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATg                                                               <  2:385947/148‑1 (MQ=255)
                                                                                         ccGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGac                                                             >  1:931195/1‑150 (MQ=255)
                                                                                           gTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACAcc                                                          >  1:874071/1‑151 (MQ=255)
                                                                                                 cataaGAAAAGGGCGATTCTGCTTGCTGGGAATCCCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCg                                                         <  2:287675/142‑1 (MQ=255)
                                                                                                     aGAAAAAGGCGATCCTGCATGCGGGCCATCCCAGCCAGAAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGtt                                                  <  1:350969/149‑1 (MQ=255)
                                                                                                       aaaaagtcatcgcTGCAGGCGGGGCATCCAAGCCAAAAACGCGAAAGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTc                                                 <  2:745359/136‑1 (MQ=255)
                                                                                                       aaTAGGGAGATTCTGCAGGCGGGCAACCCCAGCCAGAAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAg                                              <  2:547972/151‑1 (MQ=255)
                                                                                                       aaTAAGGCGATCCTGCATGCGGGGCATCCCAGCCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAg                                              <  1:45288/151‑1 (MQ=255)
                                                                                                       aaTAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTcc                                                >  2:723105/1‑149 (MQ=255)
                                                                                                            ggCTATACTGCATGCTGGACATCCCAGTCAGCAATGCTAATACCCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCt                                          <  2:921492/150‑1 (MQ=255)
                                                                                                             gCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCttt                                        >  2:1048669/1‑151 (MQ=255)
                                                                                                             gCGATACTGCATGCTGGCCATCCCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCt                                          <  1:165689/149‑1 (MQ=255)
                                                                                                               gATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTa                                       >  2:1032953/1‑150 (MQ=255)
                                                                                                                    tgcgtgcgGGGCATCCCAGTCACCAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTa                                   <  1:232748/149‑1 (MQ=255)
                                                                                                                      cgttttGGGCATACCCGTCGAAAACGATTATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTa                                   <  2:895085/142‑1 (MQ=255)
                                                                                                                      cATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGc                               >  1:133833/1‑151 (MQ=255)
                                                                                                                       ttgctgcGCACCCCAGCCAACAATGCTAATGCCCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAAtt                                    <  2:859010/144‑1 (MQ=255)
                                                                                                                            ggccATACCAGTCAGCAATGCTAATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAAcc                           <  1:766806/146‑1 (MQ=255)
                                                                                                                                 tACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACgg                     >  1:608667/1‑150 (MQ=255)
                                                                                                                                   ccAGTAAGAAATGCGGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGa                    <  1:206926/149‑1 (MQ=255)
                                                                                                                                      gtcaaaaaTGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATc               <  1:647429/145‑1 (MQ=255)
                                                                                                                                         aGCAGTCCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAAt                <  2:767415/147‑1 (MQ=255)
                                                                                                                                          gCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTg           >  2:777756/1‑151 (MQ=255)
                                                                                                                                            aaggcagatgaCCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTgg          <  1:560208/144‑1 (MQ=255)
                                                                                                                                             cttcTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGt         <  2:666970/147‑1 (MQ=255)
                                                                                                                                              tGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCg       >  1:984516/1‑151 (MQ=255)
                                                                                                                                              tGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCg       >  1:346496/1‑151 (MQ=255)
                                                                                                                                              tGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCg       >  2:891308/1‑151 (MQ=255)
                                                                                                                                                   atgaCCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGa  <  1:998761/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
AGATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAA‑CATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGA  >  NC_000913/2348995‑2349291

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