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

CGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCG  >  NC_000913/2349003‑2349290
                                                                                                                                               |                                                                                                                                                  
cGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCa                                                                                                                                           <  2:504804/151‑1 (MQ=255)
  aaCATAACAACAGCCTTGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATCCTGCATGGTGGGCATACCAGTCAGCATTGCTGAAGACCCCCGCCat                                                                                                                                          <  1:100931/150‑1 (MQ=255)
  aaCATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCat                                                                                                                                          >  2:607984/1‑150 (MQ=255)
     aaaccaccaGCCATGCTGACTATTAATGTTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTCCATGCTGGGCATCCCAGTCAGAATTGCTGATGAACCCCGCCat                                                                                                                                          <  1:452643/140‑1 (MQ=255)
     aaaacaacAGCCATGTTGACGATTAATGCTATGCATTTTCTAAAACGCCAGCCCCCCTGCTTCCAACAGGGCTTGCCGTACTTCGTAAGAATAAGGCGATCCGGCATGCTGGGCATCCCAGTCAGAAGTGCTGATGACCCCCGCCATATc                                                                                                                                       <  1:480012/148‑1 (MQ=255)
     aTAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCt                                                                                                                                      >  2:508468/1‑151 (MQ=255)
      ttacaacaGCCATGCTGACGATTAATGCTATGCATATTCTTAAATGTCAGCGCCACTGCTTTCAACATGGCTCGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGCCAGCAATGCTGATGAACCCCGCCATATCt                                                                                                                                      >  1:600595/3‑150 (MQ=255)
                       actaTTAATGCTGGTTATTTTCTAAAACCCCAGCGCCCCTGTTTACAACAGGGCTTGCCTTATTTCGTCAAAAAAAGGCGATCCTCCATGCGGGGCATCCCAGCCAGCAAGGCTGAGGCACCCCGCCATATCTTCAGTCACCAGCAACGc                                                                                                                     <  1:47939/147‑1 (MQ=255)
                          aTTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGc                                                                                                                 >  2:366904/1‑151 (MQ=255)
                           ttAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGcc                                                                                                                >  2:233265/1‑151 (MQ=255)
                            tAATGCTATGCTTTTTCTAAAACCCCTCCCCCACTTTTTTCAATATGGTTTGCTGTTCTCCGTAAAAAAAAGGCGATCCTGCATTTTGGCCATACCACTCACAAATCCTGATGACCCCCGCCATATCTTCAGTCACCAGCAAc                                                                                                                       <  1:627907/143‑1 (MQ=255)
                              aTGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCAtt                                                                                                             >  1:445447/1‑151 (MQ=255)
                                    tGCTTTTTCTAAAACGCCAGCGCCACTGCTTTCAACATGGTTGGCCGTATTTCGTCAGAATAAGGCGTTACTGCATGCGGGGCATACCAGCCAGAAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCAcc                                                                                                         <  2:637274/149‑1 (MQ=255)
                                    cgCATTTTCTTAAACGCCAGCACCACTGCTCTCAACACGGCCTGCCGTCCTTCGTCAGAATAAAACGATACCTCATGCGGGGCATACCAGTCAGCAATGCTGATGAACCCCTCCATATCGTCAGTCAACAGCAACGCTCGCCATTCACCCt                                                                                                       >  1:448277/2‑151 (MQ=255)
                                      cATTTTTTTAACCGCCAGCCCCCCTGCTTTCAACGTGGTTTCCCGTACTTCCTAAGAATAAGGCCATACTGCATTCGGGGCATACCAGCCAGTATTGCTTATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCAc                                                                                                          <  1:499089/146‑1 (MQ=255)
                                      cATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTc                                                                                                      >  2:411521/1‑150 (MQ=255)
                                           tCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAAtttt                                                                                                >  2:355148/1‑151 (MQ=255)
                                               aacccccAGCCCCCCTCCTTCCACCTGGGTTTGCCGTCCTTCGTCAAAAAAGGGCTCTCCGCCGTGCTGGGAATACCAGCCAGCAATGCTGATGCCCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAAtttt                                                                                                <  1:607495/142‑1 (MQ=255)
                                                 aCGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCACTTTTCACCGc                                                                                          >  1:570592/1‑151 (MQ=255)
                                                  ccccACCGCCTCTGCTTACACCATGTCTTCCCGTCTTCCGTAAGAAAAAGGCTATCCTCCATGCGGGCCATCCAATTCGGAAATGCTAATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCAcc                                                                                            <  1:569197/146‑1 (MQ=255)
                                                  cGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCgat                                                                                         >  1:115878/1‑149 (MQ=255)
                                                   gCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCtt                                                                                        >  2:461934/1‑151 (MQ=255)
                                                      aGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGACATTCACCCTCAATTTTCACCGCTTCAc                                                                                     >  1:414716/1‑151 (MQ=255)
                                                       gcgcCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCAcc                                                                                    >  2:546466/1‑151 (MQ=255)
                                                        cgcCACTTCTTTAAACAGGGCTGGCCGTCCTTCGTCAAAAAAAGGCGATACTGCATGCTGGGCATACCAGCCAGCAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCa                                                                                   <  2:443902/151‑1 (MQ=255)
                                                          ccACTGTTTTCAACTTGCTTTGCTTTCCTTCTTAAAAAAAAGGCTATCCTCCATGCGGGCCATACCAGTCAGCAATGCGGAGAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaa                                                                                  <  2:157686/150‑1 (MQ=255)
                                                           cctcgCTTTCAAACATGGCCTGCCCTACTCCGACAGAAAACTGCAACCCTGCATGCTGGCCACCCCAGTCAGCAATCCTAATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTGCACCGCTTCTCCa                                                                                   <  1:566927/145‑1 (MQ=255)
                                                            aCTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaaaa                                                                                >  1:109219/1‑150 (MQ=255)
                                                                      aaCTTGGTTTGCCTTTTTGCGTCAAAATAGGGCGATACTGCATGTGGGGAATACCAGACAGCAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCt                                                                         <  2:126401/148‑1 (MQ=255)
                                                                      aaCAGGGCTTCCCTTCCTCCGTAAAAAAAAGGCGATACGGCATGCGGGGCATACCAGTCAGAAATGCTGAGGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCtt                                                                        <  2:599553/149‑1 (MQ=255)
                                                                          aggCTTGCCGTACTTAGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCa                                                                   >  1:324211/2‑150 (MQ=255)
                                                                                  cGTCCTTCGTCAAAATAGGGAGATACTGCTTGCGGGGCAACCCACTCAGAAAGGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGacac                                                          <  2:254874/151‑1 (MQ=255)
                                                                                       ttCGTAAAAATAGGGCGATACTGCATGCTGGGCATCCCAGCCAGCAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCg                                                        <  2:175581/148‑1 (MQ=255)
                                                                                         cgtaaaaaaagggCGATCCTGCATGCGGGGCATCCCAGCCAGAAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGccc                                                    <  2:648997/139‑1 (MQ=255)
                                                                                         cGTAAGAATAAGGCGATACTGCATGCGGGCCATACCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCg                                                   <  2:357293/151‑1 (MQ=255)
                                                                                             aGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTcc                                               >  2:596993/1‑151 (MQ=255)
                                                                                                    ggCGTTCCTGCATGCGGGGCATCCCATTCAGAAATCTGGAGAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCtt                                        <  1:354767/151‑1 (MQ=255)
                                                                                                     gCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCttt                                       >  1:52622/1‑151 (MQ=255)
                                                                                                           cTGCATGCTGGGCATACCAGTAAGCCATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTATTTAc                                 >  1:614378/1‑151 (MQ=255)
                                                                                                                tGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCa                             >  1:450632/1‑150 (MQ=255)
                                                                                                                 gCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAAc                           >  2:592730/1‑151 (MQ=255)
                                                                                                                    gggCATCCAAGCCAGAAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTc                        <  2:639765/151‑1 (MQ=255)
                                                                                                                     gccATCCCAGTCAGCAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCa                       <  2:433843/149‑1 (MQ=255)
                                                                                                                              gccAGCAATGCTGATGCACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACaa                <  2:583739/147‑1 (MQ=255)
                                                                                                                                       gCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCgg     >  1:576578/1‑151 (MQ=255)
                                                                                                                                        cTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGATCGCCATTCACCCTCATTTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACAGCAACCTCACGGACAATCGTTGGTCGGc    >  1:271879/1‑151 (MQ=255)
                                                                                                                                           atgaACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCg  >  1:260039/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGGCCG  >  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: