Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 775,580 A→G 8.9% K401K (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,5800AG 8.9% 250.0 / 13.3 90K401K (AAA→AAGicdisocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base A (43/39);  new base G (3/5);  total (46/44)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATT  >  minE/775516‑775647
                                                                |                                                                   
tGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTGATGGAAGGCGCTAAGctgctg                                                               <  1:3745558/71‑1 (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)
        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)
             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:852117/1‑71 (MQ=255)
              aaCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGtt                                                 >  1:892116/1‑71 (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)
                   gAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGt                                             <  1:443054/70‑1 (MQ=255)
                    aaaaCCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGtt                                                 <  1:1861229/65‑1 (MQ=255)
                     aaaCCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGa                                           >  1:589595/1‑70 (MQ=255)
                      aaCCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcg                                         >  1:1380692/1‑71 (MQ=255)
                      aaCCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcg                                         >  1:3461826/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:2762455/71‑1 (MQ=255)
                                    aCTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAAc                           <  1:1865233/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:2674273/37‑1 (MQ=25)
                                      ttCGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAAt                                                           <  1:595213/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:3347349/1‑70 (MQ=255)
                                        cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                        <  1:3236137/70‑1 (MQ=255)
                                        cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                        >  1:3953634/1‑70 (MQ=255)
                                        cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                        >  1:3995278/1‑70 (MQ=255)
                                        cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                        >  1:2405560/1‑70 (MQ=255)
                                        cGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATg                        >  1:2338659/1‑70 (MQ=255)
                                        cGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATg                        >  1:964617/1‑70 (MQ=255)
                                         gaacGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTa                      <  1:3776972/68‑1 (MQ=255)
                                         gAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGt                       <  1:2778753/70‑1 (MQ=255)
                                         gAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGt                       <  1:1592153/70‑1 (MQ=255)
                                          aacGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTa                      <  1:2672830/68‑1 (MQ=255)
                                            cGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcatc                                 >  1:1855742/1‑57 (MQ=255)
                                              tCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTg                                               <  1:974736/41‑1 (MQ=255)
                                               cTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGa                               >  1:1873766/1‑56 (MQ=255)
                                               cTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCg                >  1:1897193/1‑71 (MQ=255)
                                               cTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGatcat                                  >  1:2997723/1‑53 (MQ=255)
                                                tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGat                                     <  1:633284/49‑1 (MQ=255)
                                                tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGat                                     <  1:3160492/49‑1 (MQ=255)
                                                tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATg                   >  1:1065761/1‑67 (MQ=255)
                                                tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATg                   >  1:1770847/1‑67 (MQ=255)
                                                tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCg                >  1:810503/1‑70 (MQ=255)
                                                  atggatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTaa                     >  1:309759/1‑63 (MQ=255)
                                                  atggatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAg             >  1:341330/1‑71 (MQ=255)
                                                     gatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACAt                         <  1:3061178/56‑1 (MQ=255)
                                                     gatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGttt          >  1:136183/1‑71 (MQ=255)
                                                      atggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTaa                     >  1:3669567/1‑59 (MQ=255)
                                                       tggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGt        >  1:3609091/1‑71 (MQ=255)
                                                       tggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGt        <  1:1641257/71‑1 (MQ=255)
                                                         gcgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTa      <  1:1266873/71‑1 (MQ=255)
                                                          cgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcatc                                 >  1:2916913/1‑43 (MQ=255)
                                                          cgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGtt       <  1:1202046/69‑1 (MQ=255)
                                                           gcTAAACTGCTGAAATGTTCAGAGTTTGGTGCCGCGATCATCGaaa                             >  1:4145209/1‑46 (MQ=255)
                                                           gcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTaaa    <  1:648918/71‑1 (MQ=255)
                                                            cTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaa                              <  1:1873003/44‑1 (MQ=255)
                                                            cTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAAt                    >  1:3191364/1‑54 (MQ=255)
                                                            cTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGt            <  1:331815/62‑1 (MQ=255)
                                                            cTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGt            <  1:1975351/62‑1 (MQ=255)
                                                              aaaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTaa                     >  1:1282353/1‑51 (MQ=255)
                                                              aaaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAtt  >  1:3135333/1‑70 (MQ=255)
                                                              aaaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAtt  >  1:1124964/1‑70 (MQ=255)
                                                               aaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAtt  >  1:2634726/1‑69 (MQ=255)
                                                                |                                                                   
TGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATT  >  minE/775516‑775647

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: