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

GCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGG  >  NC_000913/2349002‑2349287
                                                                                                                                                 |                                                                                                                                              
gCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGcc                                                                                                                                         >  2:38570/1‑151 (MQ=255)
                  aTGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGc                                                                                                                       >  2:855740/1‑151 (MQ=255)
                       gACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGc                                                                                                                  >  1:336215/1‑151 (MQ=255)
                         cGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTc                                                                                                                >  2:606217/1‑151 (MQ=255)
                            ttAATGTTTTGATTTTTCTAAACCGCCGCCCCCCTGGTTTTAAACTGGGTTTCCCGTCCTTCTTAAAAAAAGGGCGTTCCTGCATGCTGGCCATCCCAGACAGCATTGCTAAAGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCg                                                                                                               <  1:850393/149‑1 (MQ=255)
                            ttAATGCTATGCTTTTTCTTAACCCCCAGCGCCACTGTTTTCAACTTGGTTTGCCGTTCTTCGTCAGAATAAGGCGATCCTGCATGCTGGGCATCCCAGTCAGAAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCaa                                                                                                                     <  2:37207/143‑1 (MQ=255)
                              aaTGTTATGCTTTTTCTAAACCGCCGGCCCCCCTGTTTTAAACGTGCCTTCCCGTACTTCGTCAAAATAGGGCGATCCTGCAGGCGGGGAATCCCAGTCAGCAATGCTGAGAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGcc                                                                                                             <  2:542968/149‑1 (MQ=255)
                                  ctatttaTTTTCTAAACCCCCAGCGCCACTGCTTTCAACATGGCTTGCCGTTCTCCGTAAAAATAAGGCGATCCTGCATGCGGGGCATCCCAGACAGAATTGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCAtt                                                                                                          <  2:789035/142‑1 (MQ=255)
                                   tATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCAcc                                                                                                      >  1:476362/1‑151 (MQ=255)
                                                aacccccAGCGCCCTTCTTTTAAACATGGCTTCCCTTACTTCGCAAAAAAAGGGAATTCCTGCTTGCGGGGCAACCCAGCCAACAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCa                                                                                           <  1:19323/144‑1 (MQ=255)
                                                   cGCCAGCCCCCCTGCTTTCAACATGGCTTGCCGTACTTCGTCAAAATAAGGCGATACTCCTTGCGGGGCATCCCAGTCACCATTTCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTc                                                                                            <  1:850820/145‑1 (MQ=255)
                                                    gCCATCGCCACTGCTTTCAACATGGCTTGCCGGACTTCGTCAGAATAAGGCGATACTGCACGCTGGGCATCCCAGTTAGCAATCCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCAcc                                                                                         >  1:119167/1‑147 (MQ=255)
                                                     ccggccccccTGCTTCAAACAGGGTTTCCCTTACTTCTAAAAAAAAAGGCGATACTGCATGCTGGGCATCCCAGCCAGAAATGTTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTc                                                                                    <  2:272482/142‑1 (MQ=255)
                                                      cAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTc                                                                                    >  1:806311/1‑150 (MQ=255)
                                                         cgcCACTGTTTTCAACATGGCTTCCCGTACTTCGTAAAAATAGGGGGATACTGCATGCGGGCCATCCCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCa                                                                                <  2:582460/151‑1 (MQ=255)
                                                            cACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaaa                                                                              >  2:921731/1‑150 (MQ=255)
                                                                   ttCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCt                                                                      >  2:766969/1‑151 (MQ=255)
                                                                      aaCATGGCTTGCCGTACTTCGTCAAAATAAGGCGATACTGCATGCGGGGCATCCCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCt                                                                      <  2:702064/148‑1 (MQ=255)
                                                                       aCATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCgg                                                                  >  1:629203/1‑151 (MQ=255)
                                                                        cATGGTTCGCCGTACTTCCACAAAAGAGGGGTATCTCGCAGGCGGGTCATCCCAGAAAGAAATCCAGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCt                                                                      <  2:227622/146‑1 (MQ=255)
                                                                              ttGCCTTTCTTCGTAAAAATAAGGCGATCCTGCATGCGGGGCATACCAGTCAGCAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACg                                                              <  2:600723/148‑1 (MQ=255)
                                                                              ttGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATg                                                           >  2:76639/1‑151 (MQ=255)
                                                                                   cGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGacac                                                       >  2:726949/1‑151 (MQ=255)
                                                                                   cGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGacac                                                       >  1:280816/1‑151 (MQ=255)
                                                                                   cGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGacac                                                       >  1:357948/1‑151 (MQ=255)
                                                                                    gTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACAcc                                                      >  1:234327/1‑151 (MQ=255)
                                                                                        ttcttaaaaaTAAGGCTATACTGCGTGCTGGGCATACCAGCCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCgg                                                    <  2:337310/141‑1 (MQ=255)
                                                                                             cAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTc                                             >  2:131718/1‑151 (MQ=255)
                                                                                                aaaaaGGCTATCCGGCAGCCTGGCCATACCAGTCAGCAGTACACATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCg                                                <  1:166347/142‑1 (MQ=255)
                                                                                                aaTAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTcc                                            >  2:249504/1‑149 (MQ=255)
                                                                                                 aTAAGGCTATACTGCTGGCTGGGCTTCCCAGTCAGAATTGTGGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCa                                           <  2:254120/149‑1 (MQ=255)
                                                                                                     ggCGATGCTGCGGGCGGGGCATCCCAGACAGCAATGCGGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTc                                       <  1:411671/149‑1 (MQ=255)
                                                                                                     ggCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCtt                                     >  1:2433/1‑151 (MQ=255)
                                                                                                     gccGATCCTGCATGCGGGGCAACCCAGCCAGAATTGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCtt                                     <  2:292267/149‑1 (MQ=255)
                                                                                                         aTACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAAt                                 >  2:752251/1‑151 (MQ=255)
                                                                                                                     tGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCt                      >  1:176300/1‑151 (MQ=255)
                                                                                                                      gggCATCCCAGCCAGCATGGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGc                           <  2:339198/145‑1 (MQ=255)
                                                                                                                      gggCATCCCAGCCAGCAATGCGGAACACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTc                     <  1:594882/151‑1 (MQ=255)
                                                                                                                         cATCCCAGCCAGCAATGTTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCt                      <  1:809867/147‑1 (MQ=255)
                                                                                                                               aGTCAGCAATGCTAATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAAt            <  1:542886/151‑1 (MQ=255)
                                                                                                                                gacAGCATTGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATc           <  1:580657/149‑1 (MQ=255)
                                                                                                                                             atgcACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCgg  <  1:825578/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
GCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTGCCGTACTTCGTCAGAATAAGGCGATACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGG  >  NC_000913/2349002‑2349287

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