Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 775,535 G→C 15.9% A386A (GCG→GCC icd → isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE775,5350GC15.9% 148.1 / 19.5 63A386A (GCG→GCCicdisocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base G (29/24);  new base C (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 = 8.61e-01

CTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGT  >  minE/775467‑775604
                                                                    |                                                                     
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)
             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)
                  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:2496505/1‑56 (MQ=255)
                                                              cAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAActgc                      >  1:4079056/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:892116/1‑71 (MQ=255)
                                                               aaCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGtt      >  1:852117/1‑71 (MQ=255)
                                                                aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTg                                         <  1:1798323/35‑1 (MQ=255)
                                                                aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATg               <  1:2277803/61‑1 (MQ=255)
                                                                aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt     <  1:757741/71‑1 (MQ=255)
                                                                aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt     <  1:3611050/71‑1 (MQ=255)
                                                                aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt     <  1:2036013/71‑1 (MQ=255)
                                                                 cgcgAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTaaa                          >  1:554113/1‑49 (MQ=255)
                                                                  gcgAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTCCTGAAATGTTCAGAGTTTgg   >  1:2691290/1‑71 (MQ=255)
                                                                   cgAAAACCGTAACCTATGACTTCGAGCGTCTgatgg                                     >  1:1664229/1‑36 (MQ=255)
                                                                    gAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGt  <  1:443054/70‑1 (MQ=255)
                                                                    |                                                                     
CTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGT  >  minE/775467‑775604

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: