Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913362,7300GA47.2% 9.3 / 40.9 36G150G (GGC→GGT) lacYlactose permease
Reads supporting (aligned to +/- strand):  ref base G (9/10);  new base A (8/9);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AACATACAGTGACAAAAACCACAGTTTTGGCTGTCTGAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTT  >  NC_000913/362493‑362974
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
aaCATACAGTGACAAAAACCACAGTTTTGGCTGTCTGAACAGTTCCAGTGCCAGCTTAATGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCa                                                                                                                                                                                                                                          >  1:65322/1‑250 (MQ=255)
aaCATACAGTGACAAAAACCACAGTTTTGGCTGTCTGAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCaa                                                                                                                                                                                                                                         <  1:183008/251‑1 (MQ=255)
            cAAAAACCACAGTTTTGGCTGTCTGAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCcgcgcgc                                                                                                                                                                                                                             >  2:271931/1‑251 (MQ=255)
                aaCCACAGTTTTGGCTGTCTGAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGAcc                                                                                                                                                                                                                         <  2:216981/251‑1 (MQ=255)
                aaCCACAGTTTTGGCTGTCTGAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCGACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATGGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGAcc                                                                                                                                                                                                                         <  2:142220/251‑1 (MQ=255)
                           gggCTGTCTGAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAAt                                                                                                                                                                                                              <  2:34965/250‑1 (MQ=255)
                                    gAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGc                                                                                                                                                                                                        <  1:271931/248‑1 (MQ=255)
                                               gTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTtc                                                                                                                                                                                          >  2:146544/1‑251 (MQ=255)
                                                    aGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAAt                                                                                                                                                                                     >  2:116652/1‑251 (MQ=255)
                                                    aGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAAt                                                                                                                                                                                     >  1:220332/1‑251 (MQ=255)
                                                    aGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAAt                                                                                                                                                                                     <  1:232687/251‑1 (MQ=255)
                                                                aaTGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTAc                                                                                                                                                                         >  1:235864/1‑251 (MQ=255)
                                                                  tGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTActg                                                                                                                                                                       >  1:153443/1‑251 (MQ=255)
                                                                           gTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCAc                                                                                                                                                              <  1:5065/251‑1 (MQ=255)
                                                                              ggCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAACGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCgg                                                                                                                                                           >  1:175434/1‑251 (MQ=255)
                                                                                        gCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACa                                                                                                                                                 <  2:175434/251‑1 (MQ=255)
                                                                                            gggCAAACGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAg                                                                                                                                             <  2:243431/250‑1 (MQ=255)
                                                                                                      ggCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATaa                                                                                                                                    <  1:78011/250‑1 (MQ=255)
                                                                                                                     aTCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCctc                                                                                                                                                                            <  2:298496/195‑1 (MQ=255)
                                                                                                                     aTCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCctc                                                                                                                                                                            >  1:298496/1‑195 (MQ=255)
                                                                                                                        cGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGCGTGAACAACCAGACCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATc                                                                                                                 >  1:36072/1‑251 (MQ=255)
                                                                                                                        cGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATc                                                                                                                 <  1:121379/251‑1 (MQ=255)
                                                                                                                                  aaaaaGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAAGACCACCAACAATCGATCCTACTAAAAt                                                                                                       <  2:246313/251‑1 (MQ=255)
                                                                                                                                                ggCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCTACAAACGATCCTACTAAAATGTTGTATTGTAACa                                                                                         <  1:280092/251‑1 (MQ=255)
                                                                                                                                                   gAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTg                                                                                      <  1:15756/251‑1 (MQ=255)
                                                                                                                                                         gAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCg                                                                                 <  1:146544/250‑1 (MQ=255)
                                                                                                                                                             gCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATGGATGGTGAACACGATGCCGACAATCGAGTCACACAGCGCCCAGCCAACACAGCCAAACAGCCGCGCGCGACCAAATTCTAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGGTAAAACAAAAGCCTATATAAATACCACCAACAAGCGATCCTACTAAAATGTTGTATTGTAACAGTTGCCTGAAGa                                                                             <  2:20569/250‑1 (MQ=255)
                                                                                                                                                                     agagCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGAGACCAAATTCGACATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATa                                                                     >  1:292093/1‑250 (MQ=255)
                                                                                                                                                                            agccagAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACgg                                                             >  1:106125/1‑251 (MQ=255)
                                                                                                                                                                              ccagAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGgcg                                                           >  1:122265/1‑251 (MQ=255)
                                                                                                                                                                                                attGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCTCTAACATGc                                         <  2:165308/251‑1 (MQ=255)
                                                                                                                                                                                                             atgatgCCGACAATCGAGGCACACAGCGCCCAACCCACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCATACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATcc                             >  1:111000/1‑250 (MQ=255)
                                                                                                                                                                                                               gatgcccccAATCGAGGCACACAGCGCCCAACCCACACCGCCAAACATCCGCGCGCGACCAACTTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAACCAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCcac                           >  2:199740/9‑250 (MQ=255)
                                                                                                                                                                                                                atgCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCAcag                         <  1:237182/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTgcgc            >  1:103238/1‑250 (MQ=255)
                                                                                                                                                                                                                                       cgcCCAACCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGtt  >  2:6553/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    
AACATACAGTGACAAAAACCACAGTTTTGGCTGTCTGAACAGTTCCAGTGCCAGCTTAAGGCTAAATGCCGAATGGTTGGCACCTACCGCATTGGCAACCGTGGCAGAAGAGGGCGCATCCGTTTTGGCGAAAAAGAGTAAAACGGCGAGGATGAGTGCACAGCCAGAGCCCAGCCAGAAAACAAACTGATTATTGATGGTGAACATGATGCCGACAATCGAGGCACACAGCGCCCAGCCAACACAGCCAAACATCCGCGCGCGACCAAATTCGAAATTACTGCGACGGCTGACTTTCTCAATAAATGCCTCTACTGCTGGCGCACCGGCGTTAAAACAAAAGCCTAGATAAATACCACCAACAATCGATCCTACTAAAATGTTGTATTGTAACAGTGGCCCGAAGATAAAAATAAAGAACGGCGCAAACATCACTAACATGCCGGTAATAATCCACAGCAGGTATTTGCGCAGCCCGAGTT  >  NC_000913/362493‑362974

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