Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,302,9280TG35.5% 47.0 / 13.4 31V147G (GTG→GGG) yhbVputative peptidase YhbV
Reads supporting (aligned to +/- strand):  ref base T (0/20);  new base G (11/0);  total (11/20)
Fisher's exact test for biased strand distribution p-value = 1.18e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AATATTGCTCAAACAGGGCATGATGCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTGCTTAATCAGTGCGATGAGCTGGGCATTCGTAACCAGTTTGAAGTGGAAGTCCTGAGC  >  NC_000913/3302863‑3302992
                                                                 |                                                                
aaTATTGCTCAAACAGGGCATGATGCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGa                                                                <  1:1409337/68‑1 (MQ=255)
    ttGCTCAAACAGGGCATGATGCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGGGAATCt                                                            >  1:1205481/1‑68 (MQ=255)
     tGCTCAAACAGGGCATGATGCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTg                                                           <  1:2867112/68‑1 (MQ=255)
     tGCTCAAACAGGGCATGATGCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTg                                                           <  1:1504250/68‑1 (MQ=255)
     tGCTCAAACAGGGCATGATGCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTg                                                           <  1:15177/68‑1 (MQ=255)
     tGCTCAAACAGGGCATGATGCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTg                                                           <  1:82384/68‑1 (MQ=255)
     tGCTCAAACAGGGCATGATGCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTg                                                           <  1:1836455/68‑1 (MQ=255)
     tGCTCAAACAGGGCATGATGCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTg                                                           <  1:2616056/68‑1 (MQ=255)
     tGCTCAAACAGGGCATGATGCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTg                                                           <  1:485232/68‑1 (MQ=255)
       cTCAAACAGGGCATGATGCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTGCt                                                         <  1:3339910/68‑1 (MQ=255)
            aCAGGGCATGATGCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTGCTTAATc                                                    <  1:2292133/68‑1 (MQ=255)
                ggCATANTGCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGCATGGGGCATCTGCTTAATCAggg                                                >  1:2030616/1‑66 (MQ=255)
                      atgCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGGCTGGCTGGTGAATCTGCTTAATCAGTGCGATGa                                          <  1:2453084/68‑1 (MQ=255)
                      atgCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGGGAATCTTCTTAATCAGTGCGATGa                                          >  1:764113/1‑68 (MQ=255)
                      atgCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGGGAATCTGCTTAATCAGTGCGGTGa                                          >  1:3003067/1‑68 (MQ=255)
                      atgCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGGGAATCTGCTTAATCAGTGCGATGa                                          >  1:812075/1‑68 (MQ=255)
                      atgCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGGGAAACTTCTTAATCAGTGCGCTgg                                          >  1:2643863/1‑67 (MQ=255)
                       tgCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTGCTTAATCAGTGCGATGAg                                         <  1:1129486/68‑1 (MQ=255)
                               tgtATGCCAGTGGAGCTTTCCCGCGACTGGCTGGGGAAATTGCCTAATCAGTGGGATGAGCGGGGCAt                                 >  1:3139837/1‑68 (MQ=255)
                                gtATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTGCTTAATCAGTGCGATGAGCTGGGCAtt                                <  1:2218/68‑1 (MQ=255)
                                gtATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTGCTTAATCAGTGCGATGAGCTGGGCAtt                                <  1:2755716/68‑1 (MQ=255)
                                gtATGCCAGTGGAGCTTTCCCGCGACTGGCTGGGGAATCTGCTTAATCAGTGCGATGAGCCGGGCAtt                                >  1:770298/1‑68 (MQ=255)
                                  catccAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTGCTTAATCAGTGCGATGAGCTGGTCATTCg                              <  1:2076478/65‑1 (MQ=255)
                                    gCCAGTGGAGCTTTCCCGCGACTGGCTGGGGGATCTTCTTAATTAGTGGGATGAGCTGGGCATTTGTa                            >  1:2971831/1‑68 (MQ=255)
                                    gCCAGTGGAGCTTTCCCGCGACTGGCTGGGGAATTTGCTTAATCAGTGCGATGAGGTGGGCGTTCGTa                            >  1:1097047/1‑68 (MQ=38)
                                    gCCAGTGGAGCTTTCCCGCGACTGGCTGGGGAATCTGCTTTATTATTGCGATGAGCTGGGGATTCgga                            >  1:966509/1‑66 (MQ=255)
                                    gCCAGTGGAGCTTTCCCGCGACTGGCTGGGGAATCTGCTTAATCAGTGCGATGAGCCGGGCATTCGTa                            >  1:401055/1‑68 (MQ=255)
                                     ccAGTGGAGCTTTCCCGCGCCTGGCTGGGGAATCTGCTTAATCAGTGCGCTGAGGCGGGGATTCGTaa                           >  1:1406821/1‑68 (MQ=255)
                                     ccAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTGCTTAATCAGTGCGGTGAGCTGGGCATTCGTaa                           >  1:785702/1‑68 (MQ=255)
                                     ccAGTGGAGCTTTCCCGCGACTGGCTGGGGGATCTGCTTAATCCGTGCGATGAGCTGGGCATTTGTaa                           >  1:2759698/1‑68 (MQ=255)
                                     ccAGTGGAGCTTTCCCGCGACTGGCTGGGGGATCTGCTTAATCAGTGCGATGAGCTGGGCATTCGTaa                           >  1:2619919/1‑68 (MQ=255)
                                     ccAGTGGAGCTTTCCCGCGACTGGCTGGGGAATCTGCTTAATCAGTGTGACGAGCTGGGCATTCGTaa                           >  1:344894/1‑68 (MQ=255)
                                     ccAGTGGAGCTTTCCCGCGACTGGCTGGGGAATCTGCTTAATCAGTGCGATGAGCTGGGCATTCGTaa                           >  1:35125/1‑68 (MQ=255)
                                     ccAGTGGAGCTATCCCGCNACTGGCTGGGGAATCTGCTTTATCAGTGCGATGAGCTGGGCATTCGTaa                           >  1:2573018/1‑68 (MQ=38)
                                         tGGAGCTTTCCCGCGACTGGCTGGTGAATTTGCTTAATCAGTGCGATGAGCTGGGCATTCGTAACCAg                       >  1:3383726/1‑68 (MQ=255)
                                         tGGAGCTTTCCCGCGACTGGCTGGGGAATCTGCTTAATCAGTGCGATGAGCTGGGCATTCGTAACCAg                       >  1:1618425/1‑68 (MQ=255)
                                         tGGAGCTTTCCCGCGACTGGCTGGGGAATCTGCTTAATAAGTGCGATGGGCTGGGCATTCGTAACCAg                       >  1:170668/1‑68 (MQ=255)
                                          ggAGCTTTCCCGCGACTGGCTGGTGAATCTGGTTAATCAGTGCGAAGGGCCGGGCATTCGTCACCAGt                      >  1:1733400/1‑68 (MQ=255)
                                               tttCCCGCGACTGGCTGGGGAATCTGCTTAATTAGTGCGGTGAGGTGGGGATTTGTAAACAGTTTGaa                 >  1:2219350/1‑68 (MQ=255)
                                                       gACTGGCTGGTGAATCTGCTTAATCAGTGCGATGAGCTGGGCATTCGTAACCAGTTTGAAGTGGAAGt         <  1:3199954/68‑1 (MQ=255)
                                                       gACTGGCTGGTGAATCTGCTTAATCAGTGCGATGAGCTGGGCATTCGTAACCAGTTTGAAGTGGAAGt         <  1:1651692/68‑1 (MQ=255)
                                                       gACTGGCTGGTGAATCTGCTTAATCAGTGCGATGAGCTGGGCATTCGTAACCAGTTTGAAGTGGAAGt         <  1:1611313/68‑1 (MQ=255)
                                                            gctggTGAATCTGCTTAATCAGTGCGATGAGCTGGGCATTCGTAACCAGTTTGAAGTGGAAGTCCTgg    <  1:2429506/68‑2 (MQ=255)
                                                            gctggTGAATCTGCTTAATCAGTGCGATGAGCTGGGCATTCGTAACCAGTTTGAAGTGGAAGTCCTGa    <  1:2587874/68‑1 (MQ=255)
                                                            gctggTGAATCTGCTTAATCAGTGCGATGAGCCGGGCATTCGTAACCAGTTTGAAGTGGAAGTCCTGa    <  1:2655867/68‑1 (MQ=255)
                                                              tggTGAATCTGCTTAATCAGTGCGATGAGCTGTGCATTCGTAACCAGTTTGAAGTGGAAGTCCTGAGc  <  1:2659220/68‑1 (MQ=255)
                                                                 |                                                                
AATATTGCTCAAACAGGGCATGATGCGCTGGTGTATGCCAGTGGAGCTTTCCCGCGACTGGCTGGTGAATCTGCTTAATCAGTGCGATGAGCTGGGCATTCGTAACCAGTTTGAAGTGGAAGTCCTGAGC  >  NC_000913/3302863‑3302992

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