Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,476,2170CA73.3% 45.4 / 23.8 30G75V (GGT→GTT) yheOputative DNA‑binding transcriptional regulator YheO
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base A (6/16);  total (10/20)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATT  >  NC_000913/3476142‑3476297
                                                                           |                                                                                
tAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATc                                                                             <  2:610062/81‑1 (MQ=255)
  gCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCtt                                                                           <  2:610064/81‑1 (MQ=255)
  gCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCtt                                                                           <  2:610063/81‑1 (MQ=255)
     tttAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCg                                                                        >  1:610098/1‑81 (MQ=255)
     tttAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCg                                                                        <  2:610065/81‑1 (MQ=255)
       tAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCgg                                                                       <  2:610066/80‑1 (MQ=255)
           aaCGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGt                                                                  <  1:610067/81‑1 (MQ=255)
               ctgctATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGt                                                              >  1:610111/1‑81 (MQ=255)
               ctgctATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGt                                                              >  1:610120/1‑81 (MQ=255)
                  ctATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCa                                                           >  2:610099/1‑81 (MQ=255)
                     tCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCAcc                                                         <  2:610068/80‑1 (MQ=255)
                             cGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAAt                                                >  2:610123/1‑81 (MQ=255)
                                   atCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAAt                                          >  2:610124/1‑81 (MQ=255)
                                       tGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCt                                      <  1:610069/81‑1 (MQ=255)
                                       tGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCt                                      <  1:610070/81‑1 (MQ=255)
                                        gcagcaTACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCt                                      >  1:610113/1‑80 (MQ=255)
                                         cagcaTACGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTg                                     <  2:610071/80‑1 (MQ=255)
                                           gcaTACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAAc                                  <  1:610072/81‑1 (MQ=255)
                                                cGTAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACAtttt                             <  2:610073/81‑1 (MQ=255)
                                                  tAGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAg                           <  2:610074/81‑1 (MQ=255)
                                                   aGCGCCAGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACAttt                              >  1:610114/1‑77 (MQ=255)
                                                     cgcCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATc                        <  2:610075/81‑1 (MQ=255)
                                                         aGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGc                    <  1:610076/81‑1 (MQ=255)
                                                         aGGTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGc                    <  2:610077/81‑1 (MQ=255)
                                                          ggTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCa                   <  2:610078/81‑1 (MQ=255)
                                                           gTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCaa                  <  1:610079/81‑1 (MQ=255)
                                                           gTCAGTAATTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCa                   >  2:610117/1‑80 (MQ=255)
                                                             cAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCa                   >  1:610122/1‑78 (MQ=255)
                                                                  aTTGGCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGc           <  1:610080/81‑1 (MQ=255)
                                                                      gCGAAACAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCaaaa       <  2:610081/81‑1 (MQ=255)
                                                                           ccAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGAtt  >  2:610129/1‑81 (MQ=255)
                                                                           |                                                                                
TAGCATTTAGAAACGCTGCTATCCGCTCCCGTCATATCGTGCAGCATACGTAGCGCCAGGTCAGTAATTGGCGAACCAATCTTCCGGCCTGTATGTTCACCGTTAGCAATGCGAATGGCTGAACATTTTAGATCCTGCAAAGAGTGCAAAACGATT  >  NC_000913/3476142‑3476297

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: