Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606964,8050A→G27.4% 14.963484ycaMpredicted transporter
Reads supporting (aligned to +/- strand):  new base (17/0):  ref base (13/32):  total (31/32)

TGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCGCTCCTGTATGTGGCG  >  REL606/964770‑964840
                                   |                                   
tGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGa                                     <  2:2784379/36‑1 (MQ=255)
 gAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGAt                                    <  2:786346/36‑1 (MQ=255)
 gAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGAt                                    <  1:3733495/36‑1 (MQ=255)
  aaaTCGCAGAAAATCGTCGGTTCTGTGGCAGGGAtt                                   <  2:3265076/36‑1 (MQ=255)
    aTCGCTGAAAATCGTCGGTTCTGTGGGTGGGTTTgt                                 >  2:3670926/1‑35 (MQ=25)
    aTCGCTGAAAATCGTCGGTTCTGTGGCAGGGAtttc                                 >  2:1001320/1‑34 (MQ=255)
    aTCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGc                                 <  1:1258269/36‑1 (MQ=255)
     tCGCTGGAAAGCGTCGGTTCTGTGGCAGGGATTGCt                                <  1:290232/36‑1 (MQ=37)
     tCGCTGAAAATCGTCGGTTCTGTGGCAGGGGttatt                                >  1:604957/1‑33 (MQ=37)
      cGCTGAAAATCGTCGGTTCTGTGGCAGGGGTTgttt                               >  1:1053490/1‑33 (MQ=37)
       gCTGAAAATCGTCGGTTCTGTGGCAGGGGTTTCttt                              >  2:3728955/1‑34 (MQ=37)
       gCTGAAAATCGTCGGTTCTGTGGCAGGGGTTGCttt                              >  2:1722911/1‑34 (MQ=38)
       gCTGAAAATCGTCGGTTCTGTGGCAGGGAGTGCttt                              >  1:3690634/1‑34 (MQ=38)
        cTGAAAATCGTCGGTTCTGTGGCAGGGGTTTGTAtg                             >  1:2238449/1‑36 (MQ=25)
        cTGAAAATCGTCGGTTCTGTGGCAGGGGTTTCTAtg                             >  2:433269/1‑36 (MQ=37)
          gAAAATCGTCGGTTCTGTGGCAGGGGTTTCTAtttt                           >  1:1454435/1‑36 (MQ=25)
          gAAAATCGTCGGTTCTGTGGCAGGGGTTGCTAtgtt                           >  1:3387981/1‑36 (MQ=255)
          gAAAATCGTCGGTTCTGTGGCAGGGATTGCTAtgtt                           <  1:3182301/36‑1 (MQ=255)
           aaatTCGTTGGTTCTGTGGGAGGGATTGCTAtgttt                          <  1:2894868/36‑1 (MQ=37)
           aaaaTCGTCGGTTCTGTGGCTGGGGTTTTTTtgttt                          >  2:2230944/1‑36 (MQ=16)
           aaaaTCGTCGGTTCTGTGGCAGGGGTTGCTTttttt                          >  1:2702527/1‑36 (MQ=25)
            agaTCGTCGGTTCTGTGGCAGGGATTGCTAtgtttg                         <  1:2439375/34‑1 (MQ=255)
            aaaTCGTCGGTTCTGTGGCAGGGGCTTCTAtgtttt                         >  1:3669376/1‑35 (MQ=25)
            aaaCCGTCGGTTCTGTGGCAGGGATTGCTAtgtttg                         <  1:3631667/36‑1 (MQ=255)
             aaTCGTCGGTTCTGTGGCCGGGGTTGCTTtttttgt                        >  1:2297416/1‑36 (MQ=19)
             aaTCGTCGGTTCTGTGGCAGGGGTTTCTTttttttt                        >  2:2698974/1‑34 (MQ=16)
             aaTCGTCGGTTCTGTGGCAGGGGTTGCTAtgttttt                        >  2:2151270/1‑34 (MQ=38)
             aaTCGGCGGTTCGGGGGCAGGGATTGCTAtgtttgt                        <  1:3066341/36‑1 (MQ=37)
              aTCGTTGGTTCTGTGGCAGGGATTGCTAtgtttgtt                       <  2:1855648/36‑1 (MQ=255)
              aTCGTCGGTTCTGTGGCAGGGATTGCTAtgtttgtt                       <  1:3185971/36‑1 (MQ=255)
              aTCGTCGGTTCTGTGGCAGGGATTGCTAtgtttgtt                       <  1:1418652/36‑1 (MQ=255)
              aTCGTCGGTTCTGTGGCAGGGATTGCTAtgtttgtt                       <  1:3380639/36‑1 (MQ=255)
              aTCGGCGGTTCTGTGGCATGGATTGCTAtgtttgtt                       <  1:1428986/36‑1 (MQ=38)
               tcgtcgGTTCTGTGGCAGGGATTGCTAtgttttttt                      >  2:2131115/1‑35 (MQ=255)
               tcgtcgGTTCTGTGGCAGGGATTGCTATGTTTGTTa                      <  1:2479177/36‑1 (MQ=255)
                cgtcgGTTCTGTGGCAGGGATTGCTATGTTTGTTAt                     <  2:3083335/36‑1 (MQ=255)
                cgtcgGTTCTGTGGCAGGGATTGCTATGTTTGTTAt                     <  2:1917284/36‑1 (MQ=255)
                cgtcgGTTCTGTGGCAGGGATTGCTATGTTTGTTAt                     <  1:3502802/36‑1 (MQ=255)
                 gtcgGTTCTGTGGCAGGGATTGCTATGTTTGTTATg                    <  1:2954495/36‑1 (MQ=255)
                    ggTTCTGTGGCAGGGATGGCTATGTTTGTTATGTCg                 <  2:3449797/36‑1 (MQ=255)
                     gTTCTGTGGCGAGGATTGCTATGTTGTTTATTTCGc                <  1:144873/36‑1 (MQ=16)
                     gTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCGc                <  1:3370905/36‑1 (MQ=255)
                     gTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCGc                <  1:147076/36‑1 (MQ=255)
                      ttCTGTGGCAGGGATTGCTATGTTTGTTATATCGCt               >  2:1774287/1‑36 (MQ=255)
                      ttCTGTGGCAGGGATTGCTAAGTTTGTTATTTCGCt               >  2:2957053/1‑36 (MQ=37)
                       tCTGTGGCAGGGGTTGCTTTGTTTTTTATGTCTCtt              >  2:3532646/1‑35 (MQ=16)
                         tgtgGCAGGGGTTGGTATGTTTTTTTTTTCGCTCCt            >  2:3237564/1‑36 (MQ=16)
                          gtgGCAGGGGTTGCTATGTTTGTTATTTCGCTCCTg           >  2:2002458/1‑36 (MQ=38)
                          gtgGCAGGGATTGCTTTGTTTTTTATGTCGCTCCtt           >  2:1639124/1‑35 (MQ=37)
                          gtgGCAGGGATTGCTATGTTTTTTTTTTTGCTCCtt           >  1:1545877/1‑35 (MQ=16)
                          gtgGCAGGGATTGCTATGTTTTTTATGTCGCTCCtt           >  2:1265941/1‑35 (MQ=255)
                          gtgGCAGGGATTGCTATGTTTGTTATGTCGCTCCtt           >  1:3628290/1‑35 (MQ=255)
                          gtgGCAGGGATTGCTATGTTTGTTATGTCGCTCCTg           >  2:1648409/1‑36 (MQ=255)
                           tgGCAGGGATTGCTATGTTTGTTATGCCGCTCCttt          >  1:537320/1‑34 (MQ=37)
                           tagCAGGGATTGCTATGTTTGTTATGTCGCTCCTGt          <  2:3077527/34‑1 (MQ=255)
                            ggCAGGGATTGCTATGTTTGTTATGTCGCTCCTGTa         <  1:243024/36‑1 (MQ=255)
                             gCAGGGATTGCTATGTTTGTTATGTCGCTCCTGTAt        <  2:2954806/36‑1 (MQ=255)
                             gCAGGGATTGCTATGTTTGTTATGTCGCTCCTGTAt        <  1:2079329/36‑1 (MQ=255)
                              cAGGGATTGCTATGTTTGTTATGTCGCTCCTGTAtg       >  2:362624/1‑36 (MQ=255)
                               aGGGATTGCTATGTTTGTTATGTCGCTCCTGTAtgt      <  1:3575101/36‑1 (MQ=255)
                               aGGGATTGCTATGTTTGTTATGTCGCTCCTGTAtgt      >  1:3445778/1‑36 (MQ=255)
                               aGGGATTGCTATGTTTGTTATGTCGCTCCTGTAtgt      <  1:1882612/36‑1 (MQ=255)
                                  gATTGCTATGTTTGTTATGTCGCTCCTGTATGTGGc   <  2:2901429/36‑1 (MQ=255)
                                  gATTGCTATGTTTGTTATGTCGCTCCTGTATGTGGc   <  2:3186938/36‑1 (MQ=255)
                                   aTTGCTATGTTTGTTATGTCGCTCCTGTATGTGGCg  >  1:3192455/1‑36 (MQ=255)
                                   |                                   
TGAAATCGCTGAAAATCGTCGGTTCTGTGGCAGGGATTGCTATGTTTGTTATGTCGCTCCTGTATGTGGCG  >  REL606/964770‑964840

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