Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,014,764 A→T 22.2% S435S (TCT→TCA fumA ← aerobic Class I fumarate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,014,7640AT22.2% 71.1 / 17.6 36S435S (TCT→TCAfumAaerobic Class I fumarate hydratase
Reads supporting (aligned to +/- strand):  ref base A (11/17);  new base T (4/4);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

TCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATG  >  minE/1014701‑1014830
                                                               |                                                                  
tCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCAGGCGGTCGTTGGGCCAAGAGAAcc                                                               >  1:474120/1‑69 (MQ=255)
 ccGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCg                                                              >  1:960034/1‑69 (MQ=255)
 ccGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCg                                                              >  1:1328016/1‑69 (MQ=255)
   gCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAgg                                                          <  1:3066225/71‑1 (MQ=255)
     ccTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAg                                                           >  1:1075183/1‑68 (MQ=255)
       tGCCCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCAt                                                       >  1:463026/1‑70 (MQ=255)
         cgcTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATa                                                      <  1:1936669/69‑1 (MQ=255)
              gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCa                                                        >  1:684540/1‑62 (MQ=255)
              gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCtt                                                 <  1:2913606/69‑1 (MQ=255)
              gCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTc                                                <  1:2428859/70‑1 (MQ=255)
               cAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCg                                              <  1:3118911/71‑1 (MQ=255)
                 gTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGc                                            >  1:1241017/1‑71 (MQ=255)
                     aTCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGttt                                        <  1:3748408/71‑1 (MQ=255)
                     aTCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGttt                                        <  1:2785296/71‑1 (MQ=255)
                      tCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCt                                                  <  1:467292/60‑1 (MQ=255)
                        gACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTg                                      <  1:2334704/70‑1 (MQ=255)
                            tAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCgg                                 <  1:3038719/71‑1 (MQ=255)
                                 atccatccGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCcg                            <  1:2553805/71‑1 (MQ=255)
                                  tccatccGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCgg                                 <  1:1531451/65‑1 (MQ=255)
                                    catacGGCCTGCGGTGGTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCcgcg                          <  1:1982954/70‑1 (MQ=255)
                                     atccGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCcgcg                          >  1:2802443/1‑69 (MQ=255)
                                            cggcggTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAAt                  <  1:1287834/70‑1 (MQ=255)
                                            cggcggTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATc                 >  1:188493/1‑71 (MQ=255)
                                                ggTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGgat             >  1:3787772/1‑71 (MQ=255)
                                                 gtggTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGggtg            <  1:2527792/68‑1 (MQ=255)
                                                 gtggTTGGGCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGggtg            <  1:2136927/68‑1 (MQ=255)
                                                      tGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCtt       >  1:2650276/1‑71 (MQ=255)
                                                       gggCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCttt      >  1:543991/1‑71 (MQ=255)
                                                       gggCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCTTAGTAAATCGGATGATCtt       <  1:2035163/70‑1 (MQ=255)
                                                       gggCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCGGACCCGCGTAGTAAATCGGATGATCttt      <  1:2515173/70‑1 (MQ=255)
                                                        ggCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGggtgat          >  1:1382077/1‑66 (MQ=255)
                                                        ggCCAAGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGATCGGGTGATCTTTg     <  1:2918340/71‑1 (MQ=255)
                                                        ggCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCttt      <  1:230616/70‑1 (MQ=255)
                                                           cAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATg  <  1:198368/71‑1 (MQ=255)
                                                            aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCgtagtaga                   >  1:1545684/1‑51 (MQ=255)
                                                            aaGTGAACCTGATGGATAACCGGCAGGGGTTTTCGCCGGACCCGCGTAGTAGAt                  >  1:882680/1‑54 (MQ=255)
                                                               |                                                                  
TCCGCCCTGCGCTTGCAGTTGATCGACATAAGAATCCATCCGTCCGGCGGTCGTTGGGCCAAGAGAACCGGAGGCATAACCTTCCGGCGTTTTGGCCGGACCCGCGTAGTAAATCGGATGATCTTTGATG  >  minE/1014701‑1014830

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: