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

CCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAG  >  NC_000913/379866‑380161
                                                                                                                                                   |                                                                                                                                                    
ccGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAg                                                                                                                                                      <  1:251899/148‑1 (MQ=255)
 cGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggg                                                                                                                                                    <  4:55021/149‑1 (MQ=255)
 cGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgg                                                                                                                                                     <  8:186498/148‑1 (MQ=255)
  ggAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgggg                                                                                                                                                   <  8:205836/149‑1 (MQ=255)
   gAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggg                                                                                                                                                  <  8:351958/149‑1 (MQ=255)
         tGGGCCTTTCGCACCTCATCCTCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATACAATTTATCCGGGGAATGTGGGAGGAAAACAAAGAGGAAAGGGGGGGGGGc                                                                                                                                            >  1:195380/1‑149 (MQ=255)
               tttCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgg                                                                                                                                                     >  4:28662/1‑134 (MQ=255)
               tttCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAgg                                                                                                                                                     <  3:28662/134‑1 (MQ=255)
                     aCCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAgg                                                                                                                                 >  7:81837/1‑148 (MQ=255)
                       cTCATCATCTGCATGCAATATACTTTAGGGGGTTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAggg                                                                                                                                >  2:282375/1‑147 (MQ=255)
                       cTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggg                                                                                                                                                >  8:36811/1‑131 (MQ=255)
                       cTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggg                                                                                                                                                <  7:36811/131‑1 (MQ=255)
                                   atgcaATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGTAATGGGCAGGGATGGCCGCCCCGGa                                                                                                                  >  2:158119/1‑149 (MQ=255)
                                                   gggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATCGCGGGGGGCAGAAg                                                                                                  >  6:252606/1‑149 (MQ=255)
                                                                  gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAAACGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAAGCAGATTTTccc                                                                                   >  2:189221/1‑148 (MQ=255)
                                                                  gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAg                                                                                  >  1:315621/1‑149 (MQ=255)
                                                                     aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGATAATCGGCAGGGAAGGCCGCCCCGGAATGCGGGCGGCAGAAGGAATCATAATTTCCAgg                                                                                 >  3:278744/1‑148 (MQ=255)
                                                                      aTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCTGAAGGAATCAGAATTTCCAGGt                                                                                >  1:156599/1‑148 (MQ=255)
                                                                          aTACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGCCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTACCAGGTCTGAc                                                                           >  3:109771/1‑149 (MQ=255)
                                                                          aTACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCAAATCGGCAGGGAAGCCCCCCCCGGATAGCGGGCGGCAGAAGGAATCTCAATTTCCAGGTCAgcc                                                                           >  5:122406/1‑147 (MQ=255)
                                                                            aCCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACATAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCAAGGTCAGACgg                                                                         >  1:365057/1‑149 (MQ=255)
                                                                              cccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGc                                                                       >  3:213329/1‑149 (MQ=255)
                                                                                         atataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                                            <  2:315621/148‑1 (MQ=255)
                                                                                          taaaaCACCGGACCATAAAATTTAACCCGCGCATGGGGGCGCAAAACCAACCGGCACCGGGCGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCa                                                            <  6:122406/145‑1 (MQ=255)
                                                                                            taACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGATAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGAAGAAGGAATCAGAATTTCCAGGTCAGACGGGATGCAAGTTGCAGAc                                                         >  5:32639/1‑149 (MQ=255)
                                                                                               cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGCCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAAAAAGAATTTCCAGGTAAGACGGGCTGCAAGTTGCAgcacg                                                       >  5:25064/1‑144 (MQ=255)
                                                                                                    gagaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTaaaa                                                 >  1:120823/1‑148 (MQ=255)
                                                                                                       aaTAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCCGGGCACACGGGCTGCAAGTTGCAGACCGTTAAAatca                                              >  7:309029/1‑149 (MQ=255)
                                                                                                         tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGGTATTCGGCAGGGAAGGCCGCCCGGGATAGCGGGCGGAAGAAGGAATCAGAATTACCAGGTCAGACGGGCTGCAAGTTCCAGACCGTTAAAaacatc                                            >  6:112456/1‑149 (MQ=255)
                                                                                                                 aTCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGTCAGGGAAGTCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAg                                                                                  >  1:371527/1‑103 (MQ=255)
                                                                                                                      gtgcttgtgGGCGTAAAACAAAGAGGCAA‑GGGGGGGGTCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGt                              <  2:120823/144‑1 (MQ=255)
                                                                                                                       tGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTa                             >  3:267228/1‑149 (MQ=255)
                                                                                                                                         aaaGAGGAAAGGGGGGGGGGGTAATCGGCAGGGAAGGCCGCCACGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCg            >  3:337357/1‑149 (MQ=255)
                                                                                                                                                   ggggggggggCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCgag  <  2:365057/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAG  >  NC_000913/379866‑380161

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