Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,722,5320CG34.9% 111.7 / 33.4 63G104G (GGC→GGGsrlBglucitol/sorbitol‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base C (7/34);  new base G (22/0);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 2.96e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 191 total aligned reads displayed.

CGAACTGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTACCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTAAAAAATGAATCAGGT  >  REL606/2722457‑2722606
                                                                            |                                                                          
 gAACTGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGcc                                                                           >  2:5455350/1‑76 (MQ=255)
  aaCTGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGgcc                                                                          >  1:7206991/1‑76 (MQ=255)
  aaCTGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGgcc                                                                          >  1:3994362/1‑76 (MQ=255)
   aCTGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGgccc                                                                         >  1:5911427/1‑75 (MQ=255)
    cTGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTg                                                                        <  1:8189296/76‑1 (MQ=255)
    cTGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTg                                                                        <  2:6875800/76‑1 (MQ=255)
     tGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGt                                                                       <  2:4965079/76‑1 (MQ=255)
     tGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGt                                                                       <  1:6045197/76‑1 (MQ=255)
      gggTTATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTa                                                                      <  1:1442720/76‑1 (MQ=255)
      gggTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCCGTa                                                                      >  2:4326834/1‑76 (MQ=255)
      gggTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTa                                                                      <  1:1148733/76‑1 (MQ=255)
      gggTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTa                                                                      <  2:3404632/76‑1 (MQ=255)
      gggTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTa                                                                      <  1:2525118/76‑1 (MQ=255)
       ggTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTggat                                                                     >  1:3692903/1‑73 (MQ=255)
       ggTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTAc                                                                     <  1:7934902/76‑1 (MQ=255)
        gTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCTGTGcc                                                                    >  1:6488411/1‑76 (MQ=255)
          cATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCCGTGCCcg                                                                  >  1:7270314/1‑76 (MQ=255)
           aTGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCTGGAcccca                                                                 >  1:4126183/1‑74 (MQ=255)
           aTGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCTGGACCcgc                                                                 >  1:4924250/1‑75 (MQ=255)
            tGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCTGTGcccccc                                                                >  1:3318755/1‑73 (MQ=255)
            tGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTACCcgac                                                                <  2:2992129/76‑1 (MQ=255)
               cccccTGCTTTTCGATGGTTAAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGCCCCTGTCCCCGACGat                                                             <  2:2570816/74‑1 (MQ=255)
               cACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCTGTTCCcgccgac                                                             >  2:111769/1‑75 (MQ=255)
               cACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCCGTACCCGACGat                                                             >  1:3198022/1‑76 (MQ=255)
                aCCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCTGTGCCCGACGata                                                            >  1:7438730/1‑76 (MQ=255)
                   cTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCCGCACCCCACGATATcg                                                         >  1:6784856/1‑76 (MQ=255)
                   cTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTACCCGACGATATcg                                                         <  1:3911932/76‑1 (MQ=255)
                     gctcTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTCCCCGACGATATcgcg                                                       <  1:709858/73‑1 (MQ=255)
                        cTTCGTTGGTTAAACCGAACGGGAATTTCCGGGCCCTGTCCATGGGGCAGGCCCTGTGCCCGACGATATCGCGCCg                                                    <  2:2826720/76‑1 (MQ=255)
                        cTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCCGTACCCCACGCTAACGCGCCg                                                    >  1:4663629/1‑76 (MQ=255)
                         ttCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCTGTACCCGCCGATATCTCGcccg                                                   >  2:2330634/1‑74 (MQ=255)
                          tcaaTGGTTTAAACGAAGCGAAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTACCCGACGATATCGCGCCggg                                                  <  1:5817927/73‑1 (MQ=255)
                          tCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCCGTACCCGACGATAATGCGCCggg                                                  >  1:5851940/1‑76 (MQ=255)
                          tCGATGGTTAAAACGAAGCGAAATTTCCGGGCACTGTCGATGTGGCAGCCCCTGTACCCGACGATATCGCGCCggg                                                  <  2:5890870/76‑1 (MQ=255)
                           cGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCTGTTCCCGACGATATCGCGCCcggg                                                 >  1:5082354/1‑75 (MQ=255)
                           cGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCTGTACCCCACGATATCGCGCCcggg                                                 >  2:2963182/1‑75 (MQ=255)
                            gagggTATAAACGAGGCGAAATTTCGGGGCACTGTCCATGTGGCGGGCCCTGTACCCGACGATATCGCGCCGGGAt                                                <  2:1546329/73‑1 (MQ=255)
                            gATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTCCCCGACGATATCGCCCCGGGAt                                                <  2:7553781/76‑1 (MQ=255)
                             ttGGTTTAAAGGAAGCGGAATTTCGGGGCCCGGCCCATGTGGCAGTCCCTGTACCCGACGATATCGCGCAGGGATc                                               <  1:1240117/75‑1 (MQ=255)
                             aTGGTTTAAACGAACCGGACTTCCCGGGCACTGTCCATGTGGCAGGCCCTGTCCCCGACGATATCGCGCCGGGATc                                               <  2:7679791/76‑1 (MQ=255)
                              tGGTTTAAACAAAGCGGAATTCCCGGGCCCTGTCCATGTGGCGGCCCCGGTACCCGACGATATCGCGCCGGGATCg                                              <  1:7857452/76‑1 (MQ=255)
                              tGGTTAAAAGAAAGCGGATTTCCCGGTCACTGTCCATGTGGCAGGCCTTGCCCCCGACGATATCGCGCCGGGATCg                                              <  2:6992563/76‑1 (MQ=255)
                               gttttAAACGAAGGGGAATTTCGGGGCACTGTCCATGTGGCAGGCCCTGTCCCCGACGATATCGCGCCGGGATCgg                                             <  2:5686538/74‑1 (MQ=255)
                               ggTTTAAACGAAGCGGATTTTCGGGGCACTGCCCATGTGGCAGGCCCTGTACCCGACGATATCGCGCCGGGATCgg                                             <  2:3644520/76‑1 (MQ=255)
                                gTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCTGTACCCCACGATATTGCGCCGGGGTCggg                                            >  1:7728085/1‑75 (MQ=255)
                                 tttAAACGAGGCGAATTTTCGGGGCACTGTCCATGGGGCAGCCCCTGTCCCCGACGATATCGCGCCGGGATCGGtt                                           <  2:4377914/76‑1 (MQ=255)
                                 tttAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCTGTACCCGACGATATCGCGCCGGGGTCGGtt                                           >  2:7627298/1‑76 (MQ=255)
                                 tttAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCCGTCCCCCACGATAACGCGCCCGGGTCgggt                                           >  2:4780816/1‑74 (MQ=255)
                                 tttAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCCATACCACACGACAGCGGGCCGGGGTCggct                                           >  1:2669630/1‑74 (MQ=255)
                                 ttaaaaCGAGGCGGAATTGCGGGGCACTGCCCATGTGGCAGGCCCTGTACCCGACGATATCGCGCCGGGATCGGtt                                           <  2:1149367/73‑1 (MQ=255)
                                 ttaaaaCGAAGCGGTATTTCGGGGCCCTGTCGATGTGGCAGGCCTTGTCCCCGACGATATCGCGCCGGGATCGGtt                                           <  1:2873669/73‑1 (MQ=255)
                                  ttAAACGAAGCGGATTTTCCGGGCACTGTCCATGTGGCAGGCCCTGTACCCGACGATATCGCGCCGGGATCGGttt                                          <  2:4603725/76‑1 (MQ=255)
                                  ttAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCTGTACCCGACGATATCGCGCCGGGGCCGGttt                                          >  2:3725852/1‑76 (MQ=255)
                                  ttAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCCGTTCCCCGCGCGATCGCGCCGGGAGCggggt                                          >  1:6540651/1‑73 (MQ=255)
                                  ttAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCCGTACCCGACGATATCGGGCCGCGGCCGGttt                                          >  2:6029896/1‑76 (MQ=255)
                                  taaaaCGAAGCGGATTTCCCGGGCATTGTCCATGTGGCAGGCCCTGTCCCCGACGATATCGCGCCGGGATCGGttt                                          <  2:8141605/74‑1 (MQ=255)
                                   tAACCGAAGCGGATTTTCCGGGCACTGTCCATGTGGCAGGCCCTGTCCCCGACGATATCGCGCCGGGATCGGtttt                                         <  1:7357690/76‑1 (MQ=255)
                                   tAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCCGTACCCGACGACATCACGCCGGGGTCGGtttt                                         >  1:2088765/1‑76 (MQ=255)
                                   aaaaCGAAGCGGATTTTCCGGGCACTGTCCATGTGGCGGGCCCGGTCCCCGACGATATCGCGCCGGGATCGGtttt                                         <  1:2366022/75‑1 (MQ=255)
                                    aaaGGATGCGGATTTCCCGGGCCCTGCCAATGGGGCAGCCCCTGTCCCCGACGATATCGCGCCGGGATCGGTTTTg                                        <  1:2389350/76‑1 (MQ=255)
                                     aaCGAGGCGGAATTTCCGGGCACTGTCCATGTGGCAGCCCCTGTACCCGACGATATCGCGCCGGGATCGGTTTTGa                                       <  2:2258219/76‑1 (MQ=255)
                                       cGAACCGGATTTTCGGGGCACTGCCCATGGGGCAGGCCCTGTCCCCGACGATATCGCGCCGGGATCGGTTTTGAAg                                     <  1:5357072/76‑1 (MQ=255)
                                         aaGGGGAGTTCCCGGGCACTGTCCATGTGGCAGGCCCTGTACCCGACGATATCGCGCCGGGATCGGTTTTGAAGtt                                   <  1:5743686/76‑1 (MQ=255)
                                           gCGGATTTTCCGGGCACTGTCCATGTGGCAGGCCCTGTCCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTg                                 <  2:3316380/76‑1 (MQ=255)
                                            cGGAATTTCCGGGCACTGTCCATGTGGCAGGGCCTGTACCCGCCGATAACGCGGCCGGGTCGGGTTTTGAGtttta                                >  2:4366758/1‑74 (MQ=255)
                                             ggAATTTCCGGGCACTGTCCATGTGGCAGGGCCTGTACCCGACGATATCGCGCCGGGGGCGGGTTTGTAGTTTTaa                               >  1:6719227/1‑76 (MQ=255)
                                             cgAATTTCCCGGGCCTTGCCCGTGTGGCAGGCCCTTTACCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGa                                <  1:4882506/75‑1 (MQ=255)
                                                 tttCCGGGCACTGTCCATGTGGCAGGCCCTGTACCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTg                           <  2:1883199/76‑1 (MQ=255)
                                                  ttcccGGGCACTGTCCATGTGGCAGGCCCTGTACCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTg                           <  2:7731386/73‑1 (MQ=255)
                                                  ttCCGGGCACTGTCCATGTGGCAGGCCCCCTACCCGGCGGTATCGCGCCGGGGTGGGGTTTTGAGTTTTAATTTGt                          >  1:5888982/1‑76 (MQ=255)
                                                   tCCGGGCACTGTCCATTGGGCAGGCCCTGACCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGtt                         <  2:3175691/76‑1 (MQ=255)
                                                       gggCACTGTCCATGTGGCAGGGCCTGTAACCGACGGTATCGCGCCGGGGTCGGTTTTGAAGTTTGAATCTGTTTAg                      >  2:8038567/1‑76 (MQ=255)
                                                         gCACTGTCCATGTGGCAGGGCCCCTACCCGACGATATTCCGCCGGGGACGGGTTTGGAGTTTGTATCTGGTCAggg                    >  2:417760/1‑75 (MQ=255)
                                                         gCACTGCCCATGTGGCAGCCCCGGCCCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGa                    <  2:6769743/76‑1 (MQ=255)
                                                          cccTGTCCATGGGGCAGGCCCTGTCCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAg                   <  1:5168793/74‑1 (MQ=255)
                                                           aCTGTCCATGTGGCAGGGCCCGTACCCCCCGATATCGCGCCGGGGTCGGTTTTTTAGGTTTAATTTTTTAAGGAga                  >  2:2534007/1‑75 (MQ=255)
                                                            ctttCCATGTGGCAGCCCCTGTCCCCGACGATATCGCGCCGGGAGCGGTTTTGAAGTTTGAATCTGTTAAGGAGTa                 <  2:6462409/73‑1 (MQ=255)
                                                            cTGTCGATGGGGCAGGCCCTGTCCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTa                 <  1:2272231/76‑1 (MQ=255)
                                                            cTGTCCGTGTGGCCGCCCGTGCACCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTa                 <  2:2624090/76‑1 (MQ=255)
                                                            cTGTCCATGTGGCAGGGCCTGTTCCCGACGACATCGCCCCCGGGACGGGTTTTTAGTTTGAATCTCTTAAAGAGTa                 >  2:1623116/1‑76 (MQ=255)
                                                            cTGTCCATGTGGCAGGGCCTCTTCCCGACGCGATCACGCCGGGGTCGGGGTTGTAGTTTGAATCTCTTTAGGAgga                 >  2:3481428/1‑74 (MQ=255)
                                                            cTGTCCATGTGGCAGGCCCGGTACCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTa                 <  2:8366079/76‑1 (MQ=255)
                                                            cTGTCCAGGTGCCAGCCCCGGTACCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTa                 <  2:2216605/76‑1 (MQ=255)
                                                             tGTCCATGTGGCAGGCCCTGTCCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTaa                <  2:7998971/76‑1 (MQ=255)
                                                             tGTCCATGTGGCAGGCCCTGTCCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTaa                <  1:4464265/76‑1 (MQ=255)
                                                             tGTCCAGGTGGCAGCCCCGTTCCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTaa                <  2:3610422/76‑1 (MQ=255)
                                                               tCCATGTGGCAGGCCGTGTCCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTaaaa              <  1:6211957/76‑1 (MQ=255)
                                                               tCCATGGGGCAGGCCCTTTACCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTaaaa              <  1:6937565/76‑1 (MQ=255)
                                                               accATGGGGCGGGCCCTGTACCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTaaaa              <  1:6108893/75‑1 (MQ=255)
                                                                ctaTGTGCCAGGCTTTGTGCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTaaaaa             <  1:1913998/74‑1 (MQ=255)
                                                                ccATGTGGCAGGGCCCGTACCCGACGATATCGCGCCGGGGGCGGGTTTGAAGGTTGAATCTGTTTAGGAGGaaaaa             >  2:494786/1‑76 (MQ=255)
                                                                 cgtgtgGCAGGCCCTGTCCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTaaaaaa            <  2:7256058/74‑1 (MQ=255)
                                                                   tgtgGCAGGCCTTGTCCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTAAAAAATg          <  1:7439745/76‑1 (MQ=255)
                                                                   tgtgGCAGGCCCTGTCCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTAAAAAATg          <  1:1169451/76‑1 (MQ=255)
                                                                   tgtgGCAGGCCCTGTCCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTAAAAAATg          <  1:3475663/76‑1 (MQ=255)
                                                                       gCAGCCCCTGTCCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTAAAAAATGAATc      <  2:6375431/76‑1 (MQ=255)
                                                                        cAGGCCCTGTACCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTAAAAAATGAATCa     >  2:3494907/1‑76 (MQ=255)
                                                                          ggCCCTGTCCCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTAAAAAATGAATCAgg   <  1:3797787/76‑1 (MQ=255)
                                                                          ggCCCTGTACCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTAAAAAATGAATCAgg   <  2:6542672/76‑1 (MQ=255)
                                                                           gCCCTGTACCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTAAAAAATGAATCAGGt  >  2:4752392/1‑76 (MQ=255)
                                                                            |                                                                          
CGAACTGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTACCCGACGATATCGCGCCGGGATCGGTTTTGAAGTTTGAATCTGTTAAGGAGTAAAAAATGAATCAGGT  >  REL606/2722457‑2722606

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