Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,595 A→G 16.1% L492L (TTG→CTG)  fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,5950AG16.1% 137.0 / 18.4 62L492L (TTG→CTG) fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base G (3/7);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

TTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCT  >  minE/1014525‑1014663
                                                                      |                                                                    
ttCCAGATGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCaa                                                                      <  1:839897/70‑1 (MQ=255)
   cAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATa                                                                    <  1:4875235/70‑1 (MQ=255)
         gCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGc                                                             <  1:3358408/71‑1 (MQ=255)
         gCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGc                                                             >  1:838174/1‑71 (MQ=255)
         gCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGc                                                             >  1:2413168/1‑71 (MQ=255)
          cTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGcc                                                            <  1:1383428/71‑1 (MQ=255)
          cTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGcc                                                            <  1:1576348/71‑1 (MQ=255)
          cTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGcc                                                            >  1:836335/1‑71 (MQ=255)
            tCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGc                                                                >  1:4193030/1‑65 (MQ=255)
             ccATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGa                                                         >  1:3342986/1‑71 (MQ=255)
                  cccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCaa                                                                      >  1:4928335/1‑53 (MQ=255)
                  cccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATa                                                                    >  1:5630077/1‑55 (MQ=255)
                  cccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGaccacc                                                    >  1:1637925/1‑71 (MQ=255)
                  cccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGaccacc                                                    >  1:5442996/1‑71 (MQ=255)
                   ccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGaccacc                                                    <  1:1525995/70‑1 (MQ=255)
                    cAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCg                                                           >  1:1622604/1‑62 (MQ=255)
                      gTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGcc                                                            <  1:5281274/59‑1 (MQ=255)
                       ttCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATAc                                               <  1:607906/71‑1 (MQ=255)
                          cGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCg                                                           >  1:5482691/1‑56 (MQ=255)
                          cGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCg                                                           >  1:2233633/1‑56 (MQ=255)
                          cGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCg                                                           >  1:4276901/1‑56 (MQ=255)
                          cGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCg                                                           >  1:1588932/1‑56 (MQ=255)
                           tgATATTCAAAACATTCCAGGCTCTTAATACTTCCCTGCGCCAATAc                                                                   <  1:220283/46‑1 (MQ=255)
                           ggATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGcc                                           <  1:3217333/71‑1 (MQ=255)
                                ttCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCgg                                                          >  1:4067478/1‑51 (MQ=255)
                                   ctcacaTTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                    >  1:1360593/3‑70 (MQ=255)
                                   aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                    <  1:2520645/70‑1 (MQ=255)
                                   aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                    <  1:2641148/70‑1 (MQ=255)
                                   aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                    <  1:4134195/70‑1 (MQ=255)
                                   aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                    <  1:1980413/70‑1 (MQ=255)
                                    acacATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGac                                                        <  1:65243/49‑1 (MQ=255)
                                          tCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCaa                                         >  1:5568517/1‑58 (MQ=255)
                                          tCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAgccgcc                            >  1:924183/1‑71 (MQ=255)
                                             aGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtg                         <  1:3329862/71‑1 (MQ=255)
                                             aGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgt                          <  1:2919898/70‑1 (MQ=255)
                                               gCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAgccgcc                            <  1:1879519/66‑1 (MQ=255)
                                                 tctTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGa                                                  <  1:2680070/42‑1 (MQ=255)
                                                  ctTAATACTTCCCTGCGCCAATACAGCGGCCGGac                                                        <  1:4384918/35‑1 (MQ=255)
                                                  ctTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGttttt                    >  1:4765620/1‑71 (MQ=255)
                                                   ttAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCa                                          <  1:462459/48‑1 (MQ=255)
                                                    aaaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATa                                                >  1:4431457/2‑41 (MQ=255)
                                                     aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGAt                                                 >  1:5015040/1‑39 (MQ=255)
                                                     aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATa                                                >  1:1266360/1‑40 (MQ=255)
                                                     aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATa                                                >  1:4257964/1‑40 (MQ=255)
                                                     aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGa                                   >  1:1823041/1‑53 (MQ=255)
                                                       tACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAgg                                       <  1:3616730/47‑1 (MQ=255)
                                                          ttCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                    <  1:162245/47‑1 (MQ=255)
                                                            cccTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa            <  1:5435483/69‑1 (MQ=255)
                                                            cccTGCGCCAATACAGCGGCCGGACCACCGATACAGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa            <  1:3040774/69‑1 (MQ=255)
                                                              cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGcgt           <  1:2813926/68‑3 (MQ=255)
                                                              cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTcgg        >  1:5648506/1‑69 (MQ=255)
                                                              cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTcg         >  1:5141058/1‑69 (MQ=255)
                                                              cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTc          <  1:227697/69‑1 (MQ=255)
                                                              cTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa            >  1:5498553/1‑67 (MQ=255)
                                                              cTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCg        <  1:406601/71‑1 (MQ=255)
                                                                 cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGtaa     >  1:4676832/1‑69 (MQ=255)
                                                                 cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGta      <  1:2331422/70‑2 (MQ=255)
                                                                 cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGta      <  1:2271464/70‑2 (MQ=255)
                                                                 cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt       <  1:4717182/69‑1 (MQ=255)
                                                                 cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt       <  1:2188830/69‑1 (MQ=255)
                                                                 cgcCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGt       <  1:3096114/69‑1 (MQ=255)
                                                                 cgcCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTc      <  1:1776041/70‑1 (MQ=255)
                                                                    cAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct  >  1:2937934/1‑71 (MQ=255)
                                                                    cAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACct  >  1:5625313/1‑71 (MQ=255)
                                                                      aTACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGccgccg                           >  1:3263390/1‑44 (MQ=255)
                                                                      aTACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCAt           <  1:5582778/60‑1 (MQ=255)
                                                                      |                                                                    
TTCCAGATGGCTTCCATTCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCT  >  minE/1014525‑1014663

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: