Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,611 A→G 22.9% G486G (GGT→GGC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,6110AG22.9% 83.8 / 23.8 48G486G (GGT→GGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (20/17);  new base G (8/3);  total (28/20)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

TCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGC  >  minE/1014542‑1014678
                                                                     |                                                                   
tCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGaccac                                                                    >  1:3749888/1‑71 (MQ=255)
  ccAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTTCGCCAATACAGCGGCCGGaccacc                                                                   >  1:999837/1‑70 (MQ=255)
   cAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGa                                                                 >  1:998939/1‑71 (MQ=255)
    aGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGa                                                                 >  1:2439872/1‑70 (MQ=255)
        ccGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGa                                                                 >  1:3498711/1‑66 (MQ=255)
           gATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCa                                                         <  1:1178934/71‑1 (MQ=255)
            atatTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGcc                                                          >  1:321020/1‑69 (MQ=255)
                tCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAgg                                                      >  1:1381880/1‑69 (MQ=255)
                 cAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTa                                                    <  1:2045148/70‑1 (MQ=255)
                  aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                   <  1:3894155/70‑1 (MQ=255)
                  aaCACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAg                                                   <  1:192785/70‑1 (MQ=255)
                      caTTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAgcc                                              <  1:2294078/71‑1 (MQ=255)
                       aTTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAgccg                                             >  1:3128097/1‑71 (MQ=255)
                        ttCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAgccg                                             >  1:1118823/1‑70 (MQ=255)
                          ccAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGccgccg                                          >  1:122986/1‑71 (MQ=255)
                           cAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgt                                         <  1:731203/71‑1 (MQ=255)
                            aGGCTCTTAATACTTCCCTGCGCCATTACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtg                                        <  1:983636/71‑1 (MQ=255)
                             ggCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtgt                                       >  1:780063/1‑71 (MQ=255)
                             ggCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtgt                                       >  1:705820/1‑71 (MQ=255)
                             ggCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtgt                                       >  1:481533/1‑71 (MQ=255)
                             ggCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtgt                                       >  1:2730254/1‑71 (MQ=255)
                              gCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgtgt                                       <  1:3249765/70‑1 (MQ=255)
                                tctTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGa                                                  >  1:304612/1‑57 (MQ=255)
                                tctTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCgt                                         <  1:1814657/66‑1 (MQ=255)
                                tctTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGtttt                                    >  1:2830426/1‑71 (MQ=255)
                                   tAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTa                                 >  1:2726405/1‑71 (MQ=255)
                                    aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGtt                                      >  1:927201/1‑65 (MQ=255)
                                    aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTa                                 >  1:522316/1‑70 (MQ=255)
                                    aaTACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCAAAGGTAGAAGCCGCCGTGtttttt                                  <  1:1505069/69‑1 (MQ=255)
                                     aTACTTCCCTACGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTAc                                <  1:1032632/70‑1 (MQ=255)
                                           tgcTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTc                         >  1:689321/3‑71 (MQ=255)
                                           cccTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCa                           <  1:1169085/69‑1 (MQ=255)
                                            ccTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTa                                                    >  1:3333425/1‑43 (MQ=255)
                                             cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACacg                             >  1:2586187/1‑63 (MQ=255)
                                             cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGcg                           <  1:3625366/67‑2 (MQ=255)
                                             cTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTcgg                       >  1:3057193/1‑69 (MQ=255)
                                                   cagtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACct                 >  1:1661580/4‑71 (MQ=255)
                                                    aaTACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTa                                 >  1:1843042/1‑54 (MQ=255)
                                                     gtACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctg                >  1:283608/2‑70 (MQ=255)
                                                      tACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACctgct              >  1:3128585/1‑71 (MQ=255)
                                                        cAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACctgctg             <  1:3721172/70‑1 (MQ=255)
                                                           cGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctggctg         <  1:2730593/71‑1 (MQ=255)
                                                           cGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTgctggct          >  1:3006410/1‑70 (MQ=255)
                                                             gCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCg       >  1:1321545/1‑71 (MQ=255)
                                                              ccGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCgg      <  1:421606/71‑1 (MQ=255)
                                                               cGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGt     <  1:1790629/71‑1 (MQ=255)
                                                                ggACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGtt    <  1:367644/71‑1 (MQ=255)
                                                                ggACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGt     >  1:2382111/1‑70 (MQ=255)
                                                                 gACCACCGATACTGCCAAGGTAGAAGCCGCCGTGtt                                      <  1:3204316/36‑1 (MQ=255)
                                                                 gACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTAc                                <  1:1575686/42‑1 (MQ=255)
                                                                  accgccGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGttac  >  1:2005342/1‑69 (MQ=255)
                                                                  accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGc  <  1:2117837/71‑1 (MQ=255)
                                                                  accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGc  >  1:2014914/1‑71 (MQ=255)
                                                                  accaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGc  <  1:1255534/71‑1 (MQ=255)
                                                                   ccaccGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGt                      >  1:3467307/1‑50 (MQ=255)
                                                                     |                                                                   
TCCCAGTTCCGGATATTCAACACATTCCAGGCTCTTAATACTTCCCTGCGCCAATACAGCGGCCGGACCACCGATACTGCCAAGGTAGAAGCCGCCGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGC  >  minE/1014542‑1014678

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: