Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_tpiA 2,817,170 C→A E112* (GAA→TAA)  yqaB ← fructose‑1‑P and 6‑phosphogluconate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_tpiA2,817,1700CA96.1% 143.4 / NA 51E112* (GAA→TAA) yqaBfructose‑1‑P and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (18/31);  total (19/32)
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.98e-01

GTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATC  >  NC_000913_3_vch_tpiA/2817025‑2817311
                                                                                                                                                  |                                                                                                                                             
gTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACt                                                                                                                                            <  1:651782/149‑1 (MQ=255)
 tCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCCATGCCTCAGCGATGGCGCTTTAACTc                                                                                                                                           >  2:536493/1‑149 (MQ=255)
  cGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTcc                                                                                                                                          <  1:718429/149‑1 (MQ=255)
   ggTTGCACGCCCATGCGCTGCGCGCACAACAAAAATGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTccc                                                                                                                                         <  2:177292/149‑1 (MQ=255)
    gTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTcccc                                                                                                                                        <  1:640992/149‑1 (MQ=255)
      tGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGt                                                                                                                                      <  1:258747/149‑1 (MQ=255)
         aCGCCCATACGCTGCGCGCACAACAAAAATGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAACATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTcc                                                                                                                                   >  1:698646/1‑149 (MQ=255)
             ccATACGCTGCGCGCACAACAAAAATGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACa                                                                                                                               >  2:95845/1‑149 (MQ=255)
             ccATACCCTGCGCGCACAACAAAAAAGTGTCT‑GGCTCTGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCCTCAAAATAATGGCGTAATCCCAGGGGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACa                                                                                                                               <  1:341730/149‑1 (MQ=255)
                tACGCTGCGCGCACAACAAAAATGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAg                                                                                                                              >  2:454550/1‑147 (MQ=255)
                tACGCTGCGCGCACAACAAAAATGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAg                                                                                                                              <  1:454550/147‑1 (MQ=255)
                          gcACAACAAAACTG‑GTCTGGGCCGGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGcgac                                                                                                                  <  1:762046/149‑1 (MQ=255)
                              aacaaAAATGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGccc                                                                                                              >  2:74752/1‑149 (MQ=255)
                                caaAAATGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGCCGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCCCTTTAACTCCCCCTTCCTACAGCCATTGGGCCCCGACCac                                                                                                            >  2:439244/1‑148 (MQ=255)
                                 aaaaaTGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCCTTCCTACAGCCATTGGGCGACGACCATg                                                                                                           <  1:112972/149‑1 (MQ=255)
                                         tgttttgcgcGGGTTTATGGTGTTTGCCGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACtt                                                                                                    <  1:279789/143‑1 (MQ=255)
                                          tCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACtttt                                                                                                  >  1:730811/1‑149 (MQ=255)
                                            t‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTca                                                                                                >  2:314807/1‑149 (MQ=255)
                                            t‑GGCGCGGGTTTATGGTGTTTGACCTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTca                                                                                                <  1:319110/149‑1 (MQ=255)
                                                     gTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATca                                                                                        >  2:73661/1‑149 (MQ=255)
                                                       ttATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATcaac                                                                                      <  1:73661/149‑1 (MQ=255)
                                                       tgaTGGTGGTTGACGTGAGCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATCCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATcaac                                                                                      <  1:132662/147‑1 (MQ=255)
                                                         aTGGGGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCATTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATcaacaa                                                                                    <  1:275025/149‑1 (MQ=255)
                                                           ggTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAaga                                                                                  <  1:47153/149‑1 (MQ=255)
                                                           ggTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAaga                                                                                  <  2:34469/149‑1 (MQ=255)
                                                               tttGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGaa                                                                              >  1:279013/1‑149 (MQ=255)
                                                                ttGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAg                                                                             <  2:341136/149‑1 (MQ=255)
                                                                  gACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCg                                                                           <  2:765621/149‑1 (MQ=255)
                                                                     gTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGtt                                                                        >  1:403443/1‑149 (MQ=255)
                                                                                          gTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCcagca                                                         <  2:201184/143‑1 (MQ=255)
                                                                                          gTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCcagca                                                         >  1:201184/1‑143 (MQ=255)
                                                                                          gTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTAGACGCTATCCAGCAGCATAc                                                   >  1:199510/1‑149 (MQ=255)
                                                                                              aaaTAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATActtct                                               >  2:529897/1‑149 (MQ=255)
                                                                                                aaaatGGCGGAATCCCAGGTGCGCCAGCAACTCCGCAGCGATGGCGCTTTAACTCCCCTTGCCTACAGCCATTGGGCGACGACCATGCCAACGTTGCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTa                                             <  2:59002/147‑1 (MQ=255)
                                                                                                 taatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGAc                                                                    >  2:183695/1‑125 (MQ=255)
                                                                                                 taatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGAc                                                                    <  1:183695/125‑1 (MQ=255)
                                                                                                        gTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGATTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCt                                     >  2:293755/1‑149 (MQ=255)
                                                                                                           aTCGCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCGCGTTTCTACAGCCATTGGGCGACGACCAGGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGtt                                  <  2:164197/149‑1 (MQ=255)
                                                                                                                  gggcgcCAGCAATGCCTCAGCGATGGCGTTTTAACCCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCAc                           <  2:383416/147‑1 (MQ=255)
                                                                                                                  gctcgcCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCAc                           <  2:65094/146‑1 (MQ=255)
                                                                                                                         aGCAATGCCTCAGAGATGGCGCTTTAACGCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTaa                    <  2:379105/149‑1 (MQ=255)
                                                                                                                             aTGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCa                <  1:758994/149‑1 (MQ=255)
                                                                                                                               gCCCCAGCGATGGCGCTTTAACTCCCCGTTTCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATg              <  2:503530/149‑1 (MQ=255)
                                                                                                                                 cTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCg            <  1:625832/149‑1 (MQ=255)
                                                                                                                                 cTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCg            >  2:75201/1‑149 (MQ=255)
                                                                                                                                  tCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCgg           >  1:750393/1‑149 (MQ=255)
                                                                                                                                       gATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCAcg                          <  2:424936/129‑1 (MQ=255)
                                                                                                                                       gATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCAcg                          >  1:424936/1‑129 (MQ=255)
                                                                                                                                          ggCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGAt   <  1:5787/149‑1 (MQ=255)
                                                                                                                                          ggCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACAACAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGAt   <  1:406364/149‑1 (MQ=255)
                                                                                                                                           gcgcTTTAACTCCCCGTTCCTACAGCCAGTGGTCGACTACCATCACAACTATTCACCACATCACCAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATc  <  1:498256/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
GTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCT‑GGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATC  >  NC_000913_3_vch_tpiA/2817025‑2817311

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