Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,598 T→C 16.4% E157E (GAA→GAG fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,5980TC16.4% 134.6 / 18.4 61E157E (GAA→GAGfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base T (32/19);  new base C (4/6);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGA  >  minE/1015528‑1015665
                                                                      |                                                                   
ggAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCAct                                                                     <  1:3024302/71‑1 (MQ=255)
    tttGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCActtctt                                                                 <  1:4432141/71‑1 (MQ=255)
     ttGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCActtctt                                                                 >  1:4243759/1‑70 (MQ=255)
     ttGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCActtctt                                                                 >  1:4250465/1‑70 (MQ=255)
        tATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTtata                                                             >  1:142132/1‑71 (MQ=255)
        tATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTtat                                                              <  1:5629182/70‑1 (MQ=255)
         aTTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTtata                                                             >  1:4939793/1‑70 (MQ=255)
           tcgtcgCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACa                                                           >  1:2891749/1‑70 (MQ=255)
             gtcgCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACatat                                                        <  1:1752028/71‑1 (MQ=255)
                   aTCAACGGCATAATGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAgcg                                                  <  1:963476/71‑1 (MQ=255)
                      aaCGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAgcggcg                                               >  1:1751716/1‑71 (MQ=255)
                      aaCGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAgcggcg                                               >  1:2116629/1‑71 (MQ=255)
                      aaCGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAgcggcg                                               >  1:4497177/1‑71 (MQ=255)
                      aaCGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAgcggcg                                               >  1:4405777/1‑71 (MQ=255)
                      aaCGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAgcggc                                                >  1:2999858/1‑70 (MQ=255)
                           cATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGttt                                          <  1:4754688/71‑1 (MQ=255)
                            aTAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTta                                                               >  1:5240198/1‑49 (MQ=255)
                            aTAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTta                                                               >  1:789653/1‑49 (MQ=255)
                            aTAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGtttt                                         <  1:2929639/71‑1 (MQ=255)
                            aTAAAGATCGATCTGCGCTGGCAAGTTGGTGCCGGTATTCACCTCTTTATACATATCCAGCGGCGCGtttt                                         >  1:5242443/1‑71 (MQ=255)
                             tAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGtttt                                         >  1:1929062/1‑70 (MQ=255)
                               aaGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACatat                                                        >  1:1238145/1‑53 (MQ=255)
                               aaGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                       >  1:3881137/1‑70 (MQ=255)
                               aaGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                       >  1:1177020/1‑70 (MQ=255)
                                aGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                       >  1:497651/1‑69 (MQ=255)
                                aGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGc                                       >  1:4381132/1‑69 (MQ=255)
                                  atcgatcTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGt                                   <  1:2303055/71‑1 (MQ=255)
                                    cgatcTGCGCTGGCAGATTGGTGCCGGTATTCACTTCttt                                                                >  1:2619794/1‑40 (MQ=255)
                                    cgatcTGCGCTGGCAGATTGGTGCCGGTATTCACTTCttt                                                                >  1:443631/1‑40 (MQ=255)
                                    cgatcTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTta                                                               >  1:253394/1‑41 (MQ=255)
                                    cgatcTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTta                                                               >  1:5604258/1‑41 (MQ=255)
                                        cTGCGCTGGCAAGTTGGTGCCGGTATTCACCTCttt                                                                <  1:783197/36‑1 (MQ=25)
                                          gcgcTGGCAGATTGGTGCCGGTATTCACTTCTTtata                                                             >  1:2252904/1‑37 (MQ=255)
                                          gcgcTGGCAGATTGGTGCCGGTATTCACTTCTTtata                                                             >  1:3340704/1‑37 (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:465080/71‑1 (MQ=255)
                                            gcTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGAtt                         <  1:3908554/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:4143419/70‑1 (MQ=255)
                                                        gTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGatat              <  1:3461270/70‑1 (MQ=255)
                                                        gTGCCGGTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGata               >  1:5581302/1‑69 (MQ=255)
                                                        gTGCCGGTGGTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGatat              <  1:3454570/70‑1 (MQ=255)
                                                        gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat              <  1:3564154/70‑1 (MQ=255)
                                                        gTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGatat              <  1:1775278/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:2754096/1‑71 (MQ=255)
                                                                   cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa  >  1:609669/1‑71 (MQ=255)
                                                                   cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa  >  1:1504845/1‑71 (MQ=255)
                                                                     cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCtt                    >  1:2849303/1‑51 (MQ=255)
                                                                     cttcttTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGAtata             <  1:5037207/58‑1 (MQ=255)
                                                                      |                                                                   
GGAATTTGTATTCGTCGCCATCAACGGCATAAAGATCGATCTGCGCTGGCAGATTGGTGCCGGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGA  >  minE/1015528‑1015665

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