Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,730 A→G 27.1% I113I (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,7300AG27.1% 125.1 / 51.3 81I113I (ATT→ATCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (32/27);  new base G (11/11);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 101 total aligned reads displayed.

ATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGG  >  minE/1015661‑1015799
                                                                     |                                                                     
ataGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACcaaca                                                                       >  1:5127351/1‑70 (MQ=255)
 taGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACcaacaa                                                                      >  1:2528372/1‑70 (MQ=255)
 taGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACcaacaa                                                                      >  1:3402511/1‑70 (MQ=255)
 taGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaat                                                                     <  1:4727397/71‑1 (MQ=255)
   gACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaataa                                                                   >  1:3811583/1‑71 (MQ=255)
    acacCGCGCGCCAGCGCCGCTTCATCACCACCCCCGGTCCATACCCGCTGCCCTTTTTTACCAACaataat                                                                  >  1:520643/1‑71 (MQ=255)
     cacCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATc                                                                 <  1:3447993/71‑1 (MQ=255)
      acCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATcg                                                                <  1:2695996/71‑1 (MQ=255)
       ccGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATcgc                                                               >  1:3982361/1‑71 (MQ=255)
        cgcgcgCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATcg                                                                <  1:3658047/69‑1 (MQ=255)
         gcgcgcCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCgg                                                             >  1:232508/1‑71 (MQ=255)
                aGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCgg                                                       >  1:1551041/1‑70 (MQ=255)
                  cgccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTAt                                                    >  1:4722473/1‑71 (MQ=255)
                   gccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaat                                                                     >  1:5430502/1‑53 (MQ=255)
                   gccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACaat                                                                     >  1:1775721/1‑53 (MQ=255)
                   gccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATc                                                                 >  1:4629508/1‑57 (MQ=255)
                    ccgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATc                                                   <  1:1063355/70‑1 (MQ=255)
                     cgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCt                                                 <  1:3515801/71‑1 (MQ=255)
                     cgcTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCt                                                 <  1:3896823/71‑1 (MQ=255)
                        ttCATCACCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTg                                                <  1:1524142/69‑1 (MQ=255)
                         tcatcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCgg                                                       <  1:237579/61‑1 (MQ=255)
                          catcaCCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTggc                                              <  1:2661484/69‑3 (MQ=255)
                          catcaCCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGAGGTGCCGGTATCCTgtc                                              <  1:2363479/69‑3 (MQ=255)
                           atcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCgg                                                             >  1:2489002/1‑53 (MQ=255)
                           atcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                               >  1:3467340/1‑67 (MQ=255)
                           atcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                               >  1:5142457/1‑67 (MQ=255)
                           atcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAgg                                           <  1:770359/71‑1 (MQ=255)
                            tcaCCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGt                                          <  1:243615/71‑1 (MQ=255)
                            tcaCCGCCGCCGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGt                                          >  1:4035952/1‑71 (MQ=255)
                            tcaCCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                               <  1:646429/66‑1 (MQ=255)
                              accaccaccGGTCCATACACGCTGCCCTTTTTTACCAACaataat                                                                  >  1:4356031/1‑45 (MQ=255)
                              accaccaccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATcg                                                                >  1:1572158/1‑47 (MQ=255)
                               ccaccaccGGTCCATACACGCTGCCCTTTTTTACCAACaata                                                                    <  1:431388/42‑1 (MQ=255)
                               ccaccaccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATc                                                                 <  1:1340906/45‑1 (MQ=255)
                                  ccgccGGTCCACACGCGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAg                                    <  1:2558444/68‑1 (MQ=255)
                                  ccaccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAg                                    <  1:1060718/71‑1 (MQ=255)
                                   caccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGa                                   <  1:3155127/71‑1 (MQ=255)
                                   caccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGa                                   <  1:304152/71‑1 (MQ=255)
                                     ccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCgg                                                       >  1:3911609/1‑49 (MQ=255)
                                     ccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGaa                                  >  1:1963509/1‑70 (MQ=255)
                                       ggTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGaa                                  >  1:3296863/1‑68 (MQ=255)
                                       ggTCCACACGCGCTGGCCTTTTTTACCGACGATGATcgc                                                               <  1:5383109/39‑1 (MQ=9)
                                        gTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCgg                                                       <  1:2627401/46‑1 (MQ=255)
                                        gTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCgg                                                       <  1:5573713/46‑1 (MQ=255)
                                        gTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGcc                              >  1:4816617/1‑71 (MQ=255)
                                        gTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGc                               <  1:37240/70‑1 (MQ=255)
                                         tCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCgg                                                       <  1:4286053/45‑1 (MQ=255)
                                          ccATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTa                                                     >  1:3679966/1‑46 (MQ=255)
                                          ccATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGa                                               >  1:3388531/1‑52 (MQ=255)
                                            aTACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAg                                    <  1:1308118/61‑1 (MQ=255)
                                             tACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGc                               <  1:3271662/65‑1 (MQ=255)
                                              acgcGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggc                        >  1:3611482/4‑71 (MQ=255)
                                              acgcGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggc                        >  1:2231729/4‑71 (MQ=255)
                                              acgcGCTGGCCGTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggc                        >  1:5263149/4‑71 (MQ=255)
                                              acacGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGaa                                  >  1:3667765/1‑61 (MQ=255)
                                               cacGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCg                                                        >  1:3166586/1‑38 (MQ=255)
                                               cacGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgcc                       >  1:5557374/1‑71 (MQ=255)
                                                gcGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCtttgg                         >  1:3773496/2‑66 (MQ=255)
                                                 cGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcggc                     <  1:4973483/71‑1 (MQ=255)
                                                 cGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                     <  1:1316713/71‑1 (MQ=255)
                                                 cGCTGCCCTTTTTAACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                     <  1:1974887/71‑1 (MQ=255)
                                                 cGCTGCCCTTTTTAACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                     <  1:1075037/71‑1 (MQ=255)
                                                      gccTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcggc                     <  1:4051551/65‑1 (MQ=255)
                                                       ccTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTc               <  1:2657598/71‑1 (MQ=255)
                                                        cTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCt                             <  1:1546240/56‑1 (MQ=255)
                                                         ttttttACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcggc                     <  1:4500332/63‑1 (MQ=255)
                                                         ttttttACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcggc                     <  1:3290668/63‑1 (MQ=255)
                                                         ttttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                     <  1:3174130/63‑1 (MQ=255)
                                                         ttttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                     <  1:1661459/63‑1 (MQ=255)
                                                         ttttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGa                   >  1:4324635/1‑65 (MQ=255)
                                                         ttttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCg              >  1:1028860/1‑70 (MQ=255)
                                                         ttttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCg              >  1:2546137/1‑70 (MQ=255)
                                                         ttttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCg              >  1:150716/1‑70 (MQ=255)
                                                          tttttACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcggc                     <  1:3672295/62‑1 (MQ=255)
                                                          tttttACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCg              >  1:657737/1‑69 (MQ=255)
                                                          tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCttt                           >  1:5230059/1‑56 (MQ=255)
                                                          tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCttt                           >  1:5224497/1‑56 (MQ=255)
                                                          tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCttt                           >  1:5193701/1‑56 (MQ=255)
                                                          tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCg              >  1:2910695/1‑69 (MQ=255)
                                                          tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCg              >  1:4703573/1‑69 (MQ=255)
                                                          tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCg              >  1:4065337/1‑69 (MQ=255)
                                                           ttttACCGACGATGATCGCGGTGCCGGTATCCTGGCa                                             >  1:3142188/1‑37 (MQ=21)
                                                           ttttACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGt                                          >  1:4556733/1‑40 (MQ=25)
                                                           ttttACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGt                                          >  1:1528893/1‑40 (MQ=25)
                                                           ttttACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAg           >  1:5225673/1‑71 (MQ=255)
                                                           ttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGc                                      >  1:1806926/1‑44 (MQ=255)
                                                            tttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGt          >  1:1643407/1‑71 (MQ=255)
                                                             ttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGt          <  1:2777173/70‑1 (MQ=255)
                                                               aCCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGa                   >  1:5517858/1‑59 (MQ=255)
                                                               aCCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGa                   >  1:3033088/1‑59 (MQ=255)
                                                                   acaaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                     <  1:2455454/53‑1 (MQ=255)
                                                                   acaaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                     <  1:5245579/53‑1 (MQ=255)
                                                                   acaaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                     <  1:5196834/53‑1 (MQ=255)
                                                                   acaaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                     <  1:1375116/53‑1 (MQ=255)
                                                                    caaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAg                                    <  1:4250178/37‑1 (MQ=255)
                                                                    caaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAgg  <  1:1928952/71‑1 (MQ=255)
                                                                     aataatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGa                   >  1:3682989/1‑53 (MQ=255)
                                                                     aataatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGa                   >  1:3872794/1‑53 (MQ=255)
                                                                     aataatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGa                   >  1:3880397/1‑53 (MQ=255)
                                                                     aataatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAgg  >  1:3286898/1‑70 (MQ=255)
                                                                     |                                                                     
ATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGG  >  minE/1015661‑1015799

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: