Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,618 G→C 15.0% P151A (CCG→GCG)  fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,6180GC15.0% 175.6 / 20.4 73P151A (CCG→GCG) fumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base G (36/26);  new base C (6/5);  total (42/31)
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.99e-01

AACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTT  >  minE/1015550‑1015686
                                                                    |                                                                    
aaCGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAgcggcg                                                                    >  1:4405777/1‑71 (MQ=255)
aaCGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAgcggcg                                                                    >  1:2116629/1‑71 (MQ=255)
aaCGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAgcggcg                                                                    >  1:4497177/1‑71 (MQ=255)
aaCGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAgcggcg                                                                    >  1:1751716/1‑71 (MQ=255)
aaCGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAgcggc                                                                     >  1:2999858/1‑70 (MQ=255)
     cATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGttt                                                               <  1:4754688/71‑1 (MQ=255)
      aTAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGtttt                                                              <  1:2929639/71‑1 (MQ=255)
      aTAAAGATCGATCTGCGCTGGCAAGTTGGTGCCGGTATTCACCTCTTTATACATATCCAGCGGCGCGtttt                                                              >  1:5242443/1‑71 (MQ=255)
       tAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGtttt                                                              >  1:1929062/1‑70 (MQ=255)
         aaGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                                            >  1:3881137/1‑70 (MQ=255)
         aaGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                                            >  1:1177020/1‑70 (MQ=255)
          aGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                                            >  1:4381132/1‑69 (MQ=255)
          aGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                                            >  1:497651/1‑69 (MQ=255)
            atcgatcTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGt                                                        <  1:2303055/71‑1 (MQ=255)
                     cgcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                <  1:3757661/70‑1 (MQ=255)
                     cgcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                >  1:4642723/1‑70 (MQ=255)
                     cgcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAt                                               >  1:4709590/1‑71 (MQ=255)
                      gcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                              <  1:3908554/71‑1 (MQ=255)
                      gcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                              <  1:465080/71‑1 (MQ=255)
                      gcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                                              <  1:4079060/71‑1 (MQ=255)
                          gCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCt                                          <  1:1982053/71‑1 (MQ=255)
                            aGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGttt                                                               <  1:2697873/48‑1 (MQ=255)
                            aGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAg                                                 <  1:1692001/62‑1 (MQ=255)
                               ttGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGat                                     <  1:2594335/71‑1 (MQ=255)
                                 ggTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTg                                                             >  1:4238739/1‑45 (MQ=255)
                                 agTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                    <  1:3778340/69‑1 (MQ=255)
                                  gTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGatat                                   <  1:3461270/70‑1 (MQ=255)
                                  gTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGatat                                   <  1:4143419/70‑1 (MQ=255)
                                  gTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                    >  1:5581302/1‑69 (MQ=255)
                                  gTGCCGGTGGTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGatat                                   <  1:3454570/70‑1 (MQ=255)
                                  gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                   <  1:1775278/70‑1 (MQ=255)
                                  gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                   <  1:3564154/70‑1 (MQ=255)
                                       ggTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                    >  1:534164/1‑64 (MQ=255)
                                       ggTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                    >  1:1171213/1‑64 (MQ=255)
                                       ggTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGgtgt                             <  1:4440195/71‑1 (MQ=255)
                                        gTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTc                                        <  1:2761376/59‑1 (MQ=255)
                                        gTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                            <  1:4724312/71‑1 (MQ=255)
                                          aTTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGa                                                          >  1:1563653/1‑39 (MQ=255)
                                          aTTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTa                                             >  1:171411/1‑52 (MQ=255)
                                           ttCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTtata                         <  1:5054075/71‑1 (MQ=255)
                                            tCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                            >  1:294279/1‑67 (MQ=255)
                                             cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTc                                        >  1:5493322/1‑54 (MQ=255)
                                             cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                       >  1:609669/1‑71 (MQ=255)
                                             cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                       >  1:1504845/1‑71 (MQ=255)
                                             cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                       >  1:2754096/1‑71 (MQ=255)
                                               cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCtt                                         >  1:2849303/1‑51 (MQ=255)
                                               cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGAtata                                  <  1:5037207/58‑1 (MQ=255)
                                                  cttTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcct                  >  1:3645695/1‑70 (MQ=255)
                                                  cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                >  1:5226723/1‑41 (MQ=255)
                                                  cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                >  1:5193590/1‑41 (MQ=255)
                                                  cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                >  1:5365871/1‑41 (MQ=255)
                                                  cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                >  1:5372689/1‑41 (MQ=255)
                                                  cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGa                                                >  1:5617610/1‑41 (MQ=255)
                                                  cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTAt                                            >  1:3926092/1‑45 (MQ=255)
                                                  cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTAt                                            >  1:1387125/1‑45 (MQ=255)
                                                  cttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTAt                                            >  1:1919862/1‑45 (MQ=255)
                                                     tataCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACcgc                 <  1:4172446/69‑1 (MQ=255)
                                                     tataCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACcgc                 <  1:1219466/69‑1 (MQ=255)
                                                      gtaCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                    >  1:99921/2‑49 (MQ=255)
                                                      gtaCATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata                                    >  1:3952652/2‑49 (MQ=255)
                                                       taCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCgcgcgc             >  1:4243451/1‑71 (MQ=255)
                                                        aCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGAtata                                  <  1:301726/49‑1 (MQ=255)
                                                        aCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                            <  1:3652130/55‑1 (MQ=255)
                                                        aCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                            <  1:931394/55‑1 (MQ=255)
                                                        aCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                            <  1:5163592/55‑1 (MQ=255)
                                                        aCATATCCAGCGGCGCGTTTTGCGACTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCgcgcgc             >  1:3526384/1‑70 (MQ=255)
                                                         cATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCa           <  1:4582949/71‑1 (MQ=255)
                                                            atCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTc                                        >  1:3748015/1‑39 (MQ=255)
                                                            atCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                   <  1:4993292/44‑1 (MQ=255)
                                                            atCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat                                   <  1:1678795/44‑1 (MQ=255)
                                                            atCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAgcg        >  1:3766262/1‑71 (MQ=255)
                                                             tCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATaa                                 >  1:3536829/1‑45 (MQ=255)
                                                             tCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGaca                     >  1:3149096/1‑57 (MQ=255)
                                                              ccAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgcc      >  1:1071204/1‑71 (MQ=255)
                                                              ccAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgcc      >  1:5688093/1‑71 (MQ=255)
                                                              ccAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgcc      >  1:1919440/1‑71 (MQ=255)
                                                               cAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgccg     >  1:4531729/1‑71 (MQ=255)
                                                               cAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgccg     >  1:967218/1‑71 (MQ=255)
                                                                 gcggcgCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGac                      <  1:760103/52‑1 (MQ=255)
                                                                  cggcgCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCtt  <  1:2243750/71‑1 (MQ=255)
                                                                   ggcgCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGacac                    >  1:389997/1‑52 (MQ=255)
                                                                    |                                                                    
AACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTT  >  minE/1015550‑1015686

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