Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123062,758,9090TG25.6% 671.4 / 21.3 274G33G (GGT→GGGagaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (54/150);  new base G (70/0);  total (124/150)
Fisher's exact test for biased strand distribution p-value = 2.36e-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 359 total aligned reads displayed.

AGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGA  >  AP012306/2758761‑2759058
                                                                                                                                                    |                                                                                                                                                       
  cgcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGtttt                                                                                                                                                      <  1:1389595/150‑1 (MQ=255)
  cgcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGtttt                                                                                                                                                      <  2:1754076/150‑1 (MQ=255)
     aaCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTAcc                                                                                                                                                   <  1:212655/150‑1 (MQ=255)
        cGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCCgg                                                                                                                                                >  2:3210914/1‑150 (MQ=255)
         gCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTggg                                                                                                                                               >  1:3989427/1‑150 (MQ=255)
           aGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCg                                                                                                                                             <  2:2861743/150‑1 (MQ=255)
            gCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGa                                                                                                                                            >  1:1495808/1‑150 (MQ=255)
              ccGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGGGGGTTTATCTGGGGGATg                                                                                                                                          >  2:687987/1‑150 (MQ=255)
              ccGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTTTCTGGGCGCTg                                                                                                                                          >  2:4101443/1‑150 (MQ=255)
                gTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGtt                                                                                                                                        <  2:4756786/150‑1 (MQ=255)
                 ttGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTg                                                                                                                                       <  1:4507816/150‑1 (MQ=255)
                 ttGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTg                                                                                                                                       <  1:1591965/150‑1 (MQ=255)
                 gtGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTg                                                                                                                                       <  2:3999753/149‑1 (MQ=255)
                    ccAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCg                                                                                                                                    <  2:2393593/150‑1 (MQ=255)
                      aGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGAt                                                                                                                                  >  1:2776616/1‑150 (MQ=255)
                      aGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGAt                                                                                                                                  >  1:3560531/1‑150 (MQ=255)
                        tAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCTGGGCGGTGTTGCCGATcc                                                                                                                                >  1:1065236/1‑150 (MQ=255)
                        tAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTACTGGGCGAGGTTTCCGATcc                                                                                                                                >  2:1252404/1‑150 (MQ=255)
                         aaaaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCTAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCt                                                                                                                               <  2:3595831/150‑1 (MQ=255)
                           aaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGa                                                                                                                            >  1:4322817/1‑150 (MQ=255)
                            aaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGaa                                                                                                                           <  1:3899607/150‑1 (MQ=255)
                            aaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGaa                                                                                                                           >  2:4424154/1‑150 (MQ=255)
                             aTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGaaa                                                                                                                          <  2:4175482/150‑1 (MQ=255)
                              tGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGaaaa                                                                                                                         <  1:4057494/150‑1 (MQ=255)
                              tGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTGTTAACTGGGGGATGTTTCCGATCCTGaaaa                                                                                                                         >  1:4184433/1‑150 (MQ=255)
                                aagaagAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCTTGGCGGTGTTGCCGGTCCTGAAAAAg                                                                                                                       >  2:4786714/1‑150 (MQ=255)
                                aagaagAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCCGTGATTCCGAGCCTGaaaaaa                                                                                                                       >  2:3254646/1‑149 (MQ=255)
                                  gaagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAAAt                                                                                                                     >  2:3867355/1‑150 (MQ=255)
                                   aagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGAtt                                                                                                                    <  2:1963793/150‑1 (MQ=255)
                                    agagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTTCCTTGGCGATGTTTGCGATCCTGGAAAAAAttt                                                                                                                   >  1:939522/1‑150 (MQ=255)
                                     gagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTta                                                                                                                  >  1:3373453/1‑150 (MQ=255)
                                      agGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTtat                                                                                                                 >  1:3428323/1‑150 (MQ=255)
                                      agGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTtat                                                                                                                 <  1:980935/150‑1 (MQ=255)
                                      agGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCTGGGCGATGTTTGCGATCCTGAAAAAGATTtat                                                                                                                 >  2:3787593/1‑150 (MQ=255)
                                       ggACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTtata                                                                                                                <  1:1527460/150‑1 (MQ=255)
                                       ggACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCACTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCTCGGCCATGCTGTCCCTCCCGAAAAAGATTtata                                                                                                                >  1:2065251/1‑150 (MQ=255)
                                            aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCGGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGgac                                                                                                           >  1:4584583/1‑150 (MQ=255)
                                            aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTTCCTTGGCGGTGGTTCCGGACCTGGAAAAAGGttttaaagac                                                                                                            >  2:1147673/1‑143 (MQ=255)
                                            aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGgac                                                                                                           <  2:3525113/150‑1 (MQ=255)
                                             cAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCTGGGCGATGGTGCCGGTCCTGAAAAAAATTTATAAGgacg                                                                                                          >  2:541351/1‑150 (MQ=255)
                                             cAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGCTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGTGGTTTTTCCCTTGCGAAGTTTCCGATCCTGTATAAGGTTTATAAGgact                                                                                                          >  1:3238125/1‑149 (MQ=255)
                                              aGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAAgacga                                                                                                         >  1:2744317/1‑150 (MQ=255)
                                                cAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACa                                                                                                      <  1:1668592/150‑1 (MQ=255)
                                                cAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACa                                                                                                      <  1:1766684/150‑1 (MQ=255)
                                                cAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAGATTTTTTAAGACGACa                                                                                                      >  2:2474826/1‑150 (MQ=255)
                                                 aaTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTACTGGGGGGTGTTGGCGGGACTGAAAAAGAATTATAAGGACGGCaa                                                                                                     >  2:589369/1‑150 (MQ=255)
                                                  aTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGTCGGGTTTTCCTTGGCCGAGTTGCCCATACTGTAAAAAATTTTTAAAGACGACaaa                                                                                                    >  2:2062250/1‑150 (MQ=255)
                                                  aTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACaaa                                                                                                    <  1:2271439/150‑1 (MQ=255)
                                                   tGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAAc                                                                                                   <  2:3041455/150‑1 (MQ=255)
                                                    gggATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAAcc                                                                                                  <  1:478147/150‑1 (MQ=255)
                                                    gggATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAAcc                                                                                                  >  2:3744556/1‑150 (MQ=255)
                                                      gATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCACTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGAGGGTTTACCTGGGCGGAGTTTCCGATCCCTGAAAAGAGTTTTAAAGGCGCCaaaaccg                                                                                                >  2:4602266/1‑146 (MQ=255)
                                                       aTCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCggg                                                                                               >  2:4372201/1‑150 (MQ=255)
                                                       aTCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAAGACGACAAACCcgg                                                                                               >  1:2899211/1‑150 (MQ=255)
                                                         cTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTTCTTGGGGGTGTTTCCGATCCCGAAAAAAAATTTTAAAGGCGACAAAACCgggt                                                                                              >  2:3645900/1‑149 (MQ=255)
                                                          tGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCtt                                                                                            <  1:4241569/150‑1 (MQ=255)
                                                          tGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTAACTCGGCCAGGGTGGCGATTCTGCAAAAGATTTATAAAGACGGAAAACCCGGcct                                                                                            >  1:889834/1‑148 (MQ=255)
                                                           gAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGGGATGTTGCCGGTCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTa                                                                                           >  2:1468255/1‑150 (MQ=255)
                                                            aaaTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTaa                                                                                          <  1:1295090/150‑1 (MQ=255)
                                                             aaTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAg                                                                                         <  2:1687688/150‑1 (MQ=255)
                                                                  gTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCgcgg                                                                                    <  1:2464319/150‑1 (MQ=255)
                                                                  gTAAAAAAGATAGCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCgcgg                                                                                    <  1:1917795/150‑1 (MQ=255)
                                                                    aaaaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCgcggcg                                                                                  <  2:4474576/150‑1 (MQ=255)
                                                                     aaaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGa                                                                                 >  1:4672616/1‑150 (MQ=255)
                                                                       aaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTTCCGATCCTTTAAAAGATTTATAAAGGAGACAAACCCGGCTTAAGCGCGGCGATg                                                                               >  2:2129195/1‑150 (MQ=255)
                                                                        aaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTTTAAGGACGACAAAACGGGCTTAAGCGCGGCGGTGa                                                                              >  2:3674284/1‑150 (MQ=255)
                                                                           atatCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAg                                                                           <  1:543102/150‑1 (MQ=255)
                                                                             atCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGAt                                                                         <  1:4136084/150‑1 (MQ=255)
                                                                                 cccGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAAcc                                                                     >  1:3938680/1‑150 (MQ=255)
                                                                                   cGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCTGGGCGGTGTTGCCGGTCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGGTGAAAGAGAACCTc                                                                   >  2:2190450/1‑150 (MQ=255)
                                                                                   cGTCTGGGCTTTCGTTCGACGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCCAACCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGTAAGATAACCTc                                                                   >  2:4565229/1‑150 (MQ=255)
                                                                                     tCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGAAAAACTcca                                                                 >  2:2844657/1‑148 (MQ=255)
                                                                                       tGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTCGGCGATTTTGCCGATCCTGAAAAAAATTTATAAGGACGACAAACCCGGCTTAAACGCGGCGATGAAAAATAACCTCCAAt                                                               >  2:1549952/1‑150 (MQ=255)
                                                                                       tGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAAt                                                               >  1:3132405/1‑150 (MQ=255)
                                                                                        gggCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAAtt                                                              <  1:2631409/150‑1 (MQ=255)
                                                                                          gCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATtta                                                            >  1:1817447/1‑150 (MQ=255)
                                                                                           cTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGGAAGATAAACTCGAATttat                                                           >  2:2735307/1‑150 (MQ=255)
                                                                                              tcgttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTaa                                                        <  1:4700784/150‑1 (MQ=255)
                                                                                                    gtcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAAATTTATAAAGACGACAAACCGGGCTTAAACGCGGCGATGAAAAATAAACTCGAATTTATTAATaacca                                                  >  2:2648182/1‑150 (MQ=255)
                                                                                                           ctgcAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTTCCTGGGGGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAAACGGGGTTAAGCGCGGGGATGAAAGATAAACTCGAATTTATTAATACCCACCCGAAt                                           >  2:827262/1‑150 (MQ=255)
                                                                                                           ctgcAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAAt                                           <  1:3792435/150‑1 (MQ=255)
                                                                                                             gcAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATCCCCACCCGAATCt                                         <  2:1399065/150‑1 (MQ=255)
                                                                                                             gcAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCt                                         >  1:2654738/1‑150 (MQ=255)
                                                                                                                    aGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCTGGGGGATGTTGCCCGTTCTGGAAAAGATTTATAAGGGCGACAAACCGGGGTTTAGCGCGGGGATGAAAGATAAACTCGAATTTTTTAATAACCACCCGGATCTGGGCGGa                                  >  2:1030795/1‑150 (MQ=255)
                                                                                                                       tttAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTc                               <  2:918811/150‑1 (MQ=255)
                                                                                                                       tttAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTc                               <  2:2301966/150‑1 (MQ=255)
                                                                                                                        ttAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTcc                              <  1:2398726/150‑1 (MQ=255)
                                                                                                                          aaCTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTTCTg                            >  2:3930961/1‑150 (MQ=255)
                                                                                                                            cTACGAAAGGATGCAGGCGGGCGGTTTTTCCCGGGCGGTGTTGCCGATCCTGAAAAAGATTTTTAAGGACGACAAAACCGGCTTAAGCGCGGCGATGAAAGATAAACTCGGATTTATTAATACCCACCCCAATCTGGGCGGATTCCTGAt                          >  2:2935423/1‑150 (MQ=255)
                                                                                                                            cTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGAt                          <  1:1844655/150‑1 (MQ=255)
                                                                                                                            cTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGAt                          <  1:4863987/150‑1 (MQ=255)
                                                                                                                                         cAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGGTGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAAtt             >  2:2251106/1‑150 (MQ=255)
                                                                                                                                          aGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAAttt            >  2:1997544/1‑150 (MQ=255)
                                                                                                                                               ggcggTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATg       <  1:3320226/150‑1 (MQ=255)
                                                                                                                                                gcggTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATgg      <  1:1160735/150‑1 (MQ=255)
                                                                                                                                                gcggTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATgg      <  1:695594/150‑1 (MQ=255)
                                                                                                                                                 cggTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGGTTAAGCGCGGCGATGGAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGga     <  1:3666129/150‑1 (MQ=255)
                                                                                                                                                  ggTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaa    >  1:2888151/1‑150 (MQ=255)
                                                                                                                                                   gTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaag   >  2:947047/1‑150 (MQ=255)
                                                                                                                                                    ttttACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaaga  >  2:3385954/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
AGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGA  >  AP012306/2758761‑2759058

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