Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123062,758,9090TG22.8% 709.9 / 22.1 263G33G (GGT→GGG) agaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (43/160);  new base G (60/0);  total (103/160)
Fisher's exact test for biased strand distribution p-value = 1.55e-31
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 340 total aligned reads displayed.

AGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAA  >  AP012306/2758761‑2759056
                                                                                                                                                     |                                                                                                                                                     
     aaCCGCAGCCCGTTGCCGGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTAcc                                                                                                                                                 >  2:2703358/1‑150 (MQ=255)
     aaCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTAcc                                                                                                                                                 >  1:3492089/1‑150 (MQ=255)
      aCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCCCCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGCGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCt                                                                                                                                                <  2:538232/150‑1 (MQ=255)
      aCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCt                                                                                                                                                <  1:3712057/150‑1 (MQ=255)
      aCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCt                                                                                                                                                <  1:2552639/150‑1 (MQ=255)
       ccGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTg                                                                                                                                               <  1:362409/150‑1 (MQ=255)
        cGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTATCTgg                                                                                                                                              >  2:1982257/1‑150 (MQ=255)
          cAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTgggg                                                                                                                                            >  2:2358927/1‑149 (MQ=255)
               cGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGt                                                                                                                                       <  1:1022835/150‑1 (MQ=255)
               cGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGt                                                                                                                                       <  2:3308529/150‑1 (MQ=255)
                gTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGtt                                                                                                                                      <  1:486404/150‑1 (MQ=255)
                 ttGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTg                                                                                                                                     <  2:1403571/150‑1 (MQ=255)
                  tGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTG‑c                                                                                                                                   >  2:1860501/1‑150 (MQ=255)
                   gCCAGTAAAAAGGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑cc                                                                                                                                  <  1:1212322/150‑1 (MQ=255)
                     cAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGGGATGTTG‑CCGa                                                                                                                                >  2:2786938/1‑150 (MQ=255)
                       gTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTACTGGGGGAAGTTGTCCGAt                                                                                                                               >  2:4554850/1‑150 (MQ=255)
                        tAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTACC‑GGGGGGTGTTG‑TCGGTcc                                                                                                                             >  2:4454538/1‑150 (MQ=255)
                        tAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATcc                                                                                                                             <  1:1636103/150‑1 (MQ=255)
                        tAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATcc                                                                                                                             >  1:1888637/1‑150 (MQ=255)
                        tAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCCTGGCGATGGTG‑CCGATcc                                                                                                                             >  1:3825956/1‑150 (MQ=255)
                           aaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGa                                                                                                                          >  1:3569551/1‑150 (MQ=255)
                            aaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGaa                                                                                                                         >  1:1528955/1‑150 (MQ=255)
                             aTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGaaa                                                                                                                        <  1:953642/150‑1 (MQ=255)
                             aTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGaaa                                                                                                                        <  1:1168573/150‑1 (MQ=255)
                             aTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTG‑CCGATCCTGaaa                                                                                                                        >  2:3201320/1‑150 (MQ=255)
                              tGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTGCCTGGGCCATGTTG‑CCGATCCTGaaaa                                                                                                                       >  1:2990170/1‑150 (MQ=255)
                              tGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGATGTTG‑CCGATCCTGaaaa                                                                                                                       >  2:1763653/1‑150 (MQ=255)
                                aagaagAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAg                                                                                                                     >  2:4703630/1‑150 (MQ=255)
                                 agaagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGa                                                                                                                    <  2:1005727/150‑1 (MQ=255)
                                 agaagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGa                                                                                                                    <  1:4797088/150‑1 (MQ=255)
                                   aagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGAtt                                                                                                                  <  2:2460171/150‑1 (MQ=255)
                                    agagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGAttt                                                                                                                 <  1:65562/150‑1 (MQ=255)
                                     gagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGAATta                                                                                                                >  1:2316193/1‑150 (MQ=255)
                                     gagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTta                                                                                                                >  2:2958539/1‑150 (MQ=255)
                                         aCTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAg                                                                                                            <  1:3331727/150‑1 (MQ=255)
                                         aCTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAg                                                                                                            <  2:3460920/150‑1 (MQ=255)
                                         aCTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAg                                                                                                            <  1:4611304/150‑1 (MQ=255)
                                         aCTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGCC‑GGTTTTTCCCGGGGGTTGTTG‑CCCATCCTGTAAAAGATTTATAAg                                                                                                            >  2:2073111/1‑150 (MQ=255)
                                           tACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGga                                                                                                          <  1:3151483/150‑1 (MQ=255)
                                            aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGG‑GGGTTTTCCTGGGCGGTGTTT‑GCGCACCCGAAAAAAATTTTTTAGgac                                                                                                         >  2:749499/1‑150 (MQ=255)
                                            aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTTGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGgac                                                                                                         >  1:1426396/1‑150 (MQ=255)
                                            aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGgac                                                                                                         >  1:55978/1‑150 (MQ=255)
                                              aGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGgacga                                                                                                       <  1:2749906/150‑1 (MQ=255)
                                               gCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGgacgac                                                                                                      <  1:4712830/150‑1 (MQ=255)
                                               gCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAAgacgac                                                                                                      >  1:1314291/1‑150 (MQ=255)
                                                    gggATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAAcc                                                                                                <  1:3801922/150‑1 (MQ=255)
                                                    gggATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCCTGGGGGTGTTT‑CCGATCCTGAAAAAAATTTATAAAGACGGC‑AAAcc                                                                                                >  1:763856/1‑150 (MQ=255)
                                                      gATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCgg                                                                                              <  2:406641/150‑1 (MQ=255)
                                                      gATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑CGGTTTTCCCGGGCGATGGTG‑CCGATCCTGAAAAAGAATTATCAGGACGAC‑AAACCgg                                                                                              >  1:1820472/1‑150 (MQ=255)
                                                          tGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAAGAGGAC‑AAACCGGGCtt                                                                                          >  2:3819119/1‑150 (MQ=255)
                                                           gAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTa                                                                                         <  2:1199839/150‑1 (MQ=255)
                                                           gAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTa                                                                                         >  1:1991451/1‑150 (MQ=255)
                                                             aaTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAg                                                                                       <  2:625180/150‑1 (MQ=255)
                                                                cAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGGTGTTG‑CCGCTCCCGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAacgc                                                                                    >  2:1903221/1‑150 (MQ=255)
                                                                cAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGCTTTATCCGGGCGGTGTTG‑GCGAACCTGAAAAAGATTTATAAAGGCGAA‑CAACCGGGGGTAAacgc                                                                                    >  2:2167071/1‑150 (MQ=255)
                                                                    aaaaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTTTAAGGACGAC‑AAACCGGGCTTAAGCgcggcg                                                                                >  2:3370893/1‑150 (MQ=255)
                                                                      aaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGAt                                                                              <  1:3388888/150‑1 (MQ=255)
                                                                          gATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGATGTTG‑CCGATCCTGGAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGaaa                                                                          >  2:590999/1‑150 (MQ=255)
                                                                            tatCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGGTGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGGTTAAACGCGGCGGTGaaaaa                                                                        >  2:1116802/1‑148 (MQ=255)
                                                                             atCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGAt                                                                       <  1:360833/150‑1 (MQ=255)
                                                                               cACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTTTCCGGGGGATGTTG‑CCGGGCCCCGAAAAGATTTTTAAGGACGAC‑AAACCCGGCTTAAGCGCGGCGCGGAAAAATaa                                                                     >  2:1281081/1‑150 (MQ=255)
                                                                               cACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTACTCGGCGATGTTG‑CCGATCCTGAAAAAAATTTATAAAGACGAC‑AAAACCGGCTTAAACGCGGCGATGAAAAATaa                                                                     >  2:313525/1‑150 (MQ=255)
                                                                                ccccGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAAc                                                                    <  2:3295043/149‑1 (MQ=255)
                                                                                  ccGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCt                                                                  <  1:2041905/150‑1 (MQ=255)
                                                                                   cGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTGCCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTc                                                                 <  2:4608562/150‑1 (MQ=255)
                                                                                   cGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTc                                                                 <  2:1664834/150‑1 (MQ=255)
                                                                                    gTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTTCCTCGGCGATGCTG‑CCCATCCTGAAAAAAATTTATTAAGACCGC‑AAAACCGGCTTAAGCGCGGGGATGAAAGATaaacctg                                                                >  2:4279428/1‑145 (MQ=255)
                                                                                    gTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCg                                                                <  2:4653668/150‑1 (MQ=255)
                                                                                      cTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGaa                                                              >  1:4035750/1‑150 (MQ=255)
                                                                                       tGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAAt                                                             <  1:530154/150‑1 (MQ=255)
                                                                                         ggCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGGTGAAAGATAACCTCGAAttt                                                           >  1:3229695/1‑150 (MQ=255)
                                                                                           cTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATttat                                                         <  1:1517672/150‑1 (MQ=255)
                                                                                            tttCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATttatt                                                        <  1:2854176/150‑1 (MQ=255)
                                                                                               cgttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAAt                                                     <  1:4923569/150‑1 (MQ=255)
                                                                                                 ttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATac                                                   <  2:2643738/150‑1 (MQ=255)
                                                                                                   cgtcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATaccc                                                 <  1:2243097/150‑1 (MQ=255)
                                                                                                      cgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATacccacc                                              <  2:4513312/150‑1 (MQ=255)
                                                                                                              cAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTg                                      <  1:2519922/150‑1 (MQ=255)
                                                                                                              cAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTTGGGGGTGTTT‑CCGGTCCTGAAAAAGATTTTTAAGGACGAC‑AAAACGGGCTTTAACCCGGGGTTGTAAGATAAACTCGAATTTATTAATACCCACCCGAAACTg                                      >  2:4762492/1‑150 (MQ=255)
                                                                                                                  cGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCg                                  >  1:3954685/1‑150 (MQ=255)
                                                                                                                     gCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCCCCCGAATCTGGTCGGAt                               <  2:88536/150‑1 (MQ=255)
                                                                                                                     gCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGAt                               <  1:3048926/150‑1 (MQ=255)
                                                                                                                      cTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGAtt                              >  1:4354703/1‑150 (MQ=255)
                                                                                                                      cTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGGTGTTG‑GCGGTCCCGAAAAAAATTTTTAAAGACGACAAAACCGGGCTAAAG‑GCGGCGGTGGAAAATAAACTCGGATTTATTTATACCCACCCCCATCTGGTccggtt                              >  2:4979656/1‑145 (MQ=255)
                                                                                                                       tttAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTc                             <  1:4500090/150‑1 (MQ=255)
                                                                                                                       ttaaaCTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAAACTCGAATTTTTTAATACCCACCCGAATATGGTCGGGTTc                             >  2:3877795/4‑150 (MQ=255)
                                                                                                                           aCTACGAAAGGATGCAGGCGGGG‑GGGTTTTGCTGGGGGGTGGTG‑GGGGGGGTGGAAAAGAGTTATAAGGGGGAG‑AAAACGGGGGTAAGGGGGGGGGTGGAAGATAAACTCGGATTTATTAATAACCACCCGAATCTGGGCGGATTCCTGa                         >  2:1560157/1‑150 (MQ=255)
                                                                                                                             tACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATg                       <  2:3778080/150‑1 (MQ=255)
                                                                                                                             tACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATg                       <  1:1274986/150‑1 (MQ=255)
                                                                                                                              aCGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATgg                      <  1:1757557/150‑1 (MQ=255)
                                                                                                                              aCGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATgg                      <  1:2175325/150‑1 (MQ=255)
                                                                                                                                gAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGATGTTG‑CCGGTCCTGAAAAAGATTTATAAGGGCGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATgggg                    >  2:4168503/1‑150 (MQ=255)
                                                                                                                                gAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGGGGAGGTG‑GGGGTCCTGAAAAAGATTTATAAGGACGGG‑AAACCGGGGGTAAGCGGGGGGGTGGAAGATAAACACGCATTTATTTATAGCCCCCCGCATCGGGTTGGGTTCCTGgggggg                    >  2:2531205/1‑150 (MQ=255)
                                                                                                                                  aaGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGtt                  <  2:2222679/150‑1 (MQ=255)
                                                                                                                                         cAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAAtt           >  1:4266175/1‑150 (MQ=255)
                                                                                                                                          aGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGttatctttca          >  1:1669433/1‑144 (MQ=255)
                                                                                                                                             cgggc‑ggTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGa       >  2:3627447/1‑150 (MQ=255)
                                                                                                                                                 c‑ggTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGga   <  2:3481125/150‑1 (MQ=255)
                                                                                                                                                 c‑ggTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGga   <  2:258421/150‑1 (MQ=255)
                                                                                                                                                   ggTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaa  >  1:4128220/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTG‑CCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAA  >  AP012306/2758761‑2759056

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑