Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,204,731 A→G 11.4% E658G (GAA→GGA)  metH → B12‑dependent methionine synthase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,204,7310A→G11.4% 2.6 ‑3.7 0.0491973metHB12‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (20/22):  total (27/22)
Fisher's exact test strand distribution p-value = 2.48e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CGCTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAG  >  REL606/4204696‑4204766
                                   |                                   
cGCTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGa                                     <  1:3686070/36‑1 (MQ=255)
cGCTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGa                                     >  1:5329785/1‑36 (MQ=255)
 gCTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGaa                                    <  1:2664099/36‑1 (MQ=255)
  cTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGGAg                                   >  1:5836209/1‑36 (MQ=255)
  cTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAg                                   >  1:4111063/1‑36 (MQ=255)
  cTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAg                                   <  1:6970442/36‑1 (MQ=255)
  cTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAg                                   <  1:1820106/36‑1 (MQ=255)
   tAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGt                                  <  1:8418451/36‑1 (MQ=255)
    aaCGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTg                                 >  1:8161583/1‑36 (MQ=255)
     aCGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGa                                >  1:2189794/1‑36 (MQ=255)
     aCGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGa                                >  1:3977591/1‑36 (MQ=255)
     aCGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGa                                <  1:5575736/36‑1 (MQ=255)
     aCGCCCAGCAGGCGGAGTGGCGCTCGCGGGTAGTGa                                >  1:3896326/1‑36 (MQ=37)
      cGCCCAGCAGGCGGAGTGGCGCTCGGGGGAAGTGaa                               >  1:2059399/1‑36 (MQ=255)
        cccAGCAGGCGGAGTGGCGCTCGTGGGGAGGGGATa                             >  1:3613135/1‑36 (MQ=25)
         ccAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATaa                            >  1:7725842/1‑36 (MQ=255)
         ccAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATaa                            <  1:7953578/36‑1 (MQ=255)
          cagcagGCGGAGTGGCGCTCGTGGGAAGTGAATaaa                           <  1:1507739/36‑1 (MQ=255)
          cagcagGCGGAGTGGCGCTCGTGGGAAGTGAATaaa                           <  1:8299003/36‑1 (MQ=255)
          cagcagGCGGAGTGGCGCTCGGGGGAAGGGAATaaa                           >  1:2849332/1‑36 (MQ=38)
           agcagGCGGAGTGGCGCTCGTGGGGAGTGAATAAAc                          >  1:5605443/1‑36 (MQ=255)
           agcagGCGGAGTGGCGCTCGTGGGAAGTGAATAAAc                          >  1:336569/1‑36 (MQ=255)
           agcagGCGGAGTGGCGCTCGTGGGAAGTGAATAAAc                          <  1:2109251/36‑1 (MQ=255)
           agcagGCGGAGTGGCGCTCGTGGGAAGTGAATAAAc                          >  1:2551918/1‑36 (MQ=255)
            gcagGCGGAGTGGCGCTCGTGGGAAGTGAATAAACg                         >  1:4227274/1‑36 (MQ=255)
            gcagGCGGAGTGGCGCTCGTGGGAAGTGAATAAACg                         >  1:875324/1‑36 (MQ=255)
            gcagGCGGAGTGGCGCTCGGGGGAAGTGAATAAACg                         >  1:2726109/1‑36 (MQ=255)
             cagGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGt                        <  1:7646436/36‑1 (MQ=255)
              agGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTc                       >  1:2378656/1‑36 (MQ=255)
                gCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTg                     >  1:4960757/1‑36 (MQ=255)
                 cGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTgg                    >  1:7634123/1‑36 (MQ=255)
                 cGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTgg                    <  1:3016492/36‑1 (MQ=255)
                  ggAGTGGCGCTCGTGGGGAGTGAATAAACGTCTGGa                   >  1:434354/1‑36 (MQ=255)
                  ggAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGa                   <  1:139840/36‑1 (MQ=255)
                  ggAGTGGCGCTCGGGGGGAGTGGATAAACGGCTggg                   >  1:6936008/1‑35 (MQ=16)
                    aGTGGCGCTCGTGGGGAGTGAATAAAAGTCTGGGAt                 >  1:4665053/1‑36 (MQ=37)
                     gTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATa                >  1:2967036/1‑36 (MQ=255)
                     gTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATa                <  1:6690199/36‑1 (MQ=255)
                      tGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATAc               <  1:5990015/36‑1 (MQ=255)
                        gcgcTCGTGGGAAGTGAATAAACGTCTGGAATACTc             <  1:5998488/36‑1 (MQ=255)
                        gcgcTCGTGGGAAGTGAATAAACGTCTGGAATACTc             <  1:3652081/36‑1 (MQ=255)
                          gcTCGTGGGAAGTGAATAAACGTCTGGAATACTCGc           >  1:738577/1‑36 (MQ=255)
                           cTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCt          <  1:4092165/36‑1 (MQ=255)
                             cGTGGGAAGTGAATAAGCGTCTGGAATACTCGCTgg        <  1:6108231/36‑1 (MQ=255)
                              gTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGt       >  1:6078887/1‑36 (MQ=255)
                              gTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGt       >  1:2972960/1‑36 (MQ=255)
                               tGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTc      <  1:5632577/36‑1 (MQ=255)
                                tggAAGTGAATAAACGTCTGGAATACTCGCTGGTCa     <  1:6542361/35‑1 (MQ=255)
                                gggAAGTGAATAAACGTCTGGAATACTCGCTGGTCa     >  1:2625739/1‑36 (MQ=255)
                                 ggAAGTGAATAAACGTCTGGAATACTCGCTGGTCaa    <  1:643397/36‑1 (MQ=255)
                                  gAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa   <  1:796757/36‑1 (MQ=255)
                                  gAAGTGAATAAACGTCTGGAATACTCGCTGGTCaaa   <  1:2339996/36‑1 (MQ=255)
                                   aaGTGAATAAACGTCTGGAATACTCGCTGGTCAAAg  >  1:7251172/1‑36 (MQ=255)
                                   aaGTGAATAAACGTCTGGAATACTCGCTGGTCAAAg  <  1:5194051/36‑1 (MQ=255)
                                   aaGTGAATAAACGTCTGGAATACTCGCTGGTCAAAg  <  1:1401665/36‑1 (MQ=255)
                                   |                                   
CGCTAACGCCCAGCAGGCGGAGTGGCGCTCGTGGGAAGTGAATAAACGTCTGGAATACTCGCTGGTCAAAG  >  REL606/4204696‑4204766

Base quality scores:  ATCG  < 3 ≤  ATCG  < 14 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG