Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 163,666 C→T 42.4% E77K (GAA→AAA)  yaeI ← predicted phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE163,6660CT42.4% 32.2 / 72.9 59E77K (GAA→AAA) yaeIpredicted phosphatase
Reads supporting (aligned to +/- strand):  ref base C (19/15);  new base T (9/16);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

GCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGGAGATCGGC  >  minE/163596‑163732
                                                                      |                                                                  
gCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTc                                                                    <  1:2729485/71‑1 (MQ=255)
gCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTc                                                                    <  1:2046455/71‑1 (MQ=255)
      aaaTTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATg                                                                >  1:877810/1‑69 (MQ=255)
      aaaTTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGc                                                               >  1:3781445/1‑70 (MQ=255)
                  atatCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGc                                                        <  1:1093407/65‑1 (MQ=255)
                  atatCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCa                                                       <  1:2079240/66‑1 (MQ=255)
                     tCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGt                                                >  1:1657799/1‑70 (MQ=255)
                     tCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGt                                                >  1:2845038/1‑70 (MQ=255)
                     tCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGt                                                >  1:1801247/1‑70 (MQ=255)
                       aaaTAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCt                                                                   >  1:2721165/1‑49 (MQ=255)
                       aaaTAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCt                                                                   >  1:2449704/1‑49 (MQ=255)
                       aaaTAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCt                                                                   >  1:2340727/1‑49 (MQ=255)
                       aaaTAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCt                                                                   >  1:1601359/1‑49 (MQ=255)
                       aaaTAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCt                                                                   >  1:4237282/1‑49 (MQ=255)
                       aaaTAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCt                                                                   >  1:523374/1‑49 (MQ=255)
                         ataataCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTc                                                                    >  1:1759603/1‑46 (MQ=255)
                         ataataCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTc                                                                    >  1:3051838/1‑46 (MQ=255)
                         ataataCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTc                                                                    >  1:3047118/1‑46 (MQ=255)
                         ataataCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGa                                            >  1:3527563/1‑70 (MQ=255)
                            ataCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATc                                        <  1:3231427/71‑1 (MQ=255)
                             taCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCaa                                                            <  1:2800449/50‑1 (MQ=255)
                              aCGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAAt                                         <  1:1799712/68‑1 (MQ=255)
                                    tCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTaa                                <  1:3195542/71‑1 (MQ=255)
                                    tCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTaa                                <  1:2460132/71‑1 (MQ=255)
                                    tCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTaa                                <  1:901655/71‑1 (MQ=255)
                                               tAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAAc                     <  1:3961737/71‑1 (MQ=255)
                                               tAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAAc                     <  1:970834/71‑1 (MQ=255)
                                                 atatCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAAc                     >  1:4253260/1‑69 (MQ=255)
                                                 atatCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACga                   <  1:2894621/71‑1 (MQ=255)
                                                 atatCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACga                   <  1:1406029/71‑1 (MQ=255)
                                                 atatCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACga                   <  1:1437827/71‑1 (MQ=255)
                                                 atatCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTaaa                               >  1:3933596/1‑59 (MQ=255)
                                                 atatCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAAc                     >  1:3561723/1‑69 (MQ=255)
                                                 atatCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAAc                     >  1:254177/1‑69 (MQ=255)
                                                  tatCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACaa                        <  1:4168523/65‑1 (MQ=255)
                                                   atCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGc                                    >  1:1324614/1‑52 (MQ=255)
                                                    tCAAATCGTGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTaa                                <  1:3699094/55‑1 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAAt                                         <  1:1756047/46‑1 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAAt                                         <  1:4009647/46‑1 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAAt                                         <  1:332645/46‑1 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAAt                                         >  1:3362512/1‑46 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAAt                                         <  1:348578/46‑1 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAAt                                         <  1:1494328/46‑1 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTaaa                               >  1:3446083/1‑56 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATcgcg                                                 >  1:4180364/1‑38 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATcgcg                                                 >  1:3670435/1‑38 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGc                                    >  1:778781/1‑51 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTaa                                <  1:2619282/55‑1 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTaa                                <  1:1544940/55‑1 (MQ=255)
                                                    tCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTaa                                <  1:1604181/55‑1 (MQ=255)
                                                     cAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGagag                >  1:3876348/1‑70 (MQ=255)
                                                      aaaTCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTa              <  1:1438397/71‑1 (MQ=255)
                                                      aaaTCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTa              <  1:4217904/71‑1 (MQ=255)
                                                       aaTCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACga                   <  1:3711240/65‑1 (MQ=255)
                                                        aTCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAAt            >  1:3563221/1‑71 (MQ=255)
                                                        aTCGGGCTTTTGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAAt            >  1:2669898/1‑71 (MQ=255)
                                                           gtgCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGga         <  1:3255875/69‑1 (MQ=255)
                                                           gggCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGga         <  1:1912181/71‑1 (MQ=255)
                                                           gggCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGga         <  1:1345548/71‑1 (MQ=255)
                                                            ggCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGgag        <  1:3138406/71‑1 (MQ=255)
                                                                  tGTTTTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGGAGATCGGc  >  1:324021/1‑71 (MQ=255)
                                                                  tGTTTTATGCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGGAGATCGGc  >  1:3035844/1‑70 (MQ=255)
                                                                  tGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGGAGATCGGc  >  1:2158575/1‑71 (MQ=255)
                                                                    ttCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAAt            <  1:1361478/59‑1 (MQ=255)
                                                                      |                                                                  
GCCGAAAAATTCAGCGACATATCAAATAATACGTAATCGCCGCCCAGTAATATCAAATCGGGCTTTTGTTCTATGCCAAGAGCAATCGCGTCAGAAATCAGGCTTAAAGGAACAAAACGAGAGTAATGGAGATCGGC  >  minE/163596‑163732

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: