Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*U00096445,1480TA20.0% 100.9 / 11.1 45F156I (TTC→ATC) yajQnucleotide binding protein YajQ
Reads supporting (aligned to +/- strand):  ref base T (16/20);  new base A (9/0);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 2.49e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGAAAATCGTCAAGATGATCAAAGACAGCAAACTGAAAGTGCAGGCGCAAATTCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGTTTGTAATTCAGAAACTTTGTAGGCCTGATAAGCGTAGCGCAGCAGGCAATTCGCC  >  U00096/445003‑445294
                                                                                                                                                 |                                                                                                                                                  
gaagaaAATCGTCAAGATGATCAAAGACAGCAAACTGAAAGTGCAGGCGCAAATTCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTc                                                                                                                                                  >  1:94319/1‑148 (MQ=255)
 aagaaAATCGTCAAGATGATCAAAGACAGCAAACTGAAAGTGCAGGCGCAAATTCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTcc                                                                                                                                                 >  2:94317/1‑148 (MQ=255)
   gaaAATCGTCAAGATGATCAAAGACAGCAAACTGAAAGTGCAGGCGCAAATTCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAg                                                                                                                                               <  2:488127/148‑1 (MQ=255)
                  gatCAAAGACAGCAAACTGAAAGTGCAGGCGCAAATTCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCTATGGTACGTGGTGGCGATCTCGGTCAGCCGATACAGTTCAAAAACTTCcgc                                                                                                                                >  2:94325/1‑148 (MQ=255)
                    tCAAAGACAGCAAACTGAAAGTGCAGGCGCAAATTCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGa                                                                                                                              <  2:94308/148‑1 (MQ=255)
                      aaaGACAGCAAACTGAAAGTGCAGGCGCAAATTCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGAtt                                                                                                                            <  1:94309/148‑1 (MQ=255)
                             gCAAACTGAAAGTGCAGGCGCAAATTCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATcgcg                                                                                                                     >  1:94336/1‑148 (MQ=255)
                                  cTGAAAGTGCAGGCGCAAATTCAGGGCGATGAAATCCGCGGAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGGGGTGGCGATCTCGGTCCGCCGTTCCAGTTCAAAACCTTCCGCGATTCATCGCGACTAt                                                                                                                >  1:94330/1‑148 (MQ=255)
                                     aaaGTGCAGGCGCAAATTCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTAtgtg                                                                                                             <  1:94310/148‑1 (MQ=255)
                                       aGTGCAGGCGCAAATTCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGTAGATCTCGGTCATCCCATCCACATAAAAAACATCCGCGATTAATCGCGACTATGTGCg                                                                                                           >  2:94322/1‑148 (MQ=255)
                                            aGGCGCAAATTCAGGGCGATGAAATCCGCGTATCGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTg                                                                                                      <  2:94311/148‑1 (MQ=255)
                                               cgcAAATTCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGTCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTc                                                                                                   >  2:94324/1‑148 (MQ=255)
                                                    aTTCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTc                                                                                                                                   <  1:94312/111‑1 (MQ=255)
                                                    aTTCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGAGCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTc                                                                                              <  2:94313/148‑1 (MQ=255)
                                                      tCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGTAGATCTCGGTCATCCCATAAACATAAAAAACTACAGCGATTAATCGCGACTATGTGCGTTTTGGTCATGTCAg                                                                                            >  2:94339/1‑148 (MQ=255)
                                                       cAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGTAGATCTCGGTCATCCCATAAACATAAAAAAAATCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTc                                                                                              >  2:94326/1‑145 (MQ=255)
                                                       cAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTc                                                                                              >  2:94335/1‑145 (MQ=255)
                                                                        cgTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGTAGATCTCGGTCATCCGTTCCACTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTg                                                                          >  2:94346/1‑148 (MQ=255)
                                                                        cgTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGcc                                                                            >  1:94333/1‑146 (MQ=255)
                                                                         gTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGa                                                                         <  2:94314/148‑1 (MQ=255)
                                                                          tAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTa                                                                                                                           >  1:94329/1‑97 (MQ=255)
                                                                           aaCGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGAtt                                                                       <  2:94315/148‑1 (MQ=255)
                                                                            aCGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTc                                                                      <  1:94317/148‑1 (MQ=255)
                                                                            aCGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTc                                                                      <  1:94316/148‑1 (MQ=255)
                                                                             cGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGTAGATCTCGGTCATCCGATACAGATAAAAAACATCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCg                                                                     >  2:94332/1‑148 (MQ=255)
                                                                               ggCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGc                                                                   <  2:94318/148‑1 (MQ=255)
                                                                                   aaTCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTAc                                                               <  2:94319/148‑1 (MQ=255)
                                                                                     tctcGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACaa                                                             >  1:94323/1‑148 (MQ=255)
                                                                                      ctcGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGTCGATCTCGGTCATCCCATCAACTTAAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACaaa                                                            >  2:94343/1‑148 (MQ=255)
                                                                                        cGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACaaa                                                            >  1:94331/1‑146 (MQ=255)
                                                                                        cGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCGCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACaa                                                             <  1:94320/145‑1 (MQ=255)
                                                                                         gTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGt                                                         <  2:94321/148‑1 (MQ=255)
                                                                                          tgatgaTTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGtt                                                        >  1:94328/1‑148 (MQ=255)
                                                                                                 ttGCAGGCTGTCATGGCGATGGTACGTGGTGTAGATCTCGGTCATCCCATAAACATAAAAAAAAAACGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGTTTGTAAtt                                                 >  2:94341/1‑148 (MQ=255)
                                                                                                            cATGGCGATGGTACGTGGTGTAGATCTCGGTCATCACATAAACATAAAAAAAAAAAACGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGTTTGTAATTCAGAAACTTTg                                      >  2:94337/1‑148 (MQ=255)
                                                                                                               ggCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGTTTGTAATTCAGAAACTTTGTAg                                   >  1:94354/1‑148 (MQ=255)
                                                                                                                   aTGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGTTTGTAATTCAGAAACTTTGTAGGCCt                               >  1:94338/1‑148 (MQ=255)
                                                                                                                   aTGGTACGTGGGGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGGTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGTTTGTAATTCAGAAACTTTGTAGGCCt                               <  1:94322/148‑1 (MQ=255)
                                                                                                                     ggTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGTTTGTAATTCAGAAACTTTGTAGGCCTGa                             >  1:94344/1‑148 (MQ=255)
                                                                                                                         cGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGTTTGTAATTCAGAAACTTTGTAGGCCTGATAAg                         <  2:94323/148‑1 (MQ=255)
                                                                                                                            ggtgtagATCTCGGTCATCCCATACACTTAAAAAAAAACCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGTTTGTAATTCAGAAACTTTGTAGGCCTGATAAGCGt                      >  2:591870/7‑148 (MQ=255)
                                                                                                                                            aGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGTTTGTAATTCAGAAACTTTGTAGGCCTGATAAGCGt                      <  1:94324/132‑1 (MQ=255)
                                                                                                                                              ccGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGTTTGTAATTCAGAAACTTTGTAGGCCTGATAAGCGTAGCGCAGCAGGCAAtt      >  1:94356/1‑146 (MQ=255)
                                                                                                                                                gTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGTTTGTAATTCAGAAACTTTGTAGGCCTGATAAGCGTAgcg                  <  1:94325/132‑1 (MQ=255)
                                                                                                                                                gTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGTTTGTAATTCAGAAACTTTGTAGGCCTGATAAGCGTAGCGCAGCAGGCAATTcgcc  <  1:94326/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GAAGAAAATCGTCAAGATGATCAAAGACAGCAAACTGAAAGTGCAGGCGCAAATTCAGGGCGATGAAATCCGCGTAACGGGCAAATCTCGTGATGATTTGCAGGCTGTCATGGCGATGGTACGTGGTGGCGATCTCGGTCAGCCGTTCCAGTTCAAAAACTTCCGCGATTAATCGCGACTATGTGCGTTTTGTTCATGTCAGATGCGGCATGAACGCCTGATTCGGCCTACAAAGGTTTGTAATTCAGAAACTTTGTAGGCCTGATAAGCGTAGCGCAGCAGGCAATTCGCC  >  U00096/445003‑445294

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