Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,959,701 A→G 58.3% L124L (TTG→CTG)  fbaA ← fructose‑bisphosphate aldolase, class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,959,7010AG58.3% 24.0 / 47.1 67L124L (TTG→CTG) fbaAfructose‑bisphosphate aldolase, class II
Reads supporting (aligned to +/- strand):  ref base A (4/24);  new base G (21/18);  total (25/42)
Fisher's exact test for biased strand distribution p-value = 1.01e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

CAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGG  >  minE/1959633‑1959761
                                                                    |                                                            
cAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCaa                                                              <  1:1074329/69‑1 (MQ=255)
cAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCaa                                                              <  1:346601/69‑1 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:923993/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:865327/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:751493/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:67051/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:3421315/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:3330456/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:323625/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:2884708/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:2200534/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:2151425/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:2042699/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:1957190/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:1710557/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:1685713/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:157061/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:148948/1‑69 (MQ=255)
   gTCAATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCAGCGAAGTGTTTTTCACCCGCGTCcagcag                                                           >  1:1110419/1‑69 (MQ=255)
                         aGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGcagcag                                      >  1:2014005/1‑68 (MQ=255)
                         aGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGcagcag                                      >  1:213137/1‑68 (MQ=255)
                         aGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGt                                     <  1:3166381/69‑1 (MQ=255)
                         aGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGt                                     >  1:3210599/1‑69 (MQ=255)
                         aGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGt                                     <  1:1218713/69‑1 (MQ=255)
                         aGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGt                                     <  1:554378/69‑1 (MQ=255)
                         aGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGt                                     <  1:221088/69‑1 (MQ=255)
                         aGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGt                                     <  1:1500760/69‑1 (MQ=255)
                         aGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGt                                     <  1:1948541/69‑1 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGAcc                                                        >  1:545768/1‑38 (MQ=38)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATc                                                 >  1:275398/1‑45 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCgcgc                           >  1:661879/1‑67 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:293131/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:2962035/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:1461455/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:2971926/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:316625/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:2602726/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:1439370/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:2665900/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:3306210/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:3312311/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:2553749/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:792299/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:818848/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:874143/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:394451/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:1547576/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:1346895/1‑69 (MQ=255)
                                     gTAGCAGCGAAGTGTTTTTCACCCGCGTCCAGCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAg                         >  1:367294/1‑69 (MQ=255)
                                       aGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTg                       <  1:68388/69‑1 (MQ=255)
                                       aGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTg                       <  1:542291/69‑1 (MQ=255)
                                        gCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTg                       <  1:2724682/68‑1 (MQ=255)
                                          tGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTc                    <  1:1704185/69‑1 (MQ=255)
                                          tGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTc                    <  1:1427371/69‑1 (MQ=255)
                                          tGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTc                    <  1:518294/69‑1 (MQ=255)
                                          tGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTc                    >  1:3170259/1‑69 (MQ=255)
                                          tGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTc                    <  1:2977282/69‑1 (MQ=255)
                                          tGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTc                    <  1:2129388/69‑1 (MQ=255)
                                              aaGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAg                  <  1:2608092/67‑1 (MQ=255)
                                              aaGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAg                  <  1:2088061/67‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcaggag                                      <  1:1266574/40‑1 (MQ=38)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:637247/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:742353/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:499662/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:914374/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:1009329/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:1626979/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:3389813/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:3209256/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:3138557/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:310459/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:2860359/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:2700942/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:2182773/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:2657840/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:2521361/40‑1 (MQ=255)
                                                     tttCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGa         <  1:1419358/69‑1 (MQ=255)
                                                     tttCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGa         <  1:2113367/69‑1 (MQ=255)
                                                     tttCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGa         <  1:844849/69‑1 (MQ=255)
                                                     tttCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGa         <  1:348889/69‑1 (MQ=255)
                                                      ttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:782191/39‑1 (MQ=255)
                                                      ttCACCCGCGTCCAGCAGACCGTCGATCCACGGcagcag                                      <  1:2168416/39‑1 (MQ=255)
                                                      ttCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGa         <  1:1806081/68‑1 (MQ=255)
                                                            cgcgTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCgg  <  1:2499682/69‑1 (MQ=255)
                                                            cgcgTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCgg  <  1:2662456/69‑1 (MQ=255)
                                                            cgcgTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCgg  <  1:2764623/69‑1 (MQ=255)
                                                                    aCAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTg                       <  1:2761387/40‑1 (MQ=255)
                                                                    |                                                            
CAGGTCGATCATGTGAGAAGAGAACAGCGGCTTACCGGTAGCTGCGAAGTGTTTTTCACCCGCGTCCAACAGACCGTCGATCCACGGCAGCAGTTTCTTCGCGCAGTGGTCAGTGTGCAGGATAACCGG  >  minE/1959633‑1959761

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