Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,710 C→G 5.2% A453A (GCG→GCC) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,7100CG 5.2% 312.1 / 4.8 95A453A (GCG→GCC) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (57/33);  new base G (3/2);  total (60/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

TTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCAT  >  minE/1014641‑1014777
                                                                     |                                                                   
ttttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTgcg                                                                    >  1:3358151/1‑71 (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:3463334/1‑69 (MQ=255)
      aCAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTg                                                                >  1:1779442/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)
             tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc                                                       >  1:2316786/1‑71 (MQ=255)
             tCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc                                                       >  1:3298645/1‑71 (MQ=255)
               cGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCg                                                      <  1:3010248/70‑1 (MQ=255)
               cGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGa                                                     <  1:96148/71‑1 (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:1696675/1‑50 (MQ=255)
                          ggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc                                                               >  1:860589/1‑50 (MQ=255)
                           gctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                         <  1:3376674/55‑1 (MQ=255)
                           gctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                         <  1:3372384/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:354527/67‑1 (MQ=255)
                            ctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat                                            <  1:2274862/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:897784/59‑1 (MQ=255)
                                    accTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGAgtc                                           <  1:3541972/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:3639728/1‑51 (MQ=255)
                                        ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATaa                                                >  1:1834072/1‑51 (MQ=255)
                                        ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGt                            >  1:3312496/1‑71 (MQ=255)
                                        ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGt                            >  1:1426070/1‑71 (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)
                                         tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                         >  1:2329721/1‑41 (MQ=255)
                                         tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTc                           >  1:2218883/1‑71 (MQ=255)
                                         tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTc                           >  1:2096966/1‑71 (MQ=255)
                                         tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTc                           >  1:1988346/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:1549984/71‑1 (MQ=255)
                                             cAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                       <  1:1977890/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:1447578/1‑49 (MQ=255)
                                                          cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcgg                                >  1:1524483/1‑49 (MQ=255)
                                                          cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc                    >  1:2532796/1‑61 (MQ=255)
                                                          cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc                    >  1:2426018/1‑61 (MQ=255)
                                                          cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc                    >  1:763165/1‑61 (MQ=255)
                                                          cTTCCGCCCTGCGCTTGCAGTTGACCGACATAAGAATCCAtccgtccg                                 >  1:1929395/1‑48 (MQ=255)
                                                           ttCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCg         >  1:3365301/1‑71 (MQ=255)
                                                           ttCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCg         >  1:2218762/1‑71 (MQ=255)
                                                             ccccccTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCg         >  1:1393359/4‑69 (MQ=255)
                                                             ccGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                       >  1:1079546/1‑55 (MQ=255)
                                                             ccGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCg         >  1:2116742/1‑69 (MQ=255)
                                                              cGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                       <  1:110734/54‑1 (MQ=255)
                                                                cccTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAgg     >  1:3631389/1‑70 (MQ=255)
                                                                 ccTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCa   >  1:959778/1‑71 (MQ=255)
                                                                  cTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCAt  <  1:2673727/71‑1 (MQ=255)
                                                                   tGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCAt  >  1:1849560/1‑70 (MQ=255)
                                                                   tGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCAt  >  1:2384967/1‑70 (MQ=255)
                                                                     |                                                                   
TTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCAT  >  minE/1014641‑1014777

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: