Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,349,1450AC38.3% 66.1 / 10.3 47F110V (TTC→GTC) inaAacid‑inducible Kdo/WaaP family putative kinase
Reads supporting (aligned to +/- strand):  ref base A (24/5);  new base C (0/18);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 7.37e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCG  >  NC_000913/2349003‑2349290
                                                                                                                                               |                                                                                                                                                  
cGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCa                                                                                                                                           <  2:504804/151‑1 (MQ=255)
  aaCATAACAACAGCCTTGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATCCTGCATGGTGGGCATACCAGTCAGCATTGCTGAAGAC‑CCCCGCCat                                                                                                                                          <  1:100931/150‑1 (MQ=255)
  aaCATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCat                                                                                                                                          >  2:607984/1‑150 (MQ=255)
     aaaccaccaGCCATGCTGACTATTAATGTTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTCCATGCTGGGCATCCCAGTCAGAATTGCTGATGAA‑CCCCGCCat                                                                                                                                          <  1:452643/140‑1 (MQ=255)
     aaaacaacAGCCATGTTGACGATTAATGCTATGCATTTTCTAAAACGCCAGCCCCCCTGCTTCC‑AACAGGGCTTGCCGTACTTCGTAAGAATAAGGCGATCCGGCATGCTGGGCATCCCAGTCAGAAGTGCTGATGAC‑CCCCGCCATATc                                                                                                                                       <  1:480012/148‑1 (MQ=255)
     aTAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCt                                                                                                                                      >  2:508468/1‑151 (MQ=255)
      ttacaacaGCCATGCTGACGATTAATGCTATGCATATTCTTAAATGTCAGCGCCACTGCTTTC‑AACATGGCTCGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGCCAGCAATGCTGATGAA‑CCCCGCCATATCt                                                                                                                                      >  1:600595/3‑150 (MQ=255)
                       actaTTAATGCTGGTTATTTTCTAAAACCCCAGCGCCCCTGTTTAC‑AACAGGGCTTGCCTTATTTCGTCAAAAAAAGGCGATCCTCCATGCGGGGCATCCCAGCCAGCAAGGCTGAGGCA‑CCCCGCCATATCTTCAGTCACCAGCAACGc                                                                                                                     <  1:47939/147‑1 (MQ=255)
                          aTTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGc                                                                                                                 >  2:366904/1‑151 (MQ=255)
                           ttAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGcc                                                                                                                >  2:233265/1‑151 (MQ=255)
                            tAATGCTATGCTTTTTCTAAAACCCCTCCCCCACTTTTTTC‑AATATGGTTTGCTGTTCTCCGTAAAAAAAAGGCGATCCTGCATTTTGGCCATACCACTCACAAATCCTGATGAC‑CCCCGCCATATCTTCAGTCACCAGCAAc                                                                                                                       <  1:627907/143‑1 (MQ=255)
                              aTGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCAtt                                                                                                             >  1:445447/1‑151 (MQ=255)
                                    tGCTTTTTCTAAAACGCCAGCGCCACTGCTTTC‑AACATGGTTGGCCGTATTTCGTCAGAATAAGGCGTTACTGCATGCGGGGCATACCAGCCAGAAATGCTGATGAC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCAcc                                                                                                         <  2:637274/149‑1 (MQ=255)
                                    cgCATTTTCTTAAACGCCAGCACCACTGCTCTC‑AACACGGCCTGCCGTCCTTCGTCAGAATAAAACGATACCTCATGCGGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCTCCATATCGTCAGTCAACAGCAACGCTCGCCATTCACCCt                                                                                                       >  1:448277/2‑151 (MQ=255)
                                      cATTTTTTTAACCGCCAGCCCCCCTGCTTTC‑AACGTGGTTTCCCGTACTTCCTAAGAATAAGGCCATACTGCATTCGGGGCATACCAGCCAGTATTGCTTATGAC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCAc                                                                                                          <  1:499089/146‑1 (MQ=255)
                                      cATTTTCTTAAACGCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTc                                                                                                      >  2:411521/1‑150 (MQ=255)
                                           tCTTAAACGCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAAtttt                                                                                                >  2:355148/1‑151 (MQ=255)
                                               aacccccAGCCCCCCTCCTTCC‑ACCTGGGTTTGCCGTCCTTCGTCAAAAAAGGGCTCTCCGCCGTGCTGGGAATACCAGCCAGCAATGCTGATGCC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAAtttt                                                                                                <  1:607495/142‑1 (MQ=255)
                                                 aCGCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCACTTTTCACCGc                                                                                          >  1:570592/1‑151 (MQ=255)
                                                  ccccACCGCCTCTGCTTAC‑ACCATGTCTTCCCGTCTTCCGTAAGAAAAAGGCTATCCTCCATGCGGGCCATCCAATTCGGAAATGCTAATAAC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCAcc                                                                                            <  1:569197/146‑1 (MQ=255)
                                                  cGCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCgat                                                                                         >  1:115878/1‑149 (MQ=255)
                                                   gCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCtt                                                                                        >  2:461934/1‑151 (MQ=255)
                                                      aGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGACATTCACCCTCAATTTTCACCGCTTCAc                                                                                     >  1:414716/1‑151 (MQ=255)
                                                       gcgcCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCAcc                                                                                    >  2:546466/1‑151 (MQ=255)
                                                        cgcCACTTCTTTA‑AACAGGGCTGGCCGTCCTTCGTCAAAAAAAGGCGATACTGCATGCTGGGCATACCAGCCAGCAATGCTGATGAC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCa                                                                                   <  2:443902/151‑1 (MQ=255)
                                                          ccACTGTTTTC‑AACTTGCTTTGCTTTCCTTCTTAAAAAAAAGGCTATCCTCCATGCGGGCCATACCAGTCAGCAATGCGGAGAAC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaa                                                                                  <  2:157686/150‑1 (MQ=255)
                                                           cctcgCTTTCAAACATGGCCTGCCCTACTCCGACAGAAAACTGCAACCCTGCATGCTGGCCACCCCAGTCAGCAATCCTAATAAC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTGCACCGCTTCTCCa                                                                                   <  1:566927/145‑1 (MQ=255)
                                                            aCTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaaaa                                                                                >  1:109219/1‑150 (MQ=255)
                                                                      aaCTTGGTTTGCCTTTTTGCGTCAAAATAGGGCGATACTGCATGTGGGGAATACCAGACAGCAATGCTGATGAC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCt                                                                         <  2:126401/148‑1 (MQ=255)
                                                                      aaCAGGGCTTCCCTTCCTCCGTAAAAAAAAGGCGATACGGCATGCGGGGCATACCAGTCAGAAATGCTGAGGAC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCtt                                                                        <  2:599553/149‑1 (MQ=255)
                                                                          aggCTTGCCGTACTTAGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCa                                                                   >  1:324211/2‑150 (MQ=255)
                                                                                  cGTCCTTCGTCAAAATAGGGAGATACTGCTTGCGGGGCAACCCACTCAGAAAGGCTGATGAC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGacac                                                          <  2:254874/151‑1 (MQ=255)
                                                                                       ttCGTAAAAATAGGGCGATACTGCATGCTGGGCATCCCAGCCAGCAATGCTGATGAC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCg                                                        <  2:175581/148‑1 (MQ=255)
                                                                                         cgtaaaaaaagggCGATCCTGCATGCGGGGCATCCCAGCCAGAAATGCTGATAAC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGccc                                                    <  2:648997/139‑1 (MQ=255)
                                                                                         cGTAAGAATAAGGCGATACTGCATGCGGGCCATACCAGTCAGCAATGCTGATAAC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCg                                                   <  2:357293/151‑1 (MQ=255)
                                                                                             aGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTcc                                               >  2:596993/1‑151 (MQ=255)
                                                                                                    ggCGTTCCTGCATGCGGGGCATCCCATTCAGAAAT‑CTGGAGAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCtt                                        <  1:354767/151‑1 (MQ=255)
                                                                                                     gCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCttt                                       >  1:52622/1‑151 (MQ=255)
                                                                                                           cTGCATGCTGGGCATACCAGTAAGCCATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTATTTAc                                 >  1:614378/1‑151 (MQ=255)
                                                                                                                tGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCa                             >  1:450632/1‑150 (MQ=255)
                                                                                                                 gCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAAc                           >  2:592730/1‑151 (MQ=255)
                                                                                                                    gggCATCCAAGCCAGAAATGCTGATAAC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTc                        <  2:639765/151‑1 (MQ=255)
                                                                                                                     gccATCCCAGTCAGCAATGCTGATGAC‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCa                       <  2:433843/149‑1 (MQ=255)
                                                                                                                              gccAGCAATGCTGATGCA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACaa                <  2:583739/147‑1 (MQ=255)
                                                                                                                                       gCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCgg     >  1:576578/1‑151 (MQ=255)
                                                                                                                                        cTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGATCGCCATTCACCCTCATTTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACAGCAACCTCACGGACAATCGTTGGTCGGc    >  1:271879/1‑151 (MQ=255)
                                                                                                                                           atgaA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCg  >  1:260039/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTC‑AACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAA‑CCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCG  >  NC_000913/2349003‑2349290

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