Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,698,0930T→.65.7% 55.735coding (1994/2145 nt)nrdEribonucleotide‑diphosphate reductase alpha subunit
Reads supporting (aligned to +/- strand):  new base (15/8):  ref base (7/5):  total (22/13)

ACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATAT  >  REL606/2698059‑2698126
                                  |                                 
ccTCGCCATGTCGATCAGGGGCTGTCGCTGACGCtt                                  <  1:6494445/35‑1 (MQ=255)
aCTCGCCATGTCGATCAGGGGCTGTCGCTGACGCtt                                  <  1:5231833/36‑1 (MQ=255)
aCTCGCCAGGTCGATCAGGGGCTGTCGCAGACGCtt                                  <  1:2650639/36‑1 (MQ=25)
  tCGCCATGTCGATCAGGGGCTGTCGCTGACGCtttt                                <  1:4343970/36‑1 (MQ=255)
  tCGCCATGTCGATCAGGGGCTGTCGCTGAAGCtttt                                >  1:7215071/1‑36 (MQ=255)
   cGCCATGTCGATCAGGGGCTGTCGCTGACGCttttt                               <  1:614664/36‑1 (MQ=255)
    gCCATGTCGATCAGGGGCTGTCGCTGACGCtttttt                              >  1:7039468/1‑36 (MQ=255)
    gCCATGTCGAGCAGGGGCTGTCGCTGACGCtttttt                              <  1:1198469/36‑1 (MQ=255)
    gCCATGTAGATCAGGGGCTGTCGCTTACGCtttttt                              <  1:1206852/36‑1 (MQ=25)
     ccATGTCGATCAGGGGCTGTCGCTGACGGttttttc                             >  1:4890303/1‑35 (MQ=255)
     ccATGTCGATCAGGGGCTGTCGCTGACGCttttttc                             <  1:2402999/36‑2 (MQ=255)
      cATGTCGATCAGGGGCTGTCGCTGACGCttttttcc                            <  1:7243285/36‑3 (MQ=255)
       aTGTCGATCAGGGGCTGTCGCTGACGCttttttctc                           >  1:2247351/1‑33 (MQ=255)
       aTGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCcc                           >  1:6542078/1‑36 (MQ=255)
        tGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCccc                          <  1:267738/36‑1 (MQ=255)
         gTCGATCAGGGGCTGTCGCTGACGCTTTTTTCccct                         <  1:5368992/36‑2 (MQ=255)
           cGATCAGGGGCTGTCGCTGACGCTTTTTTTccccac                       >  1:96053/1‑34 (MQ=255)
           cGATCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGAt                      <  1:3436718/36‑1 (MQ=25)
           cGATCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGAt                      >  1:3987088/1‑36 (MQ=25)
           cGATCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGAt                      <  1:4661658/36‑1 (MQ=25)
            gATCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATa                     <  1:6605730/36‑1 (MQ=25)
            gATCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATa                     >  1:2379330/1‑36 (MQ=25)
            gATCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATa                     >  1:5125663/1‑36 (MQ=25)
            gATCAGGGGCTGTCGCTGACACTTTTTTCccccata                      >  1:5400247/1‑33 (MQ=21)
            gATCAGGGGCTGTCGCCGACGCTCTTTCCccccata                      >  1:3392068/1‑33 (MQ=12)
             aTCGGGGGCTGTCGCTGACGC.TTTTTTCCCCGATAc                    <  1:5986147/36‑1 (MQ=21)
             aTCAGGGGCTGTCGCTGACGGTTTTTTT.CCCGATAc                    >  1:1614721/1‑36 (MQ=21)
             aTCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATAc                    <  1:1700190/36‑1 (MQ=25)
              tCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATAcc                   >  1:5889565/1‑36 (MQ=25)
               cAGGGGCTGTCGCTGACGCTTTTTTT.CCCGATACCg                  >  1:470537/1‑36 (MQ=25)
               cAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATACCg                  <  1:5588489/36‑1 (MQ=25)
               cAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATACCg                  >  1:5795713/1‑36 (MQ=25)
               cAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATACCg                  >  1:7131876/1‑36 (MQ=25)
                aGGGGCTGTCGCTGACGCTTTTTTT.CCCGATACCGc                 >  1:5973490/1‑36 (MQ=25)
                 ggggCTGTCGCTGGCGC.TTTTTTCCCCGATAACggc                >  1:1842489/1‑34 (MQ=17)
                    gCTGTCGCTGACG.TTTTTTTCCCCGATACCGccacc             <  1:2420846/36‑1 (MQ=25)
                     cTGCCGCTGCCGC.TTTTTTCCCCGATACCGccacca            <  1:5301170/36‑1 (MQ=17)
                       gTCGCTGACGC.TTTTTTCCCCGATACCGCCACCACt          >  1:3830038/1‑36 (MQ=37)
                       gTCGCTGACGC.TTTTTTCCCCGATACCGCCACCACt          >  1:3547283/1‑36 (MQ=37)
                        tCGCTGACGC.TTTTTTCCCCGATACCGCCACCACTc         >  1:1047454/1‑36 (MQ=37)
                          gCTGACGC.TTTTTTCCCCGATACCGCCACCACAcgc       >  1:5915948/1‑36 (MQ=25)
                          ccTGACG.TTTTTTCCCCCGATACCGCCACCACTcgc       <  1:5817842/35‑1 (MQ=21)
                           cTGACGC.TTTTTTCCCCGATACCGCCACCACTcgcg      >  1:4433294/1‑36 (MQ=25)
                            tGACGC.TTTTTTCCCCGATACCGCCACCACTCGCGa     >  1:1892146/1‑36 (MQ=25)
                            tGACGC.TTTTTTCCCCGATACCGCCACCACTCGCGa     >  1:5400412/1‑36 (MQ=25)
                            tGACGC.TTTTTTCCCCGATACCGCCACCACTCCCGa     >  1:3816308/1‑36 (MQ=25)
                             ggagCTTTTTTTCCGGGTTACCGCCACCACTCGCGa     <  1:3486956/33‑1 (MQ=17)
                             gACGC.TTTTTTCCCCGATACCGCCACCACTCGCGat    <  1:6140538/36‑1 (MQ=25)
                               acgtttttttCCCCAATACCGCCACCACTCGCGata   <  1:1687672/33‑1 (MQ=25)
                                cgtttttttCCCCGATCCCGCCACCACTCGCGatat  <  1:3747792/34‑1 (MQ=25)
                                  |                                 
ACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATAT  >  REL606/2698059‑2698126

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 10 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG