Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,754 A→G 12.9% C105C (TGT→TGC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,7540AG12.9% 135.4 / 10.8 54C105C (TGT→TGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (25/22);  new base G (4/3);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATAT  >  minE/1015684‑1015817
                                                                      |                                                               
cTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                                                 >  1:2896941/1‑71 (MQ=255)
     tcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGt                                                            >  1:847111/1‑71 (MQ=255)
     tcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGt                                                            >  1:3633187/1‑71 (MQ=255)
       accaccaccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTg                                                          <  1:2462578/71‑1 (MQ=255)
       accaccaccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTg                                                          >  1:439802/1‑71 (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:1415530/1‑71 (MQ=255)
                  tCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCt                                               <  1:33834/71‑1 (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)
                      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)
                              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)
                                   tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAg                                                              <  1:3451798/39‑1 (MQ=255)
                                   tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCtt                                              <  1:1053496/55‑1 (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)
                                                    cgcgGTGCCGGTATCCTGACAGGTTGGCAGAACGCCt                                               <  1:2241640/37‑1 (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:3525072/1‑70 (MQ=255)
                                                           ccGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCca       >  1:294181/1‑70 (MQ=255)
                                                           ccGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCca       >  1:2654439/1‑70 (MQ=255)
                                                               tATCCTTACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTg               <  1:2764797/58‑1 (MQ=255)
                                                                aTCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACatat  >  1:1576625/1‑70 (MQ=255)
                                                                      |                                                               
CTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATAT  >  minE/1015684‑1015817

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: