Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 775,538 A→G 13.6% K387K (AAA→AAG icd → isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE775,5380AG13.6% 150.8 / 15.7 59K387K (AAA→AAGicdisocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base A (28/23);  new base G (2/6);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

CTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACG  >  minE/775467‑775608
                                                                       |                                                                      
ctgcGCCACATGGGTTGGACCGAAGC‑‑CTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa                                                                       >  1:3215267/1‑71 (MQ=255)
   cgcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGaaaa                                                                       <  1:2309365/70‑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: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)
                      aaGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGa                                                          >  1:3941952/1‑64 (MQ=255)
                      aaGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAg                                                   >  1:2557579/1‑71 (MQ=255)
                      aaGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAg                                                   >  1:1870311/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:689095/1‑42 (MQ=255)
                                                   gAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAg                                                   >  1:584306/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:2036013/71‑1 (MQ=255)
                                                                aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt         <  1:757741/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)
                                                                     aaaaCCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGtt          <  1:1861229/65‑1 (MQ=255)
                                                                      aaaCCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGa    >  1:589595/1‑70 (MQ=255)
                                                                       aaCCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcg  >  1:3461826/1‑71 (MQ=255)
                                                                       aaCCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcg  >  1:1380692/1‑71 (MQ=255)
                                                                       |                                                                      
CTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACG  >  minE/775467‑775608

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: