Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 775,559 G→A 9.9% E394E (GAG→GAA icd → isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE775,5590GA 9.9% 242.6 / 12.7 81E394E (GAG→GAAicdisocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base G (38/35);  new base A (4/4);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGT  >  minE/775489‑775626
                                                                      |                                                                   
aaGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAg                                                                     >  1:1870311/1‑71 (MQ=255)
aaGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAg                                                                     >  1:2557579/1‑71 (MQ=255)
  gCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCg                                                                   <  1:2171348/71‑1 (MQ=255)
    ggCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTc                                                                 >  1:3259805/1‑71 (MQ=255)
     gCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGGGAAAACCGTAACCTATGACTTCGAGCGTCt                                                                <  1:159289/71‑1 (MQ=255)
     gCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCt                                                                <  1:1317846/71‑1 (MQ=255)
     gCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCt                                                                <  1:467250/71‑1 (MQ=255)
       tGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTga                                                              >  1:694260/1‑71 (MQ=255)
            tAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatggat                                                         >  1:2556005/1‑71 (MQ=255)
              aTTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGatggatgg                                                       >  1:2693609/1‑71 (MQ=255)
                 gTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatgga                                                          >  1:56177/1‑65 (MQ=255)
                  ttAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCt                                                   <  1:3293128/71‑1 (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: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)
                                               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)
                                                   ccGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAAct                                              <  1:2900791/43‑1 (MQ=255)
                                                    tgTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcgcg                  <  1:3375972/69‑1 (MQ=255)
                                                     gTAACCTATGAGTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGaaa                                       <  1:1763227/48‑1 (MQ=255)
                                                      tAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAActgctg                                          <  1:1760055/44‑1 (MQ=255)
                                                      tAACCTATGACTTCGAGAGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGat                >  1:2277514/1‑70 (MQ=255)
                                                        aCCTATGACTTCGAGCGTCTGATGGATGGGGCTATACTGCTGata                                       <  1:2172875/45‑3 (MQ=255)
                                                        aCCTATGACTTCGAGCGTCTGATGGATGGCGCTaaa                                                >  1:3860236/1‑36 (MQ=255)
                                                        aCCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatca              <  1:942420/70‑1 (MQ=255)
                                                         ccTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGa                 >  1:1216677/1‑66 (MQ=255)
                                                           tATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCag                                <  1:1651079/49‑1 (MQ=255)
                                                            aTGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt                           <  1:2249452/53‑1 (MQ=255)
                                                             tGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGa                 <  1:179072/62‑1 (MQ=255)
                                                             tGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaa         <  1:3517652/70‑1 (MQ=255)
                                                             tGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGa          <  1:2779013/69‑1 (MQ=255)
                                                               aCTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatca              <  1:2323281/63‑1 (MQ=255)
                                                               aCTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatc               >  1:2573086/1‑62 (MQ=255)
                                                               aCTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaaa       >  1:1024636/1‑70 (MQ=255)
                                                               aCTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAAc      <  1:1865233/71‑1 (MQ=255)
                                                               aCTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAAc      <  1:2762455/71‑1 (MQ=255)
                                                                cTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCaga                               <  1:1650604/45‑1 (MQ=255)
                                                                 ttCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTc                                  <  1:1373448/41‑1 (MQ=255)
                                                                 ttCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTg                          >  1:3320992/1‑49 (MQ=255)
                                                                 ttCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGt                        >  1:2879708/1‑51 (MQ=255)
                                                                 ttCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAAt                                      <  1:595213/37‑1 (MQ=25)
                                                                 ttCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAAt                                      <  1:2674273/37‑1 (MQ=25)
                                                                  tCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAAt                                      >  1:3915664/1‑36 (MQ=255)
                                                                   cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcg                    >  1:648959/1‑53 (MQ=255)
                                                                   cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGat                >  1:3083934/1‑57 (MQ=255)
                                                                   cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg   >  1:2410303/1‑70 (MQ=255)
                                                                   cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg   >  1:2405560/1‑70 (MQ=255)
                                                                   cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg   <  1:3236137/70‑1 (MQ=255)
                                                                   cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg   >  1:3347349/1‑70 (MQ=255)
                                                                   cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg   >  1:3995278/1‑70 (MQ=255)
                                                                   cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg   >  1:3953634/1‑70 (MQ=255)
                                                                   cGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATg   >  1:2338659/1‑70 (MQ=255)
                                                                   cGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATg   >  1:964617/1‑70 (MQ=255)
                                                                    gAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGt  <  1:2778753/70‑1 (MQ=255)
                                                                    gAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGt  <  1:1592153/70‑1 (MQ=255)
                                                                      |                                                                   
AAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGT  >  minE/775489‑775626

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: