Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,735,6080GC52.9% 6.3 / 21.6 17A409P (GCC→CCC) ascFfused cellobiose/arbutin/salicin‑specific PTS enzymes: IIB component/IC component
Reads supporting (aligned to +/- strand):  ref base G (7/1);  new base C (0/9);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 4.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGCCTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCATCG  >  REL606/2735533‑2735683
                                                                           |                                                                             
cGCCTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTg                                                                               <  2:5014279/76‑1 (MQ=255)
 gCCTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgc                                                                             <  2:5185554/76‑1 (MQ=255)
  ccTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgcc                                                                            <  2:5553564/76‑1 (MQ=255)
   ctaaaaCCTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgcca                                                                           <  2:4567455/73‑1 (MQ=255)
   cTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgcca                                                                           <  1:1087741/76‑1 (MQ=255)
   cTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgcca                                                                           <  1:1953740/76‑1 (MQ=255)
    tGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgccag                                                                          <  2:3125856/76‑1 (MQ=255)
      aaaCGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgccagcc                                                                        <  2:2551098/76‑1 (MQ=255)
      aaaCGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgccagcc                                                                        <  1:2081515/76‑1 (MQ=255)
       aaCGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTgccagcca                                                                       <  1:4124092/76‑1 (MQ=255)
        aCGTCCTCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTTTGGCGGGGCTTCCCAGCCAc                                                                      <  2:6175801/76‑1 (MQ=255)
        aCGTCCTCTTATCCCCAGTCTTATCAGCGGTTTTATTGGCGGCGCGGTTGCCGGTATGGCGGGGCTTGCCAGCCAc                                                                      <  2:3346272/76‑1 (MQ=255)
           cccGCTTGTCGCCACGCTTCTCGGCGTCTTTTTTGGCGGGGCGGTTGCCGGTATGGGGGGGTTTCCCAGCCACTCa                                                                   <  2:445707/75‑1 (MQ=255)
             cGCTTGTCGCCATTCTTACCAGCGGCTTTTTTTGGGCGGCGATGGACGGTGTGGCGGGGCTTCCCAGCCACTCAAt                                                                 <  2:3099384/76‑1 (MQ=255)
              gtttATCGCCGTTTTTATCAGCGGCTTTTTTTGCGGCGCGGTTGCGGGTATGGGGGGGCTTCCCAGCCACTCAATg                                                                <  2:2278808/74‑1 (MQ=255)
                       cattCTTATCAGCGGCTTTATTTGCGGCGCGGTGGCCGGTAGGGGGGGGCTTCCCAGCCACTCAATGGCAGCGCCg                                                       <  1:2105578/73‑1 (MQ=255)
                        tgTCTTAGCAGCGGCTTTATGAGCGGGGCGGTCGCGGGTCGGGCGGGATTTGCCAGCCACTCAATGGCAGCGCCgg                                                      <  2:934261/75‑1 (MQ=255)
                        tgTCTTACCGGCGGTTTTATTGGCGGCGCGTTTGCGGGTAGGGGGGGGCTCCCCAGCCACTCAATGGCAGCGCCgg                                                      <  2:2213712/75‑1 (MQ=255)
                        aGTTTTATCAGGGGTTTTTTTTGGGGCGCGGTGCCCGGTAGGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCgg                                                      <  2:2352718/76‑1 (MQ=255)
                        aGTCTTATCGGCGGTTTTATTTGCGGCCGGTTTGCCGGTAGGGCGGGCCTTCCCAGCCACTCAATGGCAGCGCCgg                                                      <  2:1316216/76‑1 (MQ=255)
                        aGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCgg                                                      <  2:3100706/76‑1 (MQ=255)
                         gttttATCGGCGGTTTTTTTGGCGGCGCGGTTGCCGGTATGGGGGGGCTTCCCAGCCACTCAATGGCAGCGCCggg                                                     <  1:1666980/73‑1 (MQ=255)
                         gTCTTACCGGCGCTTTTATTTGCGGGGGGTTTGCCGGTAGGGGGGGGCTGCCCAGCCACTCAATGGCAGCGCCggg                                                     <  2:3890420/76‑1 (MQ=255)
                          tCTTATCAGGGGTTTTTTTTGCGGCGCGTTGGCCGGTAGGGGGGGGCTTCCCAGCCACTCAATGGCAGCGCCgggg                                                    <  1:5334170/76‑1 (MQ=255)
                           tttACCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTAGGGGGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGc                                                   <  1:4505537/75‑1 (MQ=255)
                            taaTCAGCGGCTTTCTTTGCGGCGCGGTGGCCGGTGGGGCGGGGTCTCCCAGCCCCTCAATGGCAGCGCCGGGGCt                                                  <  1:4737705/74‑1 (MQ=255)
                             ttgcAGCGGTTTTATTGGCGGCGGGGTGGCGGGTATGGGGGGCCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTa                                                 <  1:1344469/73‑1 (MQ=255)
                             tATCAGCGTCTTTATTTGCGGCGGGGTTGCGGGTTGGGCGGGCCTCGCCAGCCACTCAATGGCAGCGCCGGGGCTa                                                 <  2:4561824/76‑1 (MQ=255)
                              aTCAGCGGTTTTATTGGCGGCGCGGTTGCCGGTAGGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTAt                                                <  1:5356120/76‑1 (MQ=255)
                               tcggCGGTTTTATTTGCGGCGCTGTCGCCGTTAGGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTAtt                                               <  2:6155099/73‑1 (MQ=255)
                               tCAGCGGTTTTATTTGCGGCGCGGTTGCGGGTATGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTAtt                                               <  2:750449/76‑1 (MQ=255)
                                cAGCGGCTTTATTTGCGGCGGGGTTGCGGGTAGGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTAttt                                              <  1:1397922/76‑1 (MQ=255)
                                         tATTTGGGGCGCGGTTTCGGGTAGGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTg                                     <  2:3123684/76‑1 (MQ=255)
                                           ttggcggcgGGTTGGCGGGTAGGGCGGGCTTCCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCa                                  <  1:2017301/73‑1 (MQ=255)
                                               gggcgTGGTGGCCGGTTGGGGGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAAttc                              <  1:3673147/75‑1 (MQ=255)
                                               gggcgCGGTTGCGGGTAGGGGGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAAttc                              <  1:31558/75‑1 (MQ=255)
                                               cggcgCGGTTCCGGATATGGCGGGCCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAAttc                              <  1:5573638/76‑1 (MQ=255)
                                                 gcgcGGTGCCCGTTAGGGGGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAAttctt                            <  2:450213/76‑1 (MQ=255)
                                                 gcccgGTGGCCGGTATGGGGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAAttctt                            <  2:5626107/73‑1 (MQ=255)
                                                    gggTTGCCGGTAGGGCGGGCCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGa                         <  1:1112960/75‑1 (MQ=255)
                                                       ttGCCGGTATGGCGGGGCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATcc                      <  2:4645079/76‑1 (MQ=255)
                                                        tGCGGGTAGGGCGGGGCTCCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCg                     <  1:4050695/76‑1 (MQ=255)
                                                                gggCGGGTCTCCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATcc             <  2:4407047/75‑1 (MQ=255)
                                                                 gggggggCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCa            <  1:2784602/73‑1 (MQ=255)
                                                                 gggggggCTCCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCa            <  1:1208107/73‑1 (MQ=255)
                                                                 ggCGGGGCTTGCCAGCCACTCAATGGCAGGGCCGGGGGTATTTTACAACGTGCAAATTTTTCGTTCCGCGGAATca             >  2:2706160/1‑75 (MQ=255)
                                                                 ggCGGGGCTTGCCAGCCACTCAATGGCAGGGCCGGGGCTATTTACCAGCGTGGAAATCTTCCATCCCGCGAAACCa            >  2:2240595/1‑76 (MQ=255)
                                                                 ggCGGGGCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTTTTTGATCCGGCGGATccc            >  2:1184668/1‑75 (MQ=255)
                                                                 ggCGGGGCTCCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCa            <  1:1534283/76‑1 (MQ=255)
                                                                  ggggggCCTGCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCa            <  2:378937/74‑1 (MQ=255)
                                                                    gggtCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATg         <  1:2139572/76‑1 (MQ=255)
                                                                    ggggCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATg         >  1:366900/1‑76 (MQ=255)
                                                                    ggggCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCATCGTGCAATTCTTCGATCCGGCGAATCCAATg         <  2:6168669/76‑1 (MQ=255)
                                                                    ggggCTGCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATg         <  1:6066234/76‑1 (MQ=255)
                                                                    ggggCTCCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATg         <  2:6031893/76‑1 (MQ=255)
                                                                    gctgCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATg         <  2:1965515/73‑1 (MQ=255)
                                                                     gggCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGa        <  1:6189702/76‑1 (MQ=255)
                                                                     gggCTTCCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGa        <  1:4642109/76‑1 (MQ=255)
                                                                     gggCTTCCAAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGa        <  1:388963/76‑1 (MQ=255)
                                                                      ggCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAg       <  1:3152048/76‑1 (MQ=255)
                                                                       gCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGc      >  2:2272733/1‑76 (MQ=255)
                                                                       gCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCCGCGAATCCAATGAGc      >  2:367742/1‑76 (MQ=255)
                                                                         ttGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCAt    >  2:4508876/1‑76 (MQ=255)
                                                                          cgccagccaCTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCATc   <  2:1271093/75‑1 (MQ=255)
                                                                           gccagccaCTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCATCg  >  1:3989906/1‑76 (MQ=255)
                                                                           |                                                                             
CGCCTGAAACGTCCGCTTATCGCCAGTCTTATCAGCGGCTTTATTTGCGGCGCGGTTGCCGGTATGGCGGGGCTTGCCAGCCACTCAATGGCAGCGCCGGGGCTATTTACCAGCGTGCAATTCTTCGATCCGGCGAATCCAATGAGCATCG  >  REL606/2735533‑2735683

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