Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,015,655 A→G 21.8% T138T (ACT→ACC fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,015,6550AG21.8% 99.0 / 24.8 55T138T (ACT→ACCfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (28/15);  new base G (4/8);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 9.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

GGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTAC  >  minE/1015589‑1015725
                                                                  |                                                                      
ggTGTTGACCTCTTTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGgtgt                                                                    <  1:4440195/71‑1 (MQ=255)
 gTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                                                                   <  1:4724312/71‑1 (MQ=255)
    ttCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTtata                                                                <  1:5054075/71‑1 (MQ=255)
     tCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                                                                   >  1:294279/1‑67 (MQ=255)
      cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                                                              >  1:2754096/1‑71 (MQ=255)
      cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                                                              >  1:1504845/1‑71 (MQ=255)
      cACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                                                              >  1:609669/1‑71 (MQ=255)
           cttTGTACATGTCCAGCGCCGCATTCTGTGAATAGCGCAGGTTATCTTCGATATAGGTGTTATAGACGcct                                                         >  1:3645695/1‑70 (MQ=255)
              tataCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACcgc                                                        <  1:1219466/69‑1 (MQ=255)
              tataCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACcgc                                                        <  1:4172446/69‑1 (MQ=255)
                taCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCgcgcgc                                                    >  1:4243451/1‑71 (MQ=255)
                 aCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                                                                   <  1:5163592/55‑1 (MQ=255)
                 aCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                                                                   <  1:3652130/55‑1 (MQ=255)
                 aCATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGtt                                                                   <  1:931394/55‑1 (MQ=255)
                 aCATATCCAGCGGCGCGTTTTGCGACTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCgcgcgc                                                    >  1:3526384/1‑70 (MQ=255)
                  cATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCa                                                  <  1:4582949/71‑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:5688093/1‑71 (MQ=255)
                       ccAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgcc                                             >  1:1071204/1‑71 (MQ=255)
                       ccAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgcc                                             >  1:1919440/1‑71 (MQ=255)
                        cAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgccg                                            >  1:967218/1‑71 (MQ=255)
                        cAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGcgccg                                            >  1:4531729/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)
                                   tttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcacca                                 >  1:163289/1‑71 (MQ=255)
                                   tttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcacca                                 >  1:327459/1‑71 (MQ=255)
                                    ttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTta                                                                  >  1:4225252/1‑37 (MQ=255)
                                    ttGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcaccac                                >  1:4547975/1‑71 (MQ=255)
                                      gCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGa                                                              <  1:1262096/39‑1 (MQ=255)
                                         aGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAgcgc                                              <  1:3941957/52‑1 (MQ=255)
                                          gTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCgg                          >  1:1569488/1‑71 (MQ=255)
                                              cgcAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCa                      <  1:2678046/71‑1 (MQ=255)
                                                cAGGTTATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCaca                    >  1:5239332/1‑69 (MQ=255)
                                                cAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATa                    >  1:707321/1‑71 (MQ=255)
                                                  ggttATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCacac                   >  1:2203505/3‑70 (MQ=255)
                                                   aTTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTcc                       <  1:141177/65‑1 (MQ=255)
                                                    ttATCCTCGATATAAGTGTTATAGACACCGCGCGcc                                                   >  1:2181820/1‑36 (MQ=255)
                                                    ttATCCTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcacc                                  >  1:2000193/1‑53 (MQ=255)
                                                     tATCTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACgcgc               >  1:2326450/1‑71 (MQ=255)
                                                      aTCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTtca                                       >  1:4562658/1‑46 (MQ=255)
                                                       tCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTg             >  1:3654795/1‑71 (MQ=255)
                                                        cTTCGATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATcaccgc                                <  1:963907/51‑3 (MQ=255)
                                                        cTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATcaccaccac                             <  1:2301525/54‑1 (MQ=255)
                                                            gATATAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCtt        <  1:1704682/71‑1 (MQ=255)
                                                            gATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTtcat                                      >  1:4994944/1‑41 (MQ=255)
                                                            gATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTtcat                                      >  1:5581721/1‑41 (MQ=255)
                                                            gATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCAt                     >  1:1070353/1‑58 (MQ=255)
                                                             atatAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttt       <  1:1647556/71‑1 (MQ=255)
                                                             atatAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttt       <  1:1357801/71‑1 (MQ=255)
                                                             atatAGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttt       <  1:2118962/71‑1 (MQ=255)
                                                              tataGGTGTTATAGACGCCTTTCGACAGCGTTTCTTCATCACCGCCGCCGGTCCACACGCGCTGGCCttt       <  1:1981250/70‑1 (MQ=255)
                                                               ataAGTGTTATAGACACCGCGCGCCAGCGCCGCTCCATCACCACCACCGGTCCATACACGCTGCCCttttt     >  1:4085079/1‑71 (MQ=255)
                                                                taAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCa                      <  1:4237096/53‑1 (MQ=255)
                                                                  aGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTAc  >  1:2649153/1‑71 (MQ=255)
                                                                  |                                                                      
GGTATTCACTTCTTTATACATATCCAGCGGCGCGTTTTGCGAGTAGCGCAGATTATCTTCGATATAAGTGTTATAGACACCGCGCGCCAGCGCCGCTTCATCACCACCACCGGTCCATACACGCTGCCCTTTTTTAC  >  minE/1015589‑1015725

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: