Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,872 T→C 6.3% K399K (AAA→AAG fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,8720TC 6.3% 224.6 / 5.0 79K399K (AAA→AAGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (40/34);  new base C (3/2);  total (43/36)
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.85e-01

CGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACT  >  minE/1014803‑1014942
                                                                     |                                                                      
cgcgTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGtt                                                                       <  1:1717469/71‑1 (MQ=255)
cgcgTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGt                                                                        <  1:3527835/70‑1 (MQ=255)
 gcgTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGttt                                                                      <  1:3558368/71‑1 (MQ=255)
  cgTAGTAGATCGGGTGATCTTTGATGTACTGCTGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTg                                                                     >  1:1727343/1‑71 (MQ=255)
    tagtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGc                                                                   >  1:3248303/1‑71 (MQ=255)
      gtaAATCGGATGATCTTTGATGTACTGCGGCAGCTCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                                                 >  1:2768525/1‑71 (MQ=255)
      gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                                                 <  1:1093456/71‑1 (MQ=255)
      gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                                                 >  1:3375109/1‑71 (MQ=255)
      gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                                                 >  1:3070184/1‑71 (MQ=255)
      gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                                                 >  1:2522308/1‑71 (MQ=255)
      gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                                                 >  1:2888457/1‑71 (MQ=255)
      gtaAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                                                 <  1:2131266/71‑1 (MQ=255)
        aaaTCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                                                 >  1:2008715/1‑69 (MQ=255)
           tCGGGTGATCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGc                                                             >  1:101962/1‑70 (MQ=255)
           tCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGc                                                             >  1:1547852/1‑70 (MQ=255)
            cGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCaa                                                           >  1:3334761/1‑71 (MQ=255)
              gatgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAata                                                         >  1:2937026/1‑71 (MQ=255)
              gatgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAata                                                         <  1:2826406/71‑1 (MQ=255)
                tgatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGt                                                                 >  1:1553129/1‑61 (MQ=255)
                 gatCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCa                                                      >  1:2649052/1‑71 (MQ=255)
                 gatCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGc                                                                   <  1:448763/58‑1 (MQ=12)
                  atCTTTGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATctc                                                     >  1:3004617/1‑69 (MQ=255)
                     tttGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTg                                                                     <  1:1098068/52‑1 (MQ=255)
                     tttGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTg                                                                     <  1:2291074/52‑1 (MQ=255)
                       tGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgacc                                                 >  1:21887/1‑70 (MQ=255)
                       tGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgac                                                  >  1:2362979/1‑69 (MQ=255)
                       tGATGTACTGCGGAAGTTCTTTACCGGCGTCAATCAGCTCTTTCAGCTTGGCGTGTGCAATATCTcg                                                    <  1:3493042/67‑1 (MQ=255)
                          tGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgaccgac                                              <  1:2335525/70‑1 (MQ=255)
                           gTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                                          >  1:3587731/1‑57 (MQ=255)
                                 cGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATc                                       <  1:1534460/70‑1 (MQ=255)
                                 cGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATc                                       <  1:3243381/70‑1 (MQ=255)
                                 cGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCACGATAATCGt                                     >  1:1410276/1‑71 (MQ=255)
                                  ggCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGt                                     <  1:2026537/71‑1 (MQ=255)
                                    cAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGc                                   <  1:3929315/71‑1 (MQ=255)
                                     aGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGcc                                  <  1:2207572/71‑1 (MQ=255)
                                     aGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGcc                                  <  1:1662428/71‑1 (MQ=255)
                                        ccTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAata                                                         <  1:1140457/45‑1 (MQ=255)
                                        ccTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAata                                                         <  1:468557/45‑1 (MQ=255)
                                        ccTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAata                                                         <  1:2550910/45‑1 (MQ=255)
                                        ccTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGtt                               <  1:963106/71‑1 (MQ=255)
                                         cTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAatat                                                        <  1:3145898/45‑1 (MQ=255)
                                          ttCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAat                                                          <  1:470210/42‑1 (MQ=255)
                                          ttCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCg                                 >  1:1592187/1‑67 (MQ=255)
                                          ttCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGtt                               >  1:993642/1‑69 (MQ=255)
                                          ttCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGt                                >  1:2782717/1‑68 (MQ=255)
                                           tCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCa                                                            <  1:3006637/39‑1 (MQ=255)
                                            cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgac                                                  >  1:3925190/1‑48 (MQ=255)
                                            cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGc                                   >  1:371294/1‑63 (MQ=255)
                                            cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGc                           >  1:606587/1‑71 (MQ=255)
                                            cACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGc                           >  1:1501386/1‑71 (MQ=255)
                                             aCCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCa                                                            <  1:3747060/37‑1 (MQ=255)
                                              ccGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGa                         <  1:2276/71‑1 (MQ=255)
                                               cGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCAcgac                                                  >  1:3413877/1‑45 (MQ=255)
                                               cGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACcgacga                                            >  1:1945741/1‑51 (MQ=255)
                                               cGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGAt                        <  1:1538474/71‑1 (MQ=255)
                                                   atccatccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGc                    >  1:527868/1‑71 (MQ=255)
                                                   atccatccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGc                    >  1:235065/1‑71 (MQ=255)
                                                   atccatccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGc                    >  1:3432997/1‑71 (MQ=255)
                                                    tccatccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCg                   <  1:3866052/71‑1 (MQ=255)
                                                     ccatccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGt                                     >  1:3147005/1‑52 (MQ=255)
                                                      catccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCTACGATAATCGTGCCGTTAAGCGATAAGCgt                  <  1:2392487/70‑1 (MQ=255)
                                                      catccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCgt                  <  1:1629771/70‑1 (MQ=255)
                                                       atccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCgtg                 >  1:2052253/1‑70 (MQ=255)
                                                       atccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCgt                  <  1:2641121/69‑1 (MQ=255)
                                                        tccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCgtg                 >  1:595255/1‑69 (MQ=255)
                                                         ccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCgtg                 >  1:2645732/1‑68 (MQ=255)
                                                         ccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCgtg                 >  1:2318366/1‑68 (MQ=255)
                                                         ccGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCgtg                 >  1:525312/1‑68 (MQ=255)
                                                          cGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAg                            >  1:1255571/1‑56 (MQ=255)
                                                            ctctTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGc                    <  1:1072998/62‑1 (MQ=255)
                                                            ctctTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCgtg                 <  1:1867520/65‑1 (MQ=255)
                                                            ctctTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGaaa           >  1:3617309/1‑71 (MQ=255)
                                                            ctctTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGaaa           <  1:1391083/71‑1 (MQ=255)
                                                             tctTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAg                            >  1:3668642/1‑53 (MQ=255)
                                                             tctTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCAGTTAAGCGATAAGCGTGTAGAAAc          >  1:1288027/1‑71 (MQ=255)
                                                               tttCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTa                              <  1:1252951/49‑1 (MQ=255)
                                                               tttCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAAc          <  1:2670274/69‑1 (MQ=255)
                                                               tttCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAAc          <  1:3723039/69‑1 (MQ=255)
                                                                 tCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAg                            >  1:1554862/1‑49 (MQ=255)
                                                                 tCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACgg        <  1:1525099/69‑1 (MQ=255)
                                                                   aGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATa    >  1:1125357/1‑71 (MQ=255)
                                                                    gTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCg                          <  1:2713305/48‑1 (MQ=255)
                                                                     tttGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACt  <  1:1643697/71‑1 (MQ=255)
                                                                     tttGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACt  >  1:1771348/1‑71 (MQ=255)
                                                                     tttGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACt  <  1:3720232/71‑1 (MQ=255)
                                                                     |                                                                      
CGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGTTATCCATCCGCTCTTTCAGTTTGGCGTGAGCAATATCACGACCGACGATAATCGTGCCGTTAAGCGATAAGCGTGTAGAAACGGGATACT  >  minE/1014803‑1014942

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: