Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,437,4060GT48.4% 24.6 / 26.1 32D435E (GAC→GAAypdDfused predicted PTS enzymes: Hpr component/enzyme I component/enzyme IIA component
Reads supporting (aligned to +/- strand):  ref base G (12/4);  new base T (15/0);  total (28/4)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTCTTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCAT  >  REL606/2437331‑2437480
                                                                            |                                                                               
ctctTCCTGGGGAATATTCAGATAAGGAATTCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTg                                                                                 <  2:6939492/76‑1 (MQ=255)
ctctTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTg                                                                                 <  2:4855273/76‑1 (MQ=255)
ctctTCCTGGGGAATATTCAGAAAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTg                                                                                 <  2:2981453/76‑1 (MQ=255)
  ctTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTc                                                                               >  1:5230803/1‑76 (MQ=255)
     cgtGGGGAATATTAAGATAAGGAATCCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcc                                                                           <  2:4136008/74‑1 (MQ=255)
     ccTGGGGAATATTCAGATAAGGAATTCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcc                                                                           <  2:3684021/76‑1 (MQ=255)
     ccTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgcc                                                                           >  1:802715/1‑76 (MQ=255)
     ccTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtgcc                                                                           >  2:6286440/1‑76 (MQ=255)
      cTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTctccc                                                                          >  1:861608/1‑76 (MQ=255)
      cTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTcgccc                                                                          >  1:1065278/1‑76 (MQ=255)
      cTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgccc                                                                          <  2:2970035/76‑1 (MQ=255)
      cTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgccc                                                                          >  1:4244807/1‑76 (MQ=255)
       tGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTcgcccg                                                                         >  2:3785166/1‑76 (MQ=255)
        ggggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTtggccgc                                                                        >  2:432895/1‑76 (MQ=255)
         gggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTttcccgcc                                                                       >  1:4085262/1‑76 (MQ=255)
          ggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTcgctcgccc                                                                      >  2:7738527/1‑75 (MQ=255)
          ggAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTctccccccc                                                                      >  2:3606529/1‑75 (MQ=255)
            aaTATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCgccgcc                                                                    >  1:4411618/1‑76 (MQ=255)
            aaTATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCgccccc                                                                    >  1:1154516/1‑76 (MQ=255)
            aaTATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTTGCCCgccgcc                                                                    >  2:3829015/1‑76 (MQ=255)
            aaTATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTTGCCCgccccc                                                                    >  2:3390915/1‑76 (MQ=255)
             atatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCgccgccc                                                                   >  1:3221896/1‑75 (MQ=255)
             atatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCgcccccc                                                                   >  2:6743310/1‑75 (MQ=255)
             atatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGACCgccgccc                                                                   >  2:1303048/1‑75 (MQ=255)
              tatTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCGCGCCac                                                                  >  1:7198112/1‑75 (MQ=255)
               atTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCgccccccgc                                                                 >  1:5103348/1‑73 (MQ=255)
               atTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCCCCact                                                                 >  2:2414322/1‑74 (MQ=255)
               atTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCact                                                                 >  1:7377331/1‑74 (MQ=255)
               atTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCgccccccgg                                                                 >  1:2122470/1‑73 (MQ=255)
               atTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAgc                                                                 >  1:4633334/1‑75 (MQ=255)
                 taaGAAAAGGATTCTTTTTCGCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTAg                                                               <  1:7221075/74‑1 (MQ=255)
                  cAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCGCCCGCAGc                                                              >  2:2435777/1‑76 (MQ=255)
                  cAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCCGGAGc                                                              >  1:6787905/1‑76 (MQ=255)
                  cAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTACCACCCCGCCAGTAGc                                                              >  2:6345240/1‑76 (MQ=255)
                   aGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCCCCACTAGcc                                                             >  1:6900124/1‑75 (MQ=255)
                    gATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCGCCCGTAGCAc                                                            >  1:2295846/1‑76 (MQ=255)
                    gATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCCCCAGGAGGAc                                                            >  1:1865486/1‑76 (MQ=255)
                     aTAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCCCCAgcagcaac                                                           >  1:6249902/1‑74 (MQ=255)
                     aTAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCAGTAGCAcc                                                           >  2:923797/1‑76 (MQ=255)
                     aTAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTATTCCCCGCCCCCAGCAGCAcc                                                           >  2:6816830/1‑76 (MQ=255)
                       aaGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCCGTCGCAcccg                                                         >  2:4036420/1‑74 (MQ=255)
                       aaGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCCGGCGCACctc                                                         >  2:3740122/1‑75 (MQ=255)
                       aaGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCAGTAGCAcccc                                                         >  1:5560622/1‑74 (MQ=255)
                        aGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTCGCCCCCGCCCGCAGCACccgc                                                        >  1:653282/1‑76 (MQ=255)
                        aGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCTCCGCCAGTAGCACccgc                                                        >  1:3702052/1‑76 (MQ=255)
                         ggAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCCCCCGTCGCCCcagcg                                                       >  1:849374/1‑75 (MQ=255)
                         ggAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCACGAGCACccgct                                                      >  1:3489022/1‑76 (MQ=255)
                          gAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCTCCCGCCTCCCGTAGGACccgctg                                                     >  1:1876553/1‑76 (MQ=255)
                          gAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCAGTAGCACccgctg                                                     >  2:1962962/1‑76 (MQ=255)
                          gAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCCGGCGCACctgctg                                                     >  2:6148581/1‑76 (MQ=255)
                            aTGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTGCCCCCCCCCCCGTAGCACCCTCTggg                                                   >  1:4641934/1‑75 (MQ=255)
                              gCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCTCCCGTAGCACCCCCTGGGAg                                                >  2:3960745/1‑76 (MQ=255)
                               cTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCCGCCCCCAGAGCACCTTCTGGGAgg                                               >  1:4561425/1‑76 (MQ=255)
                                ttttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCCCCCGCAGCACCCGCTGGGAggg                                              >  1:6535314/1‑75 (MQ=255)
                                ttttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTGGCCCCCCCCCCGAAACACCCGCTGGTAggg                                              >  2:1289268/1‑75 (MQ=255)
                                ttttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTATCGCCCGCCGCCAGTAGCACCTGCTGGTAggg                                              >  2:2659368/1‑75 (MQ=255)
                                 tttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGGCCGCCGCCCGTAGCACCCGCTGGGAgggt                                             >  2:2101671/1‑74 (MQ=255)
                                 tttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCCCCAGCAACACCCGATGGTggggt                                             >  2:3129068/1‑74 (MQ=255)
                                 tttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCGCCAGCAGGACCCGCTGGGAgggt                                             >  1:847603/1‑74 (MQ=255)
                                 tttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCCCCACGAGCACCCCCTTGTAgggg                                             >  2:1356820/1‑74 (MQ=255)
                                 tttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCCGTAGCACCTCCTGGGAgggt                                             >  1:6057004/1‑74 (MQ=255)
                                 tttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTTTCCCCCGCCCCCCGTAGCCCCTCCTGGTAGGCt                                             >  2:7845402/1‑76 (MQ=255)
                                 tttATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGGCCGCCGCCCGGAGCACCTGCCGGTAGGCt                                             >  2:4655609/1‑76 (MQ=255)
                                     cccccgccgTTGCCCATCGTGCGGAAGAAAAGCGGCTTGCCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCaa                                       <  1:832298/73‑1 (MQ=255)
                                       gccgccGATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCAGCCGCCAGCAGCACCTCCTGGGAGGCTTCcaaa                                     >  2:778281/1‑75 (MQ=255)
                                       cccgcggATTCCCATGGGGCGAAAGAAATGCGGCTTGCCCCCCCCCGCCAGTAGCACCTGCTGGTAGGCTTCAAAt                                     <  2:6026411/75‑1 (MQ=255)
                                         cgccgATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCAGTAGCACCTGCTGGTAGGCGTCaaaatt                                   >  2:3692950/1‑73 (MQ=255)
                                         cgccgATGTCCATCGGGCGGAAGATAATCGGCTGGCCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTg                                   <  1:1904263/76‑1 (MQ=255)
                                          gccgATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCCGTAGCAcc                                                           >  2:221270/1‑55 (MQ=255)
                                          gccgATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTAGCAcc                                                           <  1:221270/55‑1 (MQ=255)
                                           ccgATGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCGCCCCCACTAGCACCTGCTGGGAGGCTTCCAAATTgt                                  >  2:6451927/1‑75 (MQ=255)
                                            cgATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCCGTAGCACCTGCTGGTAGGCTTCAAATTGCTc                                >  1:7846354/1‑76 (MQ=255)
                                             gATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCCGGAGCACCTGCTGGGAGGGGTCTAAATTGCTc                                >  2:1479445/1‑76 (MQ=255)
                                             gATGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCCGTATCCCCTGCTGGGAGGCTTCAAAATGCTcc                               >  2:1711396/1‑76 (MQ=255)
                                              aTGTCCATCGTGCGGAAGATAATCGGCTTTTTGCCCCCCGCCCCGAGAACCTGCTGGTGGGGTTCCAAATGTTCCt                              >  2:1507222/1‑76 (MQ=255)
                                              aTGTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCCCCCGTAGCACCCGCTGGGAGGGCTCAAAATTCTccc                              >  2:5013873/1‑75 (MQ=255)
                                                gcccATCGTGCGAAAGATAATCGTCTGGCCCCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGc                            <  1:1930180/74‑1 (MQ=255)
                                                gTCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCGCCCGTAGCACCCGCCGGGAGGGTTCTAAATGCTCCCGc                            >  1:476764/1‑76 (MQ=255)
                                                gTCCATCGTGCGGAAGATAATCGGCTTTTCGACCGCCCCCAGTAGCACCTCCCGGTAGGCTTCCAAATACTccctc                            >  2:6668300/1‑73 (MQ=255)
                                                 tCCATCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCCGAAGGACCCGCTGGTAGGCGTCAAATTTCTCCCGcc                           >  1:3249677/1‑75 (MQ=255)
                                                  ccATGGTGCGGAAGATAACCGGCTTGCCCCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTc                          <  1:279163/76‑1 (MQ=255)
                                                  ccATCGTGGGGAAGATATTCGGCTTGCCCCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTc                          <  1:2522153/76‑1 (MQ=255)
                                                   cATCGTGCGGAAGATAATCGGCTTTTCGCCCGCCGCCAGTAGCACCTGCTGGGAGGGTTTCAAATTctcctcctcc                         >  2:7080195/1‑75 (MQ=255)
                                                    accGTGGGGAAAATAATGGGCTGGCCCCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGt                        <  1:3253316/74‑1 (MQ=255)
                                                    aTCGTGCGGAAGATAATCGGCTTTTCGCCCCCCCCCAGTAGCACCCGCTGGGGGGCTTCAAATTTGTCCTGCTcct                        >  1:7427969/1‑74 (MQ=255)
                                                        tGCGGAAGATAATCGGCTTGTTGCCCCCCCCCCGTAGCCCCCGCTTGGTGGGCTTAAATTTCTTCTCCTCCTCCgg                    >  2:4306349/1‑76 (MQ=255)
                                                         gCGGAAGATAATCGGCTTTTCGCCCGCCCCCAGTAGCACCTGCTGGGAGGCTTCAAATTGCTCCTGCTCGTCCGgc                   >  2:5779609/1‑76 (MQ=255)
                                                                gataTCGGCGTTGCCCCCCCCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTg             <  2:1395039/73‑1 (MQ=255)
                                                                    gcGGCTTGCCCCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTcgcg        <  2:5890574/75‑1 (MQ=255)
                                                                        ttttGCCGCCCCCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATc     <  1:6881439/74‑1 (MQ=255)
                                                                        cTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATcc    >  2:2989911/1‑76 (MQ=255)
                                                                        cTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATcc    >  2:31549/1‑76 (MQ=255)
                                                                          ttGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCa   >  2:3637597/1‑76 (MQ=255)
                                                                          ttGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCa   >  1:2075194/1‑76 (MQ=255)
                                                                          ttGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCa   >  1:3722683/1‑76 (MQ=255)
                                                                          ttGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTAGCGATCCa   >  2:1639510/1‑76 (MQ=255)
                                                                          ttGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCCGTCGCGATCCa   >  2:4509242/1‑76 (MQ=255)
                                                                          tggTCCCCCCCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCa   <  2:487045/74‑1 (MQ=255)
                                                                           tGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCat  >  2:6961853/1‑76 (MQ=255)
                                                                           tGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCat  >  1:455042/1‑76 (MQ=255)
                                                                            |                                                                               
CTCTTCCTGGGGAATATTCAGATAAGGAATGCTTTTATCGCCGCCGATGTCCATCGTGCGGAAGATAATCGGCTTGTCGCCCGCCGCCAGTAGCACCTGCTGGTAGGCTTCAAATTGCTCCTGCTCGTCCGGCGCGCTGTCGCGATCCAT  >  REL606/2437331‑2437480

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: