Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,621,2740A→G26.0% 15.850302trpRTrp operon repressor
Reads supporting (aligned to +/- strand):  new base (13/0):  ref base (7/30):  total (20/30)

GAAAGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGAAGAGGTGTTGCTGAAAAGCGATTGATTTTGTAGGC  >  REL606/4621239‑4621309
                                   |                                   
gAAAGCCGCGCCGGTTGAGCTGCGCCAGTGGCTGGa                                     <  1:622302/36‑1 (MQ=255)
gAAAGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGa                                     <  1:177102/36‑1 (MQ=255)
gAAAGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGa                                     <  2:3429324/36‑1 (MQ=255)
 aaaGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGaa                                    <  2:774919/36‑1 (MQ=255)
 aaaGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGaa                                    <  2:1333281/36‑1 (MQ=255)
  aaGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGGag                                   >  2:2056749/1‑36 (MQ=255)
  aaGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGGag                                   >  2:1164111/1‑36 (MQ=255)
  aaGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGAag                                   <  2:664504/36‑1 (MQ=255)
  aaGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGAag                                   <  1:1352542/36‑1 (MQ=255)
   aGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGGagg                                  >  2:2900031/1‑35 (MQ=255)
     ccGCGCGCGTTGGGCTGCGCCAGTGGCTGGAAGAgg                                <  1:503535/36‑1 (MQ=37)
      cgctCCCGTTGAGCTGCGCCAGTGGCTGGAAGAGgt                               <  2:2466555/36‑1 (MQ=255)
      cgcgCCCGTTGAGCTGCGCCAGTGGCTGGGagaagg                               >  1:3345618/1‑33 (MQ=37)
      cgcgCCCGTTGAGCTGCGCCAGTGGCTGGAAGAGgt                               <  1:1905153/36‑1 (MQ=255)
      cgcgCCCGTTGAGCTGCGCCAGTGGCTGGAAGAGgt                               <  2:2139603/36‑1 (MQ=255)
       gcgcCCGTTGAGCTGCGCCAGTGGCTGGAAGAGgtg                              <  2:389563/36‑1 (MQ=255)
       gcgcCCGTTGAGCTGCGCCAGTGGCTGGAAGAGgtg                              <  2:2801456/36‑1 (MQ=255)
       gcgcCCGTTGAGCTGCGCCAGTGGCTGGAAGAGgtg                              <  2:1926846/36‑1 (MQ=255)
         gcCCGTTGAGCTGCGCCAGTGGCTGGAAGAGGTGtt                            <  2:1423021/36‑1 (MQ=255)
         gcCCGTTGAGCTGCGCCAGTGGCTGGAAAAGGCGtt                            >  2:2261360/1‑36 (MQ=37)
           ccGTTGAGCTGCGCAAGTGGCTGGGAGAGgtgtggg                          >  1:2926453/1‑33 (MQ=25)
             gTTGAGGGGCGCCAGTGGCTGGAAGAGGTGTTGCTg                        <  2:825027/36‑1 (MQ=25)
             gTTGAGCTGCGCCAGTGGCTGGGGGAGGTGTTGTTg                        >  1:1802749/1‑36 (MQ=25)
               tGAGCTGCGCCAGTGGCTGGAAGAGGTGTTGCTGaa                      <  1:2627413/36‑1 (MQ=255)
               tGAGCTGCGCCAGTGGCTGGAAGAGGGGTTGGTGaa                      >  2:2106232/1‑36 (MQ=37)
               tGAGCGGCGCCAGTGGCGGGAAGACGTGTTGCTGaa                      <  2:1081457/36‑1 (MQ=21)
                gAGTTGCGCCAGTGGCTGGAAGAGGTGTTGCTGaaa                     <  2:349921/36‑1 (MQ=255)
                gAGCTGCGCCAGTGGCTGGGAGAGGGGTTGGTGaga                     >  1:1131575/1‑34 (MQ=21)
                gAGCTGCGCCAGTGGCTGGAAGAGGTGTTGCTGaaa                     <  2:1960742/36‑1 (MQ=255)
                gAGCTGCGCCAGTGGCTGGAAGAGGTGTTGCTGaaa                     <  1:3467783/36‑1 (MQ=255)
                 aGCTGCGCCAGTGGCTGGAAGAGGTGTTGCTGaaaa                    <  1:3658115/36‑1 (MQ=255)
                 aGCTACGCCAGTGGCTGGGAGAGGGGCTGCTGaaaa                    >  1:2208348/1‑36 (MQ=19)
                     gcgcCAGTGGCTGGGAGGGGGGGTGCTGGAAAGCgg                >  2:387289/1‑35 (MQ=12)
                     gcgcCAGTGGCTGGAAGAGGTGTTGCTGAAAAGCga                <  1:1539456/36‑1 (MQ=255)
                     gcgcCAGTGGCTGGAAGAGGTGTTGCTGAAAAGCga                <  2:3108540/36‑1 (MQ=255)
                     gcgcCAGTGGCTGGAAGAGGTGTTGCTGAAAAGCga                <  1:2106405/36‑1 (MQ=255)
                      cgcCAGTGGCTGGGAGAGGGGTTGCTGGAAAGCgat               >  2:2035271/1‑36 (MQ=25)
                      cgcCAGTGGCTGGAAGAGGGGTTGGTTGAAAGCggt               >  2:1184376/1‑34 (MQ=16)
                       gcCAGTGGCTGGGAGAGGTGTTGCTGGAAAGAgatt              >  2:453959/1‑36 (MQ=25)
                       gcCAGTGGCTGGAAGGGGTGTTGCTGAAAAGCgatt              <  2:3551173/36‑1 (MQ=255)
                        cgaGTGGCTGGAAGAGGTGTTGCTGAAAAGCgattg             <  1:2480069/34‑1 (MQ=255)
                        ccAGTGGCTGGGAGAGGGGTTGCTGAAAAGCgattg             >  2:1860611/1‑36 (MQ=37)
                        ccAGTGGCTGGAAGAGGTGTTGCTGAAAAGCgattg             <  2:1368501/36‑1 (MQ=255)
                        ccAGTGGCTGGAAGAGGGGTTGCTGAAAAGAgtttg             >  2:2749353/1‑36 (MQ=21)
                          aGTGGCTGGGAGAGGGGTTGCTGAAAAGCggttttt           >  2:793660/1‑33 (MQ=16)
                           gTGGCTGGAAGAGGTGTTGCTGAAAAGCgattgatt          <  1:657293/36‑1 (MQ=255)
                            tgacTGGAAGAGGTGTTGCTGAAAAGCGATTGAttt         <  1:730888/33‑1 (MQ=255)
                             ggCTGGAAGAGGTGTTGCTGAAAAGCGATTGAtttt        >  2:2318268/1‑36 (MQ=255)
                             ggCTGGAAGAGGTGCTGCTGAAAAGCGATTGAtttt        <  2:3342455/36‑1 (MQ=255)
                              gCTAGAAGAGGTGTTGCTGAAAAGCGATTGATTTTg       <  2:1199125/36‑1 (MQ=255)
                               cTGGAAGAGGTGTTGCTGAAAAGCGATTGATTTTGt      >  1:2915747/1‑36 (MQ=255)
                               cTGGAAGAGGTGTTCCTGAAAAGCGATTGATTTTGt      >  1:2576541/1‑36 (MQ=255)
                                tGGAAGAGGTGTTGCTGAAAAGCGATTGATTTTGTa     <  2:1558036/36‑1 (MQ=255)
                                tGGAAGAGGTGTTGCTGAAAAGCGATTGATTTTGTa     <  1:945921/36‑1 (MQ=255)
                                 ggAAGAGGTGTTGCTGAAAAGCGATTGATTTTGTAg    <  1:3246316/36‑1 (MQ=255)
                                   aaGAGGTGTTGCTGAAAAGCGATTGATTTTGTAGGc  >  2:2552386/1‑36 (MQ=255)
                                   aaGAGGTGTTGCTGAAAAGCGATTGATTTTGTAGGc  <  2:580905/36‑1 (MQ=255)
                                   aaGAGGTGTTGCTGAAAAGCGATTGATTTTGTAGGc  <  1:847209/36‑1 (MQ=255)
                                   |                                   
GAAAGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGAAGAGGTGTTGCTGAAAAGCGATTGATTTTGTAGGC  >  REL606/4621239‑4621309

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 14 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG