Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,854,3880A→G33.3% 12.339251yeaOhypothetical protein
Reads supporting (aligned to +/- strand):  new base (13/0):  ref base (10/16):  total (23/16)

GCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTA  >  REL606/1854354‑1854421
                                  |                                 
gcgcGACACGAGCAAGATGGAGAGCGGCTGGCGGAc                                  <  1:1941746/36‑1 (MQ=37)
gcgcAACACGAGCAAGAAGGAAAGCGGCTGGCGGAc                                  <  2:629115/36‑1 (MQ=255)
 cgcAACACGAGCAAGAAGGAAAGCGGCTGGCGGACa                                 <  1:2051126/36‑1 (MQ=255)
  gcAACACGAGCAAGAAGGAAAGCGGCTGGCGGGCaa                                >  1:1121653/1‑35 (MQ=255)
   cAACACGAGCAAGAAGGAAAGCGGCTGGGGGGCAtg                               >  2:1669847/1‑35 (MQ=37)
   cAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATc                               <  2:2473239/36‑1 (MQ=255)
   cAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATc                               <  1:993970/36‑1 (MQ=255)
   cAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATc                               <  2:2600775/36‑1 (MQ=255)
   cAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATc                               <  2:812015/36‑1 (MQ=255)
      aacGAGCAAGAAGGAAAGCGGCTGGCGGACATCGcc                            <  1:390780/35‑1 (MQ=255)
       acGAGCAAGAAGGCAAGCGGCTGGCGGACATCGCCa                           <  2:2634186/36‑1 (MQ=255)
        cGAGCAAGAAGGAAAGCGGCTGGGGGGGATCGccca                          >  2:2019104/1‑34 (MQ=21)
        cGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCaa                          <  1:3336134/36‑1 (MQ=255)
        cGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCaa                          <  1:2589036/36‑1 (MQ=255)
         gAGCAAGAAGGAAAGCGGCTGGGGGGGAACGCCaca                         >  2:3428985/1‑34 (MQ=16)
         gAGCAAGAAGGAAAGCGGCTGGGGGAGACCGCCcaa                         >  2:1436380/1‑36 (MQ=16)
         gAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCaaa                         <  2:1245621/36‑1 (MQ=255)
          agaaagaagGAAAGCGGCTGGCGGACATCGCCaaaa                        <  1:2449932/33‑1 (MQ=255)
          aGCAAGAAGAAAAGCGGCTGGGGGGCATCCCCcaaa                        >  2:2863408/1‑36 (MQ=16)
           gCAAGAAGGAAAGCGGCTGGCGGACATCGCCaaaaa                       <  2:1033253/36‑1 (MQ=255)
             aagaagGAAAGCGGCTGGGGGGCATCGTCCAAAAAc                     >  1:3520629/1‑36 (MQ=19)
              agaagGAAAGCGGCTGGGCGGGATCGCCcaaaaaaa                    >  2:1806529/1‑34 (MQ=12)
               gaagGAAAGCGGCTGGGGGGCAATGCCcaaaaaaaa                   >  1:2823580/1‑33 (MQ=12)
                 agGAAAGCGGCTGGGGGACAATGCCAAAAAACaaac                 >  1:293884/1‑33 (MQ=19)
                 agGAAAGCGGCTGGCGGGCCTTGCCCAAAAACAGcc                 >  1:1811142/1‑36 (MQ=19)
                 agGAAAGCGGCTGGCGGACACCGCCAAAAAAAAGcc                 >  2:957895/1‑36 (MQ=37)
                  ggAAAGCGGCTGGCGGGCACCGCCCAAAAACAACCg                >  1:1065955/1‑36 (MQ=19)
                  aaaaaCCGGCTGCCGAAAACCGCAAAAAAACAGcca                <  2:286575/34‑2 (MQ=2)
                   gAAAGCGGCTGGGGGGGCACGCCCAAAAACAGCAGc               >  1:3499921/1‑36 (MQ=2)
                   gAAAGCGGCTGGGGGAGACCGGCAAACAACAGACGc               >  2:3357006/1‑36 (MQ=12)
                   gAAAGCGGCTGGCGGGCACCCCCAAAAAAAAACCGc               >  2:1814748/1‑36 (MQ=16)
                     aaGCGGCTGGCGGGCATCGGAAAAAAACAAGCGCTg             >  2:875998/1‑36 (MQ=16)
                     aaGCGGCTGGCGGACATCGCCAAAAAAAAACCggcg             >  1:316345/1‑33 (MQ=21)
                         ggCGGGCGAACAGCACAAAAACCCAGCCGCTGAccc         <  1:3194373/36‑1 (MQ=2)
                           cTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTg       >  2:2332605/1‑36 (MQ=255)
                           cTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTg       >  2:1198148/1‑36 (MQ=255)
                           cTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTg       >  1:3265316/1‑36 (MQ=255)
                               cGGACATCGCCAAAAAACAGCCGCTGACCCTGctct   >  2:1703772/1‑36 (MQ=255)
                                ggACATCGCCAAAAAACAGCCGCTGACCCTGCTCTa  >  2:1889490/1‑36 (MQ=255)
                                  |                                 
GCGCAACACGAGCAAGAAGGAAAGCGGCTGGCGGACATCGCCAAAAAACAGCCGCTGACCCTGCTCTA  >  REL606/1854354‑1854421

Base quality scores:  ATCG  < 3 ≤  ATCG  < 5 ≤  ATCG  < 9 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG