Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,509,9360A→G40.0% 18.940pseudogene (74/567 nt)yddEhypothetical protein; b1464_2
Reads supporting (aligned to +/- strand):  new base (16/0):  ref base (4/20):  total (20/20)

CTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT  >  REL606/1509902‑1509971
                                  |                                   
cTGCATTTACGCTTCGCTAAGATTATCGGCGGGGaa                                    >  1:2453006/1‑36 (MQ=255)
 ttcATTTGCGCTTCGCTAAGATTATCGGCGGGGaaa                                   <  1:3294102/34‑1 (MQ=255)
 tGCATTTGCGCTTCGCTAAGATTATCGGCGGGGgaa                                   >  2:3455433/1‑36 (MQ=255)
  tcATTTGCGCTTCGCTAAGATTAACGGCGGGGaaaa                                  <  2:527189/35‑1 (MQ=255)
  gtaTTTGCGCTTCGCTAAGATTATCGGCGGGGaaaa                                  <  1:2065483/34‑1 (MQ=255)
   cATTTGCGCTTCGCTAAGATTATCGGCGGGGGAaac                                 >  2:2104895/1‑36 (MQ=255)
   cATTTGCGCTTCGCTAAGATTATCGGCGGGGAAaac                                 <  2:3205120/36‑1 (MQ=255)
   cATTTGCGCTTCGCTAAGATTATCGGCGGGGAAaac                                 <  1:1605327/36‑1 (MQ=255)
    aTTTGCGCTTCGCTAAGATTATCGGCGGGGgaaaaa                                >  1:996206/1‑34 (MQ=38)
    aTTTGCGCTTCGCTAAGATTATCGGCGGGGAAaaca                                <  1:3051009/36‑1 (MQ=255)
      ttGCGCTTCGCTAAGATTATCGGCGGGGGAaacaac                              >  2:3563916/1‑36 (MQ=255)
      ttGCGCTTCGCTAAGATTATCGGCGGGGAAgacaaa                              >  2:2998400/1‑35 (MQ=255)
       tGCGCTTCGCTAAGATTATCGGCGGGGGAaaaaacc                             >  1:360254/1‑35 (MQ=37)
        gcgcTTCGCTAAGATTATCGGCGGGGGAAACAacac                            >  1:2406453/1‑36 (MQ=255)
        gcgcTTCGCTAAGATTATCGGCGGGGAAAACAacac                            <  2:3737413/36‑1 (MQ=255)
         cgcTTCGCTAAGATTATCGGCGGGGGAAaaaacaac                           >  1:2083137/1‑34 (MQ=37)
         cgcTTCGCTAAGATTATCGGCGGGGGAAACAAAAcc                           >  1:3505753/1‑36 (MQ=37)
          gcTTCGCTAAGATTATCGGCGGGGGGAAAAAAACCg                          >  2:2755736/1‑36 (MQ=19)
          gcTTCGCTAAGATTATCGGCGGGGGAAaccacaaca                          >  1:86573/1‑33 (MQ=21)
          gcTTCGCTAAGATTATCGCCGGGGAAAACAACAccc                          >  2:1397466/1‑35 (MQ=255)
           cTTCGCTAAGATTATCGGCGGGGGAAACCACAcccg                         >  2:3656795/1‑34 (MQ=21)
           cTTCGCGAAGATTATCGGCGGGGAAAACAACACCgg                         <  2:623860/36‑1 (MQ=255)
            ttCGCTAAGATTATCGGCGGGGGAAAAAAAACCgga                        >  2:2686278/1‑35 (MQ=25)
              cGCTAAGATTATCTGCGGGGAAAACAACACCGGCAg                      <  2:1695397/36‑1 (MQ=255)
              cGCTAAGATTATCGGCGGGGAAAACAACACCGGCAg                      <  1:1566393/36‑1 (MQ=255)
               gCTAAGATTATCGGCGGGGAAAACAACACCGGCAGa                     <  1:1024534/36‑1 (MQ=255)
                caaaGATTATCGGCGGGGAAAACAACACCGGCAGaa                    <  1:1701458/34‑1 (MQ=255)
                cTAAGATTATCGGCGGGGGAAACAAAACCGGCCGna                    >  1:647243/1‑36 (MQ=21)
                cTAAGATTATCGGCGGGGAAAACAACACCGGCAGaa                    <  1:3342677/36‑1 (MQ=255)
                 tAAGATTATCGGCGGGGGAAAAAACACCCGCCGAAt                   >  1:2503864/1‑36 (MQ=19)
                 tAAGATTATCGGCGGGGAAAACAACACCGGCAGAAt                   <  2:1062001/36‑1 (MQ=255)
                  aaGATTATCGGCGGGGAAAACAACACCGGCAGAAtt                  <  1:2542128/36‑1 (MQ=255)
                     aTTATCGGCGGGGGAAACAACAACGGCAGAATTGcc               >  1:444651/1‑36 (MQ=37)
                      ttATCGGCGGGGGAAACAACACCCGCCGAATTGGcg              >  1:3489974/1‑36 (MQ=16)
                       tATCGGCGGGGGAAACAAAAAACGCAGAAATGCggc             >  2:670820/1‑36 (MQ=0)
                        aacggcggGAAAAACAACACCGGCAGAATTGCcgcg            <  1:1540086/34‑1 (MQ=37)
                        aTCGGCGGGGGAAACAAAACCGGCAGCAATTCcgcg            >  1:1970788/1‑36 (MQ=12)
                          cggcggGGAAAACAACACCGGCAGAATTGCCGCGaa          <  1:3292568/36‑1 (MQ=255)
                           ggcggGGAAAACAACACCGGCAGAATTGCCGCGaaa         <  1:826468/36‑1 (MQ=255)
                            gcggGGAAAACAACACCGGCAGAATTGCCGCGAAAc        >  1:3230896/1‑36 (MQ=255)
                             cggGGAAAACAACACCGGCAGAATTGCCGCGAAACg       <  1:613008/36‑1 (MQ=255)
                             cggGGAAAACAACACCGGCAGAATTGCCGCGAAACg       <  1:1846038/36‑1 (MQ=255)
                             cggGGAAAACAACACCGGCAGAATTGCCGCGAAACg       >  2:1188822/1‑36 (MQ=255)
                               gggAAAACAACACCGGCAGAATTGCCGCGAAACGGt     <  1:1545925/36‑1 (MQ=255)
                               gggAAAACAACACCGACAGAATTGCCGCGAAACGGt     <  2:1416460/36‑1 (MQ=255)
                                gaaaaaCAACACCGGCAGAATTGCCGCGAAACGGtt    <  1:670799/34‑1 (MQ=255)
                                 gAAAACAACACCGGCAGAATTGCCGCGAAACGGTtg   <  1:426748/36‑1 (MQ=255)
                                  aaaaCAACACCGGCAGAATTGCCGCGAAACGGTtgt  <  2:3431378/36‑1 (MQ=255)
                                  aaaaCAACACCGGCAGAATTGCCGCGAAACGGTtgt  >  1:1353391/1‑36 (MQ=255)
                                  |                                   
CTGCATTTGCGCTTCGCTAAGATTATCGGCGGGGAAAACAACACCGGCAGAATTGCCGCGAAACGGTTGT  >  REL606/1509902‑1509971

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