Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,856,6450TG25.8% 60.3 / 11.7 31K54N (AAA→AACygbAprotein YgbA
Reads supporting (aligned to +/- strand):  ref base T (8/15);  new base G (8/0);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 2.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCTGTTTCATCTCTTCACGCTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAACACACATTTATCCAGCCGCTTTTGCGCGTAAACAAACAACGCCTC  >  NC_000913/2856581‑2856704
                                                                |                                                           
tCTGTTTCATCTCTTCACGCTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACGGGCTGGt                                                             >  1:1792378/1‑65 (MQ=255)
tCTGTTTCATCTCTTCACGCTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGt                                                             >  1:656873/1‑65 (MQ=255)
tCTGTTTCATCTCTTCACGCTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGt                                                             >  1:2206128/1‑65 (MQ=255)
    tttCATCTCTTCACGCTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGGTtct                                                         >  1:3016799/1‑65 (MQ=255)
         tctcTTCACGCTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCtttg                                                    >  1:2456189/1‑63 (MQ=255)
          ctctTCACGCTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGc                                                   <  1:1022938/65‑1 (MQ=255)
          ctctTCACGCTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGGTTTTTTTCGc                                                   >  1:1906140/1‑65 (MQ=255)
             ttCACGCTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGa                                                <  1:1801188/65‑1 (MQ=255)
               cACGCTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAac                                              <  1:2288583/65‑1 (MQ=255)
               cACGCTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAac                                              <  1:345831/65‑1 (MQ=255)
               cACGCTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAac                                              <  1:1290179/65‑1 (MQ=255)
                  gCTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGGTTTTCTTCGCCGGacaca                                           >  1:1703517/1‑65 (MQ=255)
                   cTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAacacac                                          <  1:2736452/65‑1 (MQ=255)
                    tttGCAGGCTGGTAGCAGTGCGCCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAAcacaca                                         <  1:1090484/65‑1 (MQ=255)
                    tttGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAAcacaca                                         <  1:307258/65‑1 (MQ=255)
                    tttGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAAcacaca                                         <  1:2873914/65‑1 (MQ=255)
                    tttGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAAcacaca                                         <  1:2744781/65‑1 (MQ=255)
                    tttGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAAcacaca                                         <  1:505168/65‑1 (MQ=255)
                    tttGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAAcacaca                                         <  1:530147/65‑1 (MQ=255)
                    tttGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAAcacaca                                         <  1:2412795/65‑1 (MQ=255)
                    tttGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAAcacaca                                         <  1:725950/65‑1 (MQ=255)
                                  aCAGTGCACCGGACACTGCTTACAGGCTGGTTTTTCTTCTCCCAACACCCACTTATCCAGCCGCt                           >  1:3194846/1‑65 (MQ=255)
                                  aCAGTGCACCGGACACTGCTTACAGGCTGGGTTTTCTTCGCCGAAACCACATTTATCCAGccccc                           >  1:2347417/1‑62 (MQ=255)
                                  aCAGTGCACCGGACACTGCTTACAGGCTGGGTTCTCTTCCCCGGACACACATTTATCCAGCCCCt                           >  1:1171416/1‑65 (MQ=255)
                                     gTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGCACACACCTTTATCCAGCCGCttct                        >  1:2083592/1‑63 (MQ=255)
                                     gTGCACCGGACACTGCTTACAGGCTGGGTTCTCTTTGCCGAACACACATTTATCCAGCCGCtttt                        >  1:1314066/1‑65 (MQ=255)
                                     gTGCACCGGACACTGCTTACAGGCTGGGTTCTCTTTGCCGAACACACATTTATCCAGCCGCtttt                        >  1:329802/1‑65 (MQ=255)
                                     gTGCACCGGACACTGCTTACAGGCTGGGTTCTCTTCGCCGAACACACATTTATCCAGCCGCtttt                        >  1:1453876/1‑65 (MQ=255)
                                         aCCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAACACACATTTATCCAACCCCTTTTtcgc                    >  1:1731742/1‑65 (MQ=255)
                                            ggACACTGCTTACAGGCTGGTTTCTCTTCGCCGAACACACCTTTATCCAGCCGCTTTTGCGCGTa                 >  1:41439/1‑65 (MQ=255)
                                              acacTGCTTACAGGCTGGTTTCTCTTCGCCGAACACACATTTATCCAGCCGCTTTTTCGCGTaaa               >  1:1978682/1‑65 (MQ=255)
                                              acacTGCTTACAGGCTGGTTTCTCTTCGCCGAACACACATTTATCCAGCCGCTTTTGCGCGTaaa               >  1:3088638/1‑65 (MQ=255)
                                              acacTGCTTACAGGCTGGGTTCTCTTCGCCGAACACACATTTATCCCTCCGCTTTTTCTCGTtaa               >  1:593389/1‑65 (MQ=255)
                                              acacTGCTTACAGGCTGGGTTCTCTTCGCCGAAAACACATTTATCCAGCCGCTTTTGCGCGTaaa               >  1:440037/1‑65 (MQ=255)
                                              acacTGCTTACAGGCTGGGTTCTCTTCGCCCAAAACACATTTATTCAGCCGCTTTTGCGCGGaaa               >  1:2263902/1‑65 (MQ=255)
                                                acTGCTTACAGGCTGGTTTCTCTTCGCCGAACACACATTTATCCAGCCGCTTTTGCGCGTaaaca             >  1:2206259/1‑65 (MQ=255)
                                                acTGCTTACAGGCTGGTTTCTCTTCCCCGAACACACATTTATCCAGCCGCCTTTGCGCCTaaaca             >  1:2242685/1‑65 (MQ=255)
                                                acTGCTTACAGGCTGGTTTCCCTTCGCCGAACACACATTTATCCAGCCGCTTTTGCCCGTcaaca             >  1:41282/1‑65 (MQ=255)
                                                acTGCTTACAGGCTGGGTTCTCTTTGCCGAACACACCTTTATCCAGCCGCTTTTTCGCGTaaaca             >  1:2680997/1‑65 (MQ=255)
                                                acTGCTTACAGGCTGGGTTCTCTTGTCCGAAACCCCATTTATCCCAGCGCTTTTTCCCGTaaaca             >  1:184163/1‑65 (MQ=255)
                                                acTGCTTACAGGCTGGGTTCTCTTCGCCGAACACACATTTATCCGGCCGCTTTTGCGCGTaaaca             >  1:1193530/1‑65 (MQ=255)
                                                acTGCTTACAGGCTGGGTTCTCTTCGCCGAACACACATTTATCCAGCCGCTTTTGCGCGTaaaca             >  1:1454680/1‑65 (MQ=255)
                                                 cTGCTTACAGGCTGGTTTCTCTTCGCCGAACACACATTTATCCAGCCGCTTTTGCGCGTaaacaa            >  1:3306684/1‑65 (MQ=255)
                                                   gCTTACAGGCTGGTTTCTCTTCGCCGAACCCACATTTATCCAGCCGCTTTTGCGCGTaaacaaac          >  1:1009162/1‑65 (MQ=255)
                                                        cAGGCTGGTTTCTCTTCGCCGAACACACATTTATCCAGCCGCTTTTGCGCGTAAACAAACAACGc     >  1:460179/1‑65 (MQ=255)
                                                        cAGGCTGGTTTCTCTTCGCCGAACACACATTTATCCAGCCGCTTTTGCGCGTAAACAAACAACGc     >  1:108914/1‑65 (MQ=255)
                                                         aGGCTGGTTTCTCTTCGCCGAACCCACATTTATCCAGCCGCTTTTGCGCGTAAACAAACAACGcc    <  1:3004906/65‑1 (MQ=255)
                                                         aGGCTGGTTTCTCTTCGCCGAACACACATTTATCCAGCCGCTTTTGCGCGTAAACAAACAACGcc    <  1:2119767/65‑1 (MQ=255)
                                                           gCTGGTTTCTCTTCGCCGAACACACATTTATCCAGCCGCTTTTGCGCGTAAACAAACAACGCCTc  <  1:569023/65‑1 (MQ=255)
                                                           gCTGGTTTCTCTTCGCCGAACACACATTTATCCAGCCGCTTTTGCGCGTAAACAAACAACGCCTc  <  1:647425/65‑1 (MQ=255)
                                                                |                                                           
TCTGTTTCATCTCTTCACGCTTTGCAGGCTGGTAACAGTGCACCGGACACTGCTTACAGGCTGGTTTCTCTTCGCCGAACACACATTTATCCAGCCGCTTTTGCGCGTAAACAAACAACGCCTC  >  NC_000913/2856581‑2856704

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