Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,488 C→G 27.0% V527V (GTG→GTC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,4880CG27.0% 102.5 / 40.1 63V527V (GTG→GTCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (27/19);  new base G (12/5);  total (39/24)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

TTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATT  >  minE/1014421‑1014558
                                                                   |                                                                      
tttATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAATAAGTCATTTCCTTTATCATCcac                                                                      >  1:4075677/1‑70 (MQ=255)
 ttATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCcac                                                                      >  1:1964758/1‑69 (MQ=255)
 ttATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCcac                                                                      <  1:4162374/69‑1 (MQ=255)
  tATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCAGAAg                                                                   >  1:651327/1‑71 (MQ=255)
  tATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAg                                                                   >  1:5554675/1‑71 (MQ=255)
  tATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAg                                                                   >  1:3445596/1‑71 (MQ=255)
   aTTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAg                                                                   >  1:1624172/1‑70 (MQ=255)
          acaGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATaa                                                              >  1:4695406/1‑68 (MQ=255)
          acaGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATaa                                                              >  1:1598061/1‑68 (MQ=255)
          acaGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATaa                                                              >  1:5647785/1‑68 (MQ=255)
              cGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCgg                                                       <  1:1554401/71‑1 (MQ=255)
                  tGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCcac                                                                      <  1:4179032/52‑1 (MQ=255)
                     aTTGTGTGAGTTGTATCTGGTGGAAGATGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCt                                                >  1:3113571/1‑71 (MQ=255)
                      ttGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAgg                                                                  >  1:1743950/1‑52 (MQ=255)
                      ttGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                               <  1:1742119/55‑1 (MQ=255)
                      ttGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATc                                                 <  1:1179722/69‑1 (MQ=255)
                      ttGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATc                                                 <  1:4372053/69‑1 (MQ=255)
                       tgtgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAAt                                                  >  1:1728731/1‑67 (MQ=255)
                       tgtgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAAt                                                  >  1:1273370/1‑67 (MQ=255)
                       tgtgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCtt                                               >  1:5460216/1‑70 (MQ=255)
                       tgtgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTc                                              >  1:4682293/1‑71 (MQ=255)
                       tgtgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTc                                              >  1:3988852/1‑71 (MQ=255)
                        gtgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTc                                              >  1:3169026/1‑70 (MQ=255)
                         tgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                               <  1:1142094/52‑1 (MQ=255)
                         tgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                               <  1:4757560/52‑1 (MQ=255)
                         tgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                               <  1:3759953/52‑1 (MQ=255)
                         tgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                               <  1:5010572/52‑1 (MQ=255)
                         tgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                               <  1:3197063/52‑1 (MQ=255)
                         tgtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTcc                                             >  1:5068504/1‑70 (MQ=255)
                          gtgAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCAc                                           >  1:4427468/1‑71 (MQ=255)
                            gAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                               <  1:1792050/49‑1 (MQ=255)
                            gAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                               <  1:1465795/49‑1 (MQ=255)
                             aGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATa                                                               <  1:813566/48‑1 (MQ=255)
                                tGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAAt                                                  >  1:4483828/1‑58 (MQ=255)
                                tGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAt                                     >  1:2839048/1‑71 (MQ=255)
                                 gTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAt                                     >  1:565999/1‑70 (MQ=255)
                                   aTCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAtttt                                  >  1:535536/1‑71 (MQ=255)
                                    tttgctgGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATc                                                 >  1:5242002/3‑55 (MQ=25)
                                     ttgctgGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGc                                                          >  1:3080222/2‑45 (MQ=19)
                                     ttgctgGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTc                                              >  1:3661546/2‑57 (MQ=25)
                                     ttgctgGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTc                                              >  1:443686/2‑57 (MQ=25)
                                        ctgGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAg                              <  1:3389495/70‑1 (MQ=255)
                                           gaagaaGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGa                             >  1:4708865/1‑68 (MQ=255)
                                           gaagaaGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGAtag                          >  1:4840571/1‑69 (MQ=255)
                                           gaagaaGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATg                           >  1:1281765/1‑70 (MQ=255)
                                           gaagaaGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATg                           <  1:3396548/70‑1 (MQ=255)
                                           caagaagTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGAtag                          >  1:4840576/2‑69 (MQ=255)
                                               aagTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTc                      >  1:3585077/1‑71 (MQ=255)
                                                agTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTc                      <  1:1851522/70‑1 (MQ=255)
                                                  tCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCAt                   >  1:1554285/1‑71 (MQ=255)
                                                   cATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCAtt                  >  1:3199116/1‑71 (MQ=255)
                                                    aTTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAAtt                                    <  1:3480619/52‑1 (MQ=255)
                                                    aTTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGAt                            >  1:5655099/1‑60 (MQ=255)
                                                     tttCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTc                                              <  1:1187718/41‑1 (MQ=255)
                                                     tttCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTc                                              <  1:3929580/41‑1 (MQ=255)
                                                     tttCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTcc                >  1:489388/1‑71 (MQ=255)
                                                      ttCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTc                                              <  1:4017569/40‑1 (MQ=255)
                                                      taccTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTc                                 <  1:384221/51‑1 (MQ=25)
                                                         cTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTcga                                      <  1:4102925/45‑3 (MQ=21)
                                                          tttGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCAt                   <  1:5383765/63‑1 (MQ=255)
                                                           ttgtcatcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTc          >  1:1254032/4‑71 (MQ=255)
                                                            tgtcatcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGAtttt                                  >  1:2888147/3‑46 (MQ=21)
                                                            tgtcatcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGAtttt                                  >  1:1219716/3‑46 (MQ=21)
                                                                atcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGa      >  1:2639127/1‑70 (MQ=255)
                                                                atcGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGa      >  1:5086069/1‑70 (MQ=255)
                                                                 tcCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTcc                                <  1:3869038/43‑1 (MQ=255)
                                                                 tcCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTcc                                <  1:2663123/43‑1 (MQ=255)
                                                                  ccACAAGGATAAACGCCGGGAAATCTTCCACTTCaa                                      <  1:995606/36‑1 (MQ=255)
                                                                  ccACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTccc               >  1:2747386/1‑59 (MQ=255)
                                                                   cacaAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGatat   >  1:1175426/1‑70 (MQ=255)
                                                                   cacaAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGatat   >  1:1079912/1‑70 (MQ=255)
                                                                   cacaAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATAtt  >  1:453674/1‑71 (MQ=255)
                                                                   |                                                                      
TTTATTTCACACAGCGGGTGCATTGTGTGAGTTGTATCTGCTGGAAGAAGTCATTTCCTTTATCATCCACAAGGATAAACGCCGGGAAATCTTCCACTTCAATTTTCCAGATGGCTTCCATTCCCAGTTCCGGATATT  >  minE/1014421‑1014558

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