Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 775,571 T→A 10.0% D398E (GAT→GAA icd → isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE775,5710TA10.0% 217.3 / 15.4 80D398E (GAT→GAAicdisocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base T (39/33);  new base A (3/5);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

TAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGT  >  minE/775501‑775641
                                                                      |                                                                      
tAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatggat                                                                        >  1:2556005/1‑71 (MQ=255)
  aTTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGatggatgg                                                                      >  1:2693609/1‑71 (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)
                 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)
                           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)
                              aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATg                                                    <  1:2277803/61‑1 (MQ=255)
                              aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt                                          <  1:757741/71‑1 (MQ=255)
                              aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt                                          <  1:2036013/71‑1 (MQ=255)
                              aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt                                          <  1:3611050/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: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:2405560/1‑70 (MQ=255)
                                                       cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                  >  1:3995278/1‑70 (MQ=255)
                                                       cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                  >  1:3953634/1‑70 (MQ=255)
                                                       cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                  >  1:2410303/1‑70 (MQ=255)
                                                       cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                  >  1:3347349/1‑70 (MQ=255)
                                                       cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                  <  1:3236137/70‑1 (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:3160492/49‑1 (MQ=255)
                                                               tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGat                               <  1:633284/49‑1 (MQ=255)
                                                               tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATg             >  1:1770847/1‑67 (MQ=255)
                                                               tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATg             >  1:1065761/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)
                                                                      |                                                                      
TAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGT  >  minE/775501‑775641

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: