Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,460 A→G 8.0% A598A (GCT→GCC)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,4600AG 8.0% 209.4 / 7.0 75A598A (GCT→GCC) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (33/36);  new base G (3/3);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTC  >  minE/608398‑608528
                                                              |                                                                    
ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACg                                                              >  1:4274323/1‑71 (MQ=255)
 tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGt                                                                  <  1:3955736/66‑1 (MQ=255)
  aCCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                                                             >  1:115970/1‑70 (MQ=255)
      gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                                                             >  1:1062633/1‑66 (MQ=255)
      gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                                                             >  1:3634530/1‑66 (MQ=255)
      gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                                                        >  1:4710746/1‑71 (MQ=255)
      gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                                                        >  1:2879597/1‑71 (MQ=255)
      gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt                                                        >  1:2541054/1‑71 (MQ=255)
      gTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGgt                                                        >  1:2102100/1‑71 (MQ=255)
       ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCAc                                                               <  1:4317281/63‑1 (MQ=255)
         ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgtgt                                                      <  1:4831545/70‑1 (MQ=255)
         ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGtgtg                                                     >  1:1910693/1‑71 (MQ=255)
         ttCTTACCATACACAACGTTAGAAGAGATGGTCAGAACAGACTGAGTCGGGATAg                                                                     >  1:3934166/1‑55 (MQ=255)
              aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCa                                                   <  1:3568046/68‑1 (MQ=255)
              aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt                                                 >  1:1920407/1‑70 (MQ=255)
              aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt                                                 <  1:4449503/70‑1 (MQ=255)
               ccATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATATCGTCACGGTAGGTGTGCAGttt                                               <  1:82175/71‑1 (MQ=255)
               ccATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAg                                                  <  1:5049511/68‑1 (MQ=255)
                 aTACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGt                                                            <  1:3255532/56‑1 (MQ=255)
                 aTACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCt                                             >  1:1650878/1‑71 (MQ=255)
                    caccacGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAgg                                                         >  1:471017/1‑56 (MQ=255)
                    cacaACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCa                                                                >  1:5395537/1‑49 (MQ=255)
                    cacaACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGaa                                          >  1:2005261/1‑71 (MQ=255)
                       aaCGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                       >  1:932886/1‑71 (MQ=255)
                       aaCGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                       >  1:1961911/1‑71 (MQ=255)
                       aaCGTTAGAAGTGATGGTCAGAACAGACTCAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                       <  1:1590691/71‑1 (MQ=255)
                        aCGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                       >  1:5092046/1‑70 (MQ=255)
                         cGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGGTTCTGAATTttc                                     <  1:4513366/71‑1 (MQ=255)
                           ttAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                       <  1:1665726/67‑1 (MQ=255)
                           ttAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                       <  1:478686/67‑1 (MQ=255)
                           ttAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTttctt                                   <  1:3191103/71‑1 (MQ=255)
                             aGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCa                                 >  1:4144598/1‑71 (MQ=255)
                             aGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCa                                 >  1:1923857/1‑71 (MQ=255)
                             aGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCa                                 >  1:2664794/1‑71 (MQ=255)
                             aGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCa                                 >  1:3691980/1‑71 (MQ=255)
                              gAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCAt                                <  1:3737269/71‑1 (MQ=255)
                              gAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTGAATTTTCTTCAt                                <  1:903448/71‑1 (MQ=255)
                                  tGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTttc                                     <  1:1842762/62‑1 (MQ=255)
                                  tGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGa                              <  1:5430284/69‑1 (MQ=255)
                                  tGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGa                              <  1:2555051/69‑1 (MQ=255)
                                   gATGGTCAGAACAGACTGAGTCGGGATAGAGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                           >  1:1607965/1‑71 (MQ=255)
                                    aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAg                                                          >  1:1619658/1‑39 (MQ=255)
                                    aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                           <  1:5169601/70‑1 (MQ=255)
                                    aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                           <  1:2557183/70‑1 (MQ=255)
                                    aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACg                          <  1:411506/71‑1 (MQ=255)
                                    aTGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACg                          >  1:4177682/1‑71 (MQ=255)
                                       gTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGGGTGCAGTTTCTGAATTTTCTTCATGAAACGtt                        <  1:470384/70‑1 (MQ=255)
                                        tCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCt                      <  1:836843/71‑1 (MQ=255)
                                         cAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                       >  1:3275993/1‑53 (MQ=255)
                                         cAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                       >  1:512016/1‑53 (MQ=255)
                                          aGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTa                     >  1:457559/1‑70 (MQ=255)
                                           gAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTAcc                   <  1:4320318/71‑1 (MQ=255)
                                            aaCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                           <  1:727378/62‑1 (MQ=255)
                                            aaCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                           <  1:3001829/62‑1 (MQ=255)
                                            aaCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                           <  1:4152402/62‑1 (MQ=255)
                                             aCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAg                                                  >  1:5564359/1‑38 (MQ=255)
                                             aCAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTg                                            >  1:2980204/1‑44 (MQ=255)
                                              cAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAgg                >  1:2523699/1‑71 (MQ=255)
                                                 aCTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAAc                           <  1:1551903/57‑1 (MQ=255)
                                                 aCTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCa             <  1:5237242/71‑1 (MQ=255)
                                                 aCTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCa             <  1:712001/71‑1 (MQ=255)
                                                 aCTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGATCTACCAGGTc              >  1:2196096/1‑70 (MQ=255)
                                                  cTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                       <  1:4296603/44‑1 (MQ=255)
                                                  cTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                       <  1:4267765/44‑1 (MQ=255)
                                                   tGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                       <  1:1389460/43‑1 (MQ=255)
                                                       tCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTc              <  1:2804514/64‑1 (MQ=255)
                                                        cGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAAttt                                       <  1:5526022/38‑1 (MQ=255)
                                                        cGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGt                         <  1:1309733/52‑1 (MQ=255)
                                                        cGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGt                         <  1:3884306/52‑1 (MQ=255)
                                                         gggATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTGTACCAGGTCAACAGCCAg     >  1:3342243/1‑71 (MQ=255)
                                                         gggATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAg     >  1:3315199/1‑71 (MQ=255)
                                                         gggATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAg     >  1:4812152/1‑71 (MQ=255)
                                                         gggATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAg     <  1:4270063/71‑1 (MQ=255)
                                                         gggATAGCGTCAAGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAg     >  1:1855319/1‑71 (MQ=255)
                                                            aTAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTc  <  1:5285724/71‑1 (MQ=255)
                                                              aGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACg                          >  1:4113556/1‑45 (MQ=255)
                                                              |                                                                    
TTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCTGAATTTTCTTCATGAAACGTTCTACCAGGTCAACAGCCAGGTC  >  minE/608398‑608528

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