Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 775,502 A→G 22.2% L375L (TTA→TTG) ‡ icd → isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE775,5020AG22.2% 101.6 / 29.0 54L375L (TTA→TTG) ‡icdisocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base A (22/20);  new base G (4/8);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

TCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGAT  >  minE/775433‑775571
                                                                     |                                                                     
tCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTaa                                                                      <  1:3705735/71‑1 (MQ=255)
      ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTa                                                                <  1:1579757/71‑1 (MQ=255)
      ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTa                                                                <  1:205607/71‑1 (MQ=255)
           ttattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTa                                                                <  1:1687272/66‑1 (MQ=255)
           ttattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                            >  1:88173/1‑70 (MQ=255)
           ttattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGt                                                           <  1:3123752/71‑1 (MQ=255)
                ctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATgg                                                       <  1:2526767/70‑1 (MQ=255)
                ctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATgg                                                       >  1:2374555/1‑70 (MQ=255)
                 tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATgg                                                       >  1:2554669/1‑69 (MQ=255)
                 tctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATgg                                                       >  1:3071760/1‑69 (MQ=255)
                  ctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAg                                                    >  1:958243/1‑71 (MQ=255)
                  ctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAg                                                    >  1:1019198/1‑71 (MQ=255)
                     cGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAgg                                                            <  1:3606138/60‑1 (MQ=255)
                     cGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGTTATGGAAGgcg                                                 >  1:3354213/1‑71 (MQ=255)
                     cGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGgc                                                  >  1:3228727/1‑70 (MQ=255)
                      gCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGgcgc                                                <  1:3728678/71‑1 (MQ=255)
                            ctaatgCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaa                                          <  1:899503/68‑1 (MQ=255)
                             taatgCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaa                                          >  1:1473075/3‑70 (MQ=255)
                              gatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcg                                        >  1:1728706/1‑71 (MQ=255)
                               atgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaaa                                                              >  1:2474042/1‑48 (MQ=255)
                                  ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                                              >  1:3047840/1‑61 (MQ=255)
                                  ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaa                                    <  1:78570/71‑1 (MQ=255)
                                  ctgcGCCACATGGGTTGGACCGAAGC‑‑CTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa                                  >  1:3215267/1‑71 (MQ=255)
                                    gcgcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaatc                                  >  1:2342608/1‑69 (MQ=255)
                                     cgcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa                                  <  1:2309365/70‑1 (MQ=255)
                                          taTGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCaa                                    <  1:789745/62‑1 (MQ=255)
                                            tGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATc                                            <  1:728853/53‑1 (MQ=38)
                                            tGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAct                          <  1:3221389/71‑2 (MQ=255)
                                              ggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                                         >  1:1775570/1‑54 (MQ=255)
                                              ggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                                         >  1:3253117/1‑54 (MQ=255)
                                               gTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTAt                       <  1:4015745/71‑1 (MQ=255)
                                               gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt                       <  1:1501705/71‑1 (MQ=255)
                                                ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa                                  >  1:3745480/1‑59 (MQ=255)
                                                ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCt                         <  1:3811425/68‑1 (MQ=255)
                                                ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg                      <  1:3761392/71‑1 (MQ=255)
                                                ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg                      <  1:176514/71‑1 (MQ=255)
                                                   gACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGgcgc                                                <  1:3182638/42‑1 (MQ=255)
                                                    aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGgcg                                                 <  1:190940/40‑1 (MQ=255)
                                                    aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACtt                  >  1:1902029/1‑71 (MQ=255)
                                                    aCCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACt                   <  1:2724325/70‑1 (MQ=255)
                                                      tgAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCaa                                    >  1:2941997/2‑51 (MQ=37)
                                                        aaGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGa                     >  1:3941952/1‑64 (MQ=255)
                                                        aaGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAg              >  1:1870311/1‑71 (MQ=255)
                                                        aaGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAg              >  1:2557579/1‑71 (MQ=255)
                                                        aaGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCg                <  1:2456686/69‑1 (MQ=255)
                                                         aGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGaac             <  1:2019518/71‑3 (MQ=255)
                                                          gCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCg            <  1:2171348/71‑1 (MQ=255)
                                                            ggCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTc          >  1:3259805/1‑71 (MQ=255)
                                                             gCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGGGAAAACCGTAACCTATGACTTCGAGCGTCt         <  1:159289/71‑1 (MQ=255)
                                                             gCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGt                              <  1:3101455/50‑1 (MQ=255)
                                                             gCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCt         <  1:1317846/71‑1 (MQ=255)
                                                             gCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCt         <  1:467250/71‑1 (MQ=255)
                                                               tGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTga       >  1:694260/1‑71 (MQ=255)
                                                                    tAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa                                  >  1:654740/1‑39 (MQ=255)
                                                                    tAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatggat  >  1:2556005/1‑71 (MQ=255)
                                                                     |                                                                     
TCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGAT  >  minE/775433‑775571

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: