Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,168,1360C→G26.5% 10.834R90G (CGA→GGA) pabC4‑amino‑4‑deoxychorismate lyase
Reads supporting (aligned to +/- strand):  new base (9/0):  ref base (3/22):  total (12/22)

TACTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCAGCAACG  >  REL606/1168101‑1168171
                                   |                                   
ttcTGAAAGTCGTGATCAGTCGCGGTAGTggcgggc                                     <  1:124019/34‑1 (MQ=255)
  cTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGa                                   <  1:4931421/36‑1 (MQ=255)
  cTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGa                                   <  1:7256507/36‑1 (MQ=255)
    ggaaGTGGTGATCAGTGGCGGTAGTGGCGGGCGAgg                                 <  1:647669/34‑1 (MQ=21)
    gAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAgg                                 >  1:3431262/1‑36 (MQ=255)
     ataGTCGTGATCAGTCGCGGTAGTGGCGGGCGAggg                                <  1:2123239/34‑1 (MQ=255)
      aaGTCGTGATCAGTCGCGGTAGTGGCGGGGGGgggg                               >  1:779311/1‑35 (MQ=21)
      aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGt                               <  1:1047390/36‑1 (MQ=255)
      aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGt                               <  1:5978608/36‑1 (MQ=255)
      aaGTCGGGATCAGTCGCGGTAGTGGCGGGCGAGGGt                               <  1:778018/36‑1 (MQ=255)
       aGTTGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTa                              <  1:4333832/36‑1 (MQ=255)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGGGGggggg                              >  1:6599725/1‑34 (MQ=21)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGGGGggggg                              >  1:6933710/1‑34 (MQ=21)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGGGGggggg                              >  1:4935304/1‑34 (MQ=21)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGGGGggggg                              >  1:7430582/1‑34 (MQ=21)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGGGCggggg                              >  1:6864613/1‑34 (MQ=21)
       aGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTa                              <  1:843775/36‑1 (MQ=255)
       aGTCGAGATCAGTCGCGGTAGTGGCGGGCGAGGGTa                              <  1:6829387/36‑1 (MQ=255)
        gTCGTGATCAGTCGCGGTAGTGGGGGGGGGgggggg                             >  1:3516655/1‑33 (MQ=18)
        gTCGTGATCAGTCGCGGTAGTGGCGGGGGGggggag                             >  1:3745477/1‑33 (MQ=21)
           gTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGc                          <  1:5363828/36‑1 (MQ=255)
            tGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGca                         <  1:157636/36‑1 (MQ=255)
               tCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACAt                      <  1:7053252/36‑1 (MQ=255)
                cAGGCGCGGTAGTGGCGGGCGAGGGTACAGCACAtt                     <  1:5599488/36‑1 (MQ=255)
                     gcgGTAGTGGCGGGCGAGGGTACAGCACATTGAACa                <  1:6162211/36‑1 (MQ=255)
                     gcgGTAGTGGCGGGCGAGGGTACAGCACATTGAACa                <  1:6227924/36‑1 (MQ=255)
                      cgGTAGTGGCGGGCGAGGGTACAGCACATTGAACAg               <  1:1779707/36‑1 (MQ=255)
                        gTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCg             <  1:5721846/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:514094/36‑1 (MQ=255)
                            tGGCGGGCGAGGGTACAGCACATTGAACAGCGGAcc         >  1:3170053/1‑36 (MQ=255)
                            tGGCGGGCGAGGGTACAGCACATTGAAAAGCGGagc         >  1:98873/1‑34 (MQ=37)
                             gacgggcgAGGGTACAGCACATTGAACAGCGGACCa        <  1:2700472/34‑1 (MQ=255)
                              gcgggcgAGGGTACAGCACATTGAACAGCGGACCAg       <  1:6127618/36‑1 (MQ=255)
                               cgggcgAGGGTACAGCACATTGAACAGCGGACCAGc      <  1:4972666/36‑1 (MQ=255)
                               cgggcgAGGGTACAGCACATTGAACAGCGGACCAGc      <  1:4257057/36‑1 (MQ=255)
                               cgggcgAGGGTACAGCACATTGAACAGCGGACCAGc      <  1:3173333/36‑1 (MQ=255)
                                gggcgAGGGTACAGCACATTGAACAGCGGACCAGCa     >  1:7214679/1‑36 (MQ=255)
                                   cgAGGGTACAGCACATTGAACAGCGGACCAGCAAcg  >  1:4716622/1‑36 (MQ=255)
                                   cgAGGGTACAGCACATTGAACAGCGGACCAGCAAcg  >  1:2574903/1‑36 (MQ=255)
                                   cgAGGGTACAGCACAATGAAAAGCGGACCAGCAAcg  <  1:3771594/36‑1 (MQ=255)
                                   |                                   
TACTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCAGCAACG  >  REL606/1168101‑1168171

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 28 ≤  ATCG  < 32 ≤  ATCG