Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913380,0130G.23.5% 79.9 / 25.5 36intergenic (‑132/+56)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (22/4);  new base . (5/3);  total (27/9)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCC  >  NC_000913/379874‑380158
                                                                                                                                            |                                                                                                                                                  
cTGGGCATTTCGCACCTCATCATATGAATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAAA‑ggggggtggg                                                                                                                                           >  5:387109/1‑149 (MQ=255)
 tGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGc                                                                                                                                          >  1:445197/1‑149 (MQ=255)
   ggCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑ggggggggtggc                                                                                                                                         >  7:349713/1‑144 (MQ=255)
               cTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTTTATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGCTAAT‑CGGCAGGGAAg                                                                                                                          >  3:264180/1‑149 (MQ=255)
               cTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGCGTAAG‑CGAAAGGGag                                                                                                                           >  6:134322/1‑148 (MQ=255)
                  atcatcTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGCTAAT‑CGGCAGGGAAGgcc                                                                                                                       >  7:409484/1‑149 (MQ=255)
                   tcatcTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGgcc                                                                                                                       >  7:370840/1‑149 (MQ=255)
                   tcatcTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGgcc                                                                                                                       >  2:333284/1‑149 (MQ=255)
                               aaTATACTATAGGGGGCTATTCTATCTCTCACTCCATACCCCCCTATAGTATATAACACCGGACAATAAAATCTCTCCCGGCAATGCGGCCGGAACACAACCCCGAAA‑GGGGGGGCGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAgc                                                                                                           <  2:333898/149‑1 (MQ=255)
                                   aaCTATAGGGGGGTATTATAAAAGTCAACGACTACCCCCCTATAGTAAATAACACTGGCCCAAAAAACTTACCCCGGGACCCCGGTCCGCAAACACACCCCCCC‑CGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                                        <  2:445197/147‑1 (MQ=255)
                                     cTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGATAAT‑AGGCAGGGAAGGCCGCCCCGGATAgcgggcgg                                                                                                     >  6:362932/1‑149 (MQ=255)
                                          ggggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAg                                                                                               <  8:409484/149‑1 (MQ=255)
                                                    ctcaatGTCAATCCCTCCCCCCCTATACTATATAACACTGGAGACTAAAATTTATCCGGGGCACGGGGCCCGACCCCCAACCCCACA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGa                                                                                       <  8:51841/144‑1 (MQ=255)
                                                             aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAC‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGCTGGCGGGCGGCAGAAGGAATCAGAATTTGCAggg                                                                             >  3:293333/1‑148 (MQ=255)
                                                             aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGTAAGCGGGGGGAAGAAGGATGCAGATTTTCCTGGt                                                                             >  1:7792/1‑149 (MQ=255)
                                                              aTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGt                                                                             >  7:76724/1‑148 (MQ=255)
                                                                    aCCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGAGAGCGGGGGGCAGAAGGATTCAGAATTTCCAGGGCAGACgg                                                                      >  1:302671/1‑149 (MQ=255)
                                                                      cccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGCTAAT‑CGGCAGGTAAGGCCGCCCCGAATAGCGGGCGGCAGAAGGAATCACAATTTCAAGGTAAGACggc                                                                     >  3:105731/1‑148 (MQ=255)
                                                                      cccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGAAAGCGGGCGGCAGAGGGAATCAGAATTTCCAGGTCAGCCggga                                                                    >  6:334691/1‑148 (MQ=255)
                                                                        cccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGGTAAT‑CGGCAGGGAAGGCCGCCCCGGTTAGCGGGCGGCAGAAGGAATCAGAATTTc                                                                                  >  3:126024/1‑133 (MQ=255)
                                                                             taGTATATAACACTGGAGAATAAAATTTATCCGGTGATTGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGATGGCCGCCCCGGATAGCGGGCGGCAGATGGATTCAGATTTTCCAGGTCAGACGGGCTGCAAGt                                                             >  3:108817/1‑149 (MQ=255)
                                                                                  taaaccccTGGACAATAAAATTTCACCCGGGAATGGGGTCGCAAAACACAGCGCAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGc                                                          <  8:76724/140‑1 (MQ=255)
                                                                                     aaCACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGAGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGAcc                                                     >  4:8372/1‑149 (MQ=255)
                                                                                          tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGGTAAG‑CGGCAGGGACGGCCGCACCGGATATCGGGCGGAAGATGGAATAAGAATTTCCCGGTAACACGGGATGAAAGTGGCAGCCCGTTa                                                 >  4:301733/1‑148 (MQ=255)
                                                                                            gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                                                                                 >  3:99933/1‑114 (MQ=255)
                                                                                            gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTaaaa                                              >  7:122072/1‑149 (MQ=255)
                                                                                            gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg                                                           >  6:94412/1‑134 (MQ=255)
                                                                                            gagaATAAAATTTATCCGGTGAATGGGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg                                                           <  5:94412/134‑1 (MQ=255)
                                                                                            gagaATAAAAATTATCCGGTGTATGTGGGCGGAAAACATAGAGG‑AT‑TGGGCGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTaaat                                              <  2:159498/148‑2 (MQ=255)
                                                                                            cagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGGAGAACGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCCTTaaaa                                              >  1:159498/2‑148 (MQ=255)
                                                                                                      tttATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGATTCAGAATTTCCAGGTCAGCCGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTg                                    >  3:284635/1‑149 (MQ=255)
                                                                                                         aTCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGCGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg                                  <  3:8372/148‑1 (MQ=255)
                                                                                                                   tgtgGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGATAAT‑CGGCCGGGAAGGCCGACCCGGATAGCGGCCGGCAGAAGGAATAAGAATTTCGAGGTCAGACGGGCTGCAAGTTGCCGACCGTTAAAATCATCGGTTGGGGTGTCGTGCca                       >  5:361836/1‑149 (MQ=255)
                                                                                                                       gTCGCAAAACAAAGAGGAAA‑‑GGGGGGGGGCAAATACGGAAGGGAAGCCCCCCCCGGAAACCGGGGACAAAAAGGAAAAAAAATTTCAGGGTCAGACGGCATGAAAGTTGCAGACCGTAAAAAGAAACGGTTGGGGAGGCGTACaaccat                   >  2:326469/1‑146 (MQ=255)
                                                                                                                          ggAAAACAAAGAGGAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTAAGACGGGCTGCAAGTTGCAGACCGTTAAAATCCTCGGTTGGGGTGTCGTACCACACTTTa                >  5:160958/1‑149 (MQ=255)
                                                                                                                                 aaaGAGGAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGGGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGCCCGTTAAAATCATCGGTTGGGGGGTCGTACCACCATTTACCTGCCg         >  2:319538/1‑149 (MQ=255)
                                                                                                                                      ggAAA‑‑GGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGcc   <  4:264180/149‑1 (MQ=255)
                                                                                                                                        aaa‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGCCCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGccc  >  7:86966/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
CTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGGCTAAT‑CGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCC  >  NC_000913/379874‑380158

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