Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6062,197,8670T→G28.3% 11.460141cirAferric iron‑catecholate outer membrane transporter
Reads supporting (aligned to +/- strand):  new base (17/0):  ref base (10/33):  total (27/33)

AGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAG.ATTTTGTTCCACGGAAGATGCAGT  >  REL606/2197832‑2197902
                                   |                                    
aGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGt                                      <  2:601341/36‑1 (MQ=255)
aGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGt                                      <  2:3735035/36‑1 (MQ=255)
 ggTCTTCCTGGGTAATGACGCTGATACTGGCAGGTg                                     <  2:214640/36‑1 (MQ=255)
  gTCTACCTAGGTAATGACGCTGCTACTGGCAGGTGc                                    <  2:787340/36‑1 (MQ=25)
   tCTTCCTGGGTAATGGCGCTGATACTGGCAGGTGCa                                   <  2:2149729/36‑1 (MQ=255)
    cTTGCTGGGTAATGACGCTGATACTGGCAGGTGCAt                                  <  2:2029783/36‑1 (MQ=255)
    cTTCCAGGGTAATGACGCTGATACTGGCAGGTGCAt                                  <  2:455104/36‑1 (MQ=255)
     ttCCTGGGTAATGACGCTGATACTGGCAGGGGcgtt                                 >  1:1254018/1‑33 (MQ=37)
      tCCTGGGTAATGACGCTGATACTGGCAGGTGCATCt                                <  1:257389/36‑1 (MQ=255)
      tCCTGGGTAATGACGCTGATACTGGCAGGTGCATCt                                <  1:3315916/36‑1 (MQ=255)
       ccTGGGTAATGACGCTGATACTGGCAGGTGCATCtt                               <  2:3397929/36‑1 (MQ=255)
        cTGGGTAATGACGCTGATACTGGCAGGTGCATCttt                              <  2:1967803/36‑1 (MQ=255)
        cTGGGTAATGACGCTGATACTGGCAGGGGGCTCttt                              >  2:832707/1‑36 (MQ=37)
         tGGGTAATGACGCTGATACTGGCAGGTGCATCTTTa                             <  1:2260218/36‑1 (MQ=255)
         tGGGTAATGACGCTGATACTGGCAGGTGCATCTTTa                             <  1:2664968/36‑1 (MQ=255)
          gggTAATGACGCTGATACTGGCAGGTGCATCTTTaa                            <  1:777899/36‑1 (MQ=255)
          gggTAATGACGCTGATACTGGCAGGTGCATCTTTaa                            <  1:2764934/36‑1 (MQ=255)
          gagTAATGACGCTGATACTGGCAGGTGCATCTTTaa                            <  1:1321948/34‑1 (MQ=255)
            gTAATGACGCTGATACTGGCAGGTGCATCTTTAAG.a                         <  1:3495045/36‑1 (MQ=255)
            gTAATGACGCTGATACTGGCAGGGGGATTTTTatat                          >  1:2365153/1‑33 (MQ=19)
             tAATGACTCTGATACTGGCAGGTGCATCTTTAAG.At                        <  1:3101331/36‑1 (MQ=255)
             tAATGACGCTGATACTGGCAGGTGCATCTTTAAG.At                        <  2:991245/36‑1 (MQ=255)
             tAATGACGCTGATACTGGCAGGGGCATCTTTATG.At                        >  2:1850666/1‑36 (MQ=37)
              aaTGACGCTGATACTGGCAGGTGCATCTTTAAG.Att                       <  1:39520/36‑1 (MQ=255)
              aaTGACGCTGATACTGGCAGGGGTCTCTTTTAT.Att                       >  1:184003/1‑36 (MQ=16)
               aTGACGCTGATACTGGCAGGTGCATCTTTAAG.Attt                      <  1:3128124/36‑1 (MQ=255)
               aTGACGCTGATACTGGCAGGTGCATCTTTAAG.Attt                      <  2:1995883/36‑1 (MQ=255)
                  aCGCTGATACTGGCCGGGGGATATTTTAGTAttttt                    >  1:1881682/1‑35 (MQ=0)
                  aCGCTGATACTGGCAGGTGGATCTTTTAG.Atttttt                   >  1:2627536/1‑34 (MQ=37)
                  aCGCTGATACTGGCAGGGGCCTCTTTTAT.AGTTTGt                   >  1:2484379/1‑36 (MQ=16)
                  aCGCTGATACTGGCAGGGGCATCTTTTAG.Atttttt                   >  1:1724734/1‑34 (MQ=25)
                  aCGCTGATACTGGCAGGGGCATCTTTAAA.ATTTTGt                   >  2:1752854/1‑36 (MQ=37)
                   cGCTGATACTGGCAGGGGGATCTTTAAG.TTTTTTtt                  >  1:2816388/1‑36 (MQ=19)
                    gCTGATACTGGCAGGTGCATCTTTAAG.ATTTTGTTc                 <  2:1147469/36‑1 (MQ=255)
                    gCTGATACTGGCAGGGGGGTCTTTTAG.ATTTTTttt                 >  2:249678/1‑35 (MQ=12)
                    gCTGATACTGGCAGGGGCATCCTCAAG.AATTTGttt                 >  1:3207773/1‑35 (MQ=16)
                     cTGATACTGGCAGGGGGATCTTTTAG.ATTTTTtttt                >  2:1039211/1‑34 (MQ=16)
                      tGATACTGGCAGGTGGTTTTTTTAG.ATTTTTtttct               >  1:2111572/1‑33 (MQ=12)
                      tGATACTGGCAGGTGCATCTTTAAG.ATTTTGTTCCa               <  1:1361662/36‑1 (MQ=255)
                       gATACTGGCAGGTGCATCTTTAAG.ATTTTGTTCCAc              <  1:1745883/36‑1 (MQ=255)
                       gATACTGGCAGGTGCATCTTTAAG.ATTTTGTTCCAc              <  2:2932898/36‑1 (MQ=255)
                       gATACTGGCAGGTGCATCTTTAAG.ATTTTGTTCCAc              <  1:2811315/36‑1 (MQ=255)
                        aTACTGGCAGGGGGAACTTTAAG.AATTTGTTTCACg             >  1:3385314/1‑36 (MQ=12)
                         tACTGGCAGGGGCCTCTTTTAG.ATTTTTTTTCAcag            >  2:126438/1‑34 (MQ=12)
                          aCTGGCAGGGGGATCTTTAAG.ATTTTTTTCCACGGa           >  2:3099899/1‑36 (MQ=37)
                           cTGGCAGGTGCATCTTTAAG.ATTTTGTTCCACGGaa          <  1:1291455/36‑1 (MQ=255)
                           cTGGCAGGTGCATCTTTAAG.ATTTTGTTCCACGGaa          <  2:2329616/36‑1 (MQ=255)
                            tGGCAGGTGCATCTTTAGG.ATTTTGTTCCACGGAAg         >  2:2605163/1‑36 (MQ=255)
                            tGGCAGGTGCATCTTTAAG.ATTTTGTTCCACGGAAg         >  2:443341/1‑36 (MQ=255)
                            tGGCAGGTGCATCTTTAAG.ATTTTGTTCCACGGAAg         <  2:599734/36‑1 (MQ=255)
                            tGGCAGGTGCATCTTTAAG.ATTTTGTTCCACGGAAg         <  1:884733/36‑1 (MQ=255)
                            tGGCAGGTGCATCTTTAAG.ATTTTGTTCCACGGAAg         <  1:510720/36‑1 (MQ=255)
                            tGGCAGGTGCATCTTTAAG.ATTTTGTTCCACGAAAg         >  2:1383853/1‑36 (MQ=255)
                            tGGCAGGGGCATCTTTAAG.ATTTTGTTACAAGGaaa         >  2:233789/1‑35 (MQ=25)
                             ggCAGGTGCATCTTTAAG.ATTTTGTTCCACGGAAGa        >  1:59546/1‑36 (MQ=255)
                             ggCAGGTGCATCTTTAAG.ATTTTGTTCCACGGAAGa        <  1:108089/36‑1 (MQ=255)
                               cAGGTGCATCTTTAAG.ATTTTGTTCCACGGAAGATg      >  2:1801747/1‑36 (MQ=255)
                               cAGGTGCATCTTTAAG.ATTTTGTTCCACGGAAGATg      >  2:1256202/1‑36 (MQ=255)
                                aGGTGCATCTTTAAG.ATTTTGTTCCACGGAAGATGc     >  2:1080539/1‑36 (MQ=255)
                                 gnTNCNTCTTTAAG.ATTTTGTTCCACGGAAGATGCa    <  2:2665380/36‑1 (MQ=255)
                                  gTGCATCTTTAAG.ATTTTGTTCCACGGAAGATGCAg   <  2:1735264/36‑1 (MQ=255)
                                  gTGCATCTTTAAG.ATTTTGTTCCACGGAAGATGCAg   >  1:3610257/1‑36 (MQ=255)
                                   tGCATCTTTAAG.ATTTTGTTCCACGGAAGATGCAGt  <  2:2429582/36‑1 (MQ=255)
                                   tGCATCTTTAAG.ATTTTGTTCCACGGAAGATGCAGt  <  2:1031914/36‑1 (MQ=255)
                                   tGCATCTTTAAG.ATTTTGTTCCACGGAAGATGCAGt  >  1:1496916/1‑36 (MQ=255)
                                   |                                    
AGGTCTTCCTGGGTAATGACGCTGATACTGGCAGGTGCATCTTTAAG.ATTTTGTTCCACGGAAGATGCAGT  >  REL606/2197832‑2197902

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