Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 262,850 G→A 31.8% L309L (CTG→CTA ribD → fused diaminohydroxyphosphoribosylaminopyrimidin e deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE262,8500GA31.8% 71.7 / 38.9 82L309L (CTG→CTAribDfused diaminohydroxyphosphoribosylaminopyrimidin e deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base G (24/31);  new base A (17/9);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 9.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACT  >  minE/262791‑262897
                                                           |                                               
aaaCAGCAAATTAACAGCATCTGTGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCa                                               <  1:1918739/62‑1 (MQ=255)
aaaCAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCa                                               <  1:2482233/62‑1 (MQ=255)
aaaCAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCa                                               <  1:1686007/62‑1 (MQ=255)
aaaCAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCa                                               <  1:3121066/62‑1 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACa                                               >  1:2794476/1‑61 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:2988858/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:2160468/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:2274235/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:2372531/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:1909205/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:265616/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:2902177/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:1608415/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:3024509/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:3340768/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:1088122/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:1108819/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:668389/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:638462/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCggg                                         >  1:143199/1‑67 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTACAGGCGGGt                                        >  1:651894/1‑68 (MQ=255)
 aaCAGCAGATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCT‑‑AGGCGGGTTTa                                     >  1:1559040/1‑69 (MQ=255)
    aGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTa                                     >  1:2955053/1‑68 (MQ=255)
    aGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTa                                     >  1:2025239/1‑68 (MQ=255)
    aGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAg                                    >  1:3114066/1‑69 (MQ=255)
    aGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAg                                    >  1:1406440/1‑69 (MQ=255)
    aGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAg                                    >  1:1273666/1‑69 (MQ=255)
    aGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAg                                    >  1:3038852/1‑69 (MQ=255)
    aGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAg                                    >  1:258275/1‑69 (MQ=255)
    aGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAg                                    >  1:1946240/1‑69 (MQ=255)
    aGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAg                                    >  1:1451544/1‑69 (MQ=255)
    aGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGAAGGCGGGTTTa                                     >  1:1194100/1‑68 (MQ=255)
    aGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGATGGCGCATTGCTGCAGGCGGGTTTAg                                    >  1:967922/1‑69 (MQ=255)
      cAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTc                                  >  1:1027071/1‑69 (MQ=255)
      cAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTc                                  >  1:2939255/1‑69 (MQ=255)
      cAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTc                                  >  1:2894821/1‑69 (MQ=255)
      cAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTc                                  >  1:2942993/1‑69 (MQ=255)
            aaCAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAg                            >  1:3362874/1‑69 (MQ=255)
            aaCAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAg                            >  1:3366739/1‑69 (MQ=255)
                     tGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGt                    >  1:2549959/1‑68 (MQ=255)
                     tGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTc                   >  1:608929/1‑69 (MQ=255)
                     tGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTc                   >  1:2905339/1‑69 (MQ=255)
                         tGGAAGCGGGGCCAACGCTCGCTGGGGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCtat                <  1:1424532/68‑1 (MQ=255)
                         tGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGTTTTAGTCGATGAGCTGATTGTCtat                <  1:2371849/68‑1 (MQ=255)
                         tGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCtat                <  1:3344401/68‑1 (MQ=255)
                         tGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCtat                <  1:1816036/68‑1 (MQ=255)
                         tGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCtat                <  1:2345647/68‑1 (MQ=255)
                         tGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCtat                <  1:947365/68‑1 (MQ=255)
                         tGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCta                 <  1:449569/67‑1 (MQ=255)
                         tGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCta                 <  1:22114/67‑1 (MQ=255)
                         tGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCta                 <  1:3331287/67‑1 (MQ=255)
                         tGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCta                 <  1:1484344/67‑1 (MQ=255)
                         tGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCta                 <  1:1532352/67‑1 (MQ=255)
                         tGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCta                 <  1:162401/67‑1 (MQ=255)
                         tGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCta                 <  1:1848436/67‑1 (MQ=255)
                         tGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCta                 <  1:1699672/67‑1 (MQ=255)
                           gAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCta                 <  1:2876171/65‑1 (MQ=255)
                           gAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCta                 <  1:1927796/65‑1 (MQ=255)
                           gAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCta                 <  1:2543174/65‑1 (MQ=255)
                           gAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCta                 <  1:277336/65‑1 (MQ=255)
                                ggggCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCAc         >  1:717251/1‑68 (MQ=255)
                                ggggCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCAc         >  1:2979676/1‑68 (MQ=255)
                                ggggCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCAc         >  1:2857128/1‑68 (MQ=255)
                                ggggCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCAc         >  1:252042/1‑68 (MQ=255)
                                     cAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAAc   <  1:2753949/69‑1 (MQ=255)
                                     cAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAAc   <  1:2214558/69‑1 (MQ=255)
                                     cAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAAc   <  1:817031/69‑1 (MQ=255)
                                     cAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAAc   <  1:1117329/69‑1 (MQ=255)
                                     cAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAAc   <  1:1167152/69‑1 (MQ=255)
                                     cAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAAc   <  1:3431146/69‑1 (MQ=255)
                                     cAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAAc   <  1:2008599/69‑1 (MQ=255)
                                     cAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAAc   <  1:1189857/69‑1 (MQ=255)
                                      aaCGCTCGCTGGCGCATTGCTACAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACt  <  1:2069247/69‑1 (MQ=255)
                                      aaCGCTCGCTGGCGCATTGCTACAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACt  <  1:3001523/69‑1 (MQ=255)
                                      aaCGCTCGCTGGCGCATTGCTACAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACt  <  1:2012635/69‑1 (MQ=255)
                                      aaCGCTCGCTGGCGCATTGCTACAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACt  <  1:570357/69‑1 (MQ=255)
                                      aaCGCTCGCTGGCGCATTGCTACAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACt  <  1:1657977/69‑1 (MQ=255)
                                      aaCGCTCGCTGGCGCATTGCTACAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACt  <  1:3138237/69‑1 (MQ=255)
                                      aaCGCTCGCTGGCGCATTGCTACAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACt  <  1:1409441/69‑1 (MQ=255)
                                      aaCGCTCGCTGGCGCATTGCTACAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACt  <  1:1466391/69‑1 (MQ=255)
                                                  cgcATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCtat                <  1:3189876/43‑1 (MQ=255)
                                                     aTTGCTACAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACt  <  1:1262691/54‑1 (MQ=255)
                                                           |                                               
AAACAGCAAATTAACAGCATCTGGGTGGAAGCGGGGCCAACGCTCGCTGGCGCATTGCTGCAGGCGGGTTTAGTCGATGAGCTGATTGTCTATATCGCACCTAAACT  >  minE/262791‑262897

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: