Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,401,4590TC34.3% 44.8 / 11.5 35D625G (GAT→GGT) csrDregulator of CsrB and CsrC decay
Reads supporting (aligned to +/- strand):  ref base T (17/6);  new base C (0/12);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 2.97e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTC  >  NC_000913/3401321‑3401590
                                                                                                                                          |                                                                                                                                   
gCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATc                                                                                                                                    <  1:63072/140‑1 (MQ=255)
   aGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATcccc                                                                                                                                 >  2:674993/1‑140 (MQ=255)
   aGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATcccc                                                                                                                                 >  1:76440/1‑140 (MQ=255)
       gACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGc                                                                                                                             <  1:708463/140‑1 (MQ=255)
        aCAAGCCCACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGcc                                                                                                                            <  1:460923/140‑1 (MQ=255)
         cAAGCACTCTCCCATGAGACGCAGCGCGATTTATTCTCGGTAAAACCGGATTAAAGGGCAGGTTAAACCGAGTATCTTTGTAAATTTTTTTTCACTTTAGTATAAGGTGGCGGTGAGGCGGCAAAAAAACCCCCTTGtcc                                                                                                                           <  2:76440/140‑1 (MQ=255)
         cAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGccc                                                                                                                           <  2:104545/140‑1 (MQ=255)
                  tCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAac                                                                                                                  >  1:737040/1‑140 (MQ=255)
                      aTGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACcg                                                                                                              <  1:674993/140‑1 (MQ=255)
                      aTGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACcg                                                                                                              >  1:353549/1‑140 (MQ=255)
                      aTGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCACCTTGCCCGCCTGTAACACcg                                                                                                              >  2:707966/1‑140 (MQ=255)
                        gagaCGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACcgcg                                                                                                            <  1:796092/140‑1 (MQ=255)
                        gagaCCCAGCGCGCTTTTTTCTACGTGAAAGCGGTTTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCCCGTTAGTATCAAGTGGTTGTGAGGACGCAAAAAAACCCCCTTGCCCGCCTGTAACACcgcg                                                                                                            <  2:737040/140‑1 (MQ=255)
                            cGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGa                                                                                                        >  2:135165/1‑140 (MQ=255)
                              cAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAAt                                                                                                      >  1:52445/1‑140 (MQ=255)
                              cAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAAt                                                                                                      >  1:585211/1‑140 (MQ=255)
                              cAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAAt                                                                                                      <  2:261886/140‑1 (MQ=255)
                              cAGCGCGCATTATTCTACGTAAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAAt                                                                                                      <  2:585211/140‑1 (MQ=255)
                                    ccATTATTCTACGTAAAAACGGATAAAACGGCAGTTTAACCCGATTATCTTTGTAAATATTTTTTCCCGTTAGTATCAAGTGGTTGTGAGGCCGCAAAAAAACCCCCTTCCCCGCCTGTAACACCGCGCTGAATCAGGGt                                                                                                <  2:399842/139‑1 (MQ=255)
                                                    aaaCGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTcgc                                                                                >  1:176792/1‑140 (MQ=255)
                                                                        tAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTa                                                            <  1:267787/140‑1 (MQ=255)
                                                                                   tCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGc                                                 <  1:9452/140‑1 (MQ=255)
                                                                                   tCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGc                                                 <  1:19007/140‑1 (MQ=255)
                                                                                         tgAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTcc                                           >  2:332071/1‑140 (MQ=255)
                                                                                         tgAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTcc                                           <  2:176792/140‑1 (MQ=255)
                                                                                               ttttttCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGc                                     >  1:443083/1‑140 (MQ=255)
                                                                                                                 aggTGGCTGTGAGGGCGCAAAAAAACCCCCTGGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCttt                   <  1:707966/138‑1 (MQ=255)
                                                                                                                      gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc              >  2:4842/1‑140 (MQ=255)
                                                                                                                      gCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAAcc              >  1:43061/1‑140 (MQ=255)
                                                                                                                       ttgggaggagGCGAAAAAACCCCCTGGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACca             <  2:443083/136‑1 (MQ=255)
                                                                                                                       cTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACca             >  2:116886/1‑140 (MQ=255)
                                                                                                                       cTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACca             >  2:263645/1‑140 (MQ=255)
                                                                                                                         gtgAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTccacca                         >  2:553053/1‑126 (MQ=255)
                                                                                                                         gtgAGGACGCGAAAAAACCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTccacca                         <  1:553053/126‑1 (MQ=255)
                                                                                                                           gAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACcagcag         >  1:442180/1‑140 (MQ=255)
                                                                                                                            gggACGCAAAAAACCCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCagcagc        <  2:321937/139‑1 (MQ=255)
                                                                                                                                cgcgAAAAACCCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGt    <  2:346854/140‑1 (MQ=255)
                                                                                                                                cgcgAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGt    >  2:567303/1‑140 (MQ=255)
                                                                                                                                  cgAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTtc  >  2:68534/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                   
GCGAGCAGACAAGCACACTCCCATGAGACGCAGCGCGCATTATTCTACGTGAAAACGGATTAAACGGCAGGTTAAACCGAGTATCTTTGTGAATATTTTTTCACGTTAGTATCAAGTGGCTGTGAGGACGCGAAAAAATCCCCTTGCCCGCCTGTAACACCGCGCTGAATCAGGGTCTGCCACTCGCTTCGCGAACGCACGCCGGTGGCGTAAACCTGGGTGCTGGTCCCGGAGCAGGCTTCCACCAGGCTTTGAACCAGCAGCTGGTTC  >  NC_000913/3401321‑3401590

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑