Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,833,262 T→G V127G (GTC→GGC)  ascF → PTS cellobiose/arbutin/salicin transporter subunit IIBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,833,2620TG98.4% 200.7 / ‑5.0 63V127G (GTC→GGC) ascFPTS cellobiose/arbutin/salicin transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (30/32);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGA  >  NZ_CP009273/2833189‑2833346
                                                                         |                                                                                    
ggTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGGAAACTGCTGGCAATg                                                                      >  1:189615/1‑90 (MQ=255)
ggTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATg                                                                      <  1:1070930/90‑1 (MQ=255)
ggTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATg                                                                      >  2:995673/1‑90 (MQ=255)
ggTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATg                                                                      <  2:1265544/90‑1 (MQ=255)
ggTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATg                                                                      <  1:973380/90‑1 (MQ=255)
ggTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATg                                                                      >  1:969988/1‑90 (MQ=255)
       ggATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCg                                                               <  2:70373/90‑1 (MQ=255)
       ggATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCACg                                                               <  1:297523/90‑1 (MQ=255)
          tCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCgaga                                                            >  2:85868/1‑90 (MQ=255)
          tCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCgaga                                                            >  2:1126218/1‑90 (MQ=255)
          tCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGAAGGGCAAACTGCTGGCAATGATCCTCgaga                                                            >  2:1503461/1‑90 (MQ=255)
           ccTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGAt                                                           >  2:270134/1‑90 (MQ=255)
           ccTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGAt                                                           >  2:1399418/1‑90 (MQ=255)
                aTGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAgcg                                                      >  2:706054/1‑90 (MQ=255)
                  gcgcTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCaa                                                                                   <  2:613752/59‑1 (MQ=255)
                    gcTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGCGATGAgcggcgg                                                  >  1:844397/1‑89 (MQ=255)
                    gcTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGt                                                  >  2:518074/1‑90 (MQ=255)
                       gATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCCGGCAATGATCCTCGAGATGAGCGGCGTGCt                                               >  1:954479/1‑90 (MQ=255)
                       gATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGGGCt                                               >  1:1102257/1‑90 (MQ=255)
                          cGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGAc                                            >  1:1297683/1‑90 (MQ=255)
                            gCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACaa                                          >  1:1084749/1‑90 (MQ=255)
                            gCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACaa                                          >  2:1548330/1‑90 (MQ=255)
                            gCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGGGCTGACaa                                          >  1:390431/1‑90 (MQ=255)
                             cACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACaaa                                         <  1:706054/90‑1 (MQ=255)
                                 aTGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGAt                                                                    >  1:1539519/1‑59 (MQ=255)
                                 aTGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGAt                                                                    <  2:1539519/59‑1 (MQ=255)
                                 aTGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGa                                     >  2:582490/1‑90 (MQ=255)
                                 aTGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGa                                     >  2:1218073/1‑90 (MQ=255)
                                 aTGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGa                                     >  2:988278/1‑90 (MQ=255)
                                 aTGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGa                                     >  1:895713/1‑90 (MQ=255)
                                 aTGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAAGAt                                    >  1:903298/1‑90 (MQ=255)
                                     cACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGGCAAAAGGATCGc                                 <  2:1540532/90‑1 (MQ=255)
                                     cACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGc                                 <  1:599313/90‑1 (MQ=255)
                                     cACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGc                                 <  2:626630/90‑1 (MQ=255)
                                     cACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGc                                 <  2:871581/90‑1 (MQ=255)
                                     cACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGc                                 <  2:271644/90‑1 (MQ=255)
                                     cACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGc                                 <  2:851323/90‑1 (MQ=255)
                                     cACCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGcc                                <  2:903298/90‑1 (MQ=255)
                                      aCCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGcc                                <  2:90602/90‑1 (MQ=255)
                                      aCCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGcc                                <  2:81520/90‑1 (MQ=255)
                                      aCCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGcc                                <  1:265928/90‑1 (MQ=255)
                                      aCCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGcc                                <  2:1540431/90‑1 (MQ=255)
                                      aCCGCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGcc                                >  2:1334484/1‑90 (MQ=255)
                                         gCTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGAc                             >  2:480851/1‑90 (MQ=255)
                                          cTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGAcc                            <  1:1437787/90‑1 (MQ=255)
                                          cTGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGAcc                            <  1:270134/90‑1 (MQ=255)
                                           tGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGc                                                <  2:1156672/69‑1 (MQ=255)
                                           tGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGc                                                >  1:1156672/1‑69 (MQ=255)
                                           tGATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGGGCTGACAAAAGGATCGCCGACCt                           >  1:1041203/1‑90 (MQ=255)
                                            gATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATg                                                                      >  2:438715/1‑46 (MQ=255)
                                            gATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATg                                                                      <  1:438715/46‑1 (MQ=255)
                                            gATCCCGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCtt                          <  1:480851/90‑1 (MQ=255)
                                                ccGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAAcc                      <  2:535321/90‑1 (MQ=255)
                                                ccGGCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAAcc                      <  1:1218073/90‑1 (MQ=255)
                                                   gCGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCAAGATGAGCGGCGTGGTGACAAAAGGATCGCCGACCTTAACCAtt                   >  1:266191/1‑90 (MQ=255)
                                                    cGATTATCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTc                  <  1:582490/90‑1 (MQ=255)
                                                        tATCGGCGGATCGATGGGCAAACTGCTGGCAATGATcc                                                                  <  1:192489/38‑1 (MQ=38)
                                                        tATCGGCGGATCGATGGGCAAACTGCTGGCAATGATcc                                                                  >  2:192489/1‑38 (MQ=255)
                                                          tCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAAtg            <  1:1503461/90‑1 (MQ=255)
                                                          tCGGCGGATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAAtg            <  1:1469028/90‑1 (MQ=255)
                                                           cggcggATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAAtgt           <  2:895713/90‑1 (MQ=255)
                                                           cggcggATCGATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAAtgt           <  2:1084749/90‑1 (MQ=255)
                                                                    gATGGGCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGa  >  1:1288359/1‑90 (MQ=255)
                                                                         |                                                                                    
GGTGCGGGGATCCTCGATGCGCTGATCGGCACCATGTCACCGCTGATCCCGGCGATTATCGGCGGATCGATGGTCAAACTGCTGGCAATGATCCTCGAGATGAGCGGCGTGCTGACAAAAGGATCGCCGACCTTAACCATTCTGAATGTGATTGGTGA  >  NZ_CP009273/2833189‑2833346

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