Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,113,332 G→T 13.3% D204Y (GAT→TAT)  metF → 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,113,3320G→T13.3% 3.2 ‑4.5 ‑0.051610metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  new base (6/1):  ref base (14/30):  total (20/31)
Fisher's exact test strand distribution p-value = 1.10e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.45e-01

TGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGCCGGTATCT  >  REL606/4113297‑4113367
                                   |                                   
tGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTg                                     <  1:504594/36‑1 (MQ=255)
tGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTg                                     <  1:3463949/36‑1 (MQ=255)
 gCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGa                                    <  1:6609378/36‑1 (MQ=255)
  cGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAt                                   <  1:4018316/36‑1 (MQ=255)
  cGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGAt                                   <  1:3753270/36‑1 (MQ=255)
     tttCGTGACCGCTGTGTATCGGCGGGCATTTAtgtg                                >  1:2424979/1‑36 (MQ=255)
     tttCGTGACCGCTGTGTATCGGCGGGCATTGAtgtg                                <  1:7167290/36‑1 (MQ=255)
      ttCGTGACCGCTGTGTATCGGCGGGCATTGATGTgg                               <  1:4883260/36‑1 (MQ=255)
       tCGTGACCGCTGTGTATCGGCGGGCATTGTTGTggg                              >  1:3835609/1‑35 (MQ=255)
       tCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGa                              <  1:1102967/36‑1 (MQ=255)
        cGTGACCGCTGTGTATCGGCGGGCATTGATGTGGaa                             <  1:7381475/36‑1 (MQ=255)
         ttGACCGCTGTGTATCGGCGGGCATTGATGTGGaaa                            <  1:6912046/35‑1 (MQ=255)
         gTGACCGCTGTGTATCGGCGGGCATTGATGTGGaaa                            <  1:5528936/36‑1 (MQ=255)
          tGACCGCTGTGTATCGGCGGGCCTTGATTTTGAAat                           >  1:2038685/1‑36 (MQ=25)
          tGACCGCTGTGTATCGGCGGGCATTGATGTGGAAat                           >  1:2271571/1‑36 (MQ=255)
          tGACCGCTGTGTATCGGCGGGCATTGATGTGGAAat                           <  1:4996099/36‑1 (MQ=255)
           gACCGCTGTGTATCGGCGGGCATTTATGTGGGAatt                          >  1:5029419/1‑36 (MQ=37)
           gACCGCTGTGTATCGGCGGGCATTGATGTGGAAatt                          <  1:6060707/36‑1 (MQ=255)
           gACCGCTGTGTATCGGCGGGCATTGATGTGGAAatt                          <  1:7420731/36‑1 (MQ=255)
            aCCGCTGTGTATCGGCGGGCATTTATGTGGAAatta                         >  1:6782431/1‑36 (MQ=255)
            aCCGCTGTGTATCGGCGGGCATTGATGTGGAAatta                         <  1:712874/36‑1 (MQ=255)
            aCCGCTGTGTATCGGCGGGCATTGATGTGGAAatta                         <  1:4205611/36‑1 (MQ=255)
             ccGCTGTGTATCGGCGGGCATTGATGTGGAAattat                        <  1:1661919/36‑1 (MQ=255)
              cGCTGTGTATCGGCGGGCATTGATGTGGAAattatt                       <  1:290601/36‑1 (MQ=255)
               gCTGTGTATCGGCGGGCATTGATGTGGAAattattt                      >  1:163746/1‑35 (MQ=255)
               gCTGTGTATCGGCGGGCATTGATGTGGAAATTATTc                      <  1:1059936/36‑1 (MQ=255)
               acTGTGTATCGGCGGTCATTGATGTGGAAATTATTc                      <  1:830695/35‑1 (MQ=255)
                 tgtgTATCTGCGTGCATTTATGTGGAAATTATTCCg                    <  1:1809396/36‑1 (MQ=25)
                 tgtgTATCGGCGGGCATTTATGTGGGAattattttt                    >  1:4782100/1‑33 (MQ=25)
                 tgtgTATCGGCGGGCATGTGGGTGGGAATTATTTCg                    >  1:6427492/1‑36 (MQ=9)
                  gtgtATCGGCGGGCATTGATGTGGAAATTATTCCgg                   <  1:5779661/36‑1 (MQ=255)
                    gtATCGGCGGGGATTGATGTGGGAATTATTACgggc                 >  1:1077609/1‑35 (MQ=21)
                    gtATCGGCGGGCATTGATGTGTAAATTATTCCGGGa                 <  1:5132458/36‑1 (MQ=255)
                    gtATCGGCGGGCATTGATGTGGAAATTATTCCGGGa                 <  1:4934293/36‑1 (MQ=255)
                    gtATCGGCGGGCATTGATGTGGAAATTATTCCGGGa                 <  1:3973436/36‑1 (MQ=255)
                     tATCGGCGGGCATTTATTTGGGAATTATTCCGGGaa                >  1:6907281/1‑36 (MQ=25)
                     tATCGGCGGGCATTGATGTGGAAATTATTCCGGGaa                <  1:5605944/36‑1 (MQ=255)
                      aTCGGCGGGCATTGATGTGGAAATTATTCCGGGAAt               <  1:2934008/36‑1 (MQ=255)
                       tCGGCGGGCATTGATGTGGAAATTATTCCGGGAAtt              >  1:3624922/1‑36 (MQ=255)
                        cggcggGCATTGATGTGGAAATTATTCCGGGAAttt             <  1:1829496/36‑1 (MQ=255)
                        cggcggGCATTGATGTGGAAATTATTCCGGGAAttt             >  1:2990009/1‑36 (MQ=255)
                         ggcggGCATTGATGTGGAAATTATTCCGGGAAtttt            >  1:613374/1‑36 (MQ=255)
                          gcggGCATTGATGTGGAAATTATTCCGGGAATTTTg           <  1:7454030/36‑1 (MQ=255)
                            gggCATTGATGTGGAAATTATTCCGGGAATTTTGcc         <  1:6256585/36‑1 (MQ=255)
                            gggCATTGATGTGGAAATTATTCCGGGAATTTTGcc         <  1:251745/36‑1 (MQ=255)
                            gggCATTGATGTGGAAATTATTCCGGGAATTTTGcc         >  1:6416048/1‑36 (MQ=255)
                             ggCATTGATGTGGAAATTATTCCGGGAATTTTGCCg        >  1:383351/1‑36 (MQ=255)
                              gCATTGATGTGGAAATTATTCCGGGAATTTTGCCgg       >  1:6316026/1‑36 (MQ=255)
                               cATTGATGTGGAAATTATTCCGGGAATTTTGCCGGt      <  1:3329801/36‑1 (MQ=255)
                               cATTGATGTGGAAATTATTCCGGGAATTTTGCCGGt      >  1:4460062/1‑36 (MQ=255)
                               cATTGATGTGGAAATTATTCCGGGAATTTTGCCGGt      >  1:3894244/1‑36 (MQ=255)
                                 ttGATGTGGAAATTATTCCGGGAATTTTGCCGGTAt    >  1:7369161/1‑36 (MQ=255)
                                  tGATGTGGAAATTATTCCGGGAATTTTGCCGGTATc   <  1:5829627/36‑1 (MQ=255)
                                   gATGTGGAAATTATTCCGGGTATTTTGCCGGTATCt  >  1:7878214/1‑36 (MQ=255)
                                   gATGTGGAAATTATTCCGGGAATTTTGCCGGTATCt  <  1:1755103/36‑1 (MQ=255)
                                   gATGTGGAAATTATTCCGGGAATTTTGCCGGTATCt  <  1:5721220/36‑1 (MQ=255)
                                   |                                   
TGCGTTTTCGTGACCGCTGTGTATCGGCGGGCATTGATGTGGAAATTATTCCGGGAATTTTGCCGGTATCT  >  REL606/4113297‑4113367

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG