Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,817 T→C 20.3% K84K (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,015,8170TC20.3% 119.0 / 25.5 59K84K (AAA→AAGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (16/31);  new base C (7/5);  total (23/36)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCA  >  minE/1015748‑1015884
                                                                     |                                                                   
aTCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACatat                                                                     >  1:1576625/1‑70 (MQ=255)
      gcAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTcgt                                                              <  1:1186741/70‑3 (MQ=255)
        aGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCAtt                                                             >  1:1658305/1‑70 (MQ=255)
           ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg                                                         <  1:476506/71‑1 (MQ=255)
           ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg                                                         <  1:3139360/71‑1 (MQ=255)
           ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg                                                         <  1:2836493/71‑1 (MQ=255)
           ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg                                                         <  1:2832034/71‑1 (MQ=255)
           ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg                                                         <  1:2719413/71‑1 (MQ=255)
           ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg                                                         <  1:2065601/71‑1 (MQ=255)
           tcggCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCg                                                         <  1:2486585/69‑1 (MQ=255)
           tcggCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCg                                                         >  1:1036365/3‑71 (MQ=255)
             ggCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTcgt                                                              >  1:2607250/1‑62 (MQ=255)
             ggCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGc                                                        <  1:1344389/70‑1 (MQ=255)
             ggCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCt                                                       >  1:3631046/1‑71 (MQ=255)
              gCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTg                                                      >  1:529851/1‑71 (MQ=255)
                 gAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGcc                                                   >  1:100230/1‑71 (MQ=255)
                  caCGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGct                                                   >  1:2195216/2‑69 (MQ=255)
                  aaCGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCAttt                                                            <  1:1878438/61‑1 (MQ=255)
                     gCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTctg                                               >  1:3846492/1‑69 (MQ=255)
                     gCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCg                                               <  1:1622153/71‑1 (MQ=255)
                      ccTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCgg                                              <  1:2916760/71‑1 (MQ=255)
                      ccTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCgg                                              <  1:3641654/71‑1 (MQ=255)
                       cTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg                                                         <  1:2452771/59‑1 (MQ=255)
                        tttCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGt                                            >  1:821981/1‑71 (MQ=255)
                          tCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTg                                                      >  1:586505/1‑59 (MQ=255)
                           cgccgcGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGc                                                        >  1:945173/1‑56 (MQ=255)
                              cgcgATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGAACTTGTCGTTTTCGCTGGCTTCTGGATcgtg                                        >  1:1904568/1‑66 (MQ=255)
                               gcgATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATcgtg                                        >  1:50103/1‑65 (MQ=255)
                               gcgATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAg                                     <  1:1130330/71‑1 (MQ=255)
                               gcgATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAg                                     <  1:284531/71‑1 (MQ=255)
                                cgATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGa                                    <  1:1570174/71‑1 (MQ=255)
                                cgATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGa                                    <  1:2946017/71‑1 (MQ=255)
                                  aTGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAAt                                  <  1:428186/71‑1 (MQ=255)
                                   tGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCggg                                             >  1:340474/1‑59 (MQ=255)
                                    gTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGc                                                        <  1:1947523/47‑1 (MQ=255)
                                    gTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTc                                                 <  1:2011663/54‑1 (MQ=255)
                                        gAGTTACGCATGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGc                            <  1:579382/71‑1 (MQ=255)
                                        gAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCg                              >  1:112155/1‑69 (MQ=255)
                                        gAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGc                            <  1:2280586/71‑1 (MQ=255)
                                        gAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGc                            >  1:2103028/1‑71 (MQ=255)
                                          gTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCa                          >  1:50623/1‑71 (MQ=255)
                                           ttACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTACGGGTCACGCAGAATGTCGGCCAc                         <  1:1390487/71‑1 (MQ=255)
                                              cGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCtg                      <  1:1286828/71‑1 (MQ=255)
                                                cAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCtgtt                    <  1:2228343/71‑1 (MQ=255)
                                                cAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCAGGTCACGCAGAATGTCGGCCAcc                        <  1:2561581/67‑1 (MQ=255)
                                                  ggAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGcc                           <  1:196257/62‑1 (MQ=255)
                                                  ggAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCtgttgt                  <  1:3369008/71‑1 (MQ=255)
                                                   gAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCtgtttct                 <  1:1757065/71‑4 (MQ=255)
                                                    aaTTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCtgttg                   <  1:477566/68‑1 (MQ=255)
                                                    aaTTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGttgttg                <  1:277602/71‑1 (MQ=255)
                                                      ttGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTgga                                             <  1:2691527/40‑2 (MQ=12)
                                                      ttGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGt              <  1:1844678/71‑1 (MQ=255)
                                                      ttGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTCTGGATCGTGAAGAATAGCCGCAACCTGTTTCTGGt              >  1:578541/1‑71 (MQ=255)
                                                      ttGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCtgttgt                  >  1:1546073/1‑67 (MQ=255)
                                                        gCAGCGCCACATATTTATCATTTTCGCTGGCCTCCgg                                              >  1:1952346/1‑37 (MQ=255)
                                                         cAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTgcg           <  1:3243852/71‑1 (MQ=255)
                                                         cAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTgcg           <  1:246284/71‑1 (MQ=255)
                                                           gcgcCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTgcgcg         >  1:3648842/1‑70 (MQ=255)
                                                               cacaTATTTATCATTTTCGCTGGCCTCGGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGAc     >  1:3078885/1‑71 (MQ=255)
                                                                acaTATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACg    >  1:3829596/1‑71 (MQ=255)
                                                                  atatTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACgca  <  1:2138846/71‑1 (MQ=255)
                                                                     tttATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTc                               >  1:1970993/1‑39 (MQ=255)
                                                                     |                                                                   
ATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGTTGGTGCGCCGGACGCA  >  minE/1015748‑1015884

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: