Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,735,6080GC45.5% 34.3 / 19.9 33A409P (GCC→CCC) ascFfused cellobiose/arbutin/salicin‑specific PTS enzymes: IIB component/IC component
Reads supporting (aligned to +/- strand):  ref base G (16/2);  new base C (0/15);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.47e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCATCG  >  REL606/2735534‑2735683
                                                                          |                                                                             
gCCTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgc                                                                             <  1:2646571/76‑1 (MQ=255)
 ccTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgcc                                                                            <  1:6427424/76‑1 (MQ=255)
 ccTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgcc                                                                            <  1:4192894/76‑1 (MQ=255)
  cTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgcca                                                                           <  2:6739776/76‑1 (MQ=255)
  cTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgcca                                                                           <  1:1419423/76‑1 (MQ=255)
  cTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgcca                                                                           <  1:3510221/76‑1 (MQ=255)
   tggaaCGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgccag                                                                         <  1:1773740/73‑1 (MQ=255)
   tGAACCGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgccag                                                                         <  2:6507858/76‑1 (MQ=255)
   tGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgccag                                                                         <  2:4204541/76‑1 (MQ=255)
   tGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgccag                                                                         <  1:71400/76‑1 (MQ=255)
   tGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgccag                                                                         <  1:2748822/76‑1 (MQ=255)
   tGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgccag                                                                         <  1:1235147/76‑1 (MQ=255)
   tGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgccag                                                                         <  1:1347919/76‑1 (MQ=255)
   tGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgccag                                                                         <  2:3693405/76‑1 (MQ=255)
      aaCGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgccagcca                                                                      <  1:5224835/76‑1 (MQ=255)
            cttttATCGCCAGTTTTATCAGGGGCTTTATTGGCGGCGGGTTTGCCGTTATGGGGGGGCTTCCCAGCCACTCAAt                                                                <  2:5988532/73‑1 (MQ=255)
             gtttATCGCCAGTCTTATCAGCGGCTTTATTGGCGGCGCGGTTGCCGGTATGGCGGGGCTTCCCAGCCACTCAATg                                                               <  2:1698120/74‑1 (MQ=255)
               ttATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTCCCAGCCACTCAATGGc                                                             <  2:2593650/76‑1 (MQ=255)
                ttccGCCATTCTTATAAGCGGCTTTAGTGGCGGGGCGGTTTCGGGGGGGGGGGGGCTGGCCAGCCACTCAATGGCa                                                            <  2:776591/73‑1 (MQ=255)
                taccGCCATTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTCCCAGCCACTCAATGGCa                                                            <  1:2550483/73‑1 (MQ=255)
                tATCGCCAGTCTTAGCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTGCCAGCCACTCAATGGCa                                                            <  1:4212410/76‑1 (MQ=255)
                 aTCGCCGGGCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTCCCAGCCACTCAATGGCAg                                                           <  1:6692774/76‑1 (MQ=255)
                 aTCGCCAGTTTTATCAGCCGCTTTATTTGCGGCGGGGTTGCGGGTATGGCGGGTCTTGCCAGCCACTCAATGGCAg                                                           <  1:4545087/76‑1 (MQ=255)
                  tccccAGTTTTACCAGCGTCTTTTTTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTCCCAGCCACTCAATGGCAgc                                                          <  2:6063656/73‑1 (MQ=255)
                  tCGCCAGTCTTACCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTAGGGGGGGGCTTCCCAGCCACTCAATGGCAgc                                                          <  1:2570634/76‑1 (MQ=255)
                  ccgcgAGTCTTACCAGGGGCTTTTTTTGGGGCGGGGTTGCGGGTAGGGGGGGGCTTCCCAGCCACTCAATGGCAgc                                                          <  2:4273306/75‑1 (MQ=255)
                   ctccAGTCTGTTCAGCGGTTTTATTTGCGGCGCGGTTCCGAGTTGGGCGGGCCTTCCCAGCCACTCAATGGCAgcg                                                         <  2:6896566/74‑1 (MQ=255)
                        gTCTTATCAGCGGCTTTATTGGCGGCGCGGTTGCCGGTATGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCggg                                                    <  1:6570170/76‑1 (MQ=255)
                         tCTTATCAGCTGCTTTTTGTGCGGCGCGTTGGCGGGTAGGGGGGGGCTGCCCAGCCACTCAATGGCAGCGCCgggg                                                   <  1:2119938/76‑1 (MQ=255)
                          cTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGGGTGGCAGGCCCCCCAATGGCcggggcgggg                                                   >  1:5306374/1‑76 (MQ=255)
                           ttATCAGCGGCTTTATTGGCGGCGCGGTGGCCGGTAGGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCt                                                 <  1:1425274/76‑1 (MQ=255)
                            tATCAGCGGCTTTATTGGCGGCGCGGTTGCCGGTATGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTa                                                <  1:3278162/76‑1 (MQ=255)
                              tccgCGGTTTTATTTGCGGCGGGGTTGGCGGTATGGCGGCGCTCCCCAGCCACTCAATGGCAGCGCCGGGGCTAtt                                              <  1:4484403/73‑1 (MQ=255)
                               cAGCGTTTTTATTTGCGGCGCGGTGGCGGGTAGGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTAttt                                             <  2:4680904/76‑1 (MQ=255)
                                aGCGGTTTTTTTGGCGGCGCGGTTGCCGGTATGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTa                                            <  1:658557/76‑1 (MQ=255)
                                 gtggCTTTTTTGGCGGCGCGGTTGCCGGTGGGGCGGGGCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTAc                                           <  2:5861890/74‑1 (MQ=255)
                                 gctgCTTTATTTGCGGCGGGGTTGCGGGTAGGGGGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTAc                                           <  2:6239019/73‑1 (MQ=255)
                                 gCGGTTTTATTGGCGGCGGGGTGGCGGGTAGGGGGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTAc                                           <  1:890221/76‑1 (MQ=255)
                                  cGGTTTTTTTTGGGCCGGGGTGGCGGGTAGGGCGGGGCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTAcc                                          <  2:2620496/76‑1 (MQ=255)
                                  cGGTTTTTTTGGGGGCGGGGTTGCCGGTGTGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTAcc                                          <  1:6324976/76‑1 (MQ=255)
                                   ggttttATGGGCGGCGCGGTGGCGGGTAGGGGGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCa                                         <  2:457818/73‑1 (MQ=255)
                                    gCTTTATTGGCGGCGCGGTTGCCGGTATGGCGGGGCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAg                                        <  1:5746961/76‑1 (MQ=255)
                                      ttttttGGCGGCGCGTTTGCCGGTATGGCGGGCCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCg                                      <  1:4018520/76‑1 (MQ=255)
                                       ttATTGGCGGCTAGATTCCCGGGATGGCGGGCCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGt                                     <  2:4334639/76‑1 (MQ=255)
                                        tATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTg                                    <  1:5350933/76‑1 (MQ=255)
                                        tATGTGCGGCGCTGTTCCCGGTATGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTg                                    <  2:1251968/76‑1 (MQ=255)
                                          ttcgcggggCGTTGGCGGTTAGGGGGGGGTTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCa                                  <  2:426442/73‑1 (MQ=255)
                                           ttGCGCCGCGTTTGCGGGTAGGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCaa                                 <  2:1705203/76‑1 (MQ=255)
                                           gtGCGGCGGTGTTGCGGGTAGGGGGGGCTTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCaa                                 <  2:6921235/75‑1 (MQ=255)
                                             ggggcgCGGTGGCGGGTAGGGGGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAAtt                               <  1:2998331/74‑1 (MQ=255)
                                             gcggcgGTGTTGCGGGTAGGGCGGGGCTCCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAAtt                               <  1:127542/76‑1 (MQ=255)
                                              cgacgcgGTTGCCGGTATGGGGGGGTTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAAttc                              <  2:2808229/73‑1 (MQ=255)
                                               ggcgGGGTTGCGGGTAGGGGGGGGCTGGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAAttct                             <  1:2537944/76‑1 (MQ=255)
                                                gcggGGTTGCGGGTGGGGCGGGGCTCCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAAttctt                            <  1:5471076/76‑1 (MQ=255)
                                                    gtttGCGGGTGGGGGGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGAt                        <  2:4621904/74‑1 (MQ=255)
                                                    ggTTGCGGGTGGGGGGGGCCTCCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGAt                        <  1:6904913/76‑1 (MQ=255)
                                                    ggTTGCCGGTATGGCGGGCCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGAt                        <  1:1919690/76‑1 (MQ=255)
                                                    ggTGGCCGGTGTGGCGGGCATGGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGAt                        <  2:7684236/76‑1 (MQ=255)
                                                        gcgggTGTGGCGGGCCTACCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCgg                    <  1:104528/73‑1 (MQ=255)
                                                            gTATGGGGGGGCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGaa                <  1:3775242/76‑1 (MQ=255)
                                                            gTATGGCGGGGCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTTTTTTCCAGCGGGCAATTCTTCGATCCGGCgga                >  1:1233241/1‑74 (MQ=255)
                                                            gTATGGCGGGGCTTGCCAGCCACTCAAAGGCAGCGCCGGGGGCCTTTTTCCGCGGGCAAATTTTCTATCCcgcgga                >  2:3182756/1‑74 (MQ=255)
                                                            gTAGGGGGGGCCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGaa                <  1:4349568/76‑1 (MQ=255)
                                                                ggggggtCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCa            <  1:447076/73‑1 (MQ=255)
                                                                ggCGGGGCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGGATCCa            >  1:5264157/1‑76 (MQ=255)
                                                                 gCGGGGCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTAACAGCGTGCAATTCTTCGATCCGGCGAATCCaa           >  1:7702079/1‑76 (MQ=255)
                                                                  cGGGGCTTGCCAGCCACTCAATGGCAGGGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCCAt          >  2:1258461/1‑76 (MQ=255)
                                                                  cGGGGCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCCGCGAATCCaca          >  2:5082554/1‑74 (MQ=255)
                                                                   ggggCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCCGCGTGCAATTCTTCGATCCCGCGGATCCAATg         >  2:5742292/1‑76 (MQ=255)
                                                                   ggggCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCCGCGTGCAATTCTTCGATCCCGCGAATCCCATg         >  2:259283/1‑76 (MQ=255)
                                                                    gggCTTGCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATg         <  1:1430561/76‑1 (MQ=255)
                                                                    gggCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGa        <  1:4571381/76‑1 (MQ=255)
                                                                    gggCCTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGa        <  2:7739045/76‑1 (MQ=255)
                                                                    ggccTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGa        <  2:5654064/73‑1 (MQ=255)
                                                                     ggCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAACGTGCAATTCTTCGATCCGGCGAATCCAATGAg       >  2:3899806/1‑76 (MQ=255)
                                                                       cTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCa     >  2:193853/1‑76 (MQ=255)
                                                                       cTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCa     >  1:4087326/1‑76 (MQ=255)
                                                                       cTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCa     >  1:2447663/1‑76 (MQ=255)
                                                                       cTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCa     >  2:596446/1‑76 (MQ=255)
                                                                        ttGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCAt    >  1:936781/1‑76 (MQ=255)
                                                                        ttGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCAt    >  2:888088/1‑76 (MQ=255)
                                                                         tGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCATc   >  2:2643225/1‑76 (MQ=255)
                                                                         tGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCATc   >  2:7432928/1‑76 (MQ=255)
                                                                          gccagccaCTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCATCg  >  2:91264/1‑76 (MQ=255)
                                                                          |                                                                             
GCCTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCATCG  >  REL606/2735534‑2735683

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: