Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,733 A→G 10.9% I112I (ATT→ATC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,7330AG10.9% 145.6 / 8.6 55I112I (ATT→ATCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base G (3/3);  total (27/28)
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.61e-01

ACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAT  >  minE/1015665‑1015802
                                                                    |                                                                     
acacCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTCTTACCAACaata                                                                       <  1:612243/69‑1 (MQ=255)
    cgcgcgCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATcgcg                                                                 <  1:1651817/71‑1 (MQ=255)
         gcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaataat                                                                     >  1:3707824/1‑62 (MQ=255)
         gcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaataat                                                                     >  1:3655628/1‑62 (MQ=255)
         gcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGcc                                                            >  1:3082898/1‑71 (MQ=255)
              cgccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTAt                                                       >  1:2402544/1‑71 (MQ=255)
              cgccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTAt                                                       >  1:3474463/1‑71 (MQ=255)
               gccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTa                                                        >  1:1417215/1‑69 (MQ=255)
                   cTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                                  >  1:2896941/1‑71 (MQ=255)
                     tcatcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGcc                                                            <  1:1444709/59‑1 (MQ=255)
                      catcaCCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTgg                                                  <  1:2032231/68‑2 (MQ=255)
                        tcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGt                                             >  1:847111/1‑71 (MQ=255)
                        tcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGt                                             >  1:3633187/1‑71 (MQ=255)
                          accaccaccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTg                                           >  1:439802/1‑71 (MQ=255)
                          accaccaccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTg                                           <  1:2462578/71‑1 (MQ=255)
                           ccaccaccGGTCCATACACGCTGCCCATTATTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTg                                           <  1:732587/70‑1 (MQ=255)
                            caccaccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTgg                                          >  1:224662/1‑70 (MQ=255)
                              ccgccGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAg                                       >  1:1350056/4‑71 (MQ=255)
                                  cGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACg                                   >  1:290796/1‑71 (MQ=255)
                                   ggTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGc                                  >  1:3045087/1‑71 (MQ=255)
                                    gTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGGGCCGGTATCCTGACAgg                                              <  1:1781090/58‑1 (MQ=255)
                                     tCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCt                                >  1:1480523/1‑71 (MQ=255)
                                     tCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCt                                <  1:33834/71‑1 (MQ=255)
                                     tCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCt                                >  1:1415530/1‑71 (MQ=255)
                                       cATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGcc                                 <  1:949855/68‑1 (MQ=255)
                                       aaTACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCa                                        >  1:3489735/2‑61 (MQ=255)
                                        aTACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTc                             >  1:3650727/1‑71 (MQ=255)
                                        aTACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTc                             <  1:1151863/71‑1 (MQ=255)
                                         tACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTa                                                        <  1:196216/43‑1 (MQ=255)
                                         tACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATc                                                      <  1:2073736/45‑1 (MQ=255)
                                         tACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCt                                                    <  1:3177596/47‑1 (MQ=255)
                                         tACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                                  <  1:2229768/49‑1 (MQ=255)
                                         tACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAg                                       <  1:921888/60‑1 (MQ=255)
                                         tACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCtt                               <  1:2996289/68‑1 (MQ=255)
                                          acgcGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggc                           >  1:1523251/4‑71 (MQ=255)
                                          acacGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGtt                                            <  1:2033882/54‑1 (MQ=255)
                                          acacGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgc                           >  1:599650/1‑71 (MQ=255)
                                            gcGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAg                                               <  1:3920308/48‑1 (MQ=25)
                                              gCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACa                                                <  1:3428337/46‑1 (MQ=255)
                                              gCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACa                                                <  1:568705/46‑1 (MQ=255)
                                               cTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTa                                                        <  1:3925324/37‑1 (MQ=255)
                                                 ggccTTTTTTACCGACGATAATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCAc                                    >  1:2483132/3‑55 (MQ=255)
                                                 gCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTCGCAGAACGcc                                 >  1:2014258/1‑58 (MQ=255)
                                                  cccTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATg                    <  1:2857399/70‑1 (MQ=255)
                                                    cTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATg                    <  1:2847556/68‑1 (MQ=255)
                                                    cTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCg                 >  1:1683019/1‑71 (MQ=255)
                                                     ttttttACCAACAATAATCGCGGTGCCGGTATCCTGACAg                                               <  1:3453674/40‑1 (MQ=255)
                                                     ttttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcg                            >  1:415978/1‑59 (MQ=255)
                                                      tttttACCGACGATGATCACGGTGCCGGTATCCTg                                                   >  1:1644390/1‑35 (MQ=19)
                                                      tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAg                                               <  1:3451798/39‑1 (MQ=255)
                                                      tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCtt                               <  1:1053496/55‑1 (MQ=255)
                                                      tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                        >  1:3572142/1‑62 (MQ=255)
                                                       ttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAg              >  1:1662617/1‑71 (MQ=255)
                                                       ttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAg              >  1:3476318/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)
                                                                    |                                                                     
ACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAT  >  minE/1015665‑1015802

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: