Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,637,2350TG21.4% 51.1 / 12.9 28K41T (AAG→ACG) derGTPase; multicopy suppressor of ftsJ
Reads supporting (aligned to +/- strand):  ref base T (11/11);  new base G (3/3);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TAAACAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAAT  >  NC_000913/2637095‑2637380
                                                                                                                                            |                                                                                                                                                 
tAAACAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCa                                                                                                                                           >  1:160229/1‑149 (MQ=255)
  aaCAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCAc                                                                                                                                          <  2:160229/148‑1 (MQ=255)
     aGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGt                                                                                                                                      <  4:61789/149‑1 (MQ=255)
     aGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGt                                                                                                                                      <  4:59549/149‑1 (MQ=255)
                                              cGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGa                                                                                              <  8:147656/148‑1 (MQ=255)
                                                    gCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTGCGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTa                                                                                          >  3:37119/1‑146 (MQ=255)
                                                    gCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTCCGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGTGTGCGAGTTAGa                                                                                        <  4:37119/148‑1 (MQ=255)
                                                            cTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTAAGTACGGTcgtg                                                                                                                                         <  6:43823/91‑4 (MQ=255)
                                                            cTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTAAGTACGGTcgtg                                                                                                                                         >  5:43823/1‑88 (MQ=255)
                                                               cacCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATcc                                                                                                                      <  4:61596/107‑1 (MQ=255)
                                                               cacCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTAAGGTCACGAGTCAGACCCGGGAAATcc                                                                                                                      >  3:61596/1‑107 (MQ=255)
                                                                                  aaTCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGt                                                            <  3:55240/146‑1 (MQ=255)
                                                                                  aaTCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGt                                                            >  4:55240/1‑146 (MQ=255)
                                                                                              aTCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGt                                                                                    >  4:147360/1‑110 (MQ=255)
                                                                                                                aCGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTAcc                            >  6:124180/1‑148 (MQ=255)
                                                                                                                      ttCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGctt                      >  2:6439/1‑148 (MQ=255)
                                                                                                                       tCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGtt                                                                                           >  6:164320/1‑78 (MQ=255)
                                                                                                                       tCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGtt                                                                                           <  5:164320/78‑1 (MQ=255)
                                                                                                                       tCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTa                     >  2:120895/1‑148 (MQ=255)
                                                                                                                          aTTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGc                 <  5:124180/149‑1 (MQ=255)
                                                                                                                                tcgtgACCGTAAGTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCAt                          <  1:89051/130‑1 (MQ=255)
                                                                                                                                tcgtgACCGTAAGTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCAt                          >  2:89051/5‑134 (MQ=255)
                                                                                                                                tcgtgACCGTAAGTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGt                        >  7:79631/5‑136 (MQ=255)
                                                                                                                                tcgtgACCGTAAGTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGt                        <  8:79631/132‑1 (MQ=255)
                                                                                                                                gCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGt                               <  1:5608/129‑1 (MQ=255)
                                                                                                                                gCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAAACCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGt                               >  2:5608/1‑129 (MQ=255)
                                                                                                                                    gACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAAt                                                                               >  1:153925/1‑77 (MQ=255)
                                                                                                                                    gACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAAt                                                                               <  2:153925/77‑1 (MQ=255)
                                                                                                                                         tatttACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAat  >  8:138497/4‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
TAAACAGTACGACGTCCGCTTCTTCAATCGCCAGCAGCGACTGTTCCGCCATGCGGGTTTCTACACCGTCTTCTGTGCCATCAATCCCGCCGGTATCAATACAGATAAACTCACGGCCTTCAATTTCCGCACGACCGTACTTACGGTCACGAGTCAGACCCGGGAAATCCGCAACCAGCGCATCTCGGGTGCGAGTTAGACGGTTAAATAACGTGGATTTTCCTACGTTAGGGCGCCCGACAAGCGCGACCACAGGTACCATGTTTAAAGCCTCATTTTTATAAAT  >  NC_000913/2637095‑2637380

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