Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,775 C→G 18.4% A98A (GCG→GCC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,7750CG18.4% 103.9 / 17.2 49A98A (GCG→GCCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (21/19);  new base G (5/4);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

ATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGT  >  minE/1015705‑1015842
                                                                      |                                                                   
aTACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTc                                                                     <  1:1151863/71‑1 (MQ=255)
aTACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTc                                                                     >  1:3650727/1‑71 (MQ=255)
  acgcGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggc                                                                   >  1:1523251/4‑71 (MQ=255)
  acacGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgc                                                                   >  1:599650/1‑71 (MQ=255)
          cccTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATg                                                            <  1:2857399/70‑1 (MQ=255)
            cTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATg                                                            <  1:2847556/68‑1 (MQ=255)
            cTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCg                                                         >  1:1683019/1‑71 (MQ=255)
             ttttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcg                                                                    >  1:415978/1‑59 (MQ=255)
              tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                                                                >  1:3572142/1‑62 (MQ=255)
               ttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAg                                                      >  1:3476318/1‑71 (MQ=255)
               ttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAg                                                      >  1:1662617/1‑71 (MQ=255)
                tttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGt                                                     >  1:3416392/1‑71 (MQ=255)
                 ttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGt                                                     <  1:1854993/70‑1 (MQ=255)
                  tACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGtt                                                    <  1:2008371/70‑1 (MQ=255)
                  tACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCgg                                                        <  1:1229421/66‑1 (MQ=255)
                           taatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAt                                          >  1:3900731/1‑71 (MQ=255)
                             atCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTg                                        >  1:3847882/1‑71 (MQ=255)
                              tCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTa                                                   >  1:2662165/1‑59 (MQ=255)
                               cgcgGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAg                                                      <  1:1979499/55‑1 (MQ=255)
                                gcgGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTg                                        >  1:79115/1‑68 (MQ=255)
                                gcgGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCa                                      <  1:361761/70‑1 (MQ=255)
                                  ggTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAgcg                                   <  1:800964/71‑1 (MQ=255)
                                      ccGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCca                                >  1:294181/1‑70 (MQ=255)
                                      ccGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCca                                >  1:2654439/1‑70 (MQ=255)
                                      ccGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCca                                >  1:3525072/1‑70 (MQ=255)
                                          tATCCTTACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTg                                        <  1:2764797/58‑1 (MQ=255)
                                           aTCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACatat                           >  1:1576625/1‑70 (MQ=255)
                                                 gcAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTcgt                    <  1:1186741/70‑3 (MQ=255)
                                                   aGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGaa                                           >  1:923204/1‑46 (MQ=255)
                                                   aGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCAtt                   >  1:1658305/1‑70 (MQ=255)
                                                      ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg               <  1:2836493/71‑1 (MQ=255)
                                                      ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg               <  1:2832034/71‑1 (MQ=255)
                                                      ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg               <  1:2719413/71‑1 (MQ=255)
                                                      ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg               <  1:2065601/71‑1 (MQ=255)
                                                      ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg               <  1:3139360/71‑1 (MQ=255)
                                                      ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg               <  1:476506/71‑1 (MQ=255)
                                                      tcggCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCg               >  1:1036365/3‑71 (MQ=255)
                                                      tcggCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCg               <  1:2486585/69‑1 (MQ=255)
                                                       tGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGaa                                           >  1:3794479/1‑42 (MQ=255)
                                                        ggCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTcgt                    >  1:2607250/1‑62 (MQ=255)
                                                        ggCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGc              <  1:1344389/70‑1 (MQ=255)
                                                        ggCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCt             >  1:3631046/1‑71 (MQ=255)
                                                         gCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTg            >  1:529851/1‑71 (MQ=255)
                                                            gAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGcc         >  1:100230/1‑71 (MQ=255)
                                                             caCGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGct         >  1:2195216/2‑69 (MQ=255)
                                                             aaCGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCAttt                  <  1:1878438/61‑1 (MQ=255)
                                                               cGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAt                                          >  1:1677593/1‑35 (MQ=255)
                                                               cGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAt                                          >  1:1439749/1‑35 (MQ=255)
                                                                gCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTctg     >  1:3846492/1‑69 (MQ=255)
                                                                gCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCg     <  1:1622153/71‑1 (MQ=255)
                                                                 ccTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCgg    <  1:2916760/71‑1 (MQ=255)
                                                                 ccTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCgg    <  1:3641654/71‑1 (MQ=255)
                                                                  cTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg               <  1:2452771/59‑1 (MQ=255)
                                                                   tttCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGt  >  1:821981/1‑71 (MQ=255)
                                                                     tCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTg            >  1:586505/1‑59 (MQ=255)
                                                                      cgccgcGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGc              >  1:945173/1‑56 (MQ=255)
                                                                      |                                                                   
ATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGT  >  minE/1015705‑1015842

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: