Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913821,4400CT28.9% 67.6 / 18.8 45F216F (TTC→TTT) ybhMBax1‑I family protein YbhM
Reads supporting (aligned to +/- strand):  ref base C (2/30);  new base T (13/0);  total (15/30)
Fisher's exact test for biased strand distribution p-value = 1.44e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTT  >  NC_000913/821377‑821495
                                                               |                                                       
gtgGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGcttct                                                        <  1:872497/65‑1 (MQ=255)
gtgGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGcttct                                                        <  1:2619975/65‑1 (MQ=255)
  ggAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCttt                                                      <  1:2100532/65‑1 (MQ=255)
  ggAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCttt                                                      >  1:3153827/1‑65 (MQ=255)
    aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCtttttttt                                                    >  1:2498422/1‑63 (MQ=255)
    aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCtttttttg                                                    >  1:556353/1‑63 (MQ=255)
    aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCtttttttg                                                    >  1:2316484/1‑63 (MQ=255)
    aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTgt                                                    >  1:716960/1‑64 (MQ=255)
    aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTgg                                                    >  1:1202054/1‑65 (MQ=255)
       ccGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTggggt                                                 >  1:1926859/1‑63 (MQ=255)
       ccGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGAt                                                 <  1:1489712/65‑1 (MQ=255)
        cGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATa                                                <  1:2067266/65‑1 (MQ=255)
          gTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACt                                              <  1:3178543/65‑1 (MQ=255)
            tCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTTGGATACttt                                            >  1:1885842/1‑64 (MQ=255)
            tCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTg                                            <  1:2768297/65‑1 (MQ=255)
               ttGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCcg                                         <  1:1616484/65‑1 (MQ=255)
                     ttttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAg                                   <  1:1060960/65‑1 (MQ=255)
                      tttACAGGTGCATTAACCCTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc                                  <  1:406968/65‑1 (MQ=255)
                      tttACAGGTGCATTAACACTCTATTCCTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc                                  <  1:615598/65‑1 (MQ=255)
                      tttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc                                  <  1:3147785/65‑1 (MQ=255)
                      tttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc                                  <  1:1073557/65‑1 (MQ=255)
                      tttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc                                  <  1:1972483/65‑1 (MQ=255)
                      tttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc                                  <  1:173478/65‑1 (MQ=255)
                      tttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc                                  <  1:1680546/65‑1 (MQ=255)
                      tttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc                                  <  1:1720445/65‑1 (MQ=255)
                       ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                 <  1:3045281/65‑1 (MQ=255)
                       ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                 <  1:112123/65‑1 (MQ=255)
                       ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                 <  1:761312/65‑1 (MQ=255)
                       ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                 <  1:557809/65‑1 (MQ=255)
                       ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                 <  1:1495726/65‑1 (MQ=255)
                       ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                 <  1:373342/65‑1 (MQ=255)
                       ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                 <  1:1539334/65‑1 (MQ=255)
                       ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                 <  1:1999342/65‑1 (MQ=255)
                       ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                 <  1:3034601/65‑1 (MQ=255)
                       ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                 <  1:1743078/65‑1 (MQ=255)
                       ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                 <  1:2296241/65‑1 (MQ=255)
                       ctACAGGTGGATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt                                 <  1:498948/64‑1 (MQ=255)
                         aCAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTTTGGTACCTGCCGCGCTAGCtat                               >  1:2012836/1‑65 (MQ=255)
                            ggTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAAc                            <  1:1609297/65‑1 (MQ=255)
                             gTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACg                           <  1:838614/65‑1 (MQ=255)
                               gCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTACCGCGATAGCTATAACGCt                         <  1:843973/65‑1 (MQ=255)
                                 aTTAACACTCTATTACTATTTTATCGGCTTTTTTTGGATACTTGCCCCGGTAGCTATAACGCTTg                       >  1:3040114/1‑65 (MQ=255)
                                 aTTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTg                       <  1:2580982/65‑1 (MQ=255)
                                   tAACACTCTATTACTATTTTATCGGCTTTTTTTGGGTACTTGCCGCGATAGCTATAACGCTTGTc                     >  1:3206241/1‑65 (MQ=255)
                                    aaCACTCTATTACTATTTTATCGGCTTTTTTGGGGTACTTGCCGCGATAACTATAACGCTTGtgt                    >  1:245687/1‑63 (MQ=255)
                                       acTCTATTACTATTTTATCGGCTTTTTTTGGGTACTTGCCGCGATAGCTATAACACTTTTCTggg                 >  1:1342926/1‑64 (MQ=255)
                                       acTCTATTACTATTTTATCGGCTTTTTTTGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGc                 >  1:235648/1‑65 (MQ=255)
                                       acTCTATTACTATTTTATCGGCTTTTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGc                 >  1:269166/1‑65 (MQ=255)
                                       acTCTATTACTATTTTATCGGCTTTTTTGGGATACTTGCCGCGATAGCTATAAAGCTTGTCTGGc                 >  1:1219595/1‑65 (MQ=255)
                                       acTCTATTACTATTTTATCGGCTTTTTTGGGATACTTGCCCCGCTAGCTCTAAAGCTTGCCTggg                 >  1:1632080/1‑64 (MQ=255)
                                       acTCTATTACTATTTTATCGGCTTCTTTTGGGTACTTGGCGCGGATGGTCTAACGCCTGTCTGGc                 >  1:3110302/1‑65 (MQ=255)
                                       acTCTATTACTATTTTATCGGCTTCTTTGGGGTACTTGCCGCGATAGCTATAACGCTTGTTTggg                 >  1:2766585/1‑64 (MQ=255)
                                        ctctATTACTATTTTATCGGCTTTTTTTGGATACTTTCCGCGATAGCTATAACGCTTGTCTGGCa                >  1:2443534/1‑65 (MQ=255)
                                         tctATTACTATTTTATCGGCTTTTTTGGGAGACTTGCCCCGATAGCTATTACGCCTGTCTGGCaa               >  1:436281/1‑65 (MQ=255)
                                           tATTACTATTTTATCGGCTTTTTTTGGGGACTTGCCGCGATAGCTATAACGCTTGTTTGGCAAAg             >  1:3027542/1‑65 (MQ=255)
                                           tATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAg             <  1:1417868/65‑1 (MQ=255)
                                             ttgcTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGc           <  1:2580063/62‑1 (MQ=255)
                                               aCTATTTTATCGGCTTTTTTTGGATACTTGCCGCGGTAGCTATAACGCTTGTCTTGCCAAGGcct         >  1:475423/1‑63 (MQ=255)
                                                      tATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGtt  <  1:1735261/65‑1 (MQ=255)
                                                               |                                                       
GTGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTT  >  NC_000913/821377‑821495

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