Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,205,690 A→T 19.7% K978* (AAA→TAA)  metH → B12‑dependent methionine synthase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,205,6900A→T19.7% 4.5 ‑4.0 0.0442932metHB12‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  new base (10/0):  ref base (20/14):  total (30/14)
Fisher's exact test strand distribution p-value = 1.83e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCC  >  REL606/4205657‑4205725
                                 |                                   
gAAGTGGTGGGCGTTGAGGCGCAGCGGCTGGTTtaa                                   >  1:4450455/1‑36 (MQ=37)
 aaGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTTAAg                                  >  1:3839472/1‑36 (MQ=255)
 aaGTGGTGGGCGCTGAGGCGCAGCGGCTGTTTTAAg                                  >  1:2470299/1‑36 (MQ=37)
  aGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTTAAGa                                 >  1:6637204/1‑36 (MQ=255)
  aGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGa                                 <  1:2822434/36‑1 (MQ=255)
   gtggtgGGCGTTGAGGCGCAGCGGCTGTTTAAAGAc                                <  1:6491202/36‑1 (MQ=255)
    tggtggGCGTTGAGGCGCAGCGGCTGGTTTAAGACg                               >  1:5451045/1‑36 (MQ=37)
      gtggGCGTTGAGGCGCAGCGGCTGTTTTAAGACGcc                             >  1:3039258/1‑36 (MQ=255)
       tggGCGTTGAGGCGCAGCGGCTGGTTAAAGACGCCa                            >  1:5707693/1‑36 (MQ=255)
        gggCGTTGAGGCGCAGCGGCTGTTTTAAGACGccca                           >  1:6843613/1‑34 (MQ=38)
            gTTGAGGCGCAGCGGCTGTTTTAAGACGCCAACGAc                       >  1:772314/1‑36 (MQ=255)
            gTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGAc                       >  1:7814829/1‑36 (MQ=255)
            gTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGAc                       <  1:4567650/36‑1 (MQ=255)
             ttGAGGCGCAGCGGCTGTTTAAAAACGCCACAGACa                      >  1:2506952/1‑36 (MQ=25)
              tGAGGCGCAGCGGCTGTTTTAAGACGCCAACGACaa                     >  1:5773873/1‑35 (MQ=255)
              tGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACAt                     <  1:7174368/36‑1 (MQ=255)
               gAGGCGCAGCGGCTGTTTTAAGACGCCAACGACATg                    >  1:4065445/1‑36 (MQ=255)
                 ggCGCAGCGGCTGTTTAAAGACGCCAACGACATGCt                  <  1:3581156/36‑1 (MQ=255)
                  gcgcAGCGGCTGTTTAAAGACGCCAACGAAATGCTg                 >  1:7205354/1‑36 (MQ=255)
                     cAGCGGCTGTTTAAAGACGCCAACGACATGCTGGAt              <  1:1290354/36‑1 (MQ=255)
                      aGCGGCTGGTTTAAGACGCCAACGACAAGCTGGGTa             >  1:7101608/1‑36 (MQ=19)
                        cGGCTGTTTAAAGACGCCAACGACATGCTGGATaaa           >  1:4209781/1‑36 (MQ=255)
                         ggCTGTTTAAAGACGCCAACGACATGCTGGATAAAt          >  1:3777719/1‑36 (MQ=255)
                         ggCTGTATAAACACGCCAACGACATGCTGGATAAAt          <  1:3342469/36‑1 (MQ=37)
                          gCTGTTTAAAGACGCCATCGACATGCTGGATAAAtt         <  1:2781516/36‑1 (MQ=255)
                          gCTGTTTAAAGACGCCAACGACATGCTGGATAAAtt         <  1:7016546/36‑1 (MQ=255)
                          gCTGTTTAAAGACGCCAACGACATGCTGGATAAAtt         >  1:4611501/1‑36 (MQ=255)
                          gCTGTTTAAAGACGCCAACGACATGCTGGATAAAtt         >  1:2584636/1‑36 (MQ=255)
                            taattAAAGACACCAGCGACATGCTGGATAAATTaa       <  1:6110316/33‑1 (MQ=21)
                            tGTTTAAAGACGCCAACGACATGCTGGATAAATTaa       >  1:4041013/1‑36 (MQ=255)
                            tGTTTAAAGACGCCAACGACATGCTGGATAAATTaa       >  1:1957987/1‑36 (MQ=255)
                            tGTTTAAAGACGCCAACGACATGCTGGATAAATTaa       <  1:5705613/36‑1 (MQ=255)
                            tGTTTAAAGACGACAACGACAGGCTGGATAAATTaa       >  1:1240813/1‑36 (MQ=38)
                             gTTTAAAGACGCCAACGACATGCTGGATAAATTAAg      >  1:2610120/1‑36 (MQ=255)
                             gTTTAAAGACGCCAACGACATGCTGGATAAATTAAg      >  1:7662062/1‑36 (MQ=255)
                              tttAAAGACGCCAACGACATGCTGGATAAATTAAgc     <  1:6161229/36‑1 (MQ=255)
                              tttAAAGACGCCAACGACATGCTGGATAAATTAAgc     >  1:559952/1‑36 (MQ=255)
                              tttAAAGACGCCAACGACATGCTGGATAAATTAAgc     <  1:4847964/36‑1 (MQ=255)
                              tttAAAGACGCCAACGACATGCTGGATAAATTAAgc     <  1:3278088/36‑1 (MQ=255)
                              tttAAAGACGCCAACGAAATGCTGGATAAATTAAgc     >  1:4303989/1‑36 (MQ=255)
                               ttAAAGACGCCAACGACATGCTGGATAAATTAAgcg    >  1:6428434/1‑36 (MQ=255)
                               ttAAAGACGCCAAAGACATGCTGGATAAATTAAgcg    >  1:6912601/1‑36 (MQ=255)
                                tAAAGACGCCAACGACATGCTGGATAAATTAAgcgc   >  1:2844932/1‑36 (MQ=255)
                                tAAAGACGCCAACGACATGCTGGATAAATTAAgcgc   >  1:7456140/1‑36 (MQ=255)
                                tAAAGACGCCAACGACATGCTGGATAAATTAAgcgc   >  1:1501407/1‑36 (MQ=255)
                                 aaaGACGCCAACGACATGCTGGATAAGTTAAGCGcc  >  1:999438/1‑36 (MQ=255)
                                 aaaGACGCCAACGACATGCTGGATAAATTAAGCGcc  <  1:7819854/36‑1 (MQ=255)
                                 |                                   
GAAGTGGTGGGCGTTGAGGCGCAGCGGCTGTTTAAAGACGCCAACGACATGCTGGATAAATTAAGCGCC  >  REL606/4205657‑4205725

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 16 ≤  ATCG  < 27 ≤  ATCG  < 37 ≤  ATCG