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

AGATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGA  >  NC_000913/2348995‑2349291
                                                                                                                                                       |                                                                                                                                                  
aGATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGaa                                                                                                                                                    <  2:834098/151‑1 (MQ=255)
 gATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAAc                                                                                                                                                   >  1:849274/1‑151 (MQ=255)
   tatGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAAccc                                                                                                                                                 >  2:377399/1‑151 (MQ=255)
    atGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAAcccc                                                                                                                                                <  1:84394/151‑1 (MQ=255)
            cATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATc                                                                                                                                        >  2:350882/1‑151 (MQ=255)
             aTAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCt                                                                                                                                       >  2:998438/1‑151 (MQ=255)
                  aacaGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGt                                                                                                                                  >  1:915863/1‑151 (MQ=255)
                  aacaGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGt                                                                                                                                  >  1:885001/1‑151 (MQ=255)
                   acaGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAg                                                                                                                                   >  2:3104/1‑149 (MQ=255)
                       ccATGCTGACAATTAATGATATCAATTTTCTTAAACGCCGGCCCCACTGCTTCCAACATGGTTTGCCTTATTTCGTAAAAAAAAGGCGATCCTGAATGCTGGGCATACCAGTAAGCAATGCTGATAACCCCCGCCATATCTTCAGTCa                                                                                                                                <  2:440817/148‑1 (MQ=255)
                            ctgaataTTAATGCTATGCTTTTTCTAAAACGCCAGCCCCACTGTTTTAAACAGGGCTTCCCGTACTTCGTCAGAAAAGGGCAATACGGCATGCGGGGCATACCAGCAAGAAAGGCTAATAAACCCCGCCATATCTTCAGTCACCAGCa                                                                                                                          <  2:594526/143‑1 (MQ=255)
                            cTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAAc                                                                                                                        >  1:638479/1‑151 (MQ=255)
                                 gATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCg                                                                                                                   >  2:206883/1‑151 (MQ=255)
                                 gATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCg                                                                                                                   >  1:422134/1‑151 (MQ=255)
                                   ttattccTATGTATTTTCTAAACCGCCCGCCCCCGTGTTTACACCTTGGCTTGCCGTCCTTCTTCAAAATAGGGCGTTACTGTTTGCTGGGCAAACCAGTCAGAAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCa                                                                                                                          <  2:640969/136‑1 (MQ=255)
                                          tATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCAcc                                                                                                          >  2:1004346/1‑151 (MQ=255)
                                           aTGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCAccc                                                                                                         >  1:287742/1‑151 (MQ=255)
                                            tattttttCTAAAACCCCAGCCCCCTTTTTTTAAACATGGTTTGCCGTCCTTCGAAAAAAAAGGGCGATCCGGCGTGCGGGCCATCCCAGTCAGAAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCAtt                                                                                                              <  2:25364/141‑1 (MQ=255)
                                                ttttCTAAACCGCCCGCCCCACTGCTTTCAACGGGGCTTGCCGTACTTCGACAAAAAAAGGCGATACTGCATGCTGGGCATCCCAGCCACCAAGGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCaa                                                                                                     <  2:566849/150‑1 (MQ=255)
                                                       aaccGCCAGCGCCCCTGTTTTCAACATGGCTTGCCGTACTTCGTCAAAATAGGGCGATCCTGCATGCGGGGCATCCCAGTCAGCAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTc                                                                                                <  1:357265/145‑1 (MQ=255)
                                                        aaCGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCg                                                                                            >  1:793070/1‑151 (MQ=255)
                                                                 gcCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCa                                                                                    >  1:1014085/1‑150 (MQ=255)
                                                                   cACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaaa                                                                                  >  2:232702/1‑150 (MQ=255)
                                                                      ttttttCACCAGGGTTTCCCTTTCTTCGAAAAAAAAAGGCGATCCTGCTTGCGGGACATACCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaaaaa                                                                                <  1:76538/146‑1 (MQ=255)
                                                                         tttCAACATGGTTGCCCTTACTTCGTAAAAAAAGGGCGATCCTGCATGCGGGGCATCCCAGTCAGAAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGAt                                                                            <  1:978064/150‑1 (MQ=255)
                                                                         tttCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATc                                                                           >  1:333680/1‑151 (MQ=255)
                                                                          ttCAACCATGGCTTCCCGTCCTTCGCAAAAAAAGGGCGATCCTGCAGGCGGGCCACCCCAGCCAACAAGGCTGAGAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACg                                                                              <  2:56870/148‑1 (MQ=255)
                                                                          tcaaaCATGTCTGGCCGTACTTCGTAAAAATAAGGCGATACTGCATGCTGGACACACCAGTCATCAATGCTGAAAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGAt                                                                            <  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)
                                                                                                                                                       |                                                                                                                                                  
AGATATGGCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCGA  >  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: