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

CTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCT  >  NC_000913/2840592‑2840689
                                                |                                                 
cTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAt                                                  <  1:1101305/50‑1 (MQ=255)
 tGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGtt                                                 >  1:6223142/1‑50 (MQ=255)
 tGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGtt                                                 >  1:442737/1‑50 (MQ=255)
   cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAttt                                                <  1:3445969/49‑1 (MQ=255)
   cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAttt                                                <  1:742752/49‑1 (MQ=255)
   cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAttt                                                <  1:4481730/49‑1 (MQ=255)
   cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAttt                                                <  1:3007698/49‑1 (MQ=255)
   cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGAttt                                                <  1:1692274/49‑1 (MQ=255)
   cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTc                                               <  1:4624638/50‑1 (MQ=255)
   cGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTc                                               <  1:2872431/50‑1 (MQ=255)
     tCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGttttc                                              >  1:5720336/1‑47 (MQ=39)
     tCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCg                                             <  1:2121140/50‑1 (MQ=255)
       aGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGaa                                           >  1:1178298/1‑50 (MQ=39)
       aGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGa                                            >  1:3590875/1‑49 (MQ=38)
       aGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGa                                            <  1:5330750/49‑1 (MQ=255)
        gACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGaaa                                          >  1:2540319/1‑49 (MQ=38)
           ggCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGGAcccg                                       >  1:641689/1‑48 (MQ=255)
            gcgcTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGc                                      <  1:4813691/50‑1 (MQ=255)
             cgcTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGAACCGGCg                                     >  1:3130177/1‑50 (MQ=255)
             cgcTGGCTGCGGCGGCATCAGCCATTATGGCGGGGGTTTCCCAACCGGCg                                     >  1:6673120/1‑50 (MQ=255)
             cgcTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCg                                     <  1:1343210/50‑1 (MQ=255)
              gctggctgCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGAACCcgcgg                                    >  1:1348234/1‑49 (MQ=255)
              gctggctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGt                                    <  1:2051621/50‑1 (MQ=255)
              gctggctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGt                                    <  1:1138649/50‑1 (MQ=255)
              gctggctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGt                                    <  1:2097451/50‑1 (MQ=255)
              gctggctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGt                                    <  1:778779/50‑1 (MQ=255)
               ctggctgCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGtt                                   <  1:1537177/50‑1 (MQ=255)
                 ggctgCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGAAACCGCGGta                                  >  1:318647/1‑49 (MQ=25)
                  gctgCGGCGGCATCAGCCATTATGGCGGGGGTTTTCGGAACGGCGTtata                                >  1:5933170/1‑50 (MQ=255)
                    tgCGGCGGCATCAGCCATTATGGCGGGGGGTTCCCCAACTGCGTTATTtg                              >  1:1416148/1‑50 (MQ=255)
                      cggcggCATCAGCCATTATGGCGGGGGTTTTCCAACCGGCGGCTTAtggc                            >  1:236646/1‑50 (MQ=255)
                      cggcggCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTTTAtggc                            <  1:1737056/50‑1 (MQ=255)
                      cggcggCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggc                            <  1:4954829/50‑1 (MQ=255)
                      cggcggCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtg                              <  1:5587095/48‑1 (MQ=255)
                      cggcggCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtg                              <  1:1263112/48‑1 (MQ=255)
                       ggcggcATCAGCCATTATGGCGGGGGGTTTCGGACCGGGGGTTTAtgggg                           >  1:4871568/1‑48 (MQ=255)
                       ggcggcATCAGCCATTATGGCGGGGATTTTCGAACCGGCGGTATAtggcg                           >  1:919258/1‑50 (MQ=38)
                       ggcggcATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggcg                           <  1:1781321/50‑1 (MQ=255)
                        gcggcATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggcgt                          <  1:2804132/50‑1 (MQ=255)
                            cATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATAtggcgtggcg                      <  1:5049714/50‑1 (MQ=255)
                             aTCAGCCATTATGGCGGGGGGTTTCGGACCCGCGGTTTATGGCGTGGCGa                     >  1:4753656/1‑50 (MQ=255)
                             aTCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGa                     <  1:3942698/50‑1 (MQ=255)
                              tCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGAt                    <  1:4734979/50‑1 (MQ=255)
                                aGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATcc                  <  1:1250585/50‑1 (MQ=255)
                                aGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATcc                  <  1:2442026/50‑1 (MQ=255)
                                aGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATcc                  <  1:840378/50‑1 (MQ=255)
                                  ccATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGc                <  1:4741565/50‑1 (MQ=255)
                                        tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaaa          <  1:3013396/50‑1 (MQ=255)
                                        tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaaa          >  1:4603795/1‑50 (MQ=255)
                                        tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaaa          >  1:4762595/1‑50 (MQ=255)
                                        tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaaa          <  1:2101041/50‑1 (MQ=255)
                                        tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaaa          >  1:1940163/1‑50 (MQ=255)
                                        tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaa           <  1:5158075/49‑1 (MQ=255)
                                        tGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGaa           <  1:2760646/49‑1 (MQ=255)
                                          gCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACg        >  1:3156817/1‑50 (MQ=255)
                                            ggggATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTc      <  1:6390198/50‑1 (MQ=255)
                                             gggATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTcc     >  1:888308/1‑50 (MQ=255)
                                              ggATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCg    <  1:5124294/50‑1 (MQ=255)
                                              ggATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCg    <  1:1005512/50‑1 (MQ=255)
                                               gATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGc   <  1:4937278/50‑1 (MQ=255)
                                                aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt  >  1:6669204/1‑50 (MQ=255)
                                                aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt  <  1:2797897/50‑1 (MQ=255)
                                                aTTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCt  <  1:2799924/50‑1 (MQ=255)
                                                |                                                 
CTGCGTCAGACGGCGCTGGCTGCGGCGGCATCAGCCATTATGGCGGGGATTTCCGAACCGGCGTTATATGGCGTGGCGATCCGCCTGAAACGTCCGCT  >  NC_000913/2840592‑2840689

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