Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 775,532 C→T 15.9% N385N (AAC→AAT icd → isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE775,5320CT15.9% 158.4 / 20.2 63N385N (AAC→AATicdisocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base C (29/24);  new base T (3/7);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

GATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTT  >  minE/775463‑775601
                                                                     |                                                                     
gatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAcg                                                                      >  1:1728706/1‑71 (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)
              tGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAct                                                        <  1:3221389/71‑2 (MQ=255)
                ggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                                                                       >  1:3253117/1‑54 (MQ=255)
                ggTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc                                                                       >  1:1775570/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:176514/71‑1 (MQ=255)
                  ttGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATg                                                    <  1:3761392/71‑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)
                                        aTTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGatggatgg                              >  1:2693609/1‑71 (MQ=255)
                                         ttGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCg                                                             <  1:2296297/39‑1 (MQ=255)
                                         ttGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAcc                                                        <  1:3572206/44‑1 (MQ=255)
                                           gTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt                                                     >  1:3875256/1‑45 (MQ=255)
                                           gTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatgga                                 >  1:56177/1‑65 (MQ=255)
                                            ttAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCt                          <  1:3293128/71‑1 (MQ=255)
                                             tAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa                                                          >  1:2414801/1‑38 (MQ=255)
                                              aaaGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTa                         >  1:2585758/1‑70 (MQ=255)
                                                     tGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGctgctg                 <  1:3745558/71‑1 (MQ=255)
                                                     tGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatggat                                <  1:1659914/56‑1 (MQ=255)
                                                       gAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAg                                            >  1:584306/1‑42 (MQ=255)
                                                       gAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAg                                            >  1:689095/1‑42 (MQ=255)
                                                       gAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGaa               <  1:2241532/71‑1 (MQ=255)
                                                       gAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGa                >  1:3964817/1‑70 (MQ=255)
                                                        aaGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGaaa              >  1:306101/1‑71 (MQ=255)
                                                          ggCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATg            >  1:2428758/1‑71 (MQ=255)
                                                            cgcAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTga                                     >  1:4079303/1‑44 (MQ=255)
                                                             gcAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGtt          >  1:3816672/1‑70 (MQ=255)
                                                              cAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCa        <  1:603036/71‑1 (MQ=255)
                                                                aTCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCa        >  1:438164/1‑69 (MQ=255)
                                                                 tCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTaaa                       >  1:51584/1‑53 (MQ=255)
                                                                 tCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTa                         <  1:2789576/51‑1 (MQ=255)
                                                                  cAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAActgc                   >  1:4079056/1‑56 (MQ=255)
                                                                  cAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAActgc                   >  1:2496505/1‑56 (MQ=255)
                                                                  cAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGt    >  1:745115/1‑71 (MQ=255)
                                                                   aaCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGaaa              >  1:1537653/1‑60 (MQ=255)
                                                                   aaCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCag       <  1:2553595/67‑1 (MQ=255)
                                                                   aaCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGtt   >  1:852117/1‑71 (MQ=255)
                                                                   aaCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGtt   >  1:892116/1‑71 (MQ=255)
                                                                    aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTg                                      <  1:1798323/35‑1 (MQ=255)
                                                                    aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATg            <  1:2277803/61‑1 (MQ=255)
                                                                    aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt  <  1:3611050/71‑1 (MQ=255)
                                                                    aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt  <  1:757741/71‑1 (MQ=255)
                                                                    aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt  <  1:2036013/71‑1 (MQ=255)
                                                                     cgcgAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTaaa                       >  1:554113/1‑49 (MQ=255)
                                                                     |                                                                     
GATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTT  >  minE/775463‑775601

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: