Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123062,758,9090TG27.4% 591.7 / 29.2 237G33G (GGT→GGG) agaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (40/132);  new base G (65/0);  total (105/132)
Fisher's exact test for biased strand distribution p-value = 9.49e-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 296 total aligned reads displayed.

AAGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGA  >  AP012306/2758760‑2759058
                                                                                                                                                       |                                                                                                                                                      
aaGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑t                                                                                                                                                        >  1:2480644/1‑150 (MQ=255)
aaGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑t                                                                                                                                                        >  1:3535467/1‑150 (MQ=255)
 aGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑tt                                                                                                                                                       <  2:1046988/150‑1 (MQ=255)
  gcgcAATCGAGGGCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATAACACCCGTCTGGGCTTTCCTTCGTCGCTGCCGCAAGCGAGCTTTAACTAAGAAAGGATGCAGGCGCGC‑GG‑ttt                                                                                                                                                      <  2:517841/150‑1 (MQ=255)
  gcgcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTTGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑ttt                                                                                                                                                      <  2:2236766/150‑1 (MQ=255)
    gcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTa                                                                                                                                                    >  2:1657159/1‑150 (MQ=255)
     cAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑gtttta                                                                                                                                                   >  2:818007/1‑148 (MQ=255)
       aCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCt                                                                                                                                                 <  1:2230638/150‑1 (MQ=255)
       aCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCt                                                                                                                                                 <  1:703670/150‑1 (MQ=255)
         cGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTgg                                                                                                                                               >  1:1461877/1‑150 (MQ=255)
           cAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGG‑TTTTACTTggg                                                                                                                                              >  1:2044810/1‑150 (MQ=255)
           cAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGc                                                                                                                                             >  1:656241/1‑150 (MQ=255)
            aGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTCCTGGGCg                                                                                                                                            >  2:873018/1‑150 (MQ=255)
            aGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTGGGCg                                                                                                                                            >  1:971884/1‑150 (MQ=255)
              cccGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTTCTTGGCggt                                                                                                                                          >  1:3139107/1‑148 (MQ=255)
                cGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTGGGCGgtgt                                                                                                                                        >  1:3255881/1‑150 (MQ=255)
                 gTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTTCCTGGGCGGTGtt                                                                                                                                       >  2:1434885/1‑150 (MQ=255)
                 gTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGtt                                                                                                                                       <  2:2190618/150‑1 (MQ=255)
                  ttGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTg                                                                                                                                      <  1:1571290/150‑1 (MQ=255)
                   tGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGc                                                                                                                                     <  1:859546/150‑1 (MQ=255)
                   tGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGc                                                                                                                                     <  1:3926019/150‑1 (MQ=255)
                   tGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTCCTGGGCGATGTTGc                                                                                                                                     >  2:2818887/1‑150 (MQ=255)
                   tGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTGGGCGATGTTGc                                                                                                                                     >  1:1707450/1‑150 (MQ=255)
                    gCCAGTAAAAATGAATAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGCT‑GG‑TTTTGCCTGGGCGGTGTTGcc                                                                                                                                    >  2:1681369/1‑150 (MQ=255)
                    gCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGcc                                                                                                                                    <  1:103702/150‑1 (MQ=255)
                    gCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGcc                                                                                                                                    <  1:236863/150‑1 (MQ=255)
                    gCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCCCGGCGATGTTGcc                                                                                                                                    >  1:3746600/1‑150 (MQ=255)
                     ccAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCg                                                                                                                                   <  2:2275084/150‑1 (MQ=255)
                       aGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGAt                                                                                                                                 <  2:172891/150‑1 (MQ=255)
                       aGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGAt                                                                                                                                 <  1:1931335/150‑1 (MQ=255)
                         tAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTCCTGGGGGGTGTTTCCGCTcc                                                                                                                               >  2:2871779/1‑150 (MQ=255)
                           aaaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTg                                                                                                                             >  1:3485796/1‑150 (MQ=255)
                            aaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTTACCT‑GGCGATTTTGTCGCTCCTGa                                                                                                                            >  2:1833140/1‑150 (MQ=255)
                              aTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGaaa                                                                                                                          >  1:3836980/1‑150 (MQ=255)
                                gaagaaGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTCCCGGGGGGTGTTTGCGATCCTGaaaaa                                                                                                                        >  2:2846740/1‑150 (MQ=255)
                                gaagaaGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTGGGCGATGTTGCCGATCCTGaaaaa                                                                                                                        >  1:3652363/1‑150 (MQ=255)
                                gaagaaGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTAACTGGGCGATGTTGCCGATCCTGaaaaa                                                                                                                        >  1:751569/1‑150 (MQ=255)
                                   gaagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGAt                                                                                                                     >  1:2785000/1‑150 (MQ=255)
                                    aagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGAtt                                                                                                                    <  2:3905393/150‑1 (MQ=255)
                                      gagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTta                                                                                                                  <  1:1645166/150‑1 (MQ=255)
                                      gagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTCCTGGGCGATGTTGCCGGTCCTTAAAAAGATTta                                                                                                                  >  2:1594341/1‑150 (MQ=255)
                                          aCTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAg                                                                                                              <  1:1267711/150‑1 (MQ=255)
                                          aCTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTTGGCGGTGGTTGCGCTACTGAAAAAGAatttttag                                                                                                              >  1:4222376/1‑145 (MQ=255)
                                                 cAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACa                                                                                                       <  1:3671620/150‑1 (MQ=255)
                                                  aaTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTGGGCGATGTTGCCGGTCCTGAAAAAGATTTTTTAGgacgacca                                                                                                      >  1:2184684/1‑148 (MQ=255)
                                                    tGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAAGACGACA‑AAc                                                                                                   >  1:826178/1‑150 (MQ=255)
                                                     gggATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGG‑GG‑GTTTTTCTCGGCGATGTTGTCGCACCTGTAAAAGATTTATAAAGACGACa‑aaac                                                                                                  >  2:3737532/1‑148 (MQ=255)
                                                     gggATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTATCTCGGCGATGTTTCCGATCCTGAAAAAGAATTATAAGGGCGACA‑AAcc                                                                                                  >  2:423637/1‑150 (MQ=255)
                                                       gATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCCGGGCGGTGTTGCCGGTCCTGAAAAAGAATTATAAGGACGACA‑AAcccg                                                                                                >  2:1557082/1‑148 (MQ=255)
                                                       gATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTTGGGGATGTTGCCGATCCTGAAAAAAATTTATAAGGACGACA‑AAACgg                                                                                                >  2:2632200/1‑150 (MQ=255)
                                                         tCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGc                                                                                              <  2:1909152/150‑1 (MQ=255)
                                                         tCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTGGGCGATGTTTCCGATCCTGAAAAAGAATTATAAGGACGACA‑AAACGGGc                                                                                              >  2:1109215/1‑150 (MQ=255)
                                                          cTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTCCTGGGCGATGTTGCCGATCCTGGAAAAGATTTATAAAGACGACA‑AAACCGGCt                                                                                             >  1:3028228/1‑150 (MQ=255)
                                                          cTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTTGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCt                                                                                             >  1:2891490/1‑150 (MQ=255)
                                                           tGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTGGGCGGTGTTGCCGATCCTCAAAAAGATTTTTAAGGAGGACA‑AACCCGGCtt                                                                                            >  2:1311055/1‑150 (MQ=255)
                                                              aaTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTCCTGGGCGGGGTTGCCGATTCTGAAAAAGATTTATAAAGAAGACA‑AAACCGGCTTTAg                                                                                         >  2:2593520/1‑150 (MQ=255)
                                                                 cAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGGAA‑AACCGGGGTTAAgcgc                                                                                      >  2:3528665/1‑150 (MQ=255)
                                                                    tAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCgcggc                                                                                   <  1:1992944/150‑1 (MQ=255)
                                                                     aaaaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCgcggcg                                                                                  <  1:895395/150‑1 (MQ=255)
                                                                     aaaaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCgcggcg                                                                                  <  1:210469/150‑1 (MQ=255)
                                                                     aaaaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCgcggcg                                                                                  <  2:2726643/150‑1 (MQ=255)
                                                                      aaaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTTTCTGGGCGATGGTGCCGGACCCCAAAAAGAATTTTAAGGAAGGCACAACCGGGC‑TAAGCgcggcgg                                                                                 >  2:933108/1‑149 (MQ=255)
                                                                      aaaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGa                                                                                 >  1:3609078/1‑150 (MQ=255)
                                                                        aaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATg                                                                               <  2:3116932/150‑1 (MQ=255)
                                                                            atatCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAg                                                                           <  1:3903251/150‑1 (MQ=255)
                                                                               tCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGGAAAAGATTTATAAGGACGACA‑AACCGGGGTTAAGCGCGGCGATGAAAGATa                                                                        >  1:3273522/1‑150 (MQ=255)
                                                                               tCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATa                                                                        <  1:34029/150‑1 (MQ=255)
                                                                                 aCCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTTCCTTGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAAc                                                                      >  1:384622/1‑150 (MQ=255)
                                                                                 aCCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTCCTTGGGGGTGTTTTCCGTCCTTTAAAAAAATTATTAAGGCGAAA‑AACCCGGCTTAAAGGCGGGGGGGAAAGATaaa                                                                      >  2:2003882/1‑149 (MQ=255)
                                                                                 aCCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAAc                                                                      >  1:3203722/1‑150 (MQ=255)
                                                                                   ccGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGGAAGATAACCt                                                                    >  1:2641126/1‑150 (MQ=255)
                                                                                   ccGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCt                                                                    <  1:2938151/150‑1 (MQ=255)
                                                                                    cGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTGGGGGGTGTTGCCGATCCTGAAAAAGATTTATAAAGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAAATaaaatc                                                                   >  2:2676603/1‑146 (MQ=255)
                                                                                     gTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTCCTGGGCGGTGTTGCCGATCCTGAAAAAAATTTATAAGGACGACA‑AAACGGGCTTAAGGGCGGCGGTGAAAAATAACCTCg                                                                  >  2:3692234/1‑150 (MQ=255)
                                                                                        tGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAAt                                                               <  1:589345/150‑1 (MQ=255)
                                                                                            cTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATttat                                                           <  1:1338679/150‑1 (MQ=255)
                                                                                             tttCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCCGGGGGGTGTTTTCGAACCTCAAAAAGAATTATAAGGACGGCA‑AACCCGGCCTTAGCGCGGCGGTGAAAGATAAACCCGAAattatt                                                          >  2:4026299/1‑150 (MQ=255)
                                                                                             tttCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGGTTTAGCGCGGCGATGAAAGATAACCTCGAAtttttt                                                          >  2:2449591/1‑150 (MQ=255)
                                                                                                cgttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAAt                                                       <  1:4058541/150‑1 (MQ=255)
                                                                                                 gttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGGGATGGTGCCGGTCCTGAAAAAGATTTATAAGGGCGACA‑AACCGGGCCTAAGCGCCGCGGTGAAAAATAACCTCGGATTTATTAATa                                                      >  2:3804658/1‑150 (MQ=255)
                                                                                                 gttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAAGACGACA‑AACCCGGCTTAAGCGCGGCGATGAAAGATAACCTCCAATTTATTtata                                                      >  2:3027556/1‑150 (MQ=255)
                                                                                                  ttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTCCTGGGCGGTGGTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGGTGAAAGATAACCTCGAATTTATTAATaa                                                     >  2:2341096/1‑149 (MQ=255)
                                                                                                        gctgctGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTCCTGGGCGATGTTTCCGATCCCGAAAAAGATTTATAAAGGCGACA‑AACCCGGCTTAAGCGCCGCGGAGAAAGATAACCTCGAATTTATTAATacccaccc                                               >  2:2532946/1‑150 (MQ=255)
                                                                                                               cAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTg                                        <  1:2715284/150‑1 (MQ=255)
                                                                                                                 aGCGAGCTTTAACTACGAAAGGATGCAGGCGGGG‑GG‑GTTTACCTGGGCGGTGTTTCCGAACCTGGAAAAGAATTATAAAGACGCCA‑AACCGGGCGTAAACGCGGCGGTGTAAGATAACCTAGAATTTATTAATATCCACCCGAATCTggg                                      >  2:59236/1‑149 (MQ=255)
                                                                                                                 aGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTACCTGGGGGATGTTGCCGATCCTGAAAAAGATTTATAAGGGCGACA‑AACCCGGGTTAAGCGCGGGGGTGAAAGAAAACCTCCAAATTATTAATAACCACCCGAATCTggg                                      >  2:414819/1‑149 (MQ=255)
                                                                                                                   cGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTCCTGGGGGATGTTTGCGGTCCTGAAAAAAAATTTTAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCCGATCTGGTCg                                    >  2:2926273/1‑150 (MQ=255)
                                                                                                                    gAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑GTTTTCCTGGGCGGTGTTGCCGGTCCTGAAAAAGATTTTTAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTTTTAATACCCACCCGAATCTGGTCgg                                   >  2:771187/1‑150 (MQ=255)
                                                                                                                          tAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCt                             <  1:1154007/150‑1 (MQ=255)
                                                                                                                          tAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCt                             <  2:2040311/150‑1 (MQ=255)
                                                                                                                           aaCTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTg                            >  1:438643/1‑150 (MQ=255)
                                                                                                                             cTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTTGTCGGATTCCTGAt                          >  1:484050/1‑150 (MQ=255)
                                                                                                                                cGAAAGGATGCAGGCGGGC‑GG‑TTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAAGACGACA‑AACCGGGCTTAAGCGCGGCGGTGAAAGATAAACTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATggg                       >  2:40549/1‑150 (MQ=255)
                                                                                                                                    aGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGtta                   <  1:2742873/150‑1 (MQ=255)
                                                                                                                                      gATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAGAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGttatt                 <  1:4153962/150‑1 (MQ=255)
                                                                                                                                         gCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAAt              <  2:3254709/150‑1 (MQ=255)
                                                                                                                                           aGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAAttt            <  1:2068442/150‑1 (MQ=255)
                                                                                                                                               gggc‑gg‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGAt        <  2:2715262/150‑1 (MQ=255)
                                                                                                                                                 gc‑gg‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATgg      >  2:610488/1‑150 (MQ=255)
                                                                                                                                                    gg‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaa    >  1:3073600/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
AAGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GG‑TTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACA‑AACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGA  >  AP012306/2758760‑2759058

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