Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,628 C→T 6.7% S147S (TCG→TCA fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,6280CT 6.7% 240.8 / 5.4 74S147S (TCG→TCAfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (36/33);  new base T (3/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

AAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCAC  >  minE/1015559‑1015698
                                                                     |                                                                      
aaGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                                                        >  1:821474/1‑70 (MQ=255)
aaGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                                                        >  1:3794899/1‑70 (MQ=255)
aaGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                                                        >  1:2418447/1‑70 (MQ=255)
aaGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                                                        <  1:3462149/70‑1 (MQ=255)
  gatcgatcTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAg                                                                     <  1:2445750/71‑1 (MQ=255)
    tcgatcTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAg                                                                     >  1:2144326/1‑69 (MQ=255)
         cTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAgcg                                                                >  1:1017995/1‑69 (MQ=255)
            cgcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAgc                                                                 <  1:1992101/65‑1 (MQ=255)
            cgcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCa                                                              <  1:914653/68‑1 (MQ=255)
            cgcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                            >  1:1831867/1‑70 (MQ=255)
             gcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                          <  1:2949303/71‑1 (MQ=255)
             gcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                          <  1:1557821/71‑1 (MQ=255)
                ggCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATc                                                       >  1:3230298/1‑71 (MQ=255)
                         gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                               <  1:1561873/70‑1 (MQ=255)
                         gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                               <  1:2029589/70‑1 (MQ=255)
                         gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                               <  1:592169/70‑1 (MQ=255)
                         gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGAtata                                              <  1:329373/71‑1 (MQ=255)
                         gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGAtata                                              <  1:4114088/71‑1 (MQ=255)
                          tGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAt                                                           >  1:2106549/1‑57 (MQ=255)
                           gCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                            >  1:2795268/1‑55 (MQ=255)
                           gCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGa                                                  >  1:2969842/1‑65 (MQ=255)
                            ccGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTTTGAATAGCGCAGGTTATCTTCGATATAGgt                                           <  1:3211866/71‑1 (MQ=255)
                            ccGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                                >  1:891835/1‑66 (MQ=255)
                            ccGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                                >  1:2358300/1‑66 (MQ=255)
                            ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                                                        <  1:747230/42‑1 (MQ=255)
                            ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTAt                                                        >  1:1658889/1‑58 (MQ=255)
                            ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                           <  1:3399309/71‑1 (MQ=255)
                            ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                           <  1:3687280/71‑1 (MQ=255)
                            ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                           >  1:4129317/1‑71 (MQ=255)
                            ccGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgt                                           <  1:3889355/71‑1 (MQ=255)
                              ggTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgtgt                                         <  1:3894255/71‑1 (MQ=255)
                               gTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgtgt                                         >  1:2778410/1‑70 (MQ=255)
                                tATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                            >  1:1052967/1‑50 (MQ=255)
                                tATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                       <  1:3613026/71‑1 (MQ=255)
                                  ttCACTTCTTTATACATATCGAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCg                                                   <  1:756792/57‑1 (MQ=255)
                                  ttCACTTCTTTATACATATCCAGCGGCGCGTTTTGCg                                                                       <  1:33198/37‑1 (MQ=255)
                                  ttCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                            >  1:2574574/1‑48 (MQ=255)
                                  ttCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                               <  1:464675/61‑1 (MQ=255)
                                  ttCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgtgt                                         <  1:3743346/67‑1 (MQ=255)
                                  ttCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTtata                                     <  1:1575410/71‑1 (MQ=255)
                                    cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAg                                                                  >  1:1727413/1‑40 (MQ=255)
                                    cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAg                                                                  >  1:3101627/1‑40 (MQ=255)
                                    cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                       >  1:2127447/1‑67 (MQ=255)
                                    cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                       >  1:2090403/1‑67 (MQ=255)
                                    cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                                   >  1:2286823/1‑71 (MQ=255)
                                     aCTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGac                                  <  1:2048024/71‑1 (MQ=255)
                                      cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTa                                                                   >  1:2909165/1‑37 (MQ=255)
                                      cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCTATATAAgt                                           <  1:3722959/61‑1 (MQ=255)
                                      cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                       <  1:1927646/65‑1 (MQ=255)
                                       ttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAg                                                             <  1:3751995/42‑1 (MQ=255)
                                       ttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGa                                                  <  1:1069832/53‑1 (MQ=255)
                                       ttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgtgt                                         >  1:63227/1‑62 (MQ=255)
                                         cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAt                                                           >  1:653114/1‑42 (MQ=255)
                                         cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAt                                                           >  1:3850327/1‑42 (MQ=255)
                                         cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAt                                                           >  1:2928814/1‑42 (MQ=255)
                                         cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATATTCGATATAAGTGTTATAg                                    >  1:1107220/1‑65 (MQ=255)
                                          tttATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACcgc                             >  1:566649/1‑71 (MQ=255)
                                            tataCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                          >  1:2825859/1‑40 (MQ=255)
                                            tataCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTAt                                                        >  1:1710526/1‑42 (MQ=255)
                                            tataCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACcgcgc                           >  1:669124/1‑71 (MQ=255)
                                              taCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                          <  1:964407/38‑1 (MQ=255)
                                              taCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                                          <  1:1978949/38‑1 (MQ=255)
                                              taCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCgcgcgc                         <  1:1374200/71‑1 (MQ=255)
                                              taCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCgcgcgc                         <  1:2887829/71‑1 (MQ=255)
                                                  tatCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAgtgt                                         <  1:3827516/51‑1 (MQ=255)
                                                  tatCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAgc                     >  1:1916275/1‑71 (MQ=255)
                                                   atCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCtt                                                     >  1:570182/1‑38 (MQ=255)
                                                    tCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATaa                                             >  1:2666517/1‑45 (MQ=255)
                                                    tCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                       <  1:1341398/51‑1 (MQ=255)
                                                     ccAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgcc                  >  1:4017755/1‑71 (MQ=255)
                                                      cAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                                >  1:3736814/1‑40 (MQ=16)
                                                          ggcgCGTTTTGCGAGTAGCGCAGATTATCTTCGAtata                                              >  1:3346754/1‑38 (MQ=255)
                                                          ggcgCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAg                      >  1:3534261/1‑62 (MQ=255)
                                                           ccgcaTTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGtt                                        <  1:227307/42‑1 (MQ=16)
                                                            cgcgTTTTGCGAGTAGCGCAGATTATCTTCGAtata                                              <  1:2691694/36‑1 (MQ=255)
                                                             gcgTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTtcatt          <  1:3193000/71‑2 (MQ=255)
                                                                ttttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcac        >  1:898395/1‑70 (MQ=255)
                                                                 tttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                                   >  1:903594/1‑42 (MQ=255)
                                                                 tttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcacca      >  1:3572308/1‑71 (MQ=255)
                                                                     cGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcaccaccac  <  1:490723/71‑1 (MQ=255)
                                                                     |                                                                      
AAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCAC  >  minE/1015559‑1015698

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: