Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,977,9930A→G27.1% 15.760514metE5‑methyltetrahydropteroyltriglutamate‑ homocysteine methyltransferase
Reads supporting (aligned to +/- strand):  new base (15/1):  ref base (12/31):  total (27/33)

TGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAGTTTGATCGCCTGAGCCTG  >  REL606/3977958‑3978028
                                   |                                   
tgctgtTGGGGCCGGTTACCTGGCTGTGGCTGGGGa                                     <  2:2768090/36‑1 (MQ=255)
tgctgcTGGGGCCTGTTACCTGGCTGTGGCTGGGGa                                     <  2:1083768/36‑1 (MQ=255)
tgctgcTGGGGCCGGTTACCTGGCTGTGGCTGGGGa                                     <  1:939330/36‑1 (MQ=255)
tgctgcTGGGGCCGGTTACCTGGCTGTGGCTGGGGa                                     <  2:3261902/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                                   <  2:1356317/36‑1 (MQ=255)
  cgtctgGGGCCGGTTACCTGGCTGTGGCTGGGGaaa                                   <  1:3479288/33‑1 (MQ=255)
   tgctgGGGCCGGTTACCTGGCTGTGGCTGGGGGAAg                                  >  2:2446135/1‑36 (MQ=255)
   tgctgGGGCCGGTTACCTGGCTGTGGCTGGGGAGAg                                  >  1:2452591/1‑36 (MQ=255)
   tgctgGGGCCGGTTACCTGCCTGTGGCTGGGGGAAg                                  >  1:1149674/1‑36 (MQ=37)
    gctgGGGCCGGTTACCTGGCTGTGGCTGGGGGAAgg                                 >  1:2922318/1‑35 (MQ=255)
    gctgGGGCCGGTTACCTGGCTGTGGCTGGGGGAAGt                                 >  1:695094/1‑36 (MQ=255)
     ctgGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTg                                <  1:97442/36‑1 (MQ=255)
     ctgGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTg                                <  2:3709708/36‑1 (MQ=255)
      tgGGGCGGGTTACCTGGCTGTGGCTGGGGAAAGTGa                               <  1:1664287/36‑1 (MQ=255)
       ggggCCGGTTACCTGGCTGTGGCTGGGGGTAGTgtt                              >  1:1384197/1‑34 (MQ=21)
       aaggCCGGTTACCTGGCTGTGGCTGGGGAAAGTGaa                              <  1:1046200/34‑1 (MQ=255)
        gggCTGGTTACCTGGCTGTGGCTGGGGGAAGTGaaa                             >  1:1259066/1‑36 (MQ=25)
        gggCCGGTTACCTGGCTGTGGCTGGGGAAAGTGaaa                             <  1:293039/36‑1 (MQ=255)
          gCCGGTTACCTGGCTGTGGCTGGGGGAAGTGGAAgg                           >  2:366045/1‑36 (MQ=25)
          gCCGGTTACCTGGCTGTGGCTGGGGGAAGTGGAAgg                           >  1:1969595/1‑36 (MQ=25)
          gCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAgg                           <  1:588415/36‑1 (MQ=255)
           ccGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGt                          <  2:2917984/36‑1 (MQ=255)
            cGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTg                         >  1:3441203/1‑36 (MQ=255)
            cGGTTACCTGGCTGTGGCTGGGGGAAGTGAAAAGTg                         >  2:1894069/1‑36 (MQ=25)
            cGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTg                         <  2:60042/36‑1 (MQ=255)
             ggTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGTGa                        >  1:2578449/1‑36 (MQ=255)
             ggTTACCGGGCTGTGGCTGGGGAAACGGAAAGGTGa                        <  2:2790214/36‑1 (MQ=21)
              gTTACCTGGCTGTGGCTGGGGGAAGTGAAAGGGGaa                       >  2:3486081/1‑36 (MQ=25)
              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=25)
                tgccTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACa                     <  1:1865793/34‑1 (MQ=255)
                 aCCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAg                    >  2:1870541/1‑36 (MQ=255)
                 aCCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAg                    <  2:3748405/36‑1 (MQ=255)
                 aCCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAg                    <  1:3587211/36‑1 (MQ=255)
                  ccTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAGt                   <  1:1002970/36‑1 (MQ=255)
                   cTGGCTGTGGCTGGGGTAAGGGAAAGGTGAACAGtt                  <  1:1072096/36‑1 (MQ=37)
                   cTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAGtt                  <  1:1873224/36‑1 (MQ=255)
                    tgggtgtggctgGGGAAAGTGAAAGGTGAACAGttt                 <  1:3342556/36‑1 (MQ=255)
                    tggctgtggctgGGGGAAGTGAAAGGGGAACAGttt                 >  1:600998/1‑36 (MQ=37)
                    tggctgtggctgGGGAAAGTGAAAGGTGAACAGttt                 <  2:2006318/36‑1 (MQ=255)
                    tggctgtggctgGGGAAAGTGAAAGGTGAACAGttt                 <  2:2331239/36‑1 (MQ=255)
                     ggctgtggctgGGGGAAATGAAAAGGGAACAGtttt                >  1:144817/1‑35 (MQ=16)
                     ggctgtggctgGGGAAAGTGAAAGGTGAACAGTTTg                <  2:1998828/36‑1 (MQ=255)
                      gctgtggctgGGGAAAGTGAAAGGTGAACAGTTTGa               <  2:1377480/36‑1 (MQ=255)
                      gctgtggctgGGGAAAGTGAAAGGTGAACAGTTTGa               <  1:2739524/36‑1 (MQ=255)
                       gtgtggctgGGGAAAGTGAAAGGTGAACAGTTTGAt              <  1:779722/35‑1 (MQ=255)
                       ctgtggctgGGGAAAGTGAAAGGTGAACAGTTTGAt              <  2:3683342/36‑1 (MQ=255)
                         gtggctgGGGGAAGTGAAAGGGGAACAGTTTGTTCg            >  2:1636213/1‑36 (MQ=25)
                         gtggctgGGGAAAGTGAAAGGTGAACAGTTTGATCg            <  1:2163617/36‑1 (MQ=255)
                          tggctgGGGAAAGTGAAAGGTGAACAGTTTGATCGc           >  1:3615666/1‑36 (MQ=255)
                           ggctgGGGAAAGTGAAAGGTGAACAGTTTGATCGcc          >  2:1679694/1‑36 (MQ=255)
                            gctgGGGAAAGTGAAAGGTGAACAGTTTGATCGCCt         >  2:759784/1‑36 (MQ=255)
                            gctgGGGAAAGTGAAAGGTGAACAGTTTGATCGCCt         >  1:319028/1‑36 (MQ=255)
                            gctgGGGAAAGTGAAAGGTGAACAGTTTGATCGCCt         >  1:1432440/1‑36 (MQ=255)
                            gctgGGCAAAGTGAAAGGAGAACAGTTTGACCGTCt         >  2:473183/1‑36 (MQ=19)
                              tgGGGGAAGTGAAAGGTGAACAGTTTGATCGCCTGa       <  2:639723/36‑1 (MQ=255)
                               ggggAAAGTGAAAGGTGAACAGTTTGATCGCCTGAg      >  2:3669435/1‑36 (MQ=255)
                               ggggAAAGTGAAAGGTGAACAGTTTGATCGCCTGAg      <  2:3089975/36‑1 (MQ=255)
                               ggggAAAGTGAAAGGTGAACAGTTTGATCGCCTGAg      <  2:1881847/36‑1 (MQ=255)
                                gggAAAGTGAAAGGTGAACAGTTTGATCGCCTGAGc     <  2:665850/36‑1 (MQ=255)
                                  gAAAGTGAAATGTGAACAGTTTGATCGCCTGAGCct   <  1:2427491/36‑1 (MQ=255)
                                  gAAAGTGAAAGGTGAACAGTTTGATCGCCTGAGCct   <  1:1199134/36‑1 (MQ=255)
                                  gAAAGTGAAAGGTGAACAGTTTGATCGCCTGAGCct   >  1:3080095/1‑36 (MQ=255)
                                  gAAAGTGAAAGGTGAACAGTTTGATCGCCTGAGCct   <  1:1138492/36‑1 (MQ=255)
                                   aaaGTGAAAGGTGAACAGTTTGATCGCCTGAGCctg  >  1:2824570/1‑36 (MQ=255)
                                   aaaGTGAAAGGTGAACAGTTTGATCGCCTAAGCctg  <  2:2829643/36‑1 (MQ=255)
                                   |                                   
TGCTGCTGGGGCCGGTTACCTGGCTGTGGCTGGGGAAAGTGAAAGGTGAACAGTTTGATCGCCTGAGCCTG  >  REL606/3977958‑3978028

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 15 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG