Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,168,1360C→G27.8% 10.536268pabC4‑amino‑4‑deoxychorismate lyase
Reads supporting (aligned to +/- strand):  new base (10/0):  ref base (4/22):  total (14/22)

TACTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCAGCAACG  >  REL606/1168101‑1168171
                                   |                                   
ttcTGAAAGTCGTGATCAGTCGCGGTAGTggcgggc                                     <  1:124019/34‑1 (MQ=255)
  cTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGa                                   <  2:3517297/36‑1 (MQ=255)
  cTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGa                                   <  2:1192211/36‑1 (MQ=255)
    ggaaGTGGTGATCAGTGGCGGTAGTGGCGGGCGAgg                                 <  1:647669/34‑1 (MQ=37)
    gAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAgg                                 >  1:3431262/1‑36 (MQ=255)
     ataGTCGTGATCAGTCGCGGTAGTGGCGGGCGAggg                                <  1:2123239/34‑1 (MQ=255)
      aaGTCGTGATCAGTCGCGGTAGTGGCGGGGGGgggg                               >  1:779311/1‑35 (MQ=37)
      aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGt                               <  1:1047390/36‑1 (MQ=255)
      aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGt                               <  2:2239398/36‑1 (MQ=255)
      aaGTCGGGATCAGTCGCGGTAGTGGCGGGCGAGGGt                               <  1:778018/36‑1 (MQ=255)
       aGTTGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTa                              <  2:594622/36‑1 (MQ=255)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGGGGggggg                              >  2:3691372/1‑34 (MQ=25)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGGGGggggg                              >  2:3194500/1‑34 (MQ=25)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGGGGggggg                              >  2:2860515/1‑34 (MQ=25)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGGGGggggg                              >  2:1196094/1‑34 (MQ=25)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGGGCggggg                              >  2:3125403/1‑34 (MQ=25)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTa                              <  1:843775/36‑1 (MQ=255)
       aGTCGAGATCAGTCGCGGTAGTGGCGGGCGAGGGTa                              <  2:3090177/36‑1 (MQ=255)
        gTCGTGATCAGTCGCGGTAGTGGGGGGGGGgggggg                             >  1:3516655/1‑33 (MQ=16)
        gTCGTGATCAGTCGCGGTAGTGGCGGGGGGggggag                             >  2:6267/1‑33 (MQ=25)
           gTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGc                          <  2:1624618/36‑1 (MQ=255)
            tGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGca                         <  1:157636/36‑1 (MQ=255)
               tCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACAt                      <  2:3314042/36‑1 (MQ=255)
                cAGGCGCGGTAGTGGCGGGCGAGGGTACAGCACAtt                     <  2:1860278/36‑1 (MQ=255)
                    cgcgGAAGTGGCGGGCGGGGGGACAGCACATTTAAc                 >  1:1753846/1‑36 (MQ=16)
                     gcgGTAGTGGCGGGCGAGGGTACAGCACATTGAACa                <  2:2423001/36‑1 (MQ=255)
                     gcgGTAGTGGCGGGCGAGGGTACAGCACATTGAACa                <  2:2488714/36‑1 (MQ=255)
                      cgGTAGTGGCGGGCGAGGGTACAGCACATTGAACAg               <  1:1779707/36‑1 (MQ=255)
                       ggTAGTGGCGGGGGGGGGGGCAGGACATTGAAAAGc              >  2:2198882/1‑36 (MQ=12)
                        gTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCg             <  2:1982636/36‑1 (MQ=255)
                          aGTGGCGGGCGAGGGGAAAGCACATCAAACAACggc           >  1:2606059/1‑35 (MQ=12)
                           gTGGCGGGCGAGGGTACAGCGCATTGAACAGCGGAc          <  1:3145155/36‑1 (MQ=255)
                           gTGGCGGGCGAGGGTACAGCACATTGAACAGCGGAc          <  1:1632377/36‑1 (MQ=255)
                            tGGCGGGCGAGGGTACAGCACATTGAACAGCGGAcc         <  1:2651523/36‑1 (MQ=255)
                            tGGCGGGCGAGGGTACAGCACATTGAACAGCGGAcc         >  1:3170053/1‑36 (MQ=255)
                            tGGCGGGCGAGGGTACAGCACATTGAACAGCGGAcc         <  1:514094/36‑1 (MQ=255)
                            tGGCGGGCGAGGGTACAGCACATTGAAAAGCGGagc         >  1:98873/1‑34 (MQ=38)
                             gacgggcgAGGGTACAGCACATTGAACAGCGGACCa        <  1:2700472/34‑1 (MQ=255)
                              gcgggcgAGGGTACAGCACATTGAACAGCGGACCAg       <  2:2388408/36‑1 (MQ=255)
                               cgggcgAGGGTACAGCACATTGAACAGCGGACCAGc      <  2:517847/36‑1 (MQ=255)
                               cgggcgAGGGTACAGCACATTGAACAGCGGACCAGc      <  2:1233456/36‑1 (MQ=255)
                               cgggcgAGGGTACAGCACATTGAACAGCGGACCAGc      <  1:3173333/36‑1 (MQ=255)
                                gggcgAGGGTACAGCACATTGAACAGCGGACCAGCa     >  2:3475469/1‑36 (MQ=255)
                                   cgAGGGTACAGCACATTGAACAGCGGACCAGCAAcg  >  1:2574903/1‑36 (MQ=255)
                                   cgAGGGTACAGCACATTGAACAGCGGACCAGCAAcg  >  2:977412/1‑36 (MQ=255)
                                   cgAGGGTACAGCACAATGAAAAGCGGACCAGCAAcg  <  2:32384/36‑1 (MQ=37)
                                   |                                   
TACTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCAGCAACG  >  REL606/1168101‑1168171

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 13 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG