Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606165,6690T→G30.4% 11.646S238R (AGT→AGGhrpBpredicted ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  new base (14/0):  ref base (11/21):  total (25/21)

CGTGTGCAGGAACAACTGGCTTCGCGCATCGGCAGTGATGTGTTGCTCTGCCCGCTGTATGGCGCGTTGTC  >  REL606/165634‑165704
                                   |                                   
cGTGTGCAGGAACAACTGGCTTCGCGCATCGGCAGt                                     <  1:1242765/36‑1 (MQ=255)
cGTGTGCAGGAACAACTGGCTTCGCGCATCGGCAGt                                     <  1:7312237/36‑1 (MQ=255)
  tgtgCAGGAACAACTGGCTTCGCGCATCGGCAGTGa                                   <  1:1500146/36‑1 (MQ=255)
   gtgCAGGAACAACTGGCTTCGCGCATCGGCAGTGAt                                  <  1:7125094/36‑1 (MQ=255)
   gtgCAGGAACAACTGGCTTCGCGCATCGGCAGTGAt                                  <  1:6022317/36‑1 (MQ=255)
   gtgCAGGAACAACTGGCTTCGCGCATCGGCAGGGAt                                  >  1:871892/1‑36 (MQ=255)
     tcAGGAACAACTGGCTTCGCGCATCGGCAGTGAtgt                                <  1:2280991/35‑1 (MQ=255)
     gCAGGAACAACTGGCTTCGCGCATCGGCAGGGGtgt                                >  1:5325639/1‑36 (MQ=25)
     gCAGGAACAACTGGCTTCGCGCATCGGCAGGGAggt                                >  1:6632038/1‑36 (MQ=25)
      cAGGAACAACTGGCTTCGCGCATCGGCAGGGGtgtg                               >  1:3622577/1‑36 (MQ=25)
      cAGGAACAACTGGCTTCGCGCATCGGCAGGGAtgag                               >  1:2534605/1‑34 (MQ=25)
       aGGAACAACTGGCTTCGCGCATCGGCAGGGGtgtgg                              >  1:7253552/1‑35 (MQ=21)
        ggAACAACTGGCTTCGCGCATCGGCAGGGATTTGtt                             >  1:65806/1‑36 (MQ=21)
         gAACAACTGGCTTCGCGCATCGGCAGTGATGTGTTg                            <  1:5224101/36‑1 (MQ=255)
         gAACAACTGGCTTCGCGCATCGGCAGTGATGTGTTg                            <  1:1821325/36‑1 (MQ=255)
         gAACAACTGGCTNCGCGCATCGGCAGGGTTGTGttt                            >  1:3841099/1‑35 (MQ=21)
          aactacTGGCTTCGCGCATCGGCAGTGATGTGTTGc                           <  1:4170867/36‑1 (MQ=255)
          aacaacTGGCTTCGCGCATCGGCAGTGTTGTGTTgg                           >  1:5117823/1‑35 (MQ=255)
           acaacTGGCTGCGCGCATCGGCAGTGATGTGTTGct                          <  1:4432186/36‑1 (MQ=255)
            caacTGGCTTCGCGCATCGGCAGTGATGTGTTGctc                         <  1:3974019/36‑1 (MQ=255)
            caacTGGCTTCGCGCATCGGCAGTGAGGTGTTGctc                         <  1:2389826/36‑1 (MQ=255)
            caacTGGCTTCGCGCATCGGCAGGGGTGTGTTGttc                         >  1:1601908/1‑36 (MQ=21)
             aacTGGCTTCGCGCATCGGCAGGGGGGTGTTGctct                        >  1:4991936/1‑36 (MQ=18)
               cTGGCTTCGCGCATCGGCAGTGATGTGTTGCTCTGc                      <  1:2604970/36‑1 (MQ=255)
                tGGCTTCGCGCATCGGCAGGGGTGTGTTTCTCTTcc                     >  1:2180822/1‑36 (MQ=18)
                tGGCTTCGCGCATCGGCAGGGGTGTGTTGTTCTTcc                     >  1:3095861/1‑36 (MQ=18)
                    ttCGCGCATCGGCAGTGATGTGTTTGTCTGTCCCCt                 >  1:1109633/1‑36 (MQ=18)
                    ttCGCGCATCGGCAGTGATGTGTTGCTCTGCCCGCt                 <  1:2391284/36‑1 (MQ=255)
                    ttCGCGCATCGGCAGGGGTGTGTTGTTTTTCCCGCt                 >  1:2278333/1‑36 (MQ=12)
                    tgcgcgCATCGGCAGTGATGTGTTGCTCTGCCCGCt                 <  1:2505708/34‑1 (MQ=255)
                    tgcgcgCATCGGCAGTGATGTGTTGCTCTGCCCGCt                 <  1:236814/34‑1 (MQ=255)
                     tCGCGCATCGGCAGTGATGTGTTGCTCTTCCCCCtt                >  1:322727/1‑35 (MQ=25)
                     tCGCGCATCGGCAGTGATGTGTTGCTCTGCCCGccg                >  1:5404922/1‑34 (MQ=255)
                      cgcgCATCGGCAGTGATGTGTTGCTCTTCCCGGTGt               >  1:5934968/1‑36 (MQ=25)
                      cgcgCATCGGCAGTGATGTGTTGCTCTGCCCGCTGt               <  1:966567/36‑1 (MQ=255)
                      cgcgCATCGGCAGGGGTGTGTTGGTCTGCGCGATGt               >  1:5966743/1‑36 (MQ=12)
                       gcgcATCGGCAGTGATGTGTTGCTCTGCCCGCTGTa              <  1:4765900/36‑1 (MQ=255)
                       gcgcATCGGCAGTGATGAGATGCTTCGTCCGCTgaa              >  1:1396411/1‑34 (MQ=12)
                         gcATCGGCAGTGATGTGTTGCTCTGCCCGCTGTATg            <  1:6976788/36‑1 (MQ=255)
                           aTCGGCAGTGATGTGTTGCTCTGCCCGCTGTATGgc          <  1:7318811/36‑1 (MQ=255)
                           aTCGGCAGTGATGTGTTGCTCTGCCCGCTGTATGgc          <  1:5551470/36‑1 (MQ=255)
                             cGGCAGTGATGTGTTGCTCTGCCCGCTGTATGgcgc        <  1:6312896/36‑1 (MQ=255)
                              ggCAGTGATGTGTTGCTCTGCCCGCTGTATGGggcg       >  1:912301/1‑36 (MQ=255)
                              ggCAGTGATGTGTTGCTCTGCCCGCTGTATGGcgcg       >  1:5993459/1‑36 (MQ=255)
                               gCAGTGATGTGTTGCTCTGCCCGCTGTATGGCGCGt      >  1:4741666/1‑36 (MQ=255)
                                 cttGATGTGTTGCTCTGCCCGCTGTATGGCGCGTTg    <  1:6865131/34‑1 (MQ=255)
                                  gTGATGTGTTGCTCTGCCCGCTGTATGGCGCGTTGt   >  1:5485842/1‑36 (MQ=255)
                                  gTGATGTGTTGCTCTGCCCGCTGTATGGCGCGTTGt   <  1:1840364/36‑1 (MQ=255)
                                  gTGATGTGTTGCTCTGCCCGCTGTATGGCGCGTTGt   >  1:1557719/1‑36 (MQ=255)
                                   tGATGTGTTGCTCTGCCCGCTGTATGGCGCGTTGTc  >  1:6474821/1‑36 (MQ=255)
                                   tGATGTGTTGCTCTGCCCGCTGTATGGCGCGTTGTc  >  1:6765198/1‑36 (MQ=255)
                                   tGATGTGTTGCTCTGCCCGCTGTATGGCGCGTTGTc  <  1:4611110/36‑1 (MQ=255)
                                   tGATGTGTTGCTCTGCCCGCTGTATGGCGCGTTGTc  >  1:3302606/1‑36 (MQ=255)
                                   |                                   
CGTGTGCAGGAACAACTGGCTTCGCGCATCGGCAGTGATGTGTTGCTCTGCCCGCTGTATGGCGCGTTGTC  >  REL606/165634‑165704

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 12 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG