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

CCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGC  >  NC_000913_3_pae_tpiA/2817038‑2817314
                                                                                                                                    |                                                                                                                                                
ccATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCACGTTCCTAcc                                                                                                                                  >  1:605053/1‑148 (MQ=255)
ccATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGGTTGACGTGAGCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGGGCCGGCGATGCCTCAGCGAGGGCGCTTTAGCTCCCCGTTCCTACa                                                                                                                                  >  2:470666/1‑149 (MQ=255)
             gcACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCTTACGCGCGGGCTCTTTAACTCCCCGTTCCTACAGCCATTGGGcgac                                                                                                                     >  1:638171/1‑149 (MQ=255)
                        aTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTc                                                                                                                                       <  2:445050/120‑1 (MQ=255)
                        aTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTc                                                                                                                                       >  1:445050/1‑120 (MQ=255)
                                   gcgGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTAAAAATAATGGAGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGCTGGCGCTTTGACTCCCCGTTCCTACAGCCATTGGGCGACGACCCTGCCAACTTTTcaccac                                                                                               >  1:111622/1‑149 (MQ=255)
                                      ggTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATc                                                                                            >  2:561967/1‑149 (MQ=255)
                                                     atgTGATTGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCTAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCgg                                                                             <  1:32247/147‑1 (MQ=255)
                                                          aTCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTc                                                                                                    <  1:179982/121‑1 (MQ=255)
                                                          aTCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTc                                                                                                    >  2:179982/1‑121 (MQ=255)
                                                                         ggCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGTAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAgcagca                                                         <  1:54515/149‑1 (MQ=255)
                                                                              cAAAATAATGGCGTAATCCCAGGCGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATActt                                                    <  1:254683/149‑1 (MQ=255)
                                                                               aaaaTAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGAc                                                                       >  2:80328/1‑129 (MQ=255)
                                                                               aaaaTAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGAc                                                                       <  1:80328/129‑1 (MQ=255)
                                                                                  ataatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTa                                                <  1:403859/149‑1 (MQ=255)
                                                                                  ataatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGCTCCGACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTa                                                <  2:629468/149‑1 (MQ=255)
                                                                                       ggCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGCCGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCAAGCAGCATACTTCTTACTGCt                                           >  1:82773/1‑149 (MQ=255)
                                                                                             aTCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGtt                                     >  1:348695/1‑149 (MQ=255)
                                                                                                   ggTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCa                               <  2:638171/149‑1 (MQ=255)
                                                                                                          cAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTa                        >  2:522050/1‑149 (MQ=255)
                                                                                                              aaTGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATCCTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGc                    >  1:276767/1‑149 (MQ=255)
                                                                                                              aaTGCCTCAGCGATGACGCTTTAACTTCCCGCACCTCCAGCCACTGGGCGACGACCATGACACCTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTCCTGTTTTTTCCCGCGCTACCGc                    >  2:520798/1‑149 (MQ=255)
                                                                                                                    tCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCgg              <  1:179117/149‑1 (MQ=255)
                                                                                                                               gcTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGCCGAGACCgg   >  1:178273/1‑149 (MQ=255)
                                                                                                                                ttttAACCCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGc  <  1:464455/148‑1 (MQ=255)
                                                                                                                                cTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCgcc  >  1:340107/1‑147 (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 < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: