Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,785 C→A 23.3% P428P (CCG→CCT fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,7850CA23.3% 103.8 / 17.4 47P428P (CCG→CCTfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base C (22/14);  new base A (8/3);  total (30/17)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGT  >  minE/1014715‑1014852
                                                                      |                                                                   
gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTcc                                                                     >  1:2276558/1‑71 (MQ=255)
      gATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGt                                                                <  1:2935683/70‑1 (MQ=255)
        tCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGttt                                                              >  1:116620/1‑70 (MQ=255)
           aCATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGc                                                          <  1:2207995/71‑1 (MQ=255)
             aTAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCg                                                        >  1:1045288/1‑71 (MQ=255)
               aaGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGa                                                      >  1:3023545/1‑71 (MQ=255)
                   atccatccGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGAccc                                                   >  1:2789702/1‑70 (MQ=255)
                        tccgtccgGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGttt                                                              <  1:3929613/54‑1 (MQ=255)
                          cgtcccGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCgtagta                                           <  1:3608595/71‑1 (MQ=255)
                              ctgcggTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgtag                                             >  1:2669235/3‑65 (MQ=255)
                               tgcggTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgtag                                             >  1:3835326/2‑64 (MQ=255)
                               ggcggTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCg                                      <  1:1963437/71‑1 (MQ=255)
                                 cggTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCggg                                    >  1:3979551/1‑70 (MQ=255)
                                  ggTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGgat                                   <  1:2927941/71‑1 (MQ=255)
                                   gTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGgatg                                  <  1:1281739/71‑1 (MQ=255)
                                         gggCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCttt                            >  1:2622623/1‑71 (MQ=255)
                                         gggCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCtt                             >  1:82283/1‑70 (MQ=255)
                                          ggCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt                            <  1:549856/70‑1 (MQ=255)
                                          ggCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt                            <  1:3307703/70‑1 (MQ=255)
                                            ccAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTaaa                                         <  1:1415740/55‑1 (MQ=255)
                                            ccAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTaaa                                         <  1:830660/55‑1 (MQ=255)
                                             cAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATg                        <  1:5629680/71‑1 (MQ=255)
                                              aaGTGAACCTGATGGATAACCGGCAGGTGTTTTCGCCGGACCcgc                                                 <  1:1037267/45‑1 (MQ=255)
                                              aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATg                        >  1:2722137/1‑70 (MQ=255)
                                              aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATg                        >  1:2967144/1‑70 (MQ=255)
                                              aaGAGAACGGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATc                               <  1:961661/63‑1 (MQ=255)
                                                gtgaACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgtag                                             <  1:1254013/45‑1 (MQ=255)
                                                gtgaACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgta                                              >  1:448847/3‑46 (MQ=255)
                                                gagaACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGga                                    >  1:4586304/1‑56 (MQ=255)
                                                gagaACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTAc                     >  1:47173/1‑71 (MQ=255)
                                                 agaACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt                            <  1:770335/63‑1 (MQ=255)
                                                  gaACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTGATGTACTg                   >  1:5638293/1‑71 (MQ=255)
                                                   aaCCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTg                           <  1:4151197/62‑1 (MQ=255)
                                                       ggAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                          >  1:86848/1‑59 (MQ=255)
                                                       ggAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                          >  1:590103/1‑59 (MQ=255)
                                                        gAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAg             >  1:808093/1‑71 (MQ=255)
                                                             aTAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTaa                                          >  1:3987559/1‑37 (MQ=255)
                                                               aaCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGgatga                                 >  1:314085/1‑44 (MQ=255)
                                                               aaCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCg                 >  1:3455482/1‑60 (MQ=255)
                                                               aaCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCg                 >  1:5326622/1‑60 (MQ=255)
                                                                aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGa                          >  1:1864890/1‑50 (MQ=255)
                                                                aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCa      >  1:5414016/1‑70 (MQ=255)
                                                                aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCa      >  1:3814459/1‑70 (MQ=255)
                                                                aCCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAc     >  1:101591/1‑71 (MQ=255)
                                                                 ccTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAcc    >  1:4527079/1‑71 (MQ=255)
                                                                    tCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCtt                             >  1:549322/1‑43 (MQ=255)
                                                                    tCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCAcc    >  1:2011655/1‑68 (MQ=255)
                                                                    tCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGt  >  1:3026398/1‑70 (MQ=255)
                                                                      cGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGccc          <  1:2412323/60‑1 (MQ=255)
                                                                      |                                                                   
GCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATGTACTGCGGCAGCCCTTCACCGT  >  minE/1014715‑1014852

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: