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

TATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCT  >  NC_000913/379912‑380147
                                                                                                       |                                                                                                                                      
tataGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGAt                                                                                                   >  1:250177/1‑139 (MQ=255)
  taGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGGGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGCTAATCGGCAGGGGAGGCCGCCCCGGATAgc                                                                                                >  1:440531/1‑139 (MQ=255)
     gggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgg                                                                                              >  1:265744/1‑139 (MQ=255)
       gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTCTCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                              <  1:538934/121‑1 (MQ=255)
       gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                              >  2:538934/1‑121 (MQ=255)
          tATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGggcggc                                                                                         >  2:229434/1‑139 (MQ=255)
                       aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAAtaa                                                                              >  1:283356/1‑137 (MQ=255)
                           cATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAACCAAAGAGGAAA‑‑‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTc                                                                       <  2:511699/139‑1 (MQ=255)
                           cATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTc                                                                       >  1:511699/1‑139 (MQ=255)
                              aCCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                                                                      >  2:443876/1‑139 (MQ=255)
                               ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGCTAATCGGCAGGGAAGgcc                                                                                                            >  1:160101/1‑98 (MQ=255)
                               ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGCTAATCGGCAGGGAAGgcc                                                                                                            <  2:160101/98‑1 (MQ=255)
                                          tatataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGc                                                         >  1:215651/1‑139 (MQ=255)
                                            tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑GGGGGGGGGGCTAATAGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGACGGAATCAGAATTTCCAGGTCAGACggggt                                                        >  2:88129/1‑137 (MQ=255)
                                            tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGGGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTg                                                       >  1:534422/1‑139 (MQ=255)
                                              taACACTGGAGAATAAAATTTATCCGGTGAAAGTGGTCGGAAAACAAAGAGGAAAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                                      <  1:443876/139‑1 (MQ=255)
                                               aatgatGGAGAACTAAAATTAGCCGGTGACAGTGGCCGGAAAACGAAGAGGAAAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCGCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCa                                                     <  1:88129/134‑1 (MQ=255)
                                                 cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGt                                                  >  1:79478/1‑139 (MQ=255)
                                                 cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTAAGACGGGCTGCAAGt                                                  >  1:461864/1‑139 (MQ=255)
                                                 cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGtt                                                 >  2:443011/1‑139 (MQ=255)
                                                    tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTa                                                                                                                            <  1:77547/62‑1 (MQ=255)
                                                    tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTa                                                                                                                            >  2:77547/1‑62 (MQ=255)
                                                    tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                              >  2:542971/1‑139 (MQ=255)
                                                    tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                              >  2:362074/1‑139 (MQ=255)
                                                      gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGCTAATCGGCAGGAAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                            >  1:445158/1‑139 (MQ=255)
                                                                          gAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                                      <  1:653750/110‑1 (MQ=255)
                                                                          gAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                                      >  2:653750/1‑110 (MQ=255)
                                                                            aTGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgggg                      <  1:443011/139‑1 (MQ=255)
                                                                             tgtgGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgggg                      <  2:461864/139‑1 (MQ=255)
                                                                             tgtgGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgggg                      <  2:79478/139‑1 (MQ=255)
                                                                                    ggAAAACAAAGAGGAAA‑‑‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc              <  1:362074/139‑1 (MQ=255)
                                                                                       aaaCAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACca            <  2:265744/139‑1 (MQ=255)
                                                                                               agGAAAG‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTa     >  1:48886/1‑139 (MQ=255)
                                                                                                ggAAA‑‑‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCt  >  2:317122/1‑139 (MQ=255)
                                                                                                       |                                                                                                                                      
TATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCT  >  NC_000913/379912‑380147

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