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

GAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCA  >  NC_000913/379869‑380154
                                                                                                                                                |                                                                                                                                             
gAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggg                                                                                                                                           >  1:366106/1‑149 (MQ=255)
     cTGGGCATTTCCCACCTCATGATCAGCATGCAATATACTATAGGGGGGTATTCTATTTGTCAATGCATACCCCCCTATATTATAAAAAACTGGATAATAAAATTTATCCGGTTAATGTGGTCGGAAAACAAACAGGAAAtggggggggg                                                                                                                                      <  6:152149/149‑1 (MQ=255)
      tGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggg                                                                                                                                      >  3:290159/1‑148 (MQ=255)
       gggCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggg                                                                                                                                            >  4:14543/1‑141 (MQ=255)
       gggCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggg                                                                                                                                            <  3:14543/141‑1 (MQ=255)
                      catcatCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAgg                                                                                                                     >  4:224561/1‑149 (MQ=255)
                      catcatCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGATAATCGCCAGGGAAgg                                                                                                                     >  2:364215/1‑149 (MQ=255)
                         catcTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAATATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGAAGGGACGgccg                                                                                                                  >  3:315371/1‑149 (MQ=255)
                              gcatgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAAGCGGCAGGGAAGGCCGCCCCg                                                                                                             >  8:295431/1‑149 (MQ=255)
                                            ataGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGATGGCCGCCCCGGATAgcgggcgg                                                                                                 >  3:155348/1‑147 (MQ=255)
                                              aGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGAAAGCGGGCGGCAg                                                                                              >  1:54369/1‑148 (MQ=255)
                                                          tttatGTCAATGAATACCCCCCTATACTATATAACACAGGCCAATAAAATTTACCCGGGGAATGTGGGCGGAAAACAACCCGGCAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCa                                                                                     <  8:299896/143‑1 (MQ=255)
                                                          tataTGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGACGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCa                                                                                     >  7:299896/1‑145 (MQ=255)
                                                                       aTACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAAAAAGAATTTCAAGGTCAGAc                                                                    >  2:141276/1‑149 (MQ=255)
                                                                           cccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGTAATCGGCCGGGAGGGCCGCCCCGGAAGGCGGGCGGCAGACGGAATCCAAATTTCCAGGTCAGACGGGc                                                                >  3:197518/1‑149 (MQ=255)
                                                                           cccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACggga                                                                >  4:36578/1‑148 (MQ=255)
                                                                                 ataGTTTATAACACTGTATATTAAAATTTATCCGGTTATTGTGGTCGCTAAACAAAGAGGAAA‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                                    <  4:317886/106‑1 (MQ=255)
                                                                                 ataGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGCTAAACGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                                    >  3:317886/1‑106 (MQ=255)
                                                                                     tatataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGATAATCGGCAGGGAAGGCCGCCCCGGAGAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg                                                       >  5:74289/1‑148 (MQ=255)
                                                                                         taACAATGGAGAATTAAATGTCTCCGGTGACTGGGGTCCGTAAACCAAGCGTAACCGGGGGTGGGCTAATCGGCAGGGAAGGCCGGCCCGGATAGCGTGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACTGGCTGCAAGTTGCAGa                                                   <  2:132613/148‑1 (MQ=255)
                                                                                          aaCACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCg                                                >  7:258298/1‑149 (MQ=255)
                                                                                             acTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTAAGACGGGCTGCAAGTTGCAGACCGt                                               >  6:382497/1‑148 (MQ=255)
                                                                                                  agaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCTGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAat                                         >  4:292556/1‑149 (MQ=255)
                                                                                                      tAAATTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATACCGGGCGGCAGAAGGAATCAGAATGTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatcatc                                     >  4:292741/1‑149 (MQ=255)
                                                                                                        aaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatcat                                      >  1:73031/1‑146 (MQ=255)
                                                                                                                 cGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCTGAAGGAATCAGAATTTCCAGGACAGACGGGCTGCAAGTTGCAGACCGTTAAAATAATCGGTTGGGgtgg                          >  3:100288/1‑148 (MQ=255)
                                                                                                                        tgtgGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAAACGGCAGGGAAGCCCGCCCCGGATACCGGGCGGCCGAGGGAATCAGAATTTCCAGGTAACACGGGCTGCAAGTGGCGGCCCGTTAAATTCTTCGGTTGGGGTGGCGTCCca                   >  5:315918/1‑149 (MQ=255)
                                                                                                                                   aaCAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTATAATCATCGGTTGGGGTGTCGTACCACACTTTACCTg        >  7:139084/1‑149 (MQ=255)
                                                                                                                                     cAAAGAGGAAA‑GGGGGGGGGCGAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCACAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTAACTGCCg     >  7:130016/1‑149 (MQ=255)
                                                                                                                                      aaaGAGGAAA‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGt    >  8:75474/1‑149 (MQ=255)
                                                                                                                                        agagGAAA‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCa  <  8:258298/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
GAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCA  >  NC_000913/379869‑380154

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