Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,719 T→C 5.5% Q450Q (CAA→CAG) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,7190TC 5.5% 272.9 / 4.7 91Q450Q (CAA→CAG) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (51/35);  new base C (3/2);  total (54/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCG  >  minE/1014651‑1014786
                                                                    |                                                                   
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: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)
       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)
                 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:516534/1‑55 (MQ=255)
                   tgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATc                                                                >  1:2548659/1‑55 (MQ=255)
                   tgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAatccatc                                                >  1:2617295/1‑71 (MQ=255)
                       gTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgtc                                            <  1:2894135/71‑1 (MQ=255)
                       gTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCAtccgtc                                            <  1:3149736/71‑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:2405443/56‑1 (MQ=255)
                            cTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGaa                                                      <  1:262842/56‑1 (MQ=255)
                            cTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcggc                                        <  1:1962246/70‑1 (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:2633440/1‑41 (MQ=255)
                               tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                                  >  1:3512581/1‑41 (MQ=255)
                               tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                                  >  1:3466974/1‑41 (MQ=255)
                               tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                                  >  1:1579359/1‑41 (MQ=255)
                               tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                                                  >  1:2329721/1‑41 (MQ=255)
                               tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTc                                    >  1:1988346/1‑71 (MQ=255)
                               tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTc                                    >  1:2096966/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: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:1447578/1‑49 (MQ=255)
                                                cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcgg                                         >  1:1524483/1‑49 (MQ=255)
                                                cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCcgg                                         >  1:1099523/1‑49 (MQ=255)
                                                cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc                             >  1:763165/1‑61 (MQ=255)
                                                cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc                             >  1:2426018/1‑61 (MQ=255)
                                                cTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc                             >  1:2532796/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)
                                                              ttGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGa                <  1:2219880/60‑1 (MQ=255)
                                                               tGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGc                             <  1:2071370/46‑1 (MQ=255)
                                                                gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGaa                     >  1:680422/1‑53 (MQ=255)
                                                                gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGaa                     >  1:1068468/1‑53 (MQ=255)
                                                                gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTc    <  1:3379689/70‑1 (MQ=255)
                                                                gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTc    <  1:1647673/70‑1 (MQ=255)
                                                                 cAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCg  >  1:2766370/1‑71 (MQ=255)
                                                                    ttGATCGTCATAAGAATCCATCCGTCCGGCGGTCGTTggg                              <  1:3111988/40‑1 (MQ=255)
                                                                    |                                                                   
GCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCG  >  minE/1014651‑1014786

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: