Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,009,234 G→A 15.1% D190N (GAT→AAT)  folM → dihydrofolate reductase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,009,2340GA15.1% 128.4 / 13.8 53D190N (GAT→AAT) folMdihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base G (21/24);  new base A (5/3);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAG  >  minE/1009167‑1009300
                                                                   |                                                                  
ccGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACataat                                                                  >  1:953325/1‑70 (MQ=255)
 cGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatga                                                                >  1:1246864/1‑71 (MQ=255)
  gCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACaagatga                                                                >  1:3446357/1‑70 (MQ=255)
   cAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCACCGCTGATCCTGTTTAATGAACatgaagatg                                                              >  1:2007666/1‑71 (MQ=255)
      gCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatga                                                                   >  1:1261322/1‑63 (MQ=255)
          gCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCGGAata                                                       >  1:2828176/1‑71 (MQ=255)
          gCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAat                                                        <  1:3920685/70‑1 (MQ=255)
          gCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAat                                                        <  1:3494720/70‑1 (MQ=255)
            aCCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATc                                                     >  1:1448886/1‑71 (MQ=255)
             ccGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATc                                                     >  1:3151163/1‑70 (MQ=255)
               ggAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGac                                                  <  1:665425/71‑1 (MQ=255)
                gAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGac                                                  >  1:2325116/1‑70 (MQ=255)
                 aagtgaaagtgaATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCg                                                    >  1:1766628/1‑67 (MQ=255)
                 aagtgaaagtgaATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGacaa                                                >  1:3338436/1‑71 (MQ=255)
                  agtgaaagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatg                                                                 <  1:38784/53‑1 (MQ=255)
                  agtgaaagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatg                                                                 <  1:1775461/53‑1 (MQ=255)
                  agtgaaagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATAATGATGCCGAATATCGacaa                                                >  1:379600/1‑70 (MQ=255)
                      aaagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGgc                                           <  1:1304892/71‑1 (MQ=255)
                      aaagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATAATGATGCCGAATATCGACAACAGgc                                           >  1:3606797/1‑71 (MQ=255)
                       aagtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGgc                                           <  1:1030429/70‑1 (MQ=255)
                          tgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCt                                        >  1:1713574/1‑70 (MQ=255)
                           gaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTg                                       >  1:376103/1‑70 (MQ=255)
                            aaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatg                                                                 <  1:2466619/43‑1 (MQ=255)
                            aaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgat                                                                  >  1:3309896/1‑42 (MQ=255)
                            aaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgat                                                                  >  1:2817511/1‑42 (MQ=255)
                            aaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgat                                                                  >  1:3639789/1‑42 (MQ=255)
                              ttCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatgatg                                                              >  1:3711260/1‑44 (MQ=255)
                               tCTATTGCGCCATCGCTGATCCTGTTTAATGAACatg                                                                    <  1:222865/37‑1 (MQ=255)
                                cTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatga                                                                >  1:1345182/1‑40 (MQ=255)
                                     gcgcCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAgg                                            >  1:2224566/1‑55 (MQ=255)
                                      cgcCATCGCTTATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtg                           <  1:1609946/71‑1 (MQ=255)
                                      cgcCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtg                           <  1:2657536/71‑1 (MQ=255)
                                      cgcCATCGCTGACCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACt                            <  1:3301220/70‑1 (MQ=255)
                                       gcCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCt                                        <  1:2310221/57‑1 (MQ=255)
                                       gcCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATa                                   <  1:1060280/62‑1 (MQ=255)
                                       gcCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATa                                   <  1:819115/62‑1 (MQ=255)
                                        ccATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGgc                                           >  1:1392980/1‑53 (MQ=255)
                                         cATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtgatg                        <  1:66728/71‑1 (MQ=255)
                                         cATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtgatg                        <  1:3508169/71‑1 (MQ=255)
                                         cATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtgatg                        <  1:1026550/71‑1 (MQ=255)
                                           tCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaa                      <  1:1525039/71‑1 (MQ=255)
                                             gCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaaaa                    >  1:2044311/1‑71 (MQ=255)
                                             gCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaaaa                    >  1:2724608/1‑71 (MQ=255)
                                             gCTGATCCTGTTTAATGAACATAATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaaa                     <  1:2887736/70‑1 (MQ=255)
                                              cTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAAt                                <  1:3524336/58‑1 (MQ=255)
                                               tGATCCTGTTTAATGAACATAATGATGCCGAATATCGaca                                                 <  1:3109632/40‑1 (MQ=255)
                                               tGATCCTGTTTAATGAACATAATGATGCCGAATATCGaca                                                 <  1:1121874/40‑1 (MQ=255)
                                               tGATCCTGTTTAATGAACATAATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAAcc                  >  1:1472451/1‑71 (MQ=255)
                                                 aTCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCa                              <  1:3513810/57‑1 (MQ=255)
                                                 aTCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACcgc                >  1:3001943/1‑71 (MQ=255)
                                                 aTCCTGTTTAATGAACATAATGATGCCGAATATCGacaac                                               >  1:2992844/1‑40 (MQ=255)
                                                    cTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtg                           <  1:53733/57‑1 (MQ=255)
                                                      gTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTg           >  1:1981287/1‑71 (MQ=255)
                                                       tttAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTg           <  1:853823/70‑1 (MQ=255)
                                                        ttAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtgat                         >  1:3833998/1‑55 (MQ=255)
                                                         tAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtga                          >  1:1134204/1‑53 (MQ=255)
                                                         tAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtg                           >  1:2462110/1‑52 (MQ=255)
                                                          aaTGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACcgc                <  1:2592287/62‑1 (MQ=255)
                                                          aaTGAACATAATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGc         >  1:1926109/1‑69 (MQ=255)
                                                             gAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCgaga     >  1:2676728/1‑70 (MQ=255)
                                                             gAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCgaga     >  1:824279/1‑70 (MQ=255)
                                                              aaCATGATGATGCCGAATATCGACAACAGGCGCTg                                       <  1:1688726/35‑1 (MQ=255)
                                                              aaCATGATGATGCCGAATATCGACAACAGGCGCTg                                       <  1:554959/35‑1 (MQ=255)
                                                              aaCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGagaaa   <  1:411295/71‑1 (MQ=255)
                                                               aCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGc         >  1:114930/1‑64 (MQ=255)
                                                               aCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCgaga     >  1:274006/1‑68 (MQ=255)
                                                               aCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGagaaag  <  1:863792/71‑1 (MQ=255)
                                                                  tgatgatgCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGcc             >  1:3460750/1‑57 (MQ=255)
                                                                  tgatgatgCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGcc             >  1:2818314/1‑57 (MQ=255)
                                                                   gatgatgCCGAATATCGACAACAGGCGCTGAATAAATCACtgatga                       >  1:1887244/1‑46 (MQ=255)
                                                                   |                                                                  
CCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAG  >  minE/1009167‑1009300

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: