Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606506,2070T→G25.0% 11.868126gclglyoxylate carboligase
Reads supporting (aligned to +/- strand):  new base (17/0):  ref base (14/37):  total (31/37)

AATCCGTTCTACTCAGCGATGCGTAAGCACGGCGGTATTCGTCACATTCTGGCGCGTCATGTGGAAGGTG  >  REL606/506172‑506241
                                   |                                  
aatcCGTTCTACTCAGCGATGCGTACGCACGGCGGt                                    >  1:1354155/1‑36 (MQ=255)
aatcCGTTCTACTCAGCGATGCGTAAGCACGGCGGt                                    <  1:3021040/36‑1 (MQ=255)
aatcCGTTCTACTCAGCGATGCGTAAGCACGGCGGt                                    <  1:2641269/36‑1 (MQ=255)
 atcCGTTCTACTGAGCGATGCGTAAGCACGGCGGTa                                   <  1:420859/36‑1 (MQ=255)
 atcCGTTCTACTCAGCGATGCGTAAGCACGGTGGTa                                   >  1:838611/1‑36 (MQ=255)
 atcCGTTCTACTCAGCGATGCGTAAGCACGGCGGTa                                   <  2:1569511/36‑1 (MQ=255)
  tcCGTTCTACTCAGCGATGCGTAAGCACGGCGGTAt                                  <  2:1409406/36‑1 (MQ=255)
    cGTTCTACTCAGCGATGCGTAAGCACGGCGGTGttt                                >  2:1388178/1‑35 (MQ=255)
     gTTCTACTCAGCGATGCGTAAGCACGGCGGTAtttg                               >  1:1757920/1‑34 (MQ=255)
       tCTACTCAGCGATGCGTAAGCACGGCGGTATTCGTc                             <  1:236940/36‑1 (MQ=255)
       tCTACTCAGCGATGCGTAAGCACGGCAGGATACGTc                             >  1:2056787/1‑36 (MQ=37)
        cTACTCAGCGATGCGTAAGCACTGCGGTATTCGTca                            <  2:2910341/36‑1 (MQ=255)
        cTACTCAGCGATGCGTAAGCACGGCGGTATTCGTca                            <  2:1724493/36‑1 (MQ=255)
        cTACTCAGCGATGCGTAAGCACGGCGGTATTCGTca                            <  1:109156/36‑1 (MQ=255)
        cTACTCAGCGATGCGTAAGCACGGCGGGATTCGtga                            >  2:545124/1‑34 (MQ=38)
         tACTCAGCGATGCGTACGCACGGCGGGGTTCGTcac                           >  1:3304842/1‑36 (MQ=37)
         tACTCAGCGATGCGTAAGCACGGCGGTATTCGTcac                           <  1:3342863/36‑1 (MQ=255)
         tACTCAGCGATGCGTAAGCACGGCGGGATTTGttcc                           >  2:1429694/1‑33 (MQ=37)
          aCTCAGGGATGCGTAAGCACGACGGTATTCGTcaca                          <  1:540522/36‑1 (MQ=38)
          aCTCAGCGATGCGTAAGCACGGCGGTATTCGTcaca                          <  2:2260514/36‑1 (MQ=255)
          aCTCAGCGATGCGTAAGCACGGCGGTATTCGTcaca                          <  2:2102648/36‑1 (MQ=255)
          aCTCAGCGATGCGTAAGCACGGCGGGATTTGTcatt                          >  1:1199310/1‑34 (MQ=37)
          aCTCAGCGATGCGTAAGCACGGCGGGATTCGTcaca                          >  1:2647487/1‑36 (MQ=255)
          aCTCAGCGATGCGTAAGCACGGCGGGAATTGTcaga                          >  1:3505704/1‑34 (MQ=25)
           cTCAGCGATGCGTAAGCACGGCGGTGTTCGTCACAt                         >  2:1229109/1‑36 (MQ=255)
           cTCAGCGATGCGTAAGCACGGCGGGATTTGTcacct                         >  2:2180551/1‑34 (MQ=25)
            tCAGCGATGCGTAAGCACGGCGGTATTCGTCACAtt                        <  2:3455775/36‑1 (MQ=255)
            tCAGCGATGCGTAAGCACGGCGGTATTCGTCACAtt                        <  1:177796/36‑1 (MQ=255)
            tCAGCGATGCGTAAGCACGGCGGGCTTTGTCACAtt                        >  1:2664863/1‑36 (MQ=25)
            tCAGCGATGCGTAAGCACGGCGGGATTCGTCACAtt                        >  1:213414/1‑36 (MQ=255)
             cAGCGATGCGTAAGCACGGCGGGATTCGTCACATTc                       >  1:3647283/1‑36 (MQ=255)
              aGCGATGCGTAAGCACGGCGGTATTCGTCACATTCt                      <  2:1358057/36‑1 (MQ=255)
              aGCGATGCGTAAGCACGGCGGGAATCGTCACATTcc                      >  1:617882/1‑35 (MQ=37)
               gCGATGCGTAAGCACGGCGGTATTCGTCACATTCTg                     <  1:1567088/36‑1 (MQ=255)
                tgATGCGTAAGCAAGGCGGTTTTCGTTACATTTTgt                    >  2:89447/2‑35 (MQ=16)
                 gATGCGTAAGCACGGCGGTATTCGTTACATTCttcg                   >  1:2010497/1‑33 (MQ=37)
                 gATGCGTAAGCACGGCGGTATTCGTCTCATTCTTgc                   >  1:1201395/1‑36 (MQ=37)
                 gATGCGTAAGCACGGCGGTATTCGTCACATTCTGgc                   <  2:3663459/36‑1 (MQ=255)
                 gATGCGTAAGCACGGCGGTATTCGTCACATTCTGgc                   <  1:284630/36‑1 (MQ=255)
                 gATGCGTAAGCACGGCGGGATTTGTTACATTTTggg                   >  1:1462319/1‑35 (MQ=19)
                 gATGCGTAAGCACGGCGGGATTCGTCACGTTCTGgc                   >  2:1942554/1‑36 (MQ=37)
                 gATGCGAAAACACGGCGGTATTCGTCACATTCTGgc                   <  1:1722202/36‑1 (MQ=37)
                 aaTGCGTAAGCACGGCGGTATTCGTCACATTCTGgc                   <  1:1091717/35‑1 (MQ=255)
                   tGCGTAAGCACGGCGGTATTCGTCACATTCTGgcgc                 <  1:870844/36‑1 (MQ=255)
                   tGCGTAAGCACGGCGGTATTCGTCACATTCTGgcgc                 <  2:2655863/36‑1 (MQ=255)
                    gCGTAAGCACGGCGGTATTCGTCACATTCTGGcgcg                <  1:159598/36‑1 (MQ=255)
                    gCGTAAGCACGGCGGTATTCGTCACATTCTGGcgcg                >  2:456804/1‑36 (MQ=255)
                    gCGTAAGCACGGCGGTAATCGTCACATTTTTGcgcg                >  1:810599/1‑36 (MQ=25)
                     cGTAAGCACGGCGGTATTTCTCACAATCTCgcgggc               >  1:504585/1‑33 (MQ=12)
                     cGTAACCACGGCGGTATTCGTCACATTCTGGCGCGt               <  2:273001/36‑1 (MQ=255)
                      gTAAGCACGGCGGTATTCGTCACATTCTGGCGCGTc              <  2:2934440/36‑1 (MQ=255)
                      gTAAGCACGGCGGGATTTCTTACATTCTCGggcgcg              >  1:2368626/1‑34 (MQ=9)
                       tAAGCACGGCGGGATTCGTCACATTCTGGCGCGTCa             >  1:3183227/1‑36 (MQ=255)
                         aGCACGGCGGTATTCGTCACATTCTGGCGCGTCAtg           <  1:1936244/36‑1 (MQ=255)
                          gCACGGCGGTATTCGTCACATTCTGGCGCGTCAtgt          <  1:1137716/36‑1 (MQ=255)
                          gCACGGCGGTATTCGTCACATTCTGGCGCGTCAtgt          <  1:1882190/36‑1 (MQ=255)
                          gCACGGCGGTATTCGTCACATTCTGGCGCGTCAtgt          <  2:2353762/36‑1 (MQ=255)
                          gCACGGCGGTATTCGTCACATTCTGGCGCGTCAtgt          >  2:2037342/1‑36 (MQ=255)
                          gCACGGCGGTATTCGTCACATTCTGGCGCGTCAtgt          <  2:1232537/36‑1 (MQ=255)
                           cACGGCGGTATTCGTCACATTCTGGCGCGTCAtgtg         <  1:3734159/36‑1 (MQ=255)
                           cACGGCGGGATTTGTTCCATTCTGGCGCGTCCtgga         >  2:2504774/1‑34 (MQ=12)
                            aCGGCGGTATTCGTCACATTCTGGCGCGTCATGTgg        <  1:1388223/36‑1 (MQ=255)
                             cggcggTATTCGTCACATTCTGGCGCGTCATGTGGa       <  2:1108562/36‑1 (MQ=255)
                              ggcggTATTCGTCACATTCTGGCGCGTCCTGTGGaa      >  2:1950708/1‑36 (MQ=255)
                               gcggTATTCGTCACATTCTGGCGCGTCATGTGGAAg     <  1:1058911/36‑1 (MQ=255)
                                cggTATTCGTCACATTCTGGCGCGTCATGTGGAAgg    <  2:73361/36‑1 (MQ=255)
                                 ggTATTCGTCACATTCTGGCGCGTCTTGTGGAAGGt   >  1:1058065/1‑36 (MQ=255)
                                 ggTATTCGTCACATTCTGGCGCGTCATGTGGAAGGt   <  2:365099/36‑1 (MQ=255)
                                 ggTATTCGTCACATTCCGGCGCGTCATGTGGAAgct   >  1:2876398/1‑34 (MQ=37)
                                  gTATTTGTCACATTCTGGCGCGTCATGTGGAAGGTg  <  2:2899439/36‑1 (MQ=255)
                                  gTATTCGTCACATTCTGGCGCGTCATGTGGAAGGTg  <  2:11536/36‑1 (MQ=255)
                                   |                                  
AATCCGTTCTACTCAGCGATGCGTAAGCACGGCGGTATTCGTCACATTCTGGCGCGTCATGTGGAAGGTG  >  REL606/506172‑506241

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