Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913821,4400CT31.6% 53.0 / 17.3 38F216F (TTC→TTTybhMBax1‑I family protein YbhM
Reads supporting (aligned to +/- strand):  ref base C (4/22);  new base T (12/0);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 6.72e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTGTGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGT  >  NC_000913/821374‑821494
                                                                  |                                                      
gTTGTGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGcttct                                                       <  1:2080734/68‑1 (MQ=255)
  tgtgGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCttt                                                     <  1:1643921/68‑1 (MQ=255)
    tgGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTgg                                                   <  1:2644188/68‑1 (MQ=255)
     ggAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTggg                                                  <  1:2528534/68‑1 (MQ=255)
     ggAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCCTTTTTtgg                                                  >  1:1417198/1‑68 (MQ=255)
       aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTtgggt                                                >  1:3354305/1‑66 (MQ=255)
       aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTTGGAt                                                >  1:1412825/1‑68 (MQ=255)
       aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTGGGAt                                                >  1:352100/1‑68 (MQ=255)
       aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTtgggt                                                >  1:1389733/1‑66 (MQ=255)
       aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATAGGCTTTTTTtgggt                                                >  1:3017448/1‑66 (MQ=255)
          ccGGTTCGTTGCGCTTTTACAGGTGCATTAAGACTCTATTACTATTTTATCGGCTTTTTTTGGGTACt                                             >  1:1935788/1‑68 (MQ=255)
          ccGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTTTGATACt                                             >  1:2945681/1‑68 (MQ=255)
          ccGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACt                                             <  1:2616021/68‑1 (MQ=255)
          ccGGTTCCTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACt                                             <  1:2952101/68‑1 (MQ=255)
            ggTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCAGCTTCTTTGGGATACTTg                                           <  1:2175259/68‑1 (MQ=255)
                 gTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCcgcg                                      <  1:1492341/68‑1 (MQ=255)
                      gcTTTTACAGGTGCATTAACGCTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc                                 <  1:1747217/68‑1 (MQ=255)
                      gcTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc                                 <  1:715747/68‑1 (MQ=255)
                      gcTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc                                 <  1:1417730/68‑1 (MQ=255)
                       cTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                <  1:266521/68‑1 (MQ=255)
                       cTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                <  1:89457/68‑1 (MQ=255)
                       cTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                <  1:2353378/68‑1 (MQ=255)
                            aCAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAAc                           <  1:834362/68‑1 (MQ=255)
                             cAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATGCTTGCCGCGATAGCTATAACg                          <  1:1477132/68‑1 (MQ=255)
                             cAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACg                          <  1:3192759/68‑1 (MQ=255)
                              aGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTTGGGTACTTGCCCCGCTCGCTATAACGc                         >  1:2210028/1‑68 (MQ=255)
                                    aTTAACACTCTATTACTATTTTATCGGCTTCCTTGGGATAACTGTCGCGATTGCTATAAAGCTTGGCt                   >  1:1650907/1‑68 (MQ=255)
                                      tAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTgg                 <  1:1979191/68‑1 (MQ=255)
                                         cacTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCaa              <  1:1339084/68‑1 (MQ=255)
                                          acTCTATTACTATTTTATCGGCTTTTTTTGGGTACTTTCCGCGGTAGCTATAACGCTTGTCTGGGaaa             >  1:1896485/1‑68 (MQ=255)
                                          acTCTATTACTATTTTATCGGCTTTTTTTGGATACTTTCCGCGAGAGCTATAACGCTTGTCTGGCaaa             >  1:2867892/1‑68 (MQ=255)
                                          acTCTATTACTATTTTATCGGCTTCTTTTGGGTTCTTGTCGCGGTAGCTATAACGCTTGTCTGGGaaa             >  1:688561/1‑68 (MQ=255)
                                          acTCTANTACTATTTTATCGGCTTTTTTTGGATACTTGCCGCGATAGCTATAAAGCTTGTCTGGCaaa             >  1:2029583/1‑68 (MQ=255)
                                            tctATTACTATTTTATCGGCTTTTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAgg           >  1:268590/1‑68 (MQ=255)
                                            tctATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAgg           <  1:401674/68‑1 (MQ=255)
                                             ctATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGc          <  1:2799809/68‑1 (MQ=255)
                                             ctATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGc          <  1:3059496/68‑1 (MQ=255)
                                             ctATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGc          <  1:88615/68‑1 (MQ=255)
                                             ctATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGc          <  1:1185557/68‑1 (MQ=255)
                                               aTTACTATTTTATCGGCTTCTTTTGGATAATTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCAt        >  1:2840591/1‑68 (MQ=255)
                                                    tATTTTATCGGCTTTTTTTGGGTACTTTCCCCGACAGGTATAACACTTTTCTGGCAAAAGGAAAcgcg   >  1:2017991/1‑68 (MQ=255)
                                                    tATTTTATCGGCTTTTTTTGGGTACTTGCCCCGATAGGTATAAACCTTGTCTGGCCAAGGCAAAagcg   >  1:2749405/1‑68 (MQ=255)
                                                    tATTTTATCGGCTTTTTTTGGGTACTTGCCCCGATAGCTATAACGCTTGTCTGGCAAAGGCATAcgcg   >  1:1809670/1‑68 (MQ=255)
                                                    caTTTTATCGGCTTCTTTTGGATACTTTCCGCGATAGCTATAACCCTTGTCTGGCAAAGGCATAcgcg   >  1:2736264/2‑68 (MQ=255)
                                                     aTTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGt  <  1:991922/68‑1 (MQ=255)
                                                                  |                                                      
GTTGTGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGT  >  NC_000913/821374‑821494

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