Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123062,758,9090TG29.9% 509.7 / 27.6 217G33G (GGT→GGG) agaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (48/104);  new base G (65/0);  total (113/104)
Fisher's exact test for biased strand distribution p-value = 1.14e-24
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 294 total aligned reads displayed.

AAGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGA  >  AP012306/2758760‑2759058
                                                                                                                                                      |                                                                                                                                                      
    gcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTa                                                                                                                                                    <  1:2939900/150‑1 (MQ=255)
     cAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTAc                                                                                                                                                   <  1:4184334/150‑1 (MQ=255)
      aaCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTTcc                                                                                                                                                  >  2:2859820/1‑150 (MQ=255)
      aaCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTTcc                                                                                                                                                  >  2:2998504/1‑150 (MQ=255)
              cccGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTTCCTGGGCGAt                                                                                                                                          >  1:2423621/1‑150 (MQ=255)
                cGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGTAGGC‑GGGCGGTTTTACCTGGGCGATGt                                                                                                                                        <  2:73296/150‑1 (MQ=255)
                cGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGGGATGt                                                                                                                                        >  2:4049267/1‑150 (MQ=255)
                 gTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGCGATGtt                                                                                                                                       >  2:146154/1‑150 (MQ=255)
                  ttGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTg                                                                                                                                      <  2:3668047/150‑1 (MQ=255)
                      cAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTTCCTGGGCGATGTTGCCGa                                                                                                                                  >  1:2713616/1‑150 (MQ=255)
                       aGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGAt                                                                                                                                 <  2:4260753/150‑1 (MQ=255)
                         tAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATcc                                                                                                                               <  2:1741350/150‑1 (MQ=255)
                         tAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGCGATGTTGCCGATcc                                                                                                                               >  1:273175/1‑150 (MQ=255)
                         tAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCCGGGGGATGTTTCCGGTcc                                                                                                                               >  2:306511/1‑150 (MQ=255)
                          aaaaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCt                                                                                                                              >  1:71528/1‑150 (MQ=255)
                          aaaaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCt                                                                                                                              >  1:1892780/1‑150 (MQ=255)
                           aaaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑g                                                                                                                            <  1:772400/150‑1 (MQ=255)
                            aaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTTCCTGGGCGATGTTGCCGCTCCT‑Ga                                                                                                                           >  1:1768623/1‑150 (MQ=255)
                              aTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTTCCCGGGCGATGTTTCCGATCCT‑Gaaa                                                                                                                         >  2:1656705/1‑150 (MQ=255)
                               tGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGGTGTTGGCGATCCT‑Gaaaa                                                                                                                        >  2:2223885/1‑150 (MQ=255)
                                gaagaaGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑Gaaaaa                                                                                                                       <  1:2298870/150‑1 (MQ=255)
                                 aagaagAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTAACTGGGCGATGTTGGCGAGCCT‑GAAAAAg                                                                                                                      >  2:2867876/1‑150 (MQ=255)
                                  agaagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGa                                                                                                                     <  1:1314521/150‑1 (MQ=255)
                                  agaagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGa                                                                                                                     <  2:1515205/150‑1 (MQ=255)
                                   gaagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGAt                                                                                                                    <  2:511528/150‑1 (MQ=255)
                                    aagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTTCCTTGGGGAATGTGCCGATCCT‑GTAAAAGAtt                                                                                                                   >  2:708890/1‑150 (MQ=255)
                                     agagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGGGGGTTTATCTCGGGGCTGTTTGCGCTCCT‑GGAAAAGAttt                                                                                                                  >  2:1587617/1‑150 (MQ=255)
                                     agagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTTGGCGATGTTTCCGATCCT‑GAAAAAGAatt                                                                                                                  >  1:304806/1‑150 (MQ=255)
                                       agGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTTCCTGGGCGATGTTGCCGGTCCT‑GAAAAAGATTtat                                                                                                                >  1:2005707/1‑150 (MQ=255)
                                        ggACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGCGATGTTTCCGATCCT‑GGAAAAGATTtata                                                                                                               >  1:1010126/1‑150 (MQ=255)
                                         gACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATaa                                                                                                              <  2:3390468/150‑1 (MQ=255)
                                         gACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTTCCTTGGCGATGTTTCCGGTCCTGGAAAAAAATTttta                                                                                                               >  2:1960015/1‑150 (MQ=255)
                                         gACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATaa                                                                                                              >  2:768070/1‑150 (MQ=255)
                                          aCTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAg                                                                                                             >  2:2359687/1‑150 (MQ=255)
                                            tACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGga                                                                                                           <  1:927997/150‑1 (MQ=255)
                                             aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGGTCCT‑GAAAAAGATTTTTAAgggc                                                                                                          >  1:1942868/1‑148 (MQ=255)
                                             aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGgac                                                                                                          >  1:386943/1‑150 (MQ=255)
                                              cAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGgacg                                                                                                         <  1:2828461/150‑1 (MQ=255)
                                               aGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGGGGGTTTTTCTTGGGGGAGGTTTGGGTTCTT‑TAAAAAAATTTTTAAGgaga                                                                                                         >  2:1067554/1‑148 (MQ=255)
                                               aGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGgacga                                                                                                        >  1:3961387/1‑150 (MQ=255)
                                               aGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGgacga                                                                                                        <  2:2325668/150‑1 (MQ=255)
                                               aGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGgacga                                                                                                        >  1:2565680/1‑150 (MQ=255)
                                                 cAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACa                                                                                                      >  1:360039/1‑150 (MQ=255)
                                                   aTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACaaa                                                                                                    <  1:1681571/150‑1 (MQ=255)
                                                    tGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAAc                                                                                                   <  2:242033/150‑1 (MQ=255)
                                                       gATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCgg                                                                                                <  2:354813/150‑1 (MQ=255)
                                                        aTCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAGGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAAAATTATAAAGACGACAAAACggg                                                                                               >  2:3987842/1‑150 (MQ=255)
                                                        aTCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTTTAAGGACGACAAACCcgg                                                                                               >  1:3543237/1‑150 (MQ=255)
                                                        aTCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCggg                                                                                               >  2:3716233/1‑150 (MQ=255)
                                                           tGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCtt                                                                                            <  2:2597622/150‑1 (MQ=255)
                                                            gAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGGGGGCTTACCTGTGCCGTCTTGCCGATCCT‑GAAAAAAATTTTTAAAGACGACAAACCTGGCttt                                                                                           >  2:3954653/1‑149 (MQ=255)
                                                             aaaTCAGTAAAAAAGATATCCCCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTaa                                                                                          <  2:4166653/150‑1 (MQ=255)
                                                             aaaTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTaa                                                                                          <  2:2512469/150‑1 (MQ=255)
                                                               aTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTTGGGGATTTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAgc                                                                                        >  2:4127142/1‑150 (MQ=255)
                                                                    tAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCgcggc                                                                                   <  2:1266238/150‑1 (MQ=255)
                                                                    tAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGGGATGTTGCCGATCCT‑GAAAAAAATTTATAAGGACGACAAACCGGGCTTAAGCgcggg                                                                                   >  2:3485473/1‑149 (MQ=255)
                                                                      aaaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGa                                                                                 >  2:3592125/1‑150 (MQ=255)
                                                                       aaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGCGGTGTTGCCCATCCC‑GAAAAAGATTTTTTAGGACGACAAACCGGGCTTAAACgcggcgga                                                                                >  2:275832/1‑148 (MQ=255)
                                                                        aaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATg                                                                               <  1:505598/150‑1 (MQ=255)
                                                                        aaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATg                                                                               <  1:410875/150‑1 (MQ=255)
                                                                          aGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGaa                                                                             <  2:2174188/150‑1 (MQ=255)
                                                                          aGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTTCCTGGGCGGTGTTGCCCATCCT‑GAAAAAGATTTATAAAGACGACAAACCGGGCTTAAGCGCGGCGATGaa                                                                             >  2:2212227/1‑150 (MQ=255)
                                                                           gATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGGGATGTTGCCCGTCCT‑GAAAAAAATTTTTAAGGACGAAAAACCCGGCTTAAACGCGGCGATGaaa                                                                            >  2:1353854/1‑150 (MQ=255)
                                                                                cACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATaa                                                                       <  1:4285573/150‑1 (MQ=255)
                                                                                 aCCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGGTCCT‑GAAAAAGATTTATAAGGACGACAAACCCGGCTTAAGCGCGGCGATGAAAGATaaa                                                                      >  1:2175875/1‑149 (MQ=255)
                                                                                 aCCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAAc                                                                      <  1:944876/150‑1 (MQ=255)
                                                                                         gggCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGGGGTGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAAACTCGAAtt                                                              >  2:1179609/1‑150 (MQ=255)
                                                                                          ggCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAAttt                                                             <  2:594864/150‑1 (MQ=255)
                                                                                             tttCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGATTAAGCGCGGCGATGAAAGATAACCTCGAATttatt                                                          <  1:3803368/150‑1 (MQ=255)
                                                                                               tcgttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTaa                                                        <  1:207646/150‑1 (MQ=255)
                                                                                                cgttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAAt                                                       <  1:4339675/150‑1 (MQ=255)
                                                                                                  ttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTAGCTCTGCGTTGTGTTCCAGCCT‑TAAAAAGATTTATAATGAGAAAAAAACGAGCTTAAACGCGGGGTGGGAAAATAACCTCGGATttatttttat                                                     >  2:1336696/1‑145 (MQ=255)
                                                                                                  ttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATac                                                     <  1:4121892/150‑1 (MQ=255)
                                                                                                   tcgtcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATacc                                                    <  2:1839538/150‑1 (MQ=255)
                                                                                                    cgtcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATaccc                                                   <  1:2124191/150‑1 (MQ=255)
                                                                                                     gtcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGGGGGGTTAACTGGGGGGGGTTGCCGGACCC‑GAAAAAGATTTATAAGGACGACAAACCGGGGTTAAGCGCGGCGATTAAAAATAACCCTGAATTTATTAAAaacca                                                  >  2:1868915/1‑150 (MQ=255)
                                                                                                     gtcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGCTTTACCTTGGCGATGTTTCCGATCCT‑GAAAAAAATTTATAAGGACGACAAACCCGGGTTAAACGCGGCGATGAAAGATAAACTCGGATTTTTTTATaacca                                                  >  2:3146727/1‑150 (MQ=255)
                                                                                                           gctgcAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGCGGTGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGaa                                            >  2:2788976/1‑150 (MQ=255)
                                                                                                                 aGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGt                                      <  2:1320916/150‑1 (MQ=255)
                                                                                                                  gCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTc                                     <  1:2750153/150‑1 (MQ=255)
                                                                                                                     aGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGGTTTTCCTGGGGGATGTTGCCGGTCCT‑GAAAAAGATTTTTAAGGACGACAAACCGGGCTTAAGCGGGGCGATGAAAGATAACCTCGAATTTATTTATACCCACCCGAATCTGGTCGGa                                  >  2:1218394/1‑150 (MQ=255)
                                                                                                                        tttAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTc                               <  2:4343562/150‑1 (MQ=255)
                                                                                                                          tAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCt                             <  1:1624967/150‑1 (MQ=255)
                                                                                                                          tAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCt                             <  1:3074113/150‑1 (MQ=255)
                                                                                                                             cTACGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAAATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCCAATCTGGTCGGATTCCTGAt                          >  1:1070875/1‑150 (MQ=255)
                                                                                                                                cGAAAGGATGCAGGC‑GGGCGGGTTTACCTGGGCGATGTTTCCGGTCCT‑GAAAAAGAATTATAAGGACGACAAACCGGGCTTAAGCGCGGCGGTGAAAAATAACCTCGAATTTATTAATACCCACCCCAATCTTGGCGGGTTCCTGATggg                       >  2:3314943/1‑150 (MQ=255)
                                                                                                                                    aGGATGCAGGC‑GGGCGGGTTTTCCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGgta                   >  2:2253785/1‑150 (MQ=255)
                                                                                                                                     ggATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGttat                  <  1:1154072/150‑1 (MQ=255)
                                                                                                                                          cAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAAtt             >  2:3684816/1‑150 (MQ=255)
                                                                                                                                          cAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCAGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCCCCCGAATCTGGTCGGATTCCTGATGGGGTTATTAAtt             >  2:2813362/1‑150 (MQ=255)
                                                                                                                                              c‑gggcggTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGa         >  1:977136/1‑150 (MQ=255)
                                                                                                                                                gggcggTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGAt        >  1:572807/1‑150 (MQ=255)
                                                                                                                                                 ggcggTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATg       >  1:2646375/1‑150 (MQ=255)
                                                                                                                                                  gcggTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGTGGTTATTAATTTCGATgg      <  1:3683487/150‑1 (MQ=255)
                                                                                                                                                  gcggTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATgg      >  2:2152941/1‑150 (MQ=255)
                                                                                                                                                    ggTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaa    >  2:2026671/1‑150 (MQ=255)
                                                                                                                                                    ggTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaa    >  1:1813372/1‑150 (MQ=255)
                                                                                                                                                     gTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaag   <  1:62065/150‑1 (MQ=255)
                                                                                                                                                     gTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaag   >  1:3453066/1‑150 (MQ=255)
                                                                                                                                                     gTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaag   >  1:1474844/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
AAGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGGCGGTTTTACCTGGGCGATGTTGCCGATCCT‑GAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGA  >  AP012306/2758760‑2759058

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