Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,356,7340GC24.0% 48.2 / 12.2 25intergenic (+23/‑170)glpC/yfaDanaerobic sn‑glycerol‑3‑phosphate dehydrogenase, C subunit, 4Fe‑4S iron‑sulfur cluster/transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base G (11/8);  new base C (3/3);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATAGAAGAAAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGC  >  NC_000913/2356589‑2356873
                                                                                                                                                 |                                                                                                                                           
aTAGAAGAAAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAt                                                                                                                                          >  5:119446/1‑149 (MQ=255)
  agaagaAAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCacat                                                                                                                                        >  7:96153/1‑146 (MQ=255)
   gaagaaAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtgtg                                                                                                                                       <  6:119446/149‑1 (MQ=255)
        aaGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTAc                                                                                                                                  <  8:96153/149‑1 (MQ=255)
               gcAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCAt                                                                                                                           >  5:64358/1‑149 (MQ=255)
                                             tGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtgt                                                                                                                                        <  8:54521/106‑1 (MQ=255)
                                             tGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtgt                                                                                                                                        >  7:54521/1‑106 (MQ=255)
                                              gTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTAt                                                                                            >  1:8579/1‑149 (MQ=255)
                                                                        aGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCctc                                                                  <  2:8579/149‑1 (MQ=255)
                                                                                                            cAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGt                              >  2:87454/1‑149 (MQ=255)
                                                                                                                cgcTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCaaa                          <  6:64358/149‑1 (MQ=255)
                                                                                                                     gCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCt                     >  4:113294/1‑149 (MQ=255)
                                                                                                                       ttAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTgg                   >  6:105677/1‑149 (MQ=255)
                                                                                                                         aaaCTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGa                 >  1:110956/1‑149 (MQ=255)
                                                                                                                               cTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAAc             <  1:87454/147‑1 (MQ=255)
                                                                                                                                  tCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTa         <  5:105677/148‑1 (MQ=255)
                                                                                                                                  aaaaaTGCCCGGTAACCACATGCTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCATaa        >  1:144254/5‑149 (MQ=255)
                                                                                                                                    tGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGAGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATg      >  6:150447/1‑149 (MQ=255)
                                                                                                                                    aaaTGCCCGGTAACCACATGCTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTgg                   >  4:20650/3‑136 (MQ=255)
                                                                                                                                    aaaTGCCCGGTAACCACATGCTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTgg                   <  3:20650/134‑1 (MQ=255)
                                                                                                                                     gATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGGAACTgg                                     <  4:139877/117‑1 (MQ=255)
                                                                                                                                     gATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGGAACTgg                                     >  3:139877/1‑117 (MQ=255)
                                                                                                                                         cccGGTAACCACATGCTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCATAATGTcgc  <  2:144254/148‑1 (MQ=255)
                                                                                                                                            ggTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGa                 >  4:112030/1‑130 (MQ=255)
                                                                                                                                            ggTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGa                 <  3:112030/130‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
ATAGAAGAAAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGC  >  NC_000913/2356589‑2356873

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