Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,856,333 T→G 14.6% E199A (GAG→GCG) ‡ ydjH ← putative sugar kinase YdjH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,856,3330TG14.6% 107.5 / 16.3 55E199A (GAG→GCG) ‡ydjHputative sugar kinase YdjH
Reads supporting (aligned to +/- strand):  ref base T (26/21);  new base G (5/3);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.01e-02

CGCCGGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCG  >  NC_000913/1856191‑1856474
                                                                                                                                              |                                                                                                                                             
cgccgGCCCCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAAGCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCt                                                                                                                                               <  2:570457/143‑1 (MQ=255)
        ccTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTATGCTTCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCaaaa                                                                                                                                       <  1:1096113/143‑1 (MQ=255)
           tCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTa                                                                                                                                    <  1:803009/143‑1 (MQ=255)
            cATCGTCATGCCACCCCGCTCGATACAGCAGCCTGCATTACCCGTTTTAACCACCCCCGTTTTTACGCCGTACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAg                                                                                                                                   <  1:369761/143‑1 (MQ=255)
             aTCGTCATGTCACCACGCTCGATCAAGCACCCCACTCAATCTTGATGTGGCGCCACCGTTTTTACGCCGCTCGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTact                                                                                                                                  <  1:248658/143‑3 (MQ=255)
                gTCATGTCACCGCGCTTGATAAAGCAGCCGGCTTTACCCGTTTTCATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGTTTCCTAATTTTGCCGAGGCAAAATTactctc                                                                                                                               <  2:739094/143‑6 (MQ=255)
                 tCATGTCACCACGCTTGATAAAGCAGCCGTCTCGACCCGTTTTAATGACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAAc                                                                                                                              <  2:662274/143‑1 (MQ=255)
                       caccacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGCGTCATTTTGCCTCGGCAACATTAGGAAGCAGATaa                                                                                                                        >  1:623038/1‑143 (MQ=255)
                       caccacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATaa                                                                                                                        >  1:301954/1‑143 (MQ=255)
                          cgcGCTGGATATTGCAGCCGTCGTGACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCGATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCg                                                                                                                     <  2:430892/141‑1 (MQ=255)
                          cacGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCg                                                                                                                     >  1:890005/1‑143 (MQ=255)
                           acGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGa                                                                                                                    <  2:37662/143‑1 (MQ=255)
                           acGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTCCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAACTAGTAAACAGAtgaaaga                                                                                                                    >  2:79703/1‑137 (MQ=255)
                            cGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGCGCGCCAGGTTTCCTCGGCAAAATTAGGCAACACACAATCGAc                                                                                                                   >  2:486875/1‑143 (MQ=255)
                             gCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACa                                                                                                                  <  1:1030198/143‑1 (MQ=255)
                                       gcagcCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGGGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATg                                                                                                        <  1:1052352/143‑1 (MQ=255)
                                       gcagcCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAACATTAGGAAACAGATAATCGACACACCTCAATg                                                                                                        >  1:757552/1‑143 (MQ=255)
                                       gcagcCGTCTTTACCCGTTTTAAGCACCACCGTTTGTACGCCGCACGCAAGAAAGCAGTCAGCAATGTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATg                                                                                                        <  1:182531/143‑1 (MQ=255)
                                                  tACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAAATCAATGCTTGGCAAata                                                                                             >  1:718705/1‑143 (MQ=255)
                                                     ccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCCGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatca                                                                                          <  1:417819/143‑1 (MQ=255)
                                                      cGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAAGTTTGCCTCGTCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCTCAAATatcat                                                                                         >  2:302881/1‑143 (MQ=255)
                                                       gTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTa                                                                                                                                    <  2:858796/99‑1 (MQ=255)
                                                       gTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTa                                                                                                                                    >  1:858796/1‑99 (MQ=255)
                                                        ttttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCCAATATCATcc                                                                                       >  1:468041/1‑143 (MQ=255)
                                                                accaccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAACATTATGAAACAGATAATCGACATAACTCAATGCATCGCAAATATCAACCAGCGTttc                                                                               >  2:662469/1‑143 (MQ=255)
                                                                   accGGTGTTACGCCGCCCTCAAGAAATCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcatt                                                                            <  2:347000/143‑1 (MQ=255)
                                                                      gTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCaa                                                                         <  1:72193/143‑1 (MQ=255)
                                                                           tACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCCAAATTAGGAAACAGTTCATCGACACAACGCAACGCTTCGCAAATATCATCCAGCGTTTCCTTCAACcgcg                                                                    >  1:466229/1‑143 (MQ=255)
                                                                            aCGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCgg                                                                   >  2:97621/1‑143 (MQ=255)
                                                                                  cacgcaATAAAGCAGTCCGCAATTTCATCCAGTTGCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAACCAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGAt                                                             <  2:413791/143‑1 (MQ=255)
                                                                                      caAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGTAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCATTTCATTCACCCGCGGgtttctcttt                                                         >  1:773483/1‑136 (MQ=255)
                                                                                       aaGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCCACCGCGGTTTGATCagat                                                        >  1:865369/1‑143 (MQ=255)
                                                                                         gAAAGCAGTCAGCAACTTCATCCAGTGTCTCTTTCCCGCTGAGTAATTTTGCCTCCGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCg                                                      <  1:835582/143‑1 (MQ=255)
                                                                                            aGCAGTCTGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGca                                                   <  1:221209/143‑1 (MQ=255)
                                                                                             gCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGcac                                                  >  2:448989/1‑143 (MQ=255)
                                                                                                 tCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACCTAACTCAATGCTTCGCAAATCCCATCCCGCGCTTCATTGAACCGCGGTTTGCTCAAAGCGACACAGac                                              >  2:437997/1‑142 (MQ=255)
                                                                                                 tCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATCTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGATGCATTCAACCGCGGGTTGATGATATCGGCACAGAt                                              >  2:362152/1‑143 (MQ=255)
                                                                                                       aattCAGCCAGTGGCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCACTGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAatcat                                        <  2:232674/141‑1 (MQ=255)
                                                                                                          tCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTg                                     <  1:730812/143‑1 (MQ=255)
                                                                                                           cATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGc                                    <  1:476546/143‑1 (MQ=255)
                                                                                                                   gtCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCAGTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCttt                            >  1:615028/1‑143 (MQ=255)
                                                                                                                    tctcTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGAACAGATAATCATCTGCCGGGCTTTa                           >  1:681139/1‑143 (MQ=255)
                                                                                                                     ctctTTCCCGGTGAGTACTATTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAg                          <  2:1011310/143‑1 (MQ=255)
                                                                                                                     ctctTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGGCGGGCTGTAg                          >  1:32638/1‑143 (MQ=255)
                                                                                                                      tctTTCCCGGTGAGTAATTTTGCCTCGGCGAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCACAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCAAAGATAATCATCTGCCGTGGTTTAgg                         >  2:45515/1‑142 (MQ=255)
                                                                                                                        tttCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCCGACGGGCTTTAGCtt                       >  2:1019282/1‑143 (MQ=255)
                                                                                                                            ccGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGt                   >  2:72041/1‑143 (MQ=255)
                                                                                                                                   gTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCCATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCCTCTGCCGGGCCTACGCTTGCGTAAAACtt            >  2:113076/1‑143 (MQ=255)
                                                                                                                                   ctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTGTCGCTTGAGTAAAGAtt            >  1:32395/2‑143 (MQ=255)
                                                                                                                                   ctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGAGGTTTGAACATATCGGCAGAGACATTCACCGCCCGGACTTAAGCCTGCGCAACAAta            >  2:532379/2‑142 (MQ=255)
                                                                                                                                    tAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCCTTCAACCGCGGTTTGCTCCCATCGTTACAGATAAGCATCTGCCGGGCTTTAGCTTGCGTCAAAAttt           >  2:658693/1‑143 (MQ=255)
                                                                                                                                    tAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATATGCTTGGCTTTAGCTTGCGTAAAACttt           >  1:689934/1‑143 (MQ=255)
                                                                                                                                       ttttGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTg        >  1:81687/1‑143 (MQ=255)
                                                                                                                                       ttttGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGCCCCAGATAATCATCTCCCGGGCTTTATCTTGCGTAAAAATTTcat        >  2:861084/1‑141 (MQ=255)
                                                                                                                                            ccTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAgc   <  1:191681/143‑1 (MQ=255)
                                                                                                                                             cTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAgcg  <  1:97621/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
CGCCGGCACCTTCATCGTCATGTCACCACGCTTGATAAAGCAGCCGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAATCATCTGCCGGGCTTTAGCTTGCGTAAAAATTTCTGTTAGCG  >  NC_000913/1856191‑1856474

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