Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913380,0130G.23.1% 65.2 / 17.2 27intergenic (‑132/+56)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (20/0);  new base . (4/2);  total (24/3)
Fisher's exact test for biased strand distribution p-value = 4.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGT  >  NC_000913/379868‑380152
                                                                                                                                                 |                                                                                                                                           
ggAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggg                                                                                                                                          <  1:144433/149‑1 (MQ=255)
ggAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggg                                                                                                                                          <  4:29502/149‑1 (MQ=255)
 gAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggg                                                                                                                                         >  3:240376/1‑149 (MQ=255)
            aTTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTaaac                                                                                                                              >  8:320053/1‑147 (MQ=255)
            aTTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGCTATTCTATATGTCAATGCATACGCCCCTATAGTATATAACCCTGGAGAATAAAATTTATCGGGTGAATGTGGTCGGAAAACAATCAGGAAA‑GGGGGGGGGCTAATCg                                                                                                                             >  4:354063/1‑149 (MQ=255)
            aTTTCGAACCTAATAATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCAACTATAGTATATAACACGGGTGAATAAAATTTATACGGTGAATTTGGTCGGAAAACAAATAGGAAA‑GGGGGGGGGCTAATCg                                                                                                                             >  4:348397/1‑149 (MQ=255)
                  cACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGCTAATCGGCAgg                                                                                                                        >  1:317007/1‑148 (MQ=255)
                    ccTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCTGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGTGGGCTAATCGGCAggg                                                                                                                       >  1:85376/1‑148 (MQ=255)
                     cTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGGGGAAAGGGGGGGGGGGTAAACGGCCggg                                                                                                                       >  3:146727/1‑147 (MQ=255)
                                                    ggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggcta                                                                                                                                  >  6:190554/1‑102 (MQ=255)
                                                    ggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggcta                                                                                                                                  <  5:190554/105‑4 (MQ=255)
                                                       aTTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGACGGCCGCCCCGGATAGCGGGCGGCAGAAGGAAt                                                                                     >  3:213470/1‑147 (MQ=255)
                                                         tCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGTAAAACAAAGAGGAAAGGGGGGGGGGGTATTCGGAAGGGACGGCCGCCCCGGATTGCGGGCGGAAGAAGGAAGCAg                                                                                  >  1:244342/1‑148 (MQ=255)
                                                           tataTGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATTCGGGGCGGCAGAAGGAAGCAGAAt                                                                               >  3:347120/1‑149 (MQ=255)
                                                                gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                                                                           >  5:322938/1‑148 (MQ=255)
                                                                 tCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTAGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGCAATCAGAATTTCCAg                                                                         >  5:109672/1‑149 (MQ=255)
                                                                            cccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGa                                                                                       >  7:10911/1‑124 (MQ=255)
                                                                            cccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGAGAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCCGAAGGAATCAGAATTTCAAGGTCAGACgggg                                                              >  2:91873/1‑148 (MQ=255)
                                                                                 tataGTATATAACACTGGAGAATAAATTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                                        >  6:289249/1‑149 (MQ=255)
                                                                                 tataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                                           >  7:183187/1‑147 (MQ=255)
                                                                                     gTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGCATCAGAATTTCCCGGTCAGACGGGCTGCAAGTTg                                                     >  3:277931/1‑149 (MQ=255)
                                                                                         ataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAAGCAg                                                                                  >  7:92372/1‑116 (MQ=255)
                                                                                             cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGt                                             >  5:302132/1‑149 (MQ=255)
                                                                                             cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAAACGCCAGGGAAGGCCGCCCCGCATAGCGGTCCGCAGAAGGAATCAAAATTACCAGTTCAAACGAGCTGAAAGTTCCAGACCGt                                             >  6:174562/1‑149 (MQ=255)
                                                                                                tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGCCCGTTa                                           >  3:184947/1‑148 (MQ=255)
                                                                                                  gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTaaa                                         <  2:317007/147‑1 (MQ=255)
                                                                                                              tctCCGGGGAATGGGGTCGGAAAACAAAGAGG‑AACGGGGGGGGGCAAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg                            <  3:354063/146‑1 (MQ=255)
                                                                                                              tATCCGGGGAACGTGGTCGGAAAACAAATAGGAAA‑GGGGGGCGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTggg                            <  3:348397/148‑1 (MQ=255)
                                                                                                                                 gAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTAAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTAc         >  5:76988/1‑149 (MQ=255)
                                                                                                                                      cAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGGGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGcc    >  6:135975/1‑149 (MQ=255)
                                                                                                                                        aaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGt  >  3:59445/1‑149 (MQ=255)
                                                                                                                                          gagGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGCCCGTTAAAATCATCGGTTGGGGTGTCGTCCCCCCCTTTACCTGCCGt  >  2:94294/1‑147 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
GGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGT  >  NC_000913/379868‑380152

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: ‑