Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 2,817,170 C→A 100% 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,1700CA100.0% 98.0 / ‑6.7 35E112* (GAA→TAA) yqaBfructose‑1‑P and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (15/18);  total (16/19)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCC  >  NC_000913_3_pae_tpiA/2817029‑2817315
                                                                                                                                             |                                                                                                                                                 
gTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTTAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTcccc                                                                                                                                            <  1:339291/149‑1 (MQ=255)
gTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCACCGATGGCGCTTTAACTcccc                                                                                                                                            >  2:205774/1‑149 (MQ=255)
     aCGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTcc                                                                                                                                       >  1:193425/1‑149 (MQ=255)
       gCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCACGTACGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTa                                                                                                                                     >  1:557098/1‑149 (MQ=255)
         ccATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACa                                                                                                                                   <  1:296772/149‑1 (MQ=255)
          taTACGCTGGGCGTACAACGAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTTAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGATTTAACTCCCCGTTCCTACAg                                                                                                                                  <  2:126884/148‑1 (MQ=255)
                 tGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTgg                                                                                                                           >  1:641994/1‑149 (MQ=255)
                 tGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTgg                                                                                                                           >  2:627680/1‑149 (MQ=255)
                             aaaaaTGTGTCTTGGGCGGGTGTTGGCTGTTTGACGTGATCGGCATCCACCACTGCGTCAAAATAATGGCGTAATACCAGGGGCGCCAGCAATGGCTCAGCGATGGCGATTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATg                                                                                                               <  2:116455/149‑1 (MQ=255)
                              aaaaTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGTCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGc                                                                                                              >  1:577909/1‑149 (MQ=255)
                                 aTGTGTCTGGGGCGGGTTTAGGGTGTGTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCaa                                                                                                           <  1:149597/149‑1 (MQ=255)
                                 aTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCaa                                                                                                           <  1:227650/149‑1 (MQ=255)
                                     gacTGGCGCGGGTTTATGGAGTTTGACGTGATCGGCAGCTACGACGGCGGCAAAATAATGCCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGATTTAACTCCCCGTTCCTACAGCCATTGGGCGTCGCCCATGCCAACttt                                                                                                       >  2:30527/3‑149 (MQ=255)
                                           cgcgGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTcacca                                                                                                 >  1:640360/1‑149 (MQ=255)
                                                       gtgtTTGCCGCGTCCGGCAGCGAGGCCGGCGTCAAAAGAATGCCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGTTGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAagag                                                                                     <  1:227723/149‑1 (MQ=255)
                                                        tgtTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAgg                                                                                    >  1:443287/1‑149 (MQ=255)
                                                           ttGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAg                                                                                 <  1:66410/149‑1 (MQ=255)
                                                              ctgTGATCGGGAGCGACGACGTCGCCAAAATAACGGCGTAATCCCAGGGGCGCCAGCCATGCCTCAGCGATGGCGCTGTCACCCCCCCTTCCTACAGCCATTGGGCGACGACCCTGCCAACTTTTCACCACATCAACAAGAGGAAGCgg                                                                              <  2:197992/147‑1 (MQ=255)
                                                               cGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGt                                                                             >  1:3896/1‑149 (MQ=255)
                                                                    tCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGAc                                                                        <  1:244259/149‑1 (MQ=255)
                                                                       gCAGCGACGACGGCGTCAAAATGATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCt                                                                     <  2:100395/149‑1 (MQ=255)
                                                                            gacgacGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTGACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCGAACTTTTCACGCCATCAACCAGAAGAAGCGGGTCGGCCCTAGCca                                                                >  1:41074/1‑149 (MQ=255)
                                                                                           ataatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTa                                                 <  1:535724/149‑1 (MQ=255)
                                                                                           ataatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTa                                                 <  1:416640/149‑1 (MQ=255)
                                                                                                       tCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGttt                                     <  1:272087/149‑1 (MQ=255)
                                                                                                         ccAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGttttt                                   >  1:143264/1‑149 (MQ=255)
                                                                                                                 gcCAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCCTACTTCTTACTGCTTCTGTTTTTTCACgcgc                           >  1:266532/1‑149 (MQ=255)
                                                                                                                 gacAGTAATGCCTCAGCGATGGCGCTTTAACTGCCCGTTCCTACAGCCATTGTGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAATCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACgcgc                           <  2:626119/147‑1 (MQ=255)
                                                                                                                   cAGCAATGCCTCAGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGGCGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTACCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTa                         >  1:142250/1‑149 (MQ=255)
                                                                                                                     gCAATGCCTCATCGATGGCGATTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAAc                       <  2:45039/149‑1 (MQ=255)
                                                                                                                            cTCAGCGATGGCGCTTTAACCCCCCGTTCCGACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTCCTTTCGTCATCTGGGTTTTCACGCGCTAACGCATGCg                <  2:526901/149‑1 (MQ=255)
                                                                                                                               aGCGATGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGt             <  2:613185/149‑1 (MQ=255)
                                                                                                                                   aTGGCGCTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCgag         >  1:396070/1‑149 (MQ=255)
                                                                                                                                         cTTTAACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAAGTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACCCATGCGGGTCGAGATCGGc   <  2:620046/149‑1 (MQ=255)
                                                                                                                                          tttAACTCCCCGTTCCTACAACCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCCGcc  >  2:73603/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
GTTGCACGCCCATACGCTGCGCGCACAACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAACTTTTCACCACATCAACAAGAGGAAGCGGTTCGACGCTATCCAGCAGCATACTTCTTACTGCTTCTGTTTTTTCACGCGCTAACGCATGCGGGTCGAGATCGGCC  >  NC_000913_3_pae_tpiA/2817029‑2817315

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: