Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 775,547 C→T 10.6% T390T (ACC→ACT icd → isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE775,5470CT10.6% 193.6 / 11.2 66T390T (ACC→ACTicdisocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base C (32/27);  new base T (2/5);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

GTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCA  >  minE/775480‑775614
                                                                   |                                                                   
gTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTAt                                                                  <  1:4015745/71‑1 (MQ=255)
gTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTAt                                                                  <  1:1501705/71‑1 (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)
         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)
              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)
                        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)
                             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: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)
                                                   aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTg                                                   <  1:1798323/35‑1 (MQ=255)
                                                   aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATg                         <  1:2277803/61‑1 (MQ=255)
                                                   aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt               <  1:3611050/71‑1 (MQ=255)
                                                   aCGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGttt               <  1:757741/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: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)
                                                                   |                                                                   
GTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCA  >  minE/775480‑775614

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: