Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,701 T→G 5.6% G456G (GGA→GGC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,7010TG 5.6% 243.9 / 8.5 89G456G (GGA→GGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (53/31);  new base G (3/2);  total (56/33)
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.37e-01

CGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCG  >  minE/1014637‑1014770
                                                                |                                                                     
cgTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTgccgccg                                                                 <  1:3049470/71‑2 (MQ=255)
 gtgtTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTgccgccg                                                                 <  1:1691351/70‑2 (MQ=255)
    ttttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgcg                                                             >  1:3358151/1‑71 (MQ=255)
     tttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCt                                                                >  1:1178129/1‑67 (MQ=255)
     tttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCt                                                                >  1:2549758/1‑67 (MQ=255)
      ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                                                           >  1:1725149/1‑71 (MQ=255)
      ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                                                           >  1:117427/1‑71 (MQ=255)
      ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                                                           >  1:2443948/1‑71 (MQ=255)
       tttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgc                                                              >  1:66361/1‑67 (MQ=255)
       tttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCtt                                                          >  1:3504993/1‑71 (MQ=255)
        ttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCtt                                                          <  1:63543/70‑1 (MQ=255)
          aCAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg                                                         >  1:1779442/1‑69 (MQ=255)
          aCAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg                                                         >  1:3463334/1‑69 (MQ=255)
              gCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTg                                                   <  1:1115746/71‑1 (MQ=255)
                aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt                                                 >  1:256349/1‑71 (MQ=255)
                aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt                                                 >  1:688963/1‑71 (MQ=255)
                aTCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCACCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                  >  1:1528744/1‑70 (MQ=255)
                aTCCGTCACCTGCTGGCTGCGGTTGCCTTTAGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGAt                                                 >  1:1630231/1‑71 (MQ=255)
                 tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACtt                                                                       <  1:842147/48‑1 (MQ=255)
                 tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc                                                >  1:3298645/1‑71 (MQ=255)
                 tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc                                                >  1:2316786/1‑71 (MQ=255)
                   cGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCg                                               <  1:3010248/70‑1 (MQ=255)
                   cGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGa                                              <  1:96148/71‑1 (MQ=255)
                    gTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTc                                                                      >  1:3073524/1‑46 (MQ=255)
                     tCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGtt                                                    <  1:1729560/63‑1 (MQ=255)
                      cACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACa                                            >  1:1543197/1‑70 (MQ=255)
                       aCCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACAt                                           >  1:3553830/1‑70 (MQ=255)
                       aCCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCAcg                                            >  1:1696151/1‑68 (MQ=255)
                          tgctgGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGa                                              >  1:3612077/1‑64 (MQ=255)
                          tgctgGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg                                        >  1:1654041/1‑70 (MQ=255)
                          tgctgGCTGCGGTTGCCTTTCGCCAGCATGATCAAACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAg                                        >  1:2290286/1‑70 (MQ=255)
                            ctgactgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGagt                                     >  1:1372482/1‑69 (MQ=255)
                            ctgactgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGagt                                     >  1:3391991/1‑69 (MQ=255)
                              ggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc                                                        >  1:860589/1‑50 (MQ=255)
                              ggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc                                                        >  1:1696675/1‑50 (MQ=255)
                               gctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                  <  1:3372384/55‑1 (MQ=255)
                               gctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                  <  1:3376674/55‑1 (MQ=255)
                               gctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatcc                                   <  1:1438703/70‑1 (MQ=255)
                               cctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCTCCCTGCGCTTGCAGTTGATCGACATAAGAatcca                                  <  1:2763591/70‑1 (MQ=255)
                                ctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                  <  1:1423252/54‑1 (MQ=255)
                                ctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat                                     <  1:2274862/67‑1 (MQ=255)
                                ctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat                                     <  1:354527/67‑1 (MQ=255)
                                 tgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc                                                >  1:516534/1‑55 (MQ=255)
                                 tgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc                                                >  1:2548659/1‑55 (MQ=255)
                                 tgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccatc                                >  1:2617295/1‑71 (MQ=255)
                                     gTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgtc                            <  1:3149736/71‑1 (MQ=255)
                                     gTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgtc                            <  1:2894135/71‑1 (MQ=255)
                                      ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc                                                        <  1:1810521/42‑1 (MQ=255)
                                      ttGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGAc                                             <  1:666455/53‑1 (MQ=255)
                                        accTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtc                                    <  1:3541972/59‑1 (MQ=255)
                                        accTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtc                                    <  1:897784/59‑1 (MQ=255)
                                         ccTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGtt                                                    <  1:1483595/43‑1 (MQ=255)
                                         ccTTTCGCCAGCATGATCATACTTCCGCCCTGCCCTTGCAGTTGATCGACATAAGAATCCATCCGTCcgg                         <  1:3420205/70‑1 (MQ=255)
                                          cTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGaa                                      <  1:262842/56‑1 (MQ=255)
                                          cTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGaa                                      <  1:2405443/56‑1 (MQ=255)
                                          cTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcggc                        <  1:1962246/70‑1 (MQ=255)
                                            ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc                                                        >  1:2233905/1‑36 (MQ=255)
                                            ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATaa                                         >  1:1834072/1‑51 (MQ=255)
                                            ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATaa                                         >  1:3639728/1‑51 (MQ=255)
                                            ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGt                     >  1:1426070/1‑71 (MQ=255)
                                            ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGt                     >  1:3312496/1‑71 (MQ=255)
                                             tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                  >  1:2329721/1‑41 (MQ=255)
                                             tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                  >  1:1579359/1‑41 (MQ=255)
                                             tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                  >  1:3466974/1‑41 (MQ=255)
                                             tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                  >  1:3512581/1‑41 (MQ=255)
                                             tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                  >  1:2633440/1‑41 (MQ=255)
                                             tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTc                    >  1:2096966/1‑71 (MQ=255)
                                             tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTc                    >  1:1988346/1‑71 (MQ=255)
                                             tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTc                    >  1:2218883/1‑71 (MQ=255)
                                               gCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                  >  1:1852286/1‑71 (MQ=255)
                                                ccAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGtt                 <  1:1223378/71‑1 (MQ=255)
                                                 cAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                >  1:3237529/1‑71 (MQ=255)
                                                 cAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                <  1:1977890/71‑1 (MQ=255)
                                                 cAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                <  1:1549984/71‑1 (MQ=255)
                                                 cAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                <  1:849755/71‑1 (MQ=255)
                                                 cAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCAATCCGTCCGGCGGTCGTTg                >  1:3627355/1‑71 (MQ=255)
                                                   gCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                  >  1:1118202/1‑67 (MQ=255)
                                                    cATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc             <  1:765202/71‑1 (MQ=255)
                                                     aTGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGcc            <  1:2368238/71‑1 (MQ=255)
                                                       gATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGcc            >  1:2595611/1‑69 (MQ=255)
                                                       gATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGcc            >  1:1274646/1‑69 (MQ=255)
                                                        aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                  <  1:2606145/62‑1 (MQ=255)
                                                        aTCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAag         <  1:2876656/71‑1 (MQ=255)
                                                           aTACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCaa          <  1:664880/67‑1 (MQ=255)
                                                            tACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGaa     <  1:1289561/71‑1 (MQ=255)
                                                              cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcgg                         >  1:1099523/1‑49 (MQ=255)
                                                              cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcgg                         >  1:1524483/1‑49 (MQ=255)
                                                              cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcgg                         >  1:1447578/1‑49 (MQ=255)
                                                              cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc             >  1:763165/1‑61 (MQ=255)
                                                              cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc             >  1:2532796/1‑61 (MQ=255)
                                                              cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc             >  1:2426018/1‑61 (MQ=255)
                                                              cTTCCGCCCTGCGCTTGCAGTTGACCGACATAAGAATCCAtccgtccg                          >  1:1929395/1‑48 (MQ=255)
                                                               ttCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCg  >  1:2218762/1‑71 (MQ=255)
                                                               ttCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCg  >  1:3365301/1‑71 (MQ=255)
                                                                |                                                                     
CGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCG  >  minE/1014637‑1014770

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: