Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,698 A→G 8.2% S457S (AGT→AGC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,6980AG 8.2% 243.5 / 11.1 85S457S (AGT→AGCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (46/32);  new base G (3/4);  total (49/36)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAA  >  minE/1014637‑1014767
                                                             |                                                                     
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:117427/1‑71 (MQ=255)
      ttttACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCt                                                        >  1:1725149/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: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:2274862/67‑1 (MQ=255)
                                ctgCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAat                                  <  1:354527/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: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:2405443/56‑1 (MQ=255)
                                          cTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGaa                                   <  1:262842/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:3312496/1‑71 (MQ=255)
                                            ttCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGt                  >  1:1426070/1‑71 (MQ=255)
                                             tCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGa                                               >  1:1579359/1‑41 (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: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:1977890/71‑1 (MQ=255)
                                                 cAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg             <  1:849755/71‑1 (MQ=255)
                                                 cAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTg             <  1:1549984/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)
                                                             |                                                                     
CGTGTTTTTTACAGGCATCCGTCACCTGCTGGCTGCGGTTGCCTTTCGCCAGCATGATCATACTTCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAA  >  minE/1014637‑1014767

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: