Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,280,494 T→C V93A (GTG→GCG)  agaV → N‑acetyl‑D‑galactosamine specific PTS enzyme IIB component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,280,4940TC98.2% 184.9 / ‑4.1 54V93A (GTG→GCG) agaVN‑acetyl‑D‑galactosamine specific PTS enzyme IIB component
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (26/27);  minor base A (0/1);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

ATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGG  >  NC_000913/3280358‑3280627
                                                                                                                                        |                                                                                                                                     
aggaaaaTGGTACGGGAAGAGGGGACCCCCTTACTTTTCGGGCCGCGAAAAAAATTTATCAAAAACATTCACCGCCCCGCCGATCGAAAAAAAATCCTGTTGTTTTGTAAAACCCCCGCCGATTTCCTGACGCTGGCGaa                                                                                                                                    <  2:1629614/136‑1 (MQ=255)
 tGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGaaa                                                                                                                                   >  1:828471/1‑140 (MQ=255)
 tGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGaaa                                                                                                                                   >  1:869018/1‑140 (MQ=255)
      aTGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTgg                                                                                                                              >  2:780511/1‑140 (MQ=255)
      aTGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTgg                                                                                                                              >  1:426155/1‑140 (MQ=255)
        ggTCCTGGCAGAAGGGATCGCCGTACTTTTCTGGACGCTGAAAAAATTTATCGAAAACATTCACCGCGCCGCCGATCGAAAAAAAATCCTGCTGGTTTGAAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCg                                                                                                                            <  2:1323953/140‑1 (MQ=255)
        ggTACGGGAAGAGGGGATCGCCGTACTTTTCGGGACGCTAAAAAAAGTTATCAAAAACATTCACCGCCCCGCCGATCGAAAGAAAATCCTGCTGGTTTGAAAACCCCCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCg                                                                                                                            <  2:1393551/140‑1 (MQ=255)
              ggCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCgg                                                                                                                      >  2:605401/1‑140 (MQ=255)
                  gaggggATCGCCGTCCTTTTCTGGACGCTAAAAAAAGTTATCGAAAACTTTCACCGCGCCGCCGATCGACAAAAAATCCTGCTGTTTTGAAAAACCCCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGaa                                                                                                                  <  2:1014349/137‑1 (MQ=255)
                         tCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCAtt                                                                                                           >  2:92225/1‑140 (MQ=255)
                         tCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCAtt                                                                                                           >  1:1014349/1‑140 (MQ=255)
                         tCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCAtt                                                                                                           >  1:1578678/1‑140 (MQ=255)
                         ccGCCGTACGTTTCTGGCCGCTAAAAAAAGTTATCGACAACATCCACCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGAAAAACCCCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCAtt                                                                                                           <  2:195208/139‑1 (MQ=255)
                          cGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTa                                                                                                          >  1:821823/1‑140 (MQ=255)
                             cGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACg                                                                                                       <  1:92225/140‑1 (MQ=255)
                             cGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGCGGCGTTCCGGTGAATCGCATTAACg                                                                                                       >  1:1130559/1‑140 (MQ=255)
                                tcGTTTCTGGCCGCTAAAAAAATTTATCAACAACATTCACCGCCCCGCCGATCGACAAAAAATCCTGCTGGTTTGAAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTg                                                                                                    <  2:1556591/139‑1 (MQ=255)
                                ccttttCGGGACCCTGAAAAAAGTTTCCGAAAACATTCACCCCGCCGCCGATCAACAAAAAATCCTGCTGGTTTGAAAAACCCCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTg                                                                                                    <  2:869018/137‑1 (MQ=255)
                                 cttttCGGGCCGCTGAAAAAATTTACCAAAAACTTCCCCCGCGCCGCCGATCGAAAAAAAATCCTGTTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTgg                                                                                                   <  2:426155/138‑1 (MQ=255)
                                 cGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTgg                                                                                                   >  1:411531/1‑140 (MQ=255)
                                   tttCTGGCCGCTAAAAAAATTTATCAACAACTTCCCCCGCGCCGCCGATCGAAAAAAATTCCTGTGGGTTTGAAAAACCCCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCa                                                                                                 <  2:1578678/140‑1 (MQ=255)
                                   tttCGGGACGCTCAAAAATGTTCGCAAAAACTTTCAGCGCGCCGCCGACCGAAAAAAAATCCTCCTGGTTTGAAAACCCCCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCa                                                                                                 <  2:1130559/140‑1 (MQ=255)
                                              tGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTAc                                                                                      >  2:1658531/1‑140 (MQ=255)
                                              tGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTAc                                                                                      >  2:1077686/1‑140 (MQ=255)
                                                cAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGc                                                                                    >  2:1171426/1‑140 (MQ=255)
                                                        tATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCa                                                                            >  1:1059992/1‑140 (MQ=255)
                                                        tATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCa                                                                            >  2:1606229/1‑140 (MQ=255)
                                                          tCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCaaa                                                                          >  2:1297559/1‑140 (MQ=255)
                                                          tCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCaaa                                                                          >  1:412322/1‑140 (MQ=255)
                                                                 caTTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCAGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACaaa                                                                   >  2:1447206/1‑140 (MQ=255)
                                                                        cgcgCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCaa                                                            >  2:1621365/1‑140 (MQ=255)
                                                                           gccgccGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATACGCACTACGCCAATGGCAAACAACAAATCGCCAAAAc                                                         >  2:90267/1‑140 (MQ=255)
                                                                            ccgccgATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACg                                                        >  1:1418784/1‑140 (MQ=255)
                                                                            ccgccgATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACg                                                        >  2:346190/1‑140 (MQ=255)
                                                                              gccgATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGt                                                      <  1:1297559/140‑1 (MQ=255)
                                                                              gccgATCGACAAAAAATCCTGCTGTTTTGAAAACCACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGt                                                      <  2:880064/140‑1 (MQ=255)
                                                                              gccgATCGAAAAAAAATCTTGTTGGTTTGAAAAACCCCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGt                                                      <  1:1654622/140‑1 (MQ=255)
                                                                              cccgATCGACAGAAAATCCTGGTGGTTTGAAAAACCCCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGt                                                      <  2:735461/139‑1 (MQ=255)
                                                                                      aCAGAAAATCCTGCTGTTTTGAAAACCACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAGACAACAAATCGCCAAAACGGTTTCTGTgg                                              <  1:708126/140‑1 (MQ=255)
                                                                                        aaaaaaCCCTGCTGGTTTGAAAACCCCCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGAt                                            <  2:1418784/138‑1 (MQ=255)
                                                                                             acccTGCTGTTTTGAAAAACACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCggg                                       <  2:412322/138‑1 (MQ=255)
                                                                                             aTCCTGCTGTTTTGTAAAACCCCCGCCGATTTCCTGACGCTGGAGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCggg                                       <  1:103151/140‑1 (MQ=255)
                                                                                                ctgctgGTTTGAAAAACCCCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGa                                    <  1:692167/140‑1 (MQ=255)
                                                                                                     ggTTTGAAAACCCCCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCg                               <  2:1577026/140‑1 (MQ=255)
                                                                                                       tttGTAAACCCCCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCgca                             <  1:1125170/140‑1 (MQ=255)
                                                                                                        ttGAAAACCCCCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCgcag                            <  2:531973/140‑1 (MQ=255)
                                                                                                              aaCACCCGCCGATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTa                      <  1:780511/140‑1 (MQ=255)
                                                                                                                   ccgccgATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGAc                 >  1:1078903/1‑140 (MQ=255)
                                                                                                                   ccgccgATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGAc                 >  1:1414911/1‑140 (MQ=255)
                                                                                                                   ccgccgATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGAc                 <  1:483310/140‑1 (MQ=255)
                                                                                                                     cccgATTTCCTGACGCTGGCGAAAGGCGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCt               >  2:1559289/2‑140 (MQ=255)
                                                                                                                        gATTTCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGaa            <  2:903352/140‑1 (MQ=255)
                                                                                                                           ttCCTGACGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAAc         <  2:1078903/140‑1 (MQ=255)
                                                                                                                                 aCGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTgg   >  2:3718/1‑140 (MQ=255)
                                                                                                                                  cGCTGGCGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTggg  <  1:1171426/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
ATGGAAATGGTACTGGCAGAAGGGATCGCCGTACGTTTCTGGACGCTGCAAAAAGTTATCGACAACATTCATCGCGCCGCCGATCGACAGAAAATCCTGCTGGTTTGTAAAACACCCGCCGATTTCCTGACGCTGGTGAAAGGTGGCGTTCCGGTGAATCGCATTAACGTTGGCAATATGCACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGG  >  NC_000913/3280358‑3280627

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑