Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,840,6400AG21.3% 193.6 / 17.7 75I373V (ATT→GTT) ascFbeta‑glucoside specific PTS enzyme IIBC component
Reads supporting (aligned to +/- strand):  ref base A (6/53);  new base G (16/0);  total (22/53)
Fisher's exact test for biased strand distribution p-value = 8.73e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCT  >  NC_000913/2840591‑2840689
                                                  |                                                 
aCTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGa                                                   <  1:3644895/50‑1 (MQ=255)
 cTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAt                                                  <  1:5414897/50‑1 (MQ=255)
 cTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAt                                                  <  1:6129793/50‑1 (MQ=255)
 cTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAt                                                  <  1:3785668/50‑1 (MQ=255)
  tGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGtt                                                 >  1:3450140/1‑50 (MQ=255)
    cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAttt                                                <  1:6206222/49‑1 (MQ=255)
    cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAtt                                                 <  1:4109163/48‑1 (MQ=255)
    cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAtt                                                 <  1:5270294/48‑1 (MQ=255)
    cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAtt                                                 <  1:3765976/48‑1 (MQ=255)
    cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTc                                               <  1:5334505/50‑1 (MQ=255)
     gTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGTGAgtttc                                              >  1:1220311/1‑48 (MQ=38)
     gTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTcc                                              <  1:1973544/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: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:2546261/50‑1 (MQ=255)
         gACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAAc                                          <  1:5506347/50‑1 (MQ=255)
         gACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAAc                                          <  1:7432630/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:3992407/49‑1 (MQ=255)
                       cggcggCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtgg                             <  1:147284/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:1977181/50‑1 (MQ=255)
                             cATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggcgtggcg                      <  1:6190340/50‑1 (MQ=255)
                               tCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGAt                    <  1:7273943/50‑1 (MQ=255)
                                 aGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATcc                  <  1:2901980/50‑1 (MQ=255)
                                 aGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATcc                  <  1:4001160/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:3760122/50‑1 (MQ=255)
                                              gggATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTcc     <  1:4572383/50‑1 (MQ=255)
                                                gATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGc   <  1:169611/50‑1 (MQ=255)
                                                  aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt  >  1:6408167/1‑50 (MQ=255)
                                                  aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt  >  1:2447789/1‑50 (MQ=255)
                                                  aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt  <  1:2658921/50‑1 (MQ=255)
                                                  aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt  >  1:2963063/1‑50 (MQ=255)
                                                  aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt  <  1:3327128/50‑1 (MQ=255)
                                                  aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt  <  1:819463/50‑1 (MQ=255)
                                                  |                                                 
ACTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCT  >  NC_000913/2840591‑2840689

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: