Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_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_pae_tpiA2,817,1700CA96.2% 131.3 / NA 52E112* (GAA→TAA) yqaBfructose‑1‑P and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (26/24);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGC  >  NC_000913_3_pae_tpiA/2817038‑2817314
                                                                                                                                    |                                                                                                                                                
ccagccGCTGCGGGAACAACAAAAATGTGCTTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACa                                                                                                                                  <  2:448105/144‑1 (MQ=255)
ccATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCCGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCcaca                                                                                                                                  >  1:710180/1‑149 (MQ=255)
ccATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCCAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACa                                                                                                                                  >  2:399124/1‑149 (MQ=255)
ccATACGATGCGCGCACCAAAACAATGCGTCTAGCGCGTGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACa                                                                                                                                  <  2:987468/149‑1 (MQ=255)
   tACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGTCGTCAAAGTAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGCTGGCGCTTTAACTCCCCGTTCCTCCAGcc                                                                                                                               >  2:263106/1‑149 (MQ=255)
       cTGCGCGCACAACAAAAATGTGTTTGGCGCGGGTTTATGTTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTg                                                                                                                           <  1:856007/149‑1 (MQ=255)
               acaacaAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGAGCACGGCGTAAAAATAATGGGGTAATCCTAGGTGCGCCAGCAATGCCTCACACATGGCGCTGTAAATCCCCGTTCCTAAAGCCATTTGGcgcagc                                                                                                                   >  2:1042051/1‑145 (MQ=255)
                    aaaaaTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATg                                                                                                              <  1:646854/149‑1 (MQ=255)
                        aTGTGTCTGGCGAGGGGTTATGGTGTTTGAAGTGATCGCCAGAGACGACGGCTCCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATTCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCaa                                                                                                          <  2:121265/149‑1 (MQ=255)
                             tCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCCGCGATGGCGCTTTAACTCCCCGTTCCTACCGCAATTGGCCGACGACCATGCCAACtttt                                                                                                     >  1:189614/1‑149 (MQ=255)
                               tGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTca                                                                                                   <  1:916126/149‑1 (MQ=255)
                                ggCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTcac                                                                                                  >  1:104744/1‑149 (MQ=255)
                                ggCGCGGGTTTATGGTGTTAGACGTAATCGGCAACGACGACGGCGTCAAAATAATGGCGTGATCCCAGGTGCGCCAGCGATTCCTCAGAGTTGGCCCCTTAACTCCCCGTTCCTACAGCCATTGAGCGCCCACCATCACAACTTTTcgt                                                                                                  >  2:642353/1‑147 (MQ=255)
                                   gcgGGTTTATGGTGTTTGACGTGATCGGCATCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATg                                                                                                              <  2:521795/134‑1 (MQ=255)
                                   gcgGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATg                                                                                                              >  1:521795/1‑134 (MQ=255)
                                            tGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAag                                                                                      <  1:328209/149‑1 (MQ=255)
                                            gggggtTGGAAGTTATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTCGCAGTGGAGCATGCCAACTTTTCACCACATCAACAag                                                                                      <  2:1009444/145‑1 (MQ=255)
                                             ggtggtTGGCGGGAGCGGCAGGGGCGACGGCGACAAAATAATGGCGTAATCCCAGGCGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAaga                                                                                     <  2:45011/149‑1 (MQ=255)
                                              gtgtTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAagag                                                                                    <  1:398583/149‑1 (MQ=255)
                                              gtgtTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAagag                                                                                    <  1:346283/149‑1 (MQ=255)
                                                gtTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTAACCACATCAACAAGAGGa                                                                                  >  1:243575/1‑149 (MQ=255)
                                                 tttGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGaa                                                                                 >  1:437670/1‑149 (MQ=255)
                                                    gACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTAACCCCATCAACAAGAGGAAGCg                                                                              >  1:572031/1‑149 (MQ=255)
                                                     aCGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCgg                                                                             >  2:169853/1‑149 (MQ=255)
                                                       gTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGTTGGCGCTTTACCTCCCCGTTCCTACAGCCATTGGGAGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGtt                                                                           >  1:306681/1‑149 (MQ=255)
                                                        ttctcagCAGCGGCGACCGTGTCAAAATAATGGCGTAAACCCAGAAGCGACAGCAATGCCTCAGCGATGGCGCTTTACCCCCCTGTTCCTACAGCCATTGGGCCACTACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTc                                                                          <  2:58517/143‑1 (MQ=255)
                                                          aTCGGCAGCGACTACGGCTGCAAAATTATTGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGTAAGCGGTTCGa                                                                        <  1:919821/149‑1 (MQ=255)
                                                           tCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCAAGCAATGCCTCAGCGATGGAGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAAGAAGCGGTTCGAc                                                                       >  2:466481/1‑149 (MQ=255)
                                                             ggCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGc                                                                     <  2:77313/149‑1 (MQ=255)
                                                                    acgacgGCCTCAAAATAACGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACAACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCcag                                                              <  2:48021/149‑1 (MQ=255)
                                                                            ttCAAAACAATAGCGTAATCCCACCTCCACCAACAATGCCTCAGCGATTGCGCTTTAACTCCCCGTACCAACAGCAATTGGGCGACGACCATGCCAACATTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATAc                                                      <  1:221107/148‑1 (MQ=255)
                                                                                  ataatGGGGTAATCCCAGGGGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTa                                                <  2:154674/149‑1 (MQ=255)
                                                                                      tGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGAGCGTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGATTCGACGCTATCCAGCAGAATACTTCTTACGGc                                            >  2:515610/1‑149 (MQ=255)
                                                                                       ggCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCt                                           >  1:782691/1‑149 (MQ=255)
                                                                                          gTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGAGGCTCGACACTATCCAGCAGCATACTTCTTACTGCTTCt                                        >  1:809892/1‑149 (MQ=255)
                                                                                                 cAGTTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCCGCATACTTCTTACTGCTTCTGtttttt                                 <  2:662087/149‑1 (MQ=255)
                                                                                                 cAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTATTGTTTATGttatta                                 >  1:104556/1‑148 (MQ=255)
                                                                                                 cAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGtttttt                                 <  1:842436/149‑1 (MQ=255)
                                                                                                  agttGCGCCAGCAATGACTCAGCGATGGCGCTTTAACCCCCCGTTCCTACAGCCAGTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTc                                <  1:676643/146‑1 (MQ=255)
                                                                                                  aGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGGTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTc                                >  1:885876/1‑149 (MQ=255)
                                                                                                    gTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGATCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTAGCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCAc                              <  1:747315/149‑1 (MQ=255)
                                                                                                         ccAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCAc                              <  1:540406/144‑1 (MQ=255)
                                                                                                         ccAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCAc                              >  2:540406/1‑144 (MQ=255)
                                                                                                         ccAGCAATGCCTCAGCGATGGCGCTTTAACACCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAAAGGAAGCGGTTCGACCCTATCCAGCAGCATACCTCTTACTGCTTCTGTTTTTTCACGCGCt                         >  2:430783/1‑149 (MQ=255)
                                                                                                                  ccTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGc                >  1:876233/1‑149 (MQ=255)
                                                                                                                    tCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCGTTGATGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCgg              >  1:773737/1‑149 (MQ=255)
                                                                                                                         ggtGGTGCTTTCACTACCCGGTCCCACAGCCATGGGGCGACGACCATGCCAACTGTTCACCACAGCAACAAGAGGAAGCGGTTCGACGGTATCCAGCAGCCTACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCga         <  2:492736/147‑1 (MQ=255)
                                                                                                                          aTGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACAGCTACTGTTTTTTCACCCGCTAACGCATGAGGGTCgag        >  1:552851/1‑149 (MQ=255)
                                                                                                                           tGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTACCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCgaga       <  2:184761/149‑1 (MQ=255)
                                                                                                                               gcTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCCTCCTTCTTACCGCTTCTGTTTTTTCTCGCGCTCCCGCCTCCGGGTCGCGTTCgc   >  2:547032/1‑148 (MQ=255)
                                                                                                                                cTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGc  <  2:154934/149‑1 (MQ=255)
                                                                                                                                cTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAAAGCATGCGGGTCGAGATCGGc  >  1:162581/1‑149 (MQ=255)
                                                                                                                                    |                                                                                                                                                
CCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGC  >  NC_000913_3_pae_tpiA/2817038‑2817314

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