Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,798 G→A 17.4% L91L (CTG→TTG) ‡ fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,7980GA17.4% 130.2 / 18.4 58L91L (CTG→TTG) ‡fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (21/26);  new base A (7/3);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGT  >  minE/1015732‑1015868
                                                                  |                                                                      
taatCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAt                                                                    >  1:3900731/1‑71 (MQ=255)
  atCGCGGTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTg                                                                  >  1:3847882/1‑71 (MQ=255)
     gcgGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTg                                                                  >  1:79115/1‑68 (MQ=255)
     gcgGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCa                                                                <  1:361761/70‑1 (MQ=255)
       ggTGCCGGTATCCTGGCAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAgcg                                                             <  1:800964/71‑1 (MQ=255)
           ccGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCca                                                          >  1:2654439/1‑70 (MQ=255)
           ccGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCca                                                          >  1:3525072/1‑70 (MQ=255)
           ccGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCca                                                          >  1:294181/1‑70 (MQ=255)
               tATCCTTACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTg                                                                  <  1:2764797/58‑1 (MQ=255)
                aTCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACatat                                                     >  1:1576625/1‑70 (MQ=255)
                      gcAGGTCGGCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTcgt                                              <  1:1186741/70‑3 (MQ=255)
                        aGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGaa                                                                     >  1:923204/1‑46 (MQ=255)
                        aGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCAtt                                             >  1:1658305/1‑70 (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:476506/71‑1 (MQ=255)
                           ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg                                         <  1:2719413/71‑1 (MQ=255)
                           ttGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCg                                         <  1:2065601/71‑1 (MQ=255)
                           tcggCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCg                                         >  1:1036365/3‑71 (MQ=255)
                           tcggCAGCACGCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCg                                         <  1:2486585/69‑1 (MQ=255)
                            tGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGaa                                                                     >  1:3794479/1‑42 (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)
                                    cGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAt                                                                    >  1:1439749/1‑35 (MQ=255)
                                    cGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAAt                                                                    >  1:1677593/1‑35 (MQ=255)
                                     gCCTTTGGCGGCGATTTCGGAGTTTCTTAAGAATTGCAGCGCCACGTACTTGTCGTTTTCGCTGGCTTctg                               >  1:3846492/1‑69 (MQ=255)
                                     gCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCg                               <  1:1622153/71‑1 (MQ=255)
                                      ccTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCgg                              <  1:3641654/71‑1 (MQ=255)
                                      ccTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCgg                              <  1:2916760/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:284531/71‑1 (MQ=255)
                                               gcgATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAg                     <  1:1130330/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:2103028/1‑71 (MQ=255)
                                                        gAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGc            <  1:2280586/71‑1 (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)
                                                                  |                                                                      
TAATCGCGGTGCCGGTATCCTGACAGGTTGGCAGAACGCCTTTCGCCGCGATGTCGGAGTTACGCAGGAATTGCAGCGCCACATATTTATCATTTTCGCTGGCCTCCGGGTCACGCAGAATGTCGGCCACCTGTTGT  >  minE/1015732‑1015868

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: