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,1700CA97.2% 94.0 / NA 36E112* (GAA→TAA) yqaBfructose‑1‑P and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (14/21);  minor base G (1/0);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCC  >  NC_000913_3_pae_tpiA/2817024‑2817315
                                                                                                                                                  |                                                                                                                                                 
cGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAAc                                                                                                                                                 <  2:577044/149‑1 (MQ=255)
   ctgTTGCACGTCCATACGCTGCGCGCACAACACAAACGTGTCTGGCGCGGGTTTATTGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTcc                                                                                                                                              <  1:150115/147‑1 (MQ=255)
      ttGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCg                                                                                                                                           >  1:458932/1‑149 (MQ=255)
         cACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTc                                                                                                                                        <  2:504656/149‑1 (MQ=255)
         cACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTc                                                                                                                                        <  2:549855/149‑1 (MQ=255)
            gCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTa                                                                                                                                     <  2:569719/149‑1 (MQ=255)
            gCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGAGTTTATGGTGTTTGACGGGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGGACTTTAACTCCCCGTTCCTa                                                                                                                                     >  2:255807/1‑149 (MQ=255)
                aTACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGc                                                                                                                                 <  2:118386/149‑1 (MQ=255)
                    gttGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCTAGGTGCGCACGCAATGGCTCAGCGCTAGCGATTTGACTCCCGGATCACACAGcccac                                                                                                                             >  2:507810/3‑146 (MQ=255)
                     cTGCGCGCACAACAAAAATGTGTCTGGGGCGGGTTTATGGTGTTTGAAGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTg                                                                                                                            <  2:553257/149‑1 (MQ=255)
                            cacaACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGcgacg                                                                                                                     >  1:251717/1‑149 (MQ=255)
                                      agctgtCTGGCGCGGGTTTATGGTGGTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCaa                                                                                                           <  1:561780/146‑1 (MQ=255)
                                          gtCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAAc                                                                                                                                                 <  2:280410/107‑1 (MQ=255)
                                          gtCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAAc                                                                                                                                                 >  1:280410/1‑107 (MQ=255)
                                             tGGCGCGGGTTTAAGGTGTTTGACGGGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTca                                                                                                    <  2:446762/149‑1 (MQ=255)
                                                        tATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGCCCATGCCACCTGTTCACCACATcaaca                                                                                         >  2:435750/1‑149 (MQ=255)
                                                             tgtTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAgg                                                                                    <  2:278571/149‑1 (MQ=255)
                                                                    cGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCACCAAGAGGAAGCGGt                                                                             >  1:403716/1‑149 (MQ=255)
                                                                     gTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGtt                                                                            >  1:456556/1‑149 (MQ=255)
                                                                     gTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTGCCTACAGCCATTGTGCGACGCCCATGCCAACTTCTCACCACATCAACAAGAGGTAACGGtg                                                                            >  2:666914/1‑148 (MQ=255)
                                                                         tCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGAGGGCGCTTTAACTCCCCCTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGAc                                                                        <  1:631198/149‑1 (MQ=255)
                                                                                    gacgGCGTCAAACTAATGGCGAAATCCCAGGTGAGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCCACAAGCCGAAGCGGTTCGACGCTATCcagca                                                             >  2:118789/1‑149 (MQ=255)
                                                                                             aaaaTAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCCACGCTATCCAGCAGCATActtc                                                    <  1:79076/149‑1 (MQ=255)
                                                                                                ataatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTa                                                 <  2:407582/149‑1 (MQ=255)
                                                                                                ataatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTa                                                 <  1:191523/149‑1 (MQ=255)
                                                                                                 taatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCCGCATACTTCTTAc                                                >  1:450470/1‑149 (MQ=255)
                                                                                                                      gcCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACgcgc                           <  1:32872/149‑1 (MQ=255)
                                                                                                                            agtGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGc                     >  2:362632/3‑149 (MQ=255)
                                                                                                                            aaTGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACACCCATTGGGCGACGACCACGCCAACTTTTCACCACAACAACAAGAGGAAGCGATTCAACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAATGc                     >  1:373652/1‑149 (MQ=255)
                                                                                                                             aTGCATCAGGGATGGCGCGTTAATTCCCCGTTCCTACAGCCATTGGGCGAGTACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCa                    <  2:199327/149‑1 (MQ=255)
                                                                                                                                gcTCAGCGATGGCGCTTGAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGc                 <  1:675791/148‑1 (MQ=255)
                                                                                                                                ccTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGttt                                     >  1:648266/1‑129 (MQ=255)
                                                                                                                                ccTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGttt                                     <  2:648266/129‑1 (MQ=255)
                                                                                                                                ccTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGc                 <  2:651333/149‑1 (MQ=255)
                                                                                                                                ccTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGc                 <  2:651334/149‑1 (MQ=255)
                                                                                                                                               tttAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGcc  >  2:570561/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGTCGGTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCC  >  NC_000913_3_pae_tpiA/2817024‑2817315

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: