Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi380,0130G.27.8% 33.9 / 14.2 18intergenic (‑132/+56)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (13/0);  new base . (4/1);  total (17/1)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTG  >  NC_000913_3_bme_pgi/379866‑380148
                                                                                                                                                   |                                                                                                                                       
ccGGAGTACAGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggg                                                                                                                                       <  2:261942/150‑1 (MQ=255)
  ggAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTTGTCGGAAAACAAAGAGGAAAggggg                                                                                                                                     >  1:385912/1‑150 (MQ=255)
              aTTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggg                                                                                                                                    >  2:688968/1‑139 (MQ=255)
              aTTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggg                                                                                                                                    <  1:688968/139‑1 (MQ=255)
                            catcTGCATGCAATATACTATAGGGGGGTTTTCTCTATGTCAATGCCTACCCCCCTATAGTATATACCACTGGAGAATATAATTTTTCCGGTGAATGTGGGCGGAAAACATAGAGGAAAGGGGGGGGGGCTAATCGGCGGGGAAGgccgc                                                                                                           >  2:438556/1‑150 (MQ=255)
                                tgcatgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCgg                                                                                                      >  1:700558/1‑150 (MQ=255)
                                   atgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGATAAGCGGCAGGGAAGGCCGCCCCGGAt                                                                                                    >  2:762072/1‑150 (MQ=255)
                                     gcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAg                                                                                                  >  1:782672/1‑150 (MQ=255)
                                           taatataTGGGGGTATTGAATATGTCCATGCATACCCCCGAATAGTATATAACACTGGAGAATAATAATTATCCGCGGAATGTGGGCGGAAAACAAACAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcgg                                                                                           <  2:867788/147‑1 (MQ=255)
                                                                   tCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGCCAGGGAAGCCCGCCCCGGATACCGGGCGGCAGAAGGAATCGGAATTTCCAgg                                                                    >  2:889030/1‑150 (MQ=255)
                                                                             ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCGTATAGGCAGGGAAGGGCGCCCCGGATAGGGGGGGGCCGAAGGAATGAGAATTTGCAGGTCAGACggggg                                                         >  1:50089/1‑148 (MQ=255)
                                                                                   tataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGGGGAAGAAGGAATAAGAATTTCCGGGTCAGACGGCCTCCaat                                                    >  1:305892/1‑149 (MQ=255)
                                                                                   tataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                                    >  1:487480/1‑150 (MQ=255)
                                                                                   tataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCCGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCCGGGCAGCCGGGCAGCAAg                                                    >  1:906965/1‑150 (MQ=255)
                                                                                        tatataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGTGATAAACGGAAGTGAAGGCCTCCCCGGATAGAGGGCGGAAGAAGGAATAAGAATTTCAAGGACAGACGGGCTGAAAGTTgaa                                               >  2:599618/1‑148 (MQ=255)
                                                                                                 cTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCTGGGAAGGCCGCCCCGGATAGCGGGCGGCAGCAGGAATCCGAATTTCCAGGACAGACGGGCTGCAAGTTGCAGCCCGGTaa                                      >  2:775262/1‑150 (MQ=255)
                                                                                                       aaTAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGAAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAATTTGCAGACCGTTAAAatcat                                >  2:139168/1‑150 (MQ=255)
                                                                                                         tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCTGGGGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCg                              >  1:214312/1‑150 (MQ=255)
                                                                                                                      gTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGCCCGTTAAAATCATCGGTTGGGGTGTCGTa                >  1:407446/1‑150 (MQ=255)
                                                                                                                              ggTAGGAAAACAAAGAGGAAAGAGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGCTAGCGGGCGGCAGAAGCAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCTTTAAAATCATCGGTTGGGGTGTCGTACCACACt         >  2:778288/1‑150 (MQ=255)
                                                                                                                                 cGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGATTCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTa      >  1:536351/1‑150 (MQ=255)
                                                                                                                                    aaaaCAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTg  >  1:728873/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                       
CCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTG  >  NC_000913_3_bme_pgi/379866‑380148

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: