Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,766 A→C 26.4% V101V (GTT→GTG fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,7660AC26.4% 103.1 / 67.6 91V101V (GTT→GTGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (45/22);  new base C (10/14);  total (55/36)
Fisher's exact test for biased strand distribution p-value = 5.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 101 total aligned reads displayed.

CACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCT  >  minE/1015696‑1015836
                                                                      |                                                                      
caccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGa                                                                        <  1:3155127/71‑1 (MQ=255)
caccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGa                                                                        <  1:304152/71‑1 (MQ=255)
  ccGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGaa                                                                       >  1:1963509/1‑70 (MQ=255)
    ggTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGaa                                                                       >  1:3296863/1‑68 (MQ=255)
     gTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGcc                                                                   >  1:4816617/1‑71 (MQ=255)
     gTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGc                                                                    <  1:37240/70‑1 (MQ=255)
          tACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGc                                                                    <  1:3271662/65‑1 (MQ=255)
           acgcGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggc                                                             >  1:2231729/4‑71 (MQ=255)
           acgcGCTGGCCTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggc                                                             >  1:3611482/4‑71 (MQ=255)
           acgcGCTGGCCGTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggc                                                             >  1:5263149/4‑71 (MQ=255)
           acacGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGaa                                                                       >  1:3667765/1‑61 (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:1075037/71‑1 (MQ=255)
              cGCTGCCCTTTTTAACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                                                          <  1:1974887/71‑1 (MQ=255)
                   gccTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcggc                                                          <  1:4051551/65‑1 (MQ=255)
                    ccTTTTTTACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTc                                                    >  1:851234/1‑71 (MQ=255)
                    ccTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTc                                                    <  1:2657598/71‑1 (MQ=255)
                     cTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCt                                                                  <  1:1546240/56‑1 (MQ=255)
                      ttttttACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcggc                                                          <  1:3290668/63‑1 (MQ=255)
                      ttttttACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTggcggc                                                          <  1:4500332/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:4065337/1‑69 (MQ=255)
                       tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCg                                                   >  1:2910695/1‑69 (MQ=255)
                       tttttACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCg                                                   >  1:4703573/1‑69 (MQ=255)
                        ttttACCGACGATGATCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAg                                                >  1:5225673/1‑71 (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:5196834/53‑1 (MQ=255)
                                acaaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                                                          <  1:1375116/53‑1 (MQ=255)
                                acaaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                                                          <  1:2455454/53‑1 (MQ=255)
                                acaaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTcgccgc                                                          <  1:5245579/53‑1 (MQ=255)
                                 caaTAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAgg                                       <  1:1928952/71‑1 (MQ=255)
                                  aataatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGa                                                        >  1:3880397/1‑53 (MQ=255)
                                  aataatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGa                                                        >  1:3872794/1‑53 (MQ=255)
                                  aataatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGa                                                        >  1:3682989/1‑53 (MQ=255)
                                  aataatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAgg                                       >  1:3286898/1‑70 (MQ=255)
                                   ataatCGCGGTGCCGGTATCCTGACAGGTTGGCAGaa                                                                       >  1:3048017/1‑37 (MQ=255)
                                   ataatCGCGGTGCCGGTATCCTGACAGGTTGGCAGaa                                                                       >  1:3298993/1‑37 (MQ=255)
                                    taatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAt                                    >  1:2007616/1‑71 (MQ=255)
                                      atCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTg                                  <  1:791435/71‑1 (MQ=255)
                                      atCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTg                                  <  1:3132523/71‑1 (MQ=255)
                                       tCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGcc                                                                   <  1:2539707/37‑1 (MQ=25)
                                              gCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCca                          >  1:4192488/1‑71 (MQ=255)
                                               ccGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGa                                                        >  1:346272/1‑40 (MQ=255)
                                                 ggTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCcac                         >  1:1102281/1‑69 (MQ=255)
                                                   tATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACata                      <  1:3301211/70‑1 (MQ=255)
                                                    aTCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATAtt                    >  1:5109385/1‑71 (MQ=255)
                                                    aTCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATAtt                    >  1:3329695/1‑71 (MQ=255)
                                                     tCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCcacg                        <  1:5676077/66‑2 (MQ=255)
                                                     tCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACgtac                     <  1:4853166/69‑2 (MQ=255)
                                                     tCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGc                                 <  1:1596908/57‑1 (MQ=255)
                                                     tCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCca                          <  1:1046048/64‑1 (MQ=255)
                                                     tCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATAtt                    >  1:542707/1‑70 (MQ=255)
                                                     tCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATAtt                    >  1:4150612/1‑70 (MQ=255)
                                                       cTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTa                  <  1:3355582/70‑1 (MQ=255)
                                                       cTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTAt                 >  1:3885883/1‑71 (MQ=255)
                                                       cTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTAt                 >  1:1937606/1‑71 (MQ=255)
                                                        tggcAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTc                >  1:3719848/4‑71 (MQ=255)
                                                         gACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTg                                  >  1:3165474/1‑52 (MQ=255)
                                                         gACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCa                                >  1:2935923/1‑54 (MQ=255)
                                                         gACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCa                                >  1:4388660/1‑54 (MQ=255)
                                                         gACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCa                                >  1:4178052/1‑54 (MQ=255)
                                                          gcAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTcgt              >  1:198594/2‑69 (MQ=255)
                                                          gcAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTcg               <  1:3760770/69‑2 (MQ=255)
                                                          aCAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATc                <  1:1445599/69‑1 (MQ=255)
                                                          aCAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCAt              >  1:4058642/1‑71 (MQ=255)
                                                          aCAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCAt              >  1:1035501/1‑71 (MQ=255)
                                                           cAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCAtt             >  1:5096043/1‑71 (MQ=255)
                                                           cAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCAtt             <  1:4048499/71‑1 (MQ=255)
                                                              gtcggCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTc          <  1:184778/68‑1 (MQ=255)
                                                               tcggCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCg         <  1:4983230/69‑1 (MQ=255)
                                                               tcggCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCg         <  1:800685/69‑1 (MQ=255)
                                                                 ggCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGc        <  1:1965448/70‑1 (MQ=255)
                                                                 ggCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCt       >  1:3020013/1‑71 (MQ=255)
                                                                   cAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTgg     >  1:1444056/1‑71 (MQ=255)
                                                                   cAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTgg     >  1:389951/1‑71 (MQ=255)
                                                                    aGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGc    >  1:2746248/1‑71 (MQ=255)
                                                                     gAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGcc   >  1:978119/1‑71 (MQ=255)
                                                                     gAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGcc   >  1:4326027/1‑71 (MQ=255)
                                                                     gAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGcc   >  1:4096225/1‑71 (MQ=255)
                                                                     gAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGcc   >  1:321512/1‑71 (MQ=255)
                                                                     gAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGc    >  1:334766/1‑70 (MQ=255)
                                                                      aaCGCCTTTCGCCGCGATGTCGGAGTTACGCAGGaa                                     <  1:2351181/36‑1 (MQ=255)
                                                                      aaCGCCTTTCGCCGCGATGTCGGAGTTACGCAGGaa                                     <  1:2324461/36‑1 (MQ=255)
                                                                      aaCGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCt  <  1:2604459/71‑1 (MQ=255)
                                                                      |                                                                      
CACCGGTCCATACACGCTGCCCTTTTTTACCAACAATAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCT  >  minE/1015696‑1015836

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: