Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,760 T→C 14.6% P103P (CCA→CCG fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,7600TC14.6% 113.2 / 12.1 48P103P (CCA→CCGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (25/16);  new base C (4/3);  total (29/19)
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.38e-01

ACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCG  >  minE/1015691‑1015829
                                                                     |                                                                     
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)
           tCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCt                                                           >  1:1415530/1‑71 (MQ=255)
           tCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCt                                                           >  1:1480523/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:1151863/71‑1 (MQ=255)
              aTACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTc                                                        >  1:3650727/1‑71 (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)
                       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)
                           ttttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcg                                                       >  1:415978/1‑59 (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)
                                           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:2654439/1‑70 (MQ=255)
                                                    ccGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCca                   >  1:294181/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:2832034/71‑1 (MQ=255)
                                                                    ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg  <  1:3139360/71‑1 (MQ=255)
                                                                    ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg  <  1:476506/71‑1 (MQ=255)
                                                                    ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg  <  1:2719413/71‑1 (MQ=255)
                                                                    ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg  <  1:2065601/71‑1 (MQ=255)
                                                                    ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg  <  1:2836493/71‑1 (MQ=255)
                                                                     tGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGaa                              >  1:3794479/1‑42 (MQ=255)
                                                                     |                                                                     
ACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCG  >  minE/1015691‑1015829

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: