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

GAAAGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGAAGAGGTGTTGCTGAAAAGCGATTGATTTTGTAGGC  >  REL606/4621239‑4621309
                                   |                                   
gAAAGCCGCGCCGGTTGAGCTGCGCCAGTGGCTGGa                                     <  1:622302/36‑1 (MQ=255)
gAAAGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGa                                     <  1:177102/36‑1 (MQ=255)
gAAAGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGa                                     <  1:7168534/36‑1 (MQ=255)
 aaaGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGaa                                    <  1:4514129/36‑1 (MQ=255)
 aaaGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGaa                                    <  1:5072491/36‑1 (MQ=255)
  aaGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGGag                                   >  1:4903321/1‑36 (MQ=255)
  aaGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGGag                                   >  1:5795959/1‑36 (MQ=255)
  aaGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGAag                                   <  1:1352542/36‑1 (MQ=255)
  aaGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGAag                                   <  1:4403714/36‑1 (MQ=255)
   aGCCGCGCCCGTTGAGCTGCGCCAGTGGCTGGGagg                                  >  1:6639241/1‑35 (MQ=255)
     ccGCGCGCGTTGGGCTGCGCCAGTGGCTGGAAGAgg                                <  1:503535/36‑1 (MQ=25)
      cgctCCCGTTGAGCTGCGCCAGTGGCTGGAAGAGgt                               <  1:6205765/36‑1 (MQ=255)
      cgcgCCCGTTGAGCTGCGCCAGTGGCTGGGagaagg                               >  1:3345618/1‑33 (MQ=21)
      cgcgCCCGTTGAGCTGCGCCAGTGGCTGGAAGAGgt                               <  1:1905153/36‑1 (MQ=255)
      cgcgCCCGTTGAGCTGCGCCAGTGGCTGGAAGAGgt                               <  1:5878813/36‑1 (MQ=255)
       gcgcCCGTTGAGCTGCGCCAGTGGCTGGAAGAGgtg                              <  1:4128773/36‑1 (MQ=255)
       gcgcCCGTTGAGCTGCGCCAGTGGCTGGAAGAGgtg                              <  1:5666056/36‑1 (MQ=255)
       gcgcCCGTTGAGCTGCGCCAGTGGCTGGAAGAGgtg                              <  1:6540666/36‑1 (MQ=255)
         gcCCGTTGAGCTGCGCCAGTGGCTGGAAGAGGTGtt                            <  1:5162231/36‑1 (MQ=255)
         gcCCGTTGAGCTGCGCCAGTGGCTGGAAAAGGCGtt                            >  1:6000570/1‑36 (MQ=21)
           ccGTTGAGCTGCGCAAGTGGCTGGGAGAGgtgtggg                          >  1:2926453/1‑33 (MQ=21)
             gTTGAGGGGCGCCAGTGGCTGGAAGAGGTGTTGCTg                        <  1:4564237/36‑1 (MQ=21)
             gTTGAGCTGCGCCAGTGGCTGGGGGAGGTGTTGTTg                        >  1:1802749/1‑36 (MQ=21)
               tGAGCTGCGCCAGTGGCTGGAAGAGGTGTTGCTGaa                      <  1:2627413/36‑1 (MQ=255)
               tGAGCTGCGCCAGTGGCTGGAAGAGGGGTTGGTGaa                      >  1:5845442/1‑36 (MQ=25)
                gAGTTGCGCCAGTGGCTGGAAGAGGTGTTGCTGaaa                     <  1:4089131/36‑1 (MQ=255)
                gAGCTGCGCCAGTGGCTGGGAGAGGGGTTGGTGaga                     >  1:1131575/1‑34 (MQ=18)
                gAGCTGCGCCAGTGGCTGGAAGAGGTGTTGCTGaaa                     <  1:5699952/36‑1 (MQ=255)
                gAGCTGCGCCAGTGGCTGGAAGAGGTGTTGCTGaaa                     <  1:3467783/36‑1 (MQ=255)
                 aGCTGCGCCAGTGGCTGGAAGAGGTGTTGCTGaaaa                    <  1:3658115/36‑1 (MQ=255)
                 aGCTACGCCAGTGGCTGGGAGAGGGGCTGCTGaaaa                    >  1:2208348/1‑36 (MQ=18)
                     gcgcCAGTGGCTGGGAGGGGGGGTGCTGGAAAGCgg                >  1:4126499/1‑35 (MQ=12)
                     gcgcCAGTGGCTGGAAGAGGTGTTGCTGAAAAGCga                <  1:1539456/36‑1 (MQ=255)
                     gcgcCAGTGGCTGGAAGAGGTGTTGCTGAAAAGCga                <  1:6847750/36‑1 (MQ=255)
                     gcgcCAGTGGCTGGAAGAGGTGTTGCTGAAAAGCga                <  1:2106405/36‑1 (MQ=255)
                      cgcCAGTGGCTGGGAGAGGGGTTGCTGGAAAGCgat               >  1:5774481/1‑36 (MQ=21)
                      cgcCAGTGGCTGGAAGAGGGGTTGGTTGAAAGCggt               >  1:4923586/1‑34 (MQ=17)
                       gcCAGTGGCTGGGAGAGGTGTTGCTGGAAAGAgatt              >  1:4193169/1‑36 (MQ=21)
                       gcCAGTGGCTGGAAGGGGTGTTGCTGAAAAGCgatt              <  1:7290383/36‑1 (MQ=255)
                        cgaGTGGCTGGAAGAGGTGTTGCTGAAAAGCgattg             <  1:2480069/34‑1 (MQ=255)
                        ccAGTGGCTGGGAGAGGGGTTGCTGAAAAGCgattg             >  1:5599821/1‑36 (MQ=25)
                        ccAGTGGCTGGAAGAGGTGTTGCTGAAAAGCgattg             <  1:5107711/36‑1 (MQ=255)
                        ccAGTGGCTGGAAGAGGGGTTGCTGAAAAGAgtttg             >  1:6488563/1‑36 (MQ=21)
                          aGTGGCTGGGAGAGGGGTTGCTGAAAAGCggttttt           >  1:4532870/1‑33 (MQ=17)
                           gTGGCTGGAAGAGGTGTTGCTGAAAAGCgattgatt          <  1:657293/36‑1 (MQ=255)
                            tgacTGGAAGAGGTGTTGCTGAAAAGCGATTGAttt         <  1:730888/33‑1 (MQ=255)
                             ggCTGGAAGAGGTGTTGCTGAAAAGCGATTGAtttt        >  1:6057478/1‑36 (MQ=255)
                             ggCTGGAAGAGGTGCTGCTGAAAAGCGATTGAtttt        <  1:7081665/36‑1 (MQ=255)
                              gCTAGAAGAGGTGTTGCTGAAAAGCGATTGATTTTg       <  1:4938335/36‑1 (MQ=255)
                               cTGGAAGAGGTGTTGCTGAAAAGCGATTGATTTTGt      >  1:2915747/1‑36 (MQ=255)
                               cTGGAAGAGGTGTTCCTGAAAAGCGATTGATTTTGt      >  1:2576541/1‑36 (MQ=255)
                                tGGAAGAGGTGTTGCTGAAAAGCGATTGATTTTGTa     <  1:5297246/36‑1 (MQ=255)
                                tGGAAGAGGTGTTGCTGAAAAGCGATTGATTTTGTa     <  1:945921/36‑1 (MQ=255)
                                 ggAAGAGGTGTTGCTGAAAAGCGATTGATTTTGTAg    <  1:3246316/36‑1 (MQ=255)
                                   aaGAGGTGTTGCTGAAAAGCGATTGATTTTGTAGGc  >  1:6291596/1‑36 (MQ=255)
                                   aaGAGGTGTTGCTGAAAAGCGATTGATTTTGTAGGc  <  1:4320115/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