Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,840,6440CT21.7% 173.6 / 12.6 69S374F (TCC→TTC) ascFbeta‑glucoside specific PTS enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base C (9/45);  new base T (15/0);  total (24/45)
Fisher's exact test for biased strand distribution p-value = 2.31e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCTTATC  >  NC_000913/2840595‑2840693
                                                  |                                                 
cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTc                                                   <  1:5334505/50‑1 (MQ=255)
 gTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTcc                                                  <  1:850478/50‑1 (MQ=255)
 gTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTcc                                                  <  1:5040693/50‑1 (MQ=255)
 gTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTcc                                                  <  1:323032/50‑1 (MQ=255)
 gTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTcc                                                  <  1:1973544/50‑1 (MQ=255)
 gTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTcc                                                  <  1:1305512/50‑1 (MQ=255)
  tCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCg                                                 <  1:2984383/50‑1 (MQ=255)
   cAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTCCg                                                 >  1:2608539/1‑49 (MQ=255)
    aGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTTCgca                                               >  1:6143738/1‑48 (MQ=38)
    aGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTTCgc                                                >  1:1495234/1‑48 (MQ=38)
    aGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAGTTTCgca                                               >  1:3755433/1‑48 (MQ=38)
     gACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAAc                                              <  1:7432630/50‑1 (MQ=255)
     gACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAAc                                              <  1:5506347/50‑1 (MQ=255)
     gACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAAc                                              <  1:2546261/50‑1 (MQ=255)
      aCGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGAAcc                                             >  1:3514745/1‑50 (MQ=39)
      aCGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAAcc                                             <  1:3766180/50‑1 (MQ=255)
      aCGGCGCTCGCTGCGGCGGCATCGGCCATTATGGCGGGGATTTCCGAAcc                                             <  1:6158750/50‑1 (MQ=38)
       cGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGCAccc                                            >  1:4978730/1‑49 (MQ=255)
       cGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTCCGCAccc                                            >  1:5968394/1‑49 (MQ=255)
       cGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCg                                            <  1:5994525/50‑1 (MQ=255)
        ggCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAATTTCGGAAcgg                                            >  1:2159288/1‑48 (MQ=37)
          cgcTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCg                                         <  1:4819128/50‑1 (MQ=255)
           gctggctgCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGAACCGGCgg                                        >  1:6727597/1‑49 (MQ=255)
           gctggctgCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGAACCGGCgg                                        >  1:4079482/1‑49 (MQ=255)
           gctggctgCGGCGGCATCAGCCATTATGGCGGGGGGTTTCCAACCcgcgg                                        >  1:6525648/1‑49 (MQ=255)
           gctggctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGt                                        <  1:6440590/50‑1 (MQ=255)
           gctggctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGt                                        <  1:3931530/50‑1 (MQ=255)
              ggctgCGGCGGCATCAGCCATTATGGCGGGGGTTTCCGAACCCGCGGta                                      >  1:2799051/1‑49 (MQ=38)
               gctgCGGCGGCATCAGCCATTATGGCGGGGGTTTTCCGACCGGCGTtata                                    >  1:5931046/1‑50 (MQ=255)
               gctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTtata                                    <  1:2946972/50‑1 (MQ=255)
               gctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTtat                                     <  1:6789954/49‑1 (MQ=255)
               gctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTta                                      <  1:6884498/48‑1 (MQ=255)
                ctgCGGCGGCATCAGCCATTATGGCGGGGGTTTTCCAACCCGCGTTttat                                   >  1:5228994/1‑50 (MQ=255)
                ctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTatat                                   <  1:6175956/50‑1 (MQ=255)
                 tgCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGAAACGGCGGTttatt                                  >  1:2763173/1‑49 (MQ=255)
                   cggcggCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtgg                                 <  1:147284/49‑1 (MQ=255)
                   cggcggCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtgg                                 <  1:3992407/49‑1 (MQ=255)
                    ggcggcATCAGCCATTATGGCGGGGGTTTTCGAACCCGCGGTTTAtggcg                               >  1:6621193/1‑50 (MQ=255)
                    ggcggcATCAGCCATTATGGCGGGGGTTTTCCGACCGGGGGTTTTtggcg                               >  1:6938705/1‑50 (MQ=255)
                    ggcggcATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggcg                               <  1:7295645/50‑1 (MQ=255)
                     gcggcATCAGCCATTATGGCGGGGATTTTCGCAACGGCGGTATAttgcgg                              >  1:6079926/1‑49 (MQ=25)
                      cggcATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggcgtg                             <  1:1857643/50‑1 (MQ=255)
                      cggcATCAGCCATTATGGCGGGGAGTTTCGAAACGGCGGTATAtggcggg                             >  1:3493635/1‑48 (MQ=255)
                       ggcATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggcgtgg                            <  1:1772249/50‑1 (MQ=255)
                        gcATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggcgtggc                           <  1:6418500/50‑1 (MQ=255)
                        gcATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggcgtggc                           <  1:3433263/50‑1 (MQ=255)
                         cATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggcgtggcg                          <  1:6190340/50‑1 (MQ=255)
                         cATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggcgtggcg                          <  1:1977181/50‑1 (MQ=255)
                           tCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGAt                        <  1:7273943/50‑1 (MQ=255)
                             aGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATcc                      <  1:4001160/50‑1 (MQ=255)
                             aGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATcc                      <  1:2901980/50‑1 (MQ=255)
                             aGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATc                       <  1:2577869/49‑1 (MQ=255)
                             aGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATc                       <  1:2555194/49‑1 (MQ=255)
                             aGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATc                       <  1:2094978/49‑1 (MQ=255)
                               ccATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGc                    <  1:1272372/50‑1 (MQ=255)
                               ccATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGc                    <  1:4415107/50‑1 (MQ=255)
                                 aTTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCt                  <  1:3717610/50‑1 (MQ=255)
                                  ttATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTg                 <  1:955974/50‑1 (MQ=255)
                                     tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaa               <  1:3404068/49‑1 (MQ=255)
                                      ggCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAAc             >  1:6058497/1‑50 (MQ=255)
                                       gCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAAc             <  1:5218352/49‑1 (MQ=255)
                                       gCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAAc             <  1:1295146/49‑1 (MQ=255)
                                         ggggATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTc          <  1:1161426/50‑1 (MQ=255)
                                          gggATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTcc         <  1:4572383/50‑1 (MQ=255)
                                          gggATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTcc         <  1:3760122/50‑1 (MQ=255)
                                            gATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGc       <  1:169611/50‑1 (MQ=255)
                                              aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt      >  1:6408167/1‑50 (MQ=255)
                                              aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt      <  1:3327128/50‑1 (MQ=255)
                                              aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt      >  1:2963063/1‑50 (MQ=255)
                                              aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt      <  1:819463/50‑1 (MQ=255)
                                              aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt      <  1:2658921/50‑1 (MQ=255)
                                              aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt      >  1:2447789/1‑50 (MQ=255)
                                               tttCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCtt     <  1:58043/50‑1 (MQ=255)
                                               tttCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt      >  1:6867104/1‑49 (MQ=255)
                                                 tCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCTTa    <  1:3103353/49‑1 (MQ=255)
                                                  ccGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCTTAt   <  1:2675182/49‑1 (MQ=255)
                                                  ccGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCTTATc  >  1:632268/1‑50 (MQ=255)
                                                  |                                                 
CGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCTTATC  >  NC_000913/2840595‑2840693

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