Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,538,2350G→T25.9% 11.9541437doscAMP phosphodiesterase, heme‑regulated
Reads supporting (aligned to +/- strand):  new base (14/0):  ref base (10/30):  total (24/30)

AGCGGGAAGGGTTTATCGTCAATCATTATCGGCTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGG  >  REL606/1538200‑1538268
                                   |                                 
aGCGGGAAGGGTTTATCGTCAATCATTATCGGCTTg                                   <  2:2328482/36‑1 (MQ=255)
 gCGGGAGGGGTTTATCGTCAATCATTATCGGCTTGc                                  <  1:731074/36‑1 (MQ=255)
 gCGGGAAGGGTTTATCGTCAATCATTATCGGCTTGc                                  <  1:3642831/36‑1 (MQ=255)
  ggggTAGGGTTTATCGTCAATCATTATCGGCTTGCt                                 <  1:730464/35‑1 (MQ=255)
  cGGGAAGGGTTTATCGTCAATCATTATCGGCTTTCt                                 >  2:1206043/1‑36 (MQ=255)
   gggAAGGGTTTATCGTCAATCATTATCGGCTTGCTg                                <  1:2412789/36‑1 (MQ=255)
   gggAAGGGTTTATCGTCAATCATTATCGGCTTGCTg                                <  1:2818927/36‑1 (MQ=255)
    ggAAGGGTTTATCGTCAATCATTATCGGCTTTCtct                               >  1:70235/1‑34 (MQ=38)
    ggAAGGGTTTATCGTCAATCATTATCGGCTTTCTGa                               >  2:248744/1‑36 (MQ=255)
    gaaaGGGTTTATCGTCAATCATTATCGGCTTGCTGa                               <  1:1646264/34‑1 (MQ=255)
     gAAGGGTTTATCGTCAATCATTATCGGCTTGCTGAc                              <  1:902715/36‑1 (MQ=255)
       aGGGTTTATCGTCAATCATTATCGGCTTGCTGACca                            <  2:1549501/36‑1 (MQ=255)
         ggTTTATCGTCAATCATTATCGGCTTTTTGATcaca                          >  1:2840614/1‑36 (MQ=25)
         ggTTTATCGTCAATCATTATCGGCTTTCTGACccca                          >  2:2512473/1‑36 (MQ=255)
         ggTTTATCGTCAATCATTATCGGCTTGCTGACcaca                          <  1:3076765/36‑1 (MQ=255)
           tttATCGTCAATCATTATCGGCTTGCTGACCACAtt                        <  1:3422764/36‑1 (MQ=255)
           tttATCGTCAATCATTATCGGCTTGCTGACCACAtt                        <  2:2249757/36‑1 (MQ=255)
             tATCGTCAATCATTATCGGCTTTCTTTCCACCTTcc                      >  1:365882/1‑36 (MQ=19)
             tATCGTCAATCATTATCGGCTTTCTGACCCCAATcc                      >  2:2127042/1‑36 (MQ=255)
             tATCGTCAATCATTATCGGCTTGCTGACCACATTcc                      <  1:9097/36‑1 (MQ=255)
              aTCGTCAATCATTATCGGCTTTCTGACCCCATTccc                     >  1:63756/1‑35 (MQ=255)
               tCGTCAATCATTATCGGCTTGCTGACCACATTCCGt                    <  2:1067156/36‑1 (MQ=255)
                cGTCAATCATTATCGGCTTGCTGACCACATTCCGTa                   <  2:1469849/36‑1 (MQ=255)
                cGTCAATCATTATCGGCGTGCTGACCACATTCCGTa                   <  2:978575/36‑1 (MQ=255)
                 gTCAATCATTATCGGTTTGCTGACCACATTCCGTAg                  <  2:2948096/36‑1 (MQ=255)
                 gTCAATCATTATCGGCTTGTTGACCACATTCTGttg                  >  1:834096/1‑34 (MQ=25)
                 gTCAATCATTATCGGCTTGTTGACCACATTCCGCAg                  >  1:2935180/1‑36 (MQ=37)
                 gTCAATCATTATCGGCTTGCTGACCACATTCCGTAg                  <  1:2305804/36‑1 (MQ=255)
                  tCAATCATTATCGGCTTTCTGACCCTATTCCCTAGc                 >  1:3583015/1‑36 (MQ=255)
                   cAATCATTATCGGCTTGCTGACCACATTCCGTAGCt                <  1:2922203/36‑1 (MQ=255)
                    aaTCATTATCGGCTTTCTTATCACATTCTGTAGCTc               >  2:1333550/1‑36 (MQ=19)
                     aTCATTATCGGCTTGCTTACCACATTCCGTAGCTCa              >  2:1318431/1‑36 (MQ=255)
                      tCATTATCGGCTTGCTGTCCACATTCCGTAGCTCAt             <  1:2227753/36‑1 (MQ=255)
                      tCAATATCGGCTTGCTGACCACATTCCGTAGCTCAt             <  2:797581/36‑1 (MQ=255)
                       cATTATCGGCTTTCTTATCACATTCCGTAGGTCATc            >  2:547158/1‑36 (MQ=19)
                         ttATCGGCTTTCTGACCACATTCCGTAGCTCATCgg          >  1:2184785/1‑36 (MQ=255)
                         ttATCGGCTTGCTGACCACATTCCGTAGCTCATCgg          <  1:2853847/36‑1 (MQ=255)
                          tATCGGCTTTCTGACCACATTCCGTAGCTCATTGGc         >  2:602520/1‑36 (MQ=38)
                          tATCGGCTTTCTGACCACATACCGTAGCCCCTCGGc         >  2:2342191/1‑36 (MQ=19)
                          tATCGGCTTGCTGACCACATTCCGTAGCTCATCGGc         <  1:894769/36‑1 (MQ=255)
                             cGGCTTGCTGACCACATTCCGTAGCTCATCGGCGAt      <  1:1193997/36‑1 (MQ=255)
                             cGGCTTGCTGACCACATTCCGTAGCTCATCGGCGAt      >  1:1011376/1‑36 (MQ=255)
                              ggCTTGCTGACCATATTCCGTAGCTCATCGGCGAtt     >  1:2130282/1‑36 (MQ=255)
                              ggCTTGCTGACCACCTTCCGTAGCTCATCGGCGTtt     >  1:3643039/1‑36 (MQ=37)
                              ggCTTGCTGACCACATTCCGTAGCTCATCGGCGAtt     >  1:3443507/1‑36 (MQ=255)
                              ggCTTGCTGACCACATTCCGTAGCTCATCGGCGAtt     <  2:930743/36‑1 (MQ=255)
                              ggCTTGCTGACCACATTCCGTAGCTCATCCGCGAtt     >  1:1282298/1‑36 (MQ=255)
                              ggCTTGCTGACCACATTCCGTAGATCATCGGCGAtt     >  2:3543912/1‑36 (MQ=255)
                               gtttGCTGACCACATTCCGTAGCTCATCGGCGAttt    <  1:109010/34‑1 (MQ=255)
                               gCTTGCTGACCACATTCCGTAGCTCATCGGCGAttt    <  1:2875581/36‑1 (MQ=255)
                               gCTTGCTGACCACATTCCGTAGCTCATCGGCGAttt    <  2:549240/36‑1 (MQ=255)
                               gCTGGCCGACCAGATTCCGTAGCTCATCGGCGAgtt    <  2:555681/36‑1 (MQ=16)
                                cTTGCTGACCACATTCCGTAGCTCATCGGCGATTTg   <  1:3533839/36‑1 (MQ=255)
                                 ttGCTGACCACATTCCGTAGCTCATCGGCGAttttg  >  2:1819366/1‑34 (MQ=255)
                                 ttGCTGACCACATTCCGTAGCTCATCGGCGATTTgg  <  1:2575807/36‑1 (MQ=255)
                                   |                                 
AGCGGGAAGGGTTTATCGTCAATCATTATCGGCTTGCTGACCACATTCCGTAGCTCATCGGCGATTTGG  >  REL606/1538200‑1538268

Base quality scores:  ATCG  < 3 ≤  ATCG  < 6 ≤  ATCG  < 14 ≤  ATCG  < 27 ≤  ATCG  < 32 ≤  ATCG