Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_012967 1,881,721 G→T G21* (GGA→TGA)  manY → PTS mannose transporter subunit IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0129671,881,7210GT99.6% 735.0 / ‑3.3 247G21* (GGA→TGA) manYPTS mannose transporter subunit IIC
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (145/100);  total (147/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

Only 100 of 253 total aligned reads displayed.

GAAGTACAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGGGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATCGTTCTT  >  NC_012967/1881653‑1881789
                                                                    |                                                                    
gaAGTACAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGAt                                                                   >  2:582517/1‑72 (MQ=255)
gaAGTACAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGAt                                                                   >  1:5760492/1‑72 (MQ=255)
gaAGTACAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGAt                                                                   >  1:2339390/1‑72 (MQ=255)
gaAGTACAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGAAGGTATGTGAt                                                                   >  1:2387603/1‑72 (MQ=255)
 aaGTACAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTNTATCGCAGGTATGTGATc                                                                  >  2:698845/1‑72 (MQ=255)
 aaGTACAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATc                                                                  >  2:231265/1‑72 (MQ=255)
  aGTACAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCa                                                                 >  1:4862152/1‑72 (MQ=255)
  aGTACAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCa                                                                 <  2:791188/72‑1 (MQ=255)
   gTACAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCaa                                                                >  2:6072059/1‑72 (MQ=255)
    tACAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAAt                                                               >  1:8061224/1‑72 (MQ=255)
    tACAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAAt                                                               >  2:7412165/1‑72 (MQ=255)
     aCAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATc                                                              <  2:4421006/72‑1 (MQ=255)
     aCAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATc                                                              >  1:8129256/1‑72 (MQ=255)
       aaTGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCt                                                            <  1:2831946/72‑1 (MQ=255)
       aaTGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCt                                                            >  2:5405153/1‑72 (MQ=255)
        aTGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTc                                                           >  2:3622011/1‑72 (MQ=255)
         tGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCg                                                          <  1:3592831/72‑1 (MQ=255)
          ggAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGa                                                         >  1:4595066/1‑72 (MQ=255)
            agaTTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATg                                                       <  1:2123553/72‑1 (MQ=255)
             gaTTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGa                                                      >  1:3791131/1‑72 (MQ=255)
             gaTTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGGATGTGATCAATCCTCGATGa                                                      >  2:3593308/1‑72 (MQ=255)
               ttACCACTCTTCAAATTGCGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAAt                                                    >  2:2747572/1‑72 (MQ=255)
                 aCCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAAttt                                                  <  2:6731793/72‑1 (MQ=255)
                  ccACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTc                                                 >  2:2422540/1‑72 (MQ=255)
                   cACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCa                                                <  1:6013444/72‑1 (MQ=255)
                   cACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCa                                                <  1:2046411/72‑1 (MQ=255)
                   cACTCTTCAAATTGTGCTGGTATTTATCGNAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAGTTTCa                                                >  2:744950/1‑72 (MQ=255)
                    aCTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAg                                               >  2:2787857/1‑72 (MQ=255)
                     ctctTCAAATTGTGCTGGTATTTTTCGTAGCNNNNNTCGCAGGTATGTGATCAATCCTCGATGAATTTCAGt                                              >  1:6135346/1‑72 (MQ=255)
                        ttCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTc                                           >  1:6971503/1‑72 (MQ=255)
                        ttCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTc                                           <  2:5211326/72‑1 (MQ=255)
                        ttCAAATTGTGCTGGTATTTATCGTAGCCTGNNNNNCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTc                                           >  1:2735329/1‑72 (MQ=255)
                         tccaaTTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCa                                          >  2:6700712/4‑72 (MQ=255)
                         tCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCa                                          <  2:8568764/72‑1 (MQ=255)
                            aaTTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCg                                       >  2:658619/1‑72 (MQ=255)
                            aaTTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCg                                       >  1:3940175/1‑72 (MQ=255)
                            aaTTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAAGTTCAGTTTCACCg                                       <  1:3444587/72‑1 (MQ=255)
                             aTTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGt                                      >  1:1045719/1‑72 (MQ=255)
                             aTTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGt                                      >  2:5562037/1‑72 (MQ=255)
                               tgtgCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTcc                                    >  1:7247147/1‑72 (MQ=255)
                                gtgCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCGCCGTCCg                                   <  2:4270689/72‑1 (MQ=255)
                                gtgCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCg                                   >  1:7737445/1‑72 (MQ=255)
                                gtgCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCg                                   >  1:3563875/1‑72 (MQ=255)
                                gtgCTGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCg                                   <  2:8634492/72‑1 (MQ=255)
                                    tGGTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCGGCTaa                               <  1:2993837/72‑1 (MQ=255)
                                     ggTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAAt                              <  1:7288416/72‑1 (MQ=255)
                                      gTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATc                             >  2:3708788/1‑72 (MQ=255)
                                      gTATTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATc                             >  2:2937060/1‑72 (MQ=255)
                                       tATTTATCGTAGCCTGTGTCGCAGGTATGTGATCAATCCTCGATGCATTTCAGTTTCACCGTCCGCTAATcg                            >  1:3231874/1‑72 (MQ=255)
                                       tATTTATCGTAGCCTGTATCGCAGGTAGGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATcg                            <  1:7705587/72‑1 (MQ=255)
                                        aTTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATcgc                           <  1:2209519/72‑1 (MQ=255)
                                        aTTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATcgc                           <  2:5709584/72‑1 (MQ=255)
                                        aTTTATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATcgc                           >  2:7426366/1‑72 (MQ=255)
                                         tttATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATcgcg                          >  2:1970679/1‑72 (MQ=255)
                                          ttATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCgt                         >  1:1703506/1‑72 (MQ=255)
                                          ttATCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCgt                         >  1:5341024/1‑72 (MQ=255)
                                            atggTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCgtgt                       <  2:3165945/69‑1 (MQ=255)
                                            aTCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCgtgt                       <  2:2782713/72‑1 (MQ=255)
                                            aTCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCgtgt                       >  1:4944659/1‑72 (MQ=255)
                                            aTCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCgtgt                       <  2:7262443/72‑1 (MQ=255)
                                             tCGTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTa                      >  1:2012857/1‑72 (MQ=255)
                                               gTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTAcc                    <  2:4966929/72‑1 (MQ=255)
                                               gTAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTAcc                    >  2:6876013/1‑72 (MQ=255)
                                                tAGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTAccc                   >  2:4662421/1‑72 (MQ=255)
                                                 aGCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCt                  <  2:3800969/72‑1 (MQ=255)
                                                  gCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCtg                 >  2:5933771/1‑72 (MQ=255)
                                                  gCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCtg                 <  2:6321402/72‑1 (MQ=255)
                                                  gCCTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCtg                 >  2:8381771/1‑72 (MQ=255)
                                                   ccTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCtgg                >  1:1118100/1‑72 (MQ=255)
                                                    cTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCtggt               >  2:3857070/1‑72 (MQ=255)
                                                    cTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCtggt               >  1:179026/1‑72 (MQ=255)
                                                    cTGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCtggt               <  2:5053772/72‑1 (MQ=255)
                                                    cTGTATCGCAGGTANNNGATCAATCCTCGATGAANTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCtggt               <  1:6231036/72‑1 (MQ=255)
                                                     tGTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCtgggt              >  2:3084409/1‑70 (MQ=255)
                                                      gTATCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCtggtgg             >  2:1837679/1‑72 (MQ=255)
                                                        aTCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGt           <  1:8541656/72‑1 (MQ=255)
                                                        aTCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGt           <  2:1202389/72‑1 (MQ=255)
                                                         tCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCCCCGTCCGCTAATCGCGTGTACCCTGGTGGGTa          <  2:2093005/72‑1 (MQ=255)
                                                         tCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTa          <  2:1252553/72‑1 (MQ=255)
                                                         tCGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTa          >  1:7110043/1‑72 (MQ=255)
                                                          cGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTAt         >  1:8395278/1‑72 (MQ=255)
                                                          cGCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTAt         >  1:6674836/1‑72 (MQ=255)
                                                           gCATGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATc        >  1:3301226/1‑72 (MQ=255)
                                                           gCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGTTATc        >  1:6450534/1‑72 (MQ=255)
                                                           gCAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATc        <  2:5534487/72‑1 (MQ=255)
                                                            cAGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATCg       >  1:8602328/1‑72 (MQ=255)
                                                             aGGTATGTGATCAATCCTCGATGAATTTNAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGGATCGt      >  2:3211718/1‑72 (MQ=255)
                                                             aGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATCGt      >  1:4318617/1‑72 (MQ=255)
                                                             aGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATCGt      <  2:7823823/72‑1 (MQ=255)
                                                             aGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATCGt      >  2:8114196/1‑72 (MQ=255)
                                                             aGGTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATCGt      <  2:2171786/72‑1 (MQ=255)
                                                              ggTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATTGtt     >  1:5171077/1‑72 (MQ=255)
                                                              ggTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATCGtt     <  2:7560248/72‑1 (MQ=255)
                                                              ggTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGCGTACCCTGGTGGGTATCGtt     <  2:4896806/72‑1 (MQ=255)
                                                               gTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGTGTATCGTTc    >  1:4795643/1‑72 (MQ=255)
                                                               gTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATCGTTc    <  2:6439000/72‑1 (MQ=255)
                                                               gTATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATCGTTc    >  1:6832503/1‑72 (MQ=255)
                                                                tATGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATCGTTCt   <  2:4500366/72‑1 (MQ=255)
                                                                tATGCGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATCGTTCt   >  2:1531729/1‑72 (MQ=255)
                                                                 aTGTGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGTGTATCGTTCtt  <  2:1003172/72‑1 (MQ=255)
                                                                    |                                                                    
GAAGTACAATGGAGATTACCACTCTTCAAATTGTGCTGGTATTTATCGTAGCCTGTATCGCAGGTATGGGATCAATCCTCGATGAATTTCAGTTTCACCGTCCGCTAATCGCGTGTACCCTGGTGGGTATCGTTCTT  >  NC_012967/1881653‑1881789

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