Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,074,454 A→C 9.2% F181V (TTC→GTC)  tktA ← transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,074,4540AC 9.2% 167.9 / 6.6 54F181V (TTC→GTC) tktAtransketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base A (22/27);  new base C (2/3);  total (24/30)
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.96e-01

CTTCTACTGCGCGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTT  >  CP009273/3074307‑3074603
                                                                                                                                                   |                                                                                                                                                     
cttctACTGCGCGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAAt                                                                                                                                                      <  1:148671/149‑1 (MQ=255)
 ttctACTGCGCGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATg                                                                                                                                                     >  2:66243/1‑149 (MQ=255)
   ctACTGCGCGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATg                                                                                                                                                     <  1:88820/147‑1 (MQ=255)
    tACTGCGCGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAAt                                                                                                                                                 >  2:132315/1‑150 (MQ=255)
         cgcgTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAg                                                                                                                                              <  2:73984/148‑1 (MQ=255)
         cgcgTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAg                                                                                                                                              <  1:63676/148‑1 (MQ=255)
         cgcgTTTGAGAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAg                                                                                                                                              <  1:71513/148‑1 (MQ=255)
          gcgTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGttt                                                                                                                                           <  1:105085/150‑1 (MQ=255)
                 aTAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGc                                                                                                                                    <  2:36539/150‑1 (MQ=255)
                    gATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTc                                                                                                                                 <  2:118257/150‑1 (MQ=255)
                     aTGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCtg                                                                                                                                  >  1:158364/1‑147 (MQ=255)
                        ccGCGTCCTGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGc                                                                                                                              <  1:15580/149‑1 (MQ=255)
                           cgTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTAc                                                                                                                          <  1:60006/150‑1 (MQ=255)
                               aTGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCg                                                                                                                       <  1:99261/149‑1 (MQ=255)
                                tGACCGTCAATGGCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGcc                                                                                                                     <  1:50132/150‑1 (MQ=255)
                                     gTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCagaga                                                                                                                <  1:27907/150‑1 (MQ=255)
                                     gTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCagaga                                                                                                                <  1:66243/150‑1 (MQ=255)
                                     gTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCaga                                                                                                                  >  2:185967/1‑148 (MQ=255)
                                           gTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCGAGCTTCAGCGTACCCGCCAGAGAGCAAAc                                                                                                          >  2:42808/1‑150 (MQ=255)
                                                gAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGt                                                                                                     >  1:32942/1‑150 (MQ=255)
                                                       gTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGAt                                                                                              <  2:180956/150‑1 (MQ=255)
                                                          ccAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGcc                                                                                           >  2:618/1‑150 (MQ=255)
                                                           cAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCt                                                                                          <  2:182159/150‑1 (MQ=255)
                                                             gCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTCCCCAGCTTCAGCGTCCCCGCCAGAGAGCCACCTTCGTGGGAGATGCCTTc                                                                                        <  2:46248/150‑1 (MQ=255)
                                                             gCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTc                                                                                        >  1:190504/1‑150 (MQ=255)
                                                               cGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCca                                                                                      >  2:77725/1‑150 (MQ=255)
                                                                 tAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCcatc                                                                                    <  1:124506/150‑1 (MQ=255)
                                                                 tAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCcat                                                                                     >  1:81343/1‑149 (MQ=255)
                                                                         gAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGcc                                                                            >  2:83282/1‑150 (MQ=255)
                                                                                aTTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGcc                                                                      >  1:50836/1‑149 (MQ=255)
                                                                                aTTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGcc                                                                      >  1:176723/1‑149 (MQ=255)
                                                                                    cGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATg                                                                  <  1:77725/149‑1 (MQ=255)
                                                                                           gtcgGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCg                                                            >  1:132846/1‑148 (MQ=255)
                                                                                           gtcgGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTa                                                          <  2:104852/150‑1 (MQ=255)
                                                                                                tGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGgtgt                                                     <  1:185967/150‑1 (MQ=255)
                                                                                                       gCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTc                                              >  1:152882/1‑150 (MQ=255)
                                                                                                        ccTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTc                                              >  2:49841/1‑149 (MQ=255)
                                                                                                            cAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGAc                                           <  2:81343/148‑1 (MQ=255)
                                                                                                            cAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACa                                          <  1:170380/149‑1 (MQ=255)
                                                                                                             aaCGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTa                                                    >  2:115508/1‑138 (MQ=255)
                                                                                                              aCGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAg                                                   <  1:115508/138‑1 (MQ=255)
                                                                                                              aCGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATg                                       <  1:83282/150‑1 (MQ=255)
                                                                                                                   aCCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGt                                   <  2:32942/149‑1 (MQ=255)
                                                                                                                       tCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGccc                              >  1:69388/1‑150 (MQ=255)
                                                                                                                       tCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGccc                              <  2:98425/150‑1 (MQ=255)
                                                                                                                        cGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCg                             >  2:60793/1‑150 (MQ=255)
                                                                                                                           tAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGAc                          >  2:151025/1‑150 (MQ=255)
                                                                                                                            aGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACg                         >  2:73896/1‑150 (MQ=255)
                                                                                                                                       tGTCGTCGTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTggc              >  1:75445/1‑148 (MQ=255)
                                                                                                                                        gtcgtcGTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTggc              <  1:194233/149‑3 (MQ=255)
                                                                                                                                        gTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTgcgc             >  2:157359/1‑150 (MQ=255)
                                                                                                                                         tcgtCGTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTggccc            <  2:165289/147‑5 (MQ=255)
                                                                                                                                           gtCGTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTggcccg           <  1:80840/148‑6 (MQ=255)
                                                                                                                                              gTAGACGGCAATCAGTTTGCCCAGTTTCATCGAGCCTGCCAATGAGCAAGCTTCATGCGAGATGCCTTCCATTAAGCAGCCGTCGCCCAGGAAGGCGTAGGTGTAGTGGTCGACGATTTCGTGGCCGGGGCGGTTGAATTGGCCCGCCAg       >  2:176703/1‑150 (MQ=255)
                                                                                                                                                   aTGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGttt  >  2:99941/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTTCTACTGCGCGTTTGATAGATGCCGCGTCATGACCGTCGATGTCGCGAATAACGTGCCAGCCGTAAGCTTCGAAACGCATTGCGGTGTCGTCGGTGAACCAGCCTTCAACGTGACCATCGATAGAAATACCGTTGTCATCGTAGAATGCAATCAGTTTACCCAGCTTCAGCGTACCCGCCAGAGAGCAAACTTCGTGGGAGATGCCTTCCATCATGCAGCCGTCGCCCATGAAGGCGTAGGTGTAGTGGTCGACAATGTCGTGGCCCGGACGGTTAAACTGCGCCGCCAGCGTTT  >  CP009273/3074307‑3074603

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