Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,977,9930A→G26.3% 12.258K172E (AAA→GAA) metE5‑methyltetrahydropteroyltriglutamate‑ homocysteine methyltransferase
Reads supporting (aligned to +/- strand):  new base (14/1):  ref base (11/31):  total (25/33)

TGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAGTTTGATCGCCTGAGCCTG  >  REL606/3977958‑3978028
                                   |                                   
tgctgtTGGGGCCGGTTACCTGGCTGTGGCTGGGGa                                     <  1:6507300/36‑1 (MQ=255)
tgctgcTGGGGCCTGTTACCTGGCTGTGGCTGGGGa                                     <  1:4822978/36‑1 (MQ=255)
tgctgcTGGGGCCGGTTACCTGGCTGTGGCTGGGGa                                     <  1:939330/36‑1 (MQ=255)
tgctgcTGGGGCCGGTTACCTGGCTGTGGCTGGGGa                                     <  1:7001112/36‑1 (MQ=255)
 gctgctGGGGCCGGTTACCTGGCTGTGGCTGGGGag                                    >  1:2691833/1‑35 (MQ=255)
 gctgctGGGGCCGGTTACCTGGCTGTGGCTGGGGaa                                    <  1:2158660/36‑1 (MQ=255)
  ctgctgGGGCCGGTTACCTGGCTGTGGCTGGGGgaa                                   >  1:282251/1‑36 (MQ=255)
  ctgctgGGGCCGGTTACCTGGCTGTGGCTGGGGaaa                                   <  1:5095527/36‑1 (MQ=255)
  cgtctgGGGCCGGTTACCTGGCTGTGGCTGGGGaaa                                   <  1:3479288/33‑1 (MQ=255)
   tgctgGGGCCGGTTACCTGGCTGTGGCTGGGGGAAg                                  >  1:6185345/1‑36 (MQ=255)
   tgctgGGGCCGGTTACCTGGCTGTGGCTGGGGAGAg                                  >  1:2452591/1‑36 (MQ=255)
   tgctgGGGCCGGTTACCTGCCTGTGGCTGGGGGAAg                                  >  1:1149674/1‑36 (MQ=25)
    gctgGGGCCGGTTACCTGGCTGTGGCTGGGGGAAgg                                 >  1:2922318/1‑35 (MQ=255)
    gctgGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGt                                 >  1:695094/1‑36 (MQ=255)
     ctgGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTg                                <  1:97442/36‑1 (MQ=255)
     ctgGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTg                                <  1:7448918/36‑1 (MQ=255)
      tgGGGCGGGTTACCTGGCTGTGGCTGGGGAAAGTGa                               <  1:1664287/36‑1 (MQ=255)
       ggggCCGGTTACCTGGCTGTGGCTGGGGGTAGTgtt                              >  1:1384197/1‑34 (MQ=18)
       aaggCCGGTTACCTGGCTGTGGCTGGGGAAAGTGaa                              <  1:1046200/34‑1 (MQ=255)
        gggCTGGTTACCTGGCTGTGGCTGGGGGAAGTGaaa                             >  1:1259066/1‑36 (MQ=21)
        gggCCGGTTACCTGGCTGTGGCTGGGGAAAGTGaaa                             <  1:293039/36‑1 (MQ=255)
          gCCGGTTACCTGGCTGTGGCTGGGGGAAGTGGAAgg                           >  1:1969595/1‑36 (MQ=21)
          gCCGGTTACCTGGCTGTGGCTGGGGGAAGTGGAAgg                           >  1:4105255/1‑36 (MQ=21)
          gCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAgg                           <  1:588415/36‑1 (MQ=255)
           ccGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGt                          <  1:6657194/36‑1 (MQ=255)
            cGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTg                         >  1:3441203/1‑36 (MQ=255)
            cGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAAGTg                         >  1:5633279/1‑36 (MQ=21)
            cGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTg                         <  1:3799252/36‑1 (MQ=255)
             ggTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGa                        >  1:2578449/1‑36 (MQ=255)
             ggTTACCGGGCTGTGGCTGGGGAAACGGAAAGGTGa                        <  1:6529424/36‑1 (MQ=18)
              gTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGGGaa                       >  1:7225291/1‑36 (MQ=21)
              gTTACCTGGCTGTGGCTGGGGAACGTGAAAGGTGaa                       >  1:2632700/1‑36 (MQ=255)
              gTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGaa                       <  1:2446332/36‑1 (MQ=255)
               ttACCTTGCTGTGGCTGGGGAAAGTGAAAGGTGAAc                      <  1:667687/36‑1 (MQ=255)
               ttACCTGGCTGTGGCTGGGGAAAGCGAAAGGTGGAc                      >  1:111252/1‑36 (MQ=21)
                tgccTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACa                     <  1:1865793/34‑1 (MQ=255)
                 aCCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAg                    <  1:7487615/36‑1 (MQ=255)
                 aCCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAg                    <  1:3587211/36‑1 (MQ=255)
                 aCCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAg                    >  1:5609751/1‑36 (MQ=255)
                  ccTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAGt                   <  1:1002970/36‑1 (MQ=255)
                   cTGGCTGTGGCTGGGGTAAGGGAAAGGTGAACAGtt                  <  1:1072096/36‑1 (MQ=21)
                   cTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAGtt                  <  1:1873224/36‑1 (MQ=255)
                    tgggtgtggctgGGGAAAGTGAAAGGTGAACAGttt                 <  1:3342556/36‑1 (MQ=255)
                    tggctgtggctgGGGGAAGTGAAAGGGGAACAGttt                 >  1:600998/1‑36 (MQ=21)
                    tggctgtggctgGGGAAAGTGAAAGGTGAACAGttt                 <  1:6070449/36‑1 (MQ=255)
                    tggctgtggctgGGGAAAGTGAAAGGTGAACAGttt                 <  1:5745528/36‑1 (MQ=255)
                     ggctgtggctgGGGGAAATGAAAAGGGAACAGtttt                >  1:144817/1‑35 (MQ=17)
                     ggctgtggctgGGGAAAGTGAAAGGTGAACAGTTTg                <  1:5738038/36‑1 (MQ=255)
                      gctgtggctgGGGAAAGTGAAAGGTGAACAGTTTGa               <  1:2739524/36‑1 (MQ=255)
                      gctgtggctgGGGAAAGTGAAAGGTGAACAGTTTGa               <  1:5116690/36‑1 (MQ=255)
                       gtgtggctgGGGAAAGTGAAAGGTGAACAGTTTGAt              <  1:779722/35‑1 (MQ=255)
                       ctgtggctgGGGAAAGTGAAAGGTGAACAGTTTGAt              <  1:7422552/36‑1 (MQ=255)
                         gtggctgGGGAAAGTGAAAGGTGAACAGTTTGATCg            <  1:2163617/36‑1 (MQ=255)
                          tggctgGGGAAAGTGAAAGGTGAACAGTTTGATCGc           >  1:3615666/1‑36 (MQ=255)
                           ggctgGGGAAAGTGAAAGGTGAACAGTTTGATCGcc          >  1:5418904/1‑36 (MQ=255)
                            gctgGGGAAAGTGAAAGGTGAACAGTTTGATCGCCt         >  1:4498994/1‑36 (MQ=255)
                            gctgGGGAAAGTGAAAGGTGAACAGTTTGATCGCCt         >  1:1432440/1‑36 (MQ=255)
                            gctgGGGAAAGTGAAAGGTGAACAGTTTGATCGCCt         >  1:319028/1‑36 (MQ=255)
                              tgGGGGAAGTGAAAGGTGAACAGTTTGATCGCCTGa       <  1:4378933/36‑1 (MQ=255)
                               ggggAAAGTGAAAGGTGAACAGTTTGATCGCCTGAg      <  1:6829185/36‑1 (MQ=255)
                               ggggAAAGTGAAAGGTGAACAGTTTGATCGCCTGAg      >  1:7408645/1‑36 (MQ=255)
                               ggggAAAGTGAAAGGTGAACAGTTTGATCGCCTGAg      <  1:5621057/36‑1 (MQ=255)
                                gggAAAGTGAAAGGTGAACAGTTTGATCGCCTGAGc     <  1:4405060/36‑1 (MQ=255)
                                  gAAAGTGAAATGTGAACAGTTTGATCGCCTGAGCct   <  1:2427491/36‑1 (MQ=255)
                                  gAAAGTGAAAGGTGAACAGTTTGATCGCCTGAGCct   >  1:3080095/1‑36 (MQ=255)
                                  gAAAGTGAAAGGTGAACAGTTTGATCGCCTGAGCct   <  1:1199134/36‑1 (MQ=255)
                                  gAAAGTGAAAGGTGAACAGTTTGATCGCCTGAGCct   <  1:1138492/36‑1 (MQ=255)
                                   aaaGTGAAAGGTGAACAGTTTGATCGCCTGAGCctg  >  1:2824570/1‑36 (MQ=255)
                                   aaaGTGAAAGGTGAACAGTTTGATCGCCTAAGCctg  <  1:6568853/36‑1 (MQ=255)
                                   |                                   
TGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAGTTTGATCGCCTGAGCCTG  >  REL606/3977958‑3978028

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 17 ≤  ATCG  < 28 ≤  ATCG  < 32 ≤  ATCG