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

GCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGG  >  NC_000913/2349002‑2349287
                                                                                                                                                 |                                                                                                                                              
gCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGcc                                                                                                                                         >  2:38570/1‑151 (MQ=255)
                  aTGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGc                                                                                                                       >  2:855740/1‑151 (MQ=255)
                       gACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGc                                                                                                                  >  1:336215/1‑151 (MQ=255)
                         cGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTc                                                                                                                >  2:606217/1‑151 (MQ=255)
                            ttAATGTTTTGATTTTTCTAAACCGCC‑GCCCCCCTGGTTTTAAACTGGGTTTCCCGTCCTTCTTAAAAAAAGGGCGT‑TCCTGCATGCTGGCCATCCCAGACAGCATTGCTAAAGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCg                                                                                                               <  1:850393/149‑1 (MQ=255)
                            ttAATGCTATGCTTTTTCTTAACCCCCAGCGCCACTGTTTTCAACTTGGTTTG‑CCGTTCTTCGTCAGAATAAGGCGA‑TCCTGCATGCTGGGCATCCCAGTCAGAAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCaa                                                                                                                     <  2:37207/143‑1 (MQ=255)
                              aaTGTTATGCTTTTTCTAAACCGCCGGCCCCCCTGTTTTAAACGTGCCTTC‑CCGTACTTCGTCAAAATAGGGCGA‑TCCTGCAGGCGGGGAATCCCAGTCAGCAATGCTGAGAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGcc                                                                                                             <  2:542968/149‑1 (MQ=255)
                                  ctatttaTTTTCTAAACCCCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTTCTCCGTAAAAATAAGGCGA‑TCCTGCATGCGGGGCATCCCAGACAGAATTGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCAtt                                                                                                          <  2:789035/142‑1 (MQ=255)
                                   tATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCAcc                                                                                                      >  1:476362/1‑151 (MQ=255)
                                                aacccccAGCGCCCTTCTTTTAAACATGGCTTC‑CCTTACTTCG‑CAAAAAAAGGGAATTCCTGCTTGCGGGGCAACCCAGCCAACAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCa                                                                                           <  1:19323/144‑1 (MQ=255)
                                                   cGCCAGCCCCCCTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAAAATAAGGCGA‑TACTCCTTGCGGGGCATCCCAGTCACCATTTCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTc                                                                                            <  1:850820/145‑1 (MQ=255)
                                                    gCCATCGCCACTGCTTTCAACATGGCTTG‑CCGGACTTCGTCAGAATAAGGCGA‑TACTGCACGCTGGGCATCCCAGTTAGCAATCCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCAcc                                                                                         >  1:119167/1‑147 (MQ=255)
                                                     ccggccccccTGCTTCAAACAGGGTTTC‑CCTTACTTCTAAAAAAAAAGGCGA‑TACTGCATGCTGGGCATCCCAGCCAGAAATGTTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTc                                                                                    <  2:272482/142‑1 (MQ=255)
                                                      cAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTc                                                                                    >  1:806311/1‑150 (MQ=255)
                                                         cgcCACTGTTTTCAACATGGCTTC‑CCGTACTTCGTAAAAATAGGGGGA‑TACTGCATGCGGGCCATCCCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCa                                                                                <  2:582460/151‑1 (MQ=255)
                                                            cACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCaaa                                                                              >  2:921731/1‑150 (MQ=255)
                                                                   ttCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCt                                                                      >  2:766969/1‑151 (MQ=255)
                                                                      aaCATGGCTTG‑CCGTACTTCGTCAAAATAAGGCGA‑TACTGCATGCGGGGCATCCCAGTCAGCAATGCTGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCt                                                                      <  2:702064/148‑1 (MQ=255)
                                                                       aCATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCgg                                                                  >  1:629203/1‑151 (MQ=255)
                                                                        cATGGTTCG‑CCGTACTTCCACAAAAGAGGGGTA‑TCTCGCAGGCGGGTCATCCCAGAAAGAAATCCAGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCt                                                                      <  2:227622/146‑1 (MQ=255)
                                                                              ttG‑CCTTTCTTCGTAAAAATAAGGCGA‑TCCTGCATGCGGGGCATACCAGTCAGCAATGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACg                                                              <  2:600723/148‑1 (MQ=255)
                                                                              ttG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATg                                                           >  2:76639/1‑151 (MQ=255)
                                                                                   cGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGacac                                                       >  2:726949/1‑151 (MQ=255)
                                                                                   cGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGacac                                                       >  1:280816/1‑151 (MQ=255)
                                                                                   cGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGacac                                                       >  1:357948/1‑151 (MQ=255)
                                                                                    gTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACAcc                                                      >  1:234327/1‑151 (MQ=255)
                                                                                        ttcttaaaaaTAAGGCTA‑TACTGCGTGCTGGGCATACCAGCCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCgg                                                    <  2:337310/141‑1 (MQ=255)
                                                                                             cAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTc                                             >  2:131718/1‑151 (MQ=255)
                                                                                                aaaaaGGCTA‑TCCGGCAGCCTGGCCATACCAGTCAGCAGTACACATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCg                                                <  1:166347/142‑1 (MQ=255)
                                                                                                aaTAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTcc                                            >  2:249504/1‑149 (MQ=255)
                                                                                                 aTAAGGCTA‑TACTGCTGGCTGGGCTTCCCAGTCAGAATTGTGGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCa                                           <  2:254120/149‑1 (MQ=255)
                                                                                                     ggCGA‑TGCTGCGGGCGGGGCATCCCAGACAGCAATGCGGATAACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTc                                       <  1:411671/149‑1 (MQ=255)
                                                                                                     ggCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCtt                                     >  1:2433/1‑151 (MQ=255)
                                                                                                     gccGA‑TCCTGCATGCGGGGCAACCCAGCCAGAATTGCTGATGACCCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCtt                                     <  2:292267/149‑1 (MQ=255)
                                                                                                         a‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAAt                                 >  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)
                                                                                                                                                 |                                                                                                                                              
GCGAACATAACAACAGCCATGCTGACGATTAATGCTATGCATTTTCTTAAACGCCAGCGCCACTGCTTTCAACATGGCTTG‑CCGTACTTCGTCAGAATAAGGCGA‑TACTGCATGCTGGGCATACCAGTCAGCAATGCTGATGAACCCCGCCATATCTTCAGTCACCAGCAACGCTCGCCATTCACCCTCAATTTTCACCGCTTCACCAAAAACGATCTTCGGCACGATGACACCGGCCCGTTCCAGTTCTTTAATTACCGCAACCTCACGGACAATCGTTGGTCGG  >  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: ‑