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

ACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATAT  >  REL606/2698059‑2698126
                                  |                                 
ccTCGCCATGTCGATCAGGGGCTGTCGCTGACGCtt                                  <  2:2755235/35‑1 (MQ=255)
aCTCGCCATGTCGATCAGGGGCTGTCGCTGACGCtt                                  <  2:1492623/36‑1 (MQ=255)
aCTCGCCAGGTCGATCAGGGGCTGTCGCAGACGCtt                                  <  1:2650639/36‑1 (MQ=37)
  tCGCCATGTCGATCAGGGGCTGTCGCTGACGCtttt                                <  2:604760/36‑1 (MQ=255)
  tCGCCATGTCGATCAGGGGCTGTCGCTGAAGCtttt                                >  2:3475861/1‑36 (MQ=255)
   cGCCATGTCGATCAGGGGCTGTCGCTGACGCttttt                               <  1:614664/36‑1 (MQ=255)
    gCCATGTCGATCAGGGGCTGTCGCTGACGCtttttt                              >  2:3300258/1‑36 (MQ=255)
    gCCATGTCGAGCAGGGGCTGTCGCTGACGCtttttt                              <  1:1198469/36‑1 (MQ=255)
    gCCATGTAGATCAGGGGCTGTCGCTTACGCtttttt                              <  1:1206852/36‑1 (MQ=37)
     ccATGTCGATCAGGGGCTGTCGCTGACGGttttttc                             >  2:1151093/1‑35 (MQ=255)
     ccATGTCGATCAGGGGCTGTCGCTGACGCttttttc                             <  1:2402999/36‑2 (MQ=255)
      cATGTCGATCAGGGGCTGTCGCTGACGCttttttcc                            <  2:3504075/36‑3 (MQ=255)
       aTGTCGATCAGGGGCTGTCGCTGACGCttttttctc                           >  1:2247351/1‑33 (MQ=255)
       aTGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCcc                           >  2:2802868/1‑36 (MQ=255)
        tGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCccc                          <  1:267738/36‑1 (MQ=255)
         gTCGATCAGGGGCTGTCGCTGACGCTTTTTTCccct                         <  2:1629782/36‑2 (MQ=255)
           cGATCAGGGGCTGTCGCTGACGCTTTTTTTccccac                       >  1:96053/1‑34 (MQ=255)
           cGATCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGAt                      >  2:247878/1‑36 (MQ=37)
           cGATCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGAt                      <  2:922448/36‑1 (MQ=37)
           cGATCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGAt                      <  1:3436718/36‑1 (MQ=37)
            gATCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATa                     <  2:2866520/36‑1 (MQ=37)
            gATCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATa                     >  1:2379330/1‑36 (MQ=37)
            gATCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATa                     >  2:1386453/1‑36 (MQ=37)
            gATCAGGGGCTGTCGCTGACACTTTTTTCccccata                      >  2:1661037/1‑33 (MQ=25)
            gATCAGGGGCTGTCGCCGACGCTCTTTCCccccata                      >  1:3392068/1‑33 (MQ=12)
             aTCGGGGGCTGTCGCTGACGC.TTTTTTCCCCGATAc                    <  2:2246937/36‑1 (MQ=37)
             aTCAGGGGCTGTCGCTGACGGTTTTTTT.CCCGATAc                    >  1:1614721/1‑36 (MQ=37)
             aTCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATAc                    <  1:1700190/36‑1 (MQ=38)
              tCAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATAcc                   >  2:2150355/1‑36 (MQ=38)
               cAGGGGCTGTCGCTGACGCTTTTTTT.CCCGATACCg                  >  1:470537/1‑36 (MQ=37)
               cAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATACCg                  >  2:2056503/1‑36 (MQ=38)
               cAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATACCg                  <  2:1849279/36‑1 (MQ=38)
               cAGGGGCTGTCGCTGACGC.TTTTTTCCCCGATACCg                  >  2:3392666/1‑36 (MQ=38)
                aGGGGCTGTCGCTGACGCTTTTTTT.CCCGATACCGc                 >  2:2234280/1‑36 (MQ=37)
                 ggggCTGTCGCTGGCGC.TTTTTTCCCCGATAACggc                >  1:1842489/1‑34 (MQ=16)
                    gCTGTCGCTGACG.TTTTTTTCCCCGATACCGccacc             <  1:2420846/36‑1 (MQ=38)
                    gCTCTCCCTGACGC.TTTTTCCCCCAATACCGccacc             <  2:3069526/36‑1 (MQ=12)
                     cTGCCGCTGCCGC.TTTTTTCCCCGATACCGccacca            <  2:1561960/36‑1 (MQ=18)
                       gTCGCTGACGC.TTTTTTCCCCGATACCGCCACCACt          >  2:90828/1‑36 (MQ=38)
                       gTCGCTGACGC.TTTTTTCCCCGATACCGCCACCACt          >  1:3547283/1‑36 (MQ=38)
                        tCGCTGACGC.TTTTTTCCCCGATACCGCCACCACTc         >  1:1047454/1‑36 (MQ=38)
                          gCTGACGC.TTTTTTCCCCGATACCGCCACCACAcgc       >  2:2176738/1‑36 (MQ=37)
                          ccTGACG.TTTTTTCCCCCGATACCGCCACCACTcgc       <  2:2078632/35‑1 (MQ=25)
                           cTGACGC.TTTTTTCCCCGATACCGCCACCACTcgcg      >  2:694084/1‑36 (MQ=38)
                            tGACGC.TTTTTTCCCCGATACCGCCACCACTCGCGa     >  1:1892146/1‑36 (MQ=37)
                            tGACGC.TTTTTTCCCCGATACCGCCACCACTCGCGa     >  2:1661202/1‑36 (MQ=37)
                            tGACGC.TTTTTTCCCCGATACCGCCACCACTCCCGa     >  2:77098/1‑36 (MQ=37)
                             ggagCTTTTTTTCCGGGTTACCGCCACCACTCGCGa     <  1:3486956/33‑1 (MQ=16)
                             gACGC.TTTTTTCCCCGATACCGCCACCACTCGCGat    <  2:2401328/36‑1 (MQ=37)
                               acgtttttttCCCCAATACCGCCACCACTCGCGata   <  1:1687672/33‑1 (MQ=37)
                                cgtttttttCCCCGATCCCGCCACCACTCGCGatat  <  2:8582/34‑1 (MQ=37)
                                  |                                 
ACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATAT  >  REL606/2698059‑2698126

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