Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,683,257 G→A 29.5% D190N (GAT→AAT)  folM → dihydrofolate reductase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,683,2570GA29.5% 63.4 / 32.3 44D190N (GAT→AAT) folMdihydrofolate reductase isozyme
Reads supporting (aligned to +/- strand):  ref base G (15/16);  new base A (5/8);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

CCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAG  >  W3110S.gb/1683189‑1683323
                                                                    |                                                                  
cccGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgat                                                                  <  1:3271642/71‑1 (MQ=255)
   gCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatgat                                                               >  1:1382298/1‑71 (MQ=255)
   gCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatgat                                                               >  1:2335466/1‑71 (MQ=255)
   gCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatgat                                                               <  1:4942226/71‑1 (MQ=255)
   gCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACataatgat                                                               <  1:1062652/71‑1 (MQ=255)
   gCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACataatga                                                                >  1:1118709/1‑70 (MQ=255)
       gCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCg                                                           <  1:2725494/71‑1 (MQ=255)
       gCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCg                                                           <  1:1504597/71‑1 (MQ=255)
       gCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCg                                                           <  1:535972/71‑1 (MQ=255)
            cACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATAATGATGCCGAatat                                                      <  1:1115231/71‑1 (MQ=255)
                ggAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgatga                                                                >  1:4885046/1‑57 (MQ=255)
                  aagtgaaagtgaATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGacaa                                                <  1:3747253/71‑1 (MQ=255)
                         agtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGc                                                             >  1:3438422/1‑51 (MQ=255)
                         agtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGc                                                             >  1:1669992/1‑51 (MQ=255)
                         agtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGc                                                             >  1:3800053/1‑51 (MQ=255)
                         agtgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGc                                                             >  1:2731055/1‑51 (MQ=255)
                           tgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACataat                                                                  <  1:2043664/44‑1 (MQ=255)
                           tgaaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACataat                                                                  <  1:697132/44‑1 (MQ=255)
                             aaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAAGCGCTGa                                      <  1:4553030/70‑1 (MQ=255)
                             aaTTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATAGACAACAGGCGCTGa                                      <  1:1853646/70‑1 (MQ=255)
                               ttCTATTGCGCCATCGCTGATCCTGTTTAATGAACatgat                                                                  <  1:897371/40‑1 (MQ=255)
                                 cTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGacaa                                                >  1:3242532/1‑56 (MQ=255)
                                 cTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGa                                                   >  1:30091/1‑53 (MQ=255)
                                       cgcCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACt                            <  1:3547233/70‑1 (MQ=255)
                                           aTCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtgatga                       <  1:3070648/71‑1 (MQ=255)
                                            tCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATc                                                     <  1:4503652/40‑1 (MQ=255)
                                            tCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATc                                                     <  1:3393965/40‑1 (MQ=255)
                                            tCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaa                      <  1:3008958/71‑1 (MQ=255)
                                            tCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaa                      <  1:412165/71‑1 (MQ=255)
                                            tCGCTGATCCTGTTTAATGAACATAATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaa                      <  1:3224990/71‑1 (MQ=255)
                                            tCGCTCATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaa                      <  1:517179/71‑1 (MQ=255)
                                              gCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGaca                                                 >  1:2882914/1‑42 (MQ=255)
                                              gCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGgcg                                          >  1:96933/1‑49 (MQ=255)
                                              gCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGgcg                                          >  1:1418060/1‑49 (MQ=255)
                                              gCTGATCCTGTTTAATGAACATAATGATGCCGAATATCGacaac                                               >  1:2683589/1‑44 (MQ=255)
                                              gCTGATCCTGTTTAATGAACATAATGATGCCGAATATCGacaac                                               >  1:908409/1‑44 (MQ=255)
                                              gCTGATCCTGTTTAATGAACATAATGATGCCGAATATCGacaac                                               >  1:4763396/1‑44 (MQ=255)
                                              gCTGATCCTGTTTAATGAACATAATGATGCCGAATATCGacaac                                               >  1:144703/1‑44 (MQ=255)
                                              gCTGATCCTGTTTAATGAACATAATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaaa                     <  1:962387/70‑1 (MQ=255)
                                              gCTGATCCTGTTTAATGAACATAATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaaa                     <  1:787681/70‑1 (MQ=255)
                                                         ttAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGc         >  1:1382647/1‑71 (MQ=255)
                                                           aaTGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGc         >  1:3281016/1‑69 (MQ=255)
                                                            aTGAACATAATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACcgc                >  1:5081862/1‑61 (MQ=255)
                                                             tGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaa                      <  1:1638577/54‑1 (MQ=255)
                                                             tGAACATAATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtg                           <  1:1790583/49‑1 (MQ=255)
                                                             tGAACATAATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACtg                           <  1:4109219/49‑1 (MQ=255)
                                                               aaCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaaaa                    >  1:2382844/1‑54 (MQ=255)
                                                               aaCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaaaa                    >  1:733256/1‑54 (MQ=255)
                                                               aaCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaaaa                    >  1:4806384/1‑54 (MQ=255)
                                                               aaCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCgcgc              >  1:2469144/1‑60 (MQ=255)
                                                               aaCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGagaa    <  1:4756960/70‑1 (MQ=255)
                                                                aCATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGagaaag  <  1:3303457/71‑1 (MQ=255)
                                                                aCATAATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCTTGGCGagaaag  <  1:4488386/71‑1 (MQ=255)
                                                                aCATAATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGagaaag  <  1:5113700/71‑1 (MQ=255)
                                                                   tgatgatgCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGaaa                     <  1:1949486/49‑1 (MQ=255)
                                                                    |                                                                  
CCCGCAAGCTGGCACCGGAAGTGAAAGTGAATTCTATTGCGCCATCGCTGATCCTGTTTAATGAACATGATGATGCCGAATATCGACAACAGGCGCTGAATAAATCACTGATGAAAACCGCGCCTGGCGAGAAAG  >  W3110S.gb/1683189‑1683323

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: