Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606100,4080T→G29.0% 11.263518murDUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑ glutamatesynthetase
Reads supporting (aligned to +/- strand):  new base (18/0):  ref base (14/30):  total (33/30)

CTTCCAGCTGGAAACCACCTCCAGCTTACAGGCGGTAGCAGCGACCATTCTGAACGTGACTGAAGATCATA  >  REL606/100373‑100443
                                   |                                   
ctcccAGCTGGAAACCACCTCCAGCTTACAGGCGGt                                     <  2:2352064/33‑1 (MQ=255)
 ttCCAACAGAAAACCACCTCCAGCTTACAGGCGGTa                                    <  2:3110761/36‑1 (MQ=25)
  tCCAGCTGGAAACCACCTCCAGCTTACAGGCGGTag                                   <  2:2826609/36‑1 (MQ=255)
  tCCAGCTGGAAACCACCTCCAGCTTACAGGCGGGag                                   >  2:1934643/1‑36 (MQ=255)
  tCCAGCTGGAAACCACCTCCAGCTTACAGGCGGGag                                   >  1:1336595/1‑36 (MQ=255)
  tCCACCTGGAAACCACCTCCAGCTTACAGGCGGTag                                   <  1:1205046/36‑1 (MQ=255)
   ccAGCTGGAAACCACCTCCAGCTTACAGGCGGTagc                                  <  2:1988229/36‑1 (MQ=255)
    aggCTGGAAACCACCTCCAGCTTACAGGCGGTagca                                 <  1:67485/34‑1 (MQ=255)
     aGCTGGAAACCACCTCCAGCTTACAGGCGGTagcag                                <  1:53309/36‑1 (MQ=255)
     aGCTGGAAACCACCTCCAGCTTACAGGCGGTagcaa                                >  1:1265440/1‑35 (MQ=255)
     aGCTGGAAACCACCTCCAGCTTACAGGCGGGaggtg                                >  1:72389/1‑33 (MQ=37)
     aGCTGGAAACCACCTCCAGCTTACAGGCGGGagcag                                >  2:2900438/1‑36 (MQ=255)
     aGCTGGAAACCACCTCCAGCTTACAGGCGGGagcag                                >  2:2008842/1‑36 (MQ=255)
        tgaaaaCCACCCCCAGCTTACAGGCGGTAGCAGCGa                             <  1:3248658/33‑1 (MQ=37)
        tGGAAACCACCTCCAGCTTACAGGCGGGagcagggg                             >  2:2336697/1‑33 (MQ=37)
        tGGAAACCACCTCCAGCTCACAGGCGGTAACCGCGa                             >  2:2833690/1‑36 (MQ=25)
          gAAACCACCTCCAGCTTACAGGCGGTAGCAGCGAcc                           <  2:1332351/36‑1 (MQ=255)
            aaCCACCTCCAGCTTACAGGCGGGAGCAGCCAcccc                         >  1:1119271/1‑34 (MQ=37)
             accaccTCCAGCTTACAGGCGGTAGGAGCGAGCact                        >  1:258606/1‑34 (MQ=37)
             accaccTCCAGCTTACAGGCGGTAGCCGCCAcccct                        >  1:1837271/1‑33 (MQ=25)
             accaccTCCAGCTTACAGGCGGGAGCAGGGACCCtt                        >  2:3526841/1‑36 (MQ=37)
               ccccTCGAGCTTACAGGCGGTAGCAGCGACCATTCt                      <  2:909496/34‑1 (MQ=38)
               caccTCCAGCTTACAGGCGGTAGCAGCGACCAtttt                      >  2:2976777/1‑34 (MQ=255)
               caccTCCAGCTTACAGGCGGTAGCAGCGACCATTCt                      <  1:1283865/36‑1 (MQ=255)
               caccTCCAGCTTACAGGCGGGAGGAGCGGCCAtttt                      >  1:1528072/1‑34 (MQ=21)
                accTCCAGCTTACAGGCGGTAGCAGCGACCATTCTg                     <  1:2292728/36‑1 (MQ=255)
                accCCCAGCTTACAGGCGGTAGCAGCGACCATTCTg                     <  2:791839/36‑1 (MQ=255)
                 ccTCCAGCTTACAGGCGGTAGCAGCGACCATTCTGa                    <  1:835503/36‑1 (MQ=255)
                 ccTCCAGCTTACAGGCGGTAGCAGCGACCATTCTGa                    <  2:521761/36‑1 (MQ=255)
                 ccTCCAGCTTACAGGCGGGTGCCGGGACCGTTCTGa                    >  2:1201959/1‑36 (MQ=12)
                  cTCCAGCTTACAGGCGGTAGCAGCGACCATTCTGaa                   <  2:2819137/36‑1 (MQ=255)
                   tCCAGTTTACAGGCGGTAGCAGCGACCATTCTGAAc                  <  1:3023636/36‑1 (MQ=255)
                   tCCAGCTTACAGGCGGTGGCAGCGGGCCTTCTggaa                  >  2:3695705/1‑33 (MQ=12)
                   tCCAGCTTACAGGCGGTAGCAGCGACCATTCTGAAc                  <  2:3040978/36‑1 (MQ=255)
                   tCCAGCTTACAGGCGGTAGCAGCGACCATTCTGAAc                  <  2:1630121/36‑1 (MQ=255)
                   tCCAGCTTACAGGCGGTAGCAGCGACCATTCTGAAc                  <  2:3567478/36‑1 (MQ=255)
                    ccAGCTTACAGGCGGGAGCAACGACCCTACCGAACg                 >  2:1973606/1‑36 (MQ=16)
                     cAGCTTACAGGCGGTAGCAGCGACCATTCTGAACGt                <  1:569302/36‑1 (MQ=255)
                     cAGCTTACAGGCGGGAGCAGCGACCACTCCTAACGt                >  2:3307820/1‑36 (MQ=16)
                      aGCTTATGGGCGGTAGCAGCGACTATTCTGAACGTg               <  2:1190995/36‑1 (MQ=25)
                      aGCTTACAGGCGGTAGGAGGGACAATTCTTTACGTg               >  1:895617/1‑36 (MQ=16)
                      aGCTTACAGGCGGGAGCAGCGACCATTTTGAAAGTg               >  2:63929/1‑36 (MQ=21)
                      aGCGTACAGGCGGTAGCAGCGACCATTCTGAACGTg               <  2:3693215/36‑1 (MQ=255)
                       gCTTACAGGCGGTAGCAGCGACCATTCTGAACGTGa              <  1:1282017/36‑1 (MQ=255)
                       gCTTACAGGCGGTAGCAGCGACCATTCTGAACGTGa              <  1:553718/36‑1 (MQ=255)
                        cTTACAGGCGGGAGGAGGGGCCATTTTTGACGTggc             >  2:384854/1‑34 (MQ=2)
                        cTTACAGGCGGGAGGAGCGGCCCTTCTGAAAGCGAc             >  2:3070626/1‑36 (MQ=255)
                        cTTACAGGCGGGAGGAGCGGCCATTTTGGAGGTggc             >  2:1799649/1‑34 (MQ=255)
                        cTTACAGGCGGGAGGAGCGACCATTCTGTATGTGAc             >  2:1619341/1‑36 (MQ=19)
                         ttACAGGCGGGAGCAGCGACCATTTTTTACGTGACt            >  2:1145193/1‑36 (MQ=16)
                          tACAGGCGGTAGCAGCGACCATTCTGAACGTGACTg           <  1:318900/36‑1 (MQ=255)
                           aCAGGCGGTAGCAGCGACCATTCTGAACGTGACTGa          <  1:1943719/36‑1 (MQ=255)
                            cAGGCGGTAGCAGCGACCATTCTGAACGTGACTgcg         >  1:877988/1‑34 (MQ=255)
                            cAGGCGGTAGCAGCGACCATTCTGAACGTGACTGaa         >  1:2852622/1‑36 (MQ=255)
                            cAGGCGGTAGCAGCGACCATTCTGAACGTGACCGaa         >  2:2079052/1‑36 (MQ=255)
                             aGGCGGTAGCAGCGACCATTCTGAACGTGACTGAAg        >  1:2252901/1‑36 (MQ=255)
                             aGGCGGTAGCAGCGACCATTCTGAACGTGACTGAAg        <  2:1547110/36‑1 (MQ=255)
                             aGGCGATAGCAGCGACCATTCTGAACGTGACTGAAg        <  1:3025918/36‑1 (MQ=255)
                              ggCGGTAGCAGCGACCATTCTGAACGTGACTGAAGa       <  2:3294723/36‑1 (MQ=255)
                              ggCGGTAGCAGCGACCATTCTGAACGTGACTGAAGa       >  1:2714360/1‑36 (MQ=255)
                               gCGGTAGCAGCGACCATTCTGAACGTGACTGAAGAt      <  1:961450/36‑1 (MQ=255)
                                cGGTAGCAGCGACCATTCTGAACGTGACTGAAGATc     >  2:3497234/1‑36 (MQ=255)
                                cGGAAGCAGCGACCATTCTGAACGCGACTGCAGATc     >  2:1216999/1‑36 (MQ=21)
                                 ggTAGCAGCGACCATTCTGAACGTGACTGAAGATCa    >  2:3510223/1‑36 (MQ=255)
                                   tAGCAGCGACCATTCTGAACGTGACTGAAGATCata  <  1:538197/36‑1 (MQ=255)
                                   tAGCAGCGACCATTCTGAACGTGACTGAAGATCata  >  2:1734252/1‑36 (MQ=255)
                                   |                                   
CTTCCAGCTGGAAACCACCTCCAGCTTACAGGCGGTAGCAGCGACCATTCTGAACGTGACTGAAGATCATA  >  REL606/100373‑100443

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 10 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG