Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,712 C→A 5.2% A453S (GCG→TCG) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,7120CA 5.2% 329.1 / 2.4 95A453S (GCG→TCG) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (57/33);  new base A (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 = 9.10e-01

TTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCT  >  minE/1014643‑1014782
                                                                     |                                                                      
ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                                                                       >  1:2443948/1‑71 (MQ=255)
ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                                                                       >  1:117427/1‑71 (MQ=255)
ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                                                                       >  1:1725149/1‑71 (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)
           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:3391991/1‑69 (MQ=255)
                      ctgactgCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGATCCACGTAGGagt                                                 >  1:1372482/1‑69 (MQ=255)
                        ggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc                                                                    >  1:1696675/1‑50 (MQ=255)
                        ggctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGc                                                                    >  1:860589/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:354527/67‑1 (MQ=255)
                          ctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat                                                 <  1:2274862/67‑1 (MQ=255)
                           tgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc                                                            >  1:2548659/1‑55 (MQ=255)
                           tgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc                                                            >  1:516534/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:1426070/1‑71 (MQ=255)
                                      ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGt                                 >  1:3312496/1‑71 (MQ=255)
                                       tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                              >  1:3466974/1‑41 (MQ=255)
                                       tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                              >  1:3512581/1‑41 (MQ=255)
                                       tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                              >  1:2329721/1‑41 (MQ=255)
                                       tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                              >  1:1579359/1‑41 (MQ=255)
                                       tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                              >  1:2633440/1‑41 (MQ=255)
                                       tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTc                                >  1:2218883/1‑71 (MQ=255)
                                       tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTc                                >  1:1988346/1‑71 (MQ=255)
                                       tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTc                                >  1:2096966/1‑71 (MQ=255)
                                         gCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGt                              >  1:1852286/1‑71 (MQ=255)
                                          ccAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGtt                             <  1:1223378/71‑1 (MQ=255)
                                           cAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                            <  1:1549984/71‑1 (MQ=255)
                                           cAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                            <  1:849755/71‑1 (MQ=255)
                                           cAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                            <  1:1977890/71‑1 (MQ=255)
                                           cAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg                            >  1:3237529/1‑71 (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:1274646/1‑69 (MQ=255)
                                                 gATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGcc                        >  1:2595611/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:1524483/1‑49 (MQ=255)
                                                        cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcgg                                     >  1:1447578/1‑49 (MQ=255)
                                                        cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcgg                                     >  1:1099523/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: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)
                                                                     cTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCt  <  1:1257919/71‑1 (MQ=255)
                                                                     cTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCt  <  1:3012474/71‑1 (MQ=255)
                                                                     |                                                                      
TTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCT  >  minE/1014643‑1014782

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: