Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913380,0130G.47.8% 92.0 / 95.4 46intergenic (‑132/+56)frmR/yaiODNA‑binding transcriptional repressor FrmR/outer membrane protein YaiO
Reads supporting (aligned to +/- strand):  ref base G (11/13);  new base . (10/12);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGG  >  NC_000913/379966‑380059
                                               |                                              
gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgg                                               <  1:4483455/49‑1 (MQ=255)
gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgg                                               <  1:5415636/49‑1 (MQ=255)
  gaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggg                                            <  1:5513855/50‑1 (MQ=255)
   aaTAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggg                                           >  1:871448/1‑50 (MQ=255)
   aaTAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggg                                           >  1:5352962/1‑50 (MQ=255)
    aTAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggg                                          >  1:962722/1‑50 (MQ=255)
    aTAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggg                                          <  1:3734162/50‑1 (MQ=255)
     tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggg                                         <  1:500277/50‑1 (MQ=255)
      aaaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggggg                                        >  1:3712840/1‑50 (MQ=255)
       aaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggggggggg                                       >  1:698100/1‑50 (MQ=255)
        aaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGc                                      <  1:1490714/50‑1 (MQ=255)
         aTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCt                                     <  1:533317/50‑1 (MQ=255)
         aTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCt                                     <  1:14463/50‑1 (MQ=255)
            tATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATc                                 >  1:4892649/1‑50 (MQ=255)
             aTCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAAt                                  <  1:4292362/49‑1 (MQ=255)
             aTCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATc                                 <  1:4077656/49‑1 (MQ=255)
              tCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATc                                 <  1:2161130/49‑1 (MQ=255)
              tCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCg                                <  1:1857350/50‑1 (MQ=255)
               ccGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCgg                               >  1:4946048/1‑50 (MQ=255)
                cGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCgg                               <  1:3312566/49‑1 (MQ=255)
                cGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGc                              <  1:4150983/49‑1 (MQ=255)
                 ggTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCgg                               >  1:1623018/1‑48 (MQ=255)
                 ggTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCa                             <  1:896415/50‑1 (MQ=255)
                 ggTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAg                            >  1:3800213/1‑50 (MQ=255)
                  gTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAg                            >  1:1745190/1‑50 (MQ=255)
                    gAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAgg                           >  1:4900305/1‑49 (MQ=255)
                     aaTGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAggg                          >  1:1575104/1‑49 (MQ=255)
                     aaTGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAggg                          >  1:1286540/1‑48 (MQ=255)
                      aTGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAg                       <  1:1704439/50‑1 (MQ=255)
                       tgtgGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAgg                      <  1:568547/50‑1 (MQ=255)
                        gtgGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAgg                      >  1:3415754/1‑50 (MQ=255)
                         tgGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGgc                     >  1:327706/1‑50 (MQ=255)
                         tgGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGgc                     >  1:3174778/1‑50 (MQ=255)
                         tgGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgcc                    <  1:5457140/50‑1 (MQ=255)
                         tgGTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgcc                    <  1:5201543/50‑1 (MQ=255)
                         tgGTCGGAAAACAAAGAGGAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGgcc                    <  1:1904458/49‑1 (MQ=255)
                          ggTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgccg                   >  1:3398096/1‑50 (MQ=255)
                          ggTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgccg                   >  1:829873/1‑50 (MQ=255)
                           gTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgccg                   <  1:5180360/49‑1 (MQ=255)
                           gTCGGAAAACAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGgccg                   <  1:427568/49‑1 (MQ=255)
                            tCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGgccgc                  <  1:4440139/50‑1 (MQ=255)
                              gaaaaaCAAAGAGAAAA‑‑GGGGGGGGCTAATCGGCAGGGAAGGCCGcccc               <  1:4847461/47‑1 (MQ=255)
                                 aaaCAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCgg             <  1:5270612/50‑1 (MQ=255)
                                 aaaCAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGa            >  1:4794099/1‑50 (MQ=255)
                                    cAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATa          <  1:766579/50‑1 (MQ=255)
                                    cAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAg         <  1:3763205/50‑1 (MQ=255)
                                    cAAAGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAg         <  1:1527215/50‑1 (MQ=255)
                                     aaaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAg         <  1:1868616/50‑1 (MQ=255)
                                     aaaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAg         >  1:4660631/1‑50 (MQ=255)
                                     aaaGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgc        >  1:133488/1‑50 (MQ=255)
                                     aaaGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgc        >  1:3753124/1‑50 (MQ=255)
                                     aaaGAGGAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgc        >  1:4135206/1‑50 (MQ=255)
                                          ggAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggc    <  1:4184484/50‑1 (MQ=255)
                                          ggAAAGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcg   >  1:2038432/1‑50 (MQ=255)
                                             aaGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcgg  >  1:3030694/1‑49 (MQ=255)
                                              aGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggcgg  >  1:3265069/1‑48 (MQ=255)
                                               |                                              
GAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGG  >  NC_000913/379966‑380059

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: