Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123062,758,9090TG31.1% 539.4 / 33.5 235G33G (GGT→GGG) agaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (42/120);  new base G (73/0);  total (115/120)
Fisher's exact test for biased strand distribution p-value = 4.81e-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 309 total aligned reads displayed.

AAGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGA  >  AP012306/2758760‑2759058
                                                                                                                                                       |                                                                                                                                                       
      aaCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTAcc                                                                                                                                                   <  1:1274201/150‑1 (MQ=255)
       aCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCt                                                                                                                                                  <  1:2121745/150‑1 (MQ=255)
            aGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCg                                                                                                                                             >  1:1128278/1‑150 (MQ=255)
               ccGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATg                                                                                                                                          >  2:1143127/1‑150 (MQ=255)
                 gTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGtt                                                                                                                                        <  1:2995643/150‑1 (MQ=255)
                  ttGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTg                                                                                                                                       <  2:526762/150‑1 (MQ=255)
                  ttGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTg                                                                                                                                       >  1:2469435/1‑150 (MQ=255)
                    gCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGcc                                                                                                                                     <  2:678535/150‑1 (MQ=255)
                     ccAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTTCCTGGGCGATGTTGCCg                                                                                                                                    >  1:333524/1‑150 (MQ=255)
                     ccAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTTcca                                                                                                                                    >  2:61773/1‑149 (MQ=255)
                      cAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTTACTTGGCGATGTTTCCGa                                                                                                                                   >  1:3251414/1‑150 (MQ=255)
                        gTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTTGGCGATGTTGCCGATc                                                                                                                                 >  1:1750222/1‑150 (MQ=255)
                          aaaaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCt                                                                                                                               <  1:1905384/150‑1 (MQ=255)
                             aaTGAAGAAGAGGACTACAGGAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTCCCTGGGCGATGTTGCCGATTCTGaa                                                                                                                            <  2:2527338/150‑1 (MQ=255)
                             aaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTGACTACGAAAGGATGCAGGCGG‑AC‑GGTTTTACCTGGGCGATGTTGGCGGTCCTTaa                                                                                                                            >  1:3484414/1‑150 (MQ=255)
                              aTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGaaa                                                                                                                           <  1:633215/150‑1 (MQ=255)
                              aTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCGGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGaaa                                                                                                                           <  1:2831101/150‑1 (MQ=255)
                                gaagaaGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGaaaaa                                                                                                                         >  1:922856/1‑150 (MQ=255)
                                gaagaaGAGGACTACAGCAATGGGATCGGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGaaaaa                                                                                                                         >  1:2676437/1‑150 (MQ=255)
                                    aagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGAtt                                                                                                                     >  1:528536/1‑150 (MQ=255)
                                     agagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGAttt                                                                                                                    >  1:240232/1‑150 (MQ=255)
                                      gagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTta                                                                                                                   <  1:2027781/150‑1 (MQ=255)
                                      gagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTTACTGGGCGAGGTTGTCGAATCTGTAAAAGATTta                                                                                                                   >  1:3082291/1‑150 (MQ=255)
                                       agGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTtat                                                                                                                  <  1:3623465/150‑1 (MQ=255)
                                       agGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTtat                                                                                                                  <  2:752562/150‑1 (MQ=255)
                                        ggACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGACAGGATGCAGGCGG‑GC‑CGCTTTAACCGGGCGAGGTTTCCGCGCCTCCAAAAAATTttta                                                                                                                 >  2:3003764/1‑150 (MQ=255)
                                           cTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GCGGGTTTTACC‑GGGCGATGGTGCCGGTCCTTAAAAAGAATTTTTAgg                                                                                                              >  1:1023561/1‑150 (MQ=255)
                                            tACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTTCCTTGGCGATGTTTCCGGTCCTGAAAAAAAATTATaaagg                                                                                                             >  1:4347488/1‑147 (MQ=255)
                                            tACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTTCCTTGGCGATGTTGCCGATCCTGAAAAAGATTTATAAAga                                                                                                             >  2:2849642/1‑150 (MQ=255)
                                             aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCCCCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGgac                                                                                                            <  1:1802946/150‑1 (MQ=255)
                                              cAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGgacg                                                                                                           <  1:2965498/150‑1 (MQ=255)
                                                 cAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTTACTGGGCGATGTTGCCGATCCTGAAAAAGATTTTTAAAgacgaaa                                                                                                        >  2:1239636/1‑148 (MQ=255)
                                                 cAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑a                                                                                                       >  2:3388824/1‑150 (MQ=255)
                                                  aaTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑aa                                                                                                      >  2:2161473/1‑150 (MQ=255)
                                                   aTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑aaa                                                                                                     >  1:1377186/1‑150 (MQ=255)
                                                   aTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTATCTCGGCGATGTTGTCGCACCTGTAAAAAATTTTTAAGGAAGAC‑caa                                                                                                     >  1:1000727/1‑150 (MQ=255)
                                                      ggATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGCGG‑GGTTTTACTTGGG‑GAATTTGTCGATGCCGGAAAAGATTTATAAAAACGAG‑aaaaag                                                                                                  >  2:2362938/1‑147 (MQ=255)
                                                      ggATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTTCCGGTCCTGAAAAAAATTTATAAAGAAGAC‑AAAccc                                                                                                  >  2:4191496/1‑149 (MQ=255)
                                                       gATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTTCCTGGGCGATGTTGCCCATCCTGAAAAAGATTTATAAGCACGAC‑AAACCgc                                                                                                 >  2:348968/1‑149 (MQ=255)
                                                       gATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCgg                                                                                                 <  1:1724628/150‑1 (MQ=255)
                                                       gATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCgg                                                                                                 <  1:3764072/150‑1 (MQ=255)
                                                       gATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGGGATGTTGCCGGTCCTGAAAAAAATTTTTAAGGACGAC‑aaaaccg                                                                                                 >  2:1173221/1‑146 (MQ=255)
                                                        aTCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCcgg                                                                                                >  1:1232937/1‑150 (MQ=255)
                                                         tCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGc                                                                                               <  1:3694627/150‑1 (MQ=255)
                                                           tGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAAGACGAC‑AAACCCGGCtt                                                                                             >  2:2260835/1‑150 (MQ=255)
                                                            gAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGACCGTCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTa                                                                                            >  2:3797632/1‑150 (MQ=255)
                                                              aaTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAg                                                                                          <  1:318885/150‑1 (MQ=255)
                                                               aTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCCGGCTTAAgc                                                                                         >  1:803404/1‑150 (MQ=255)
                                                                tCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAgcg                                                                                        <  1:1288052/150‑1 (MQ=255)
                                                                tCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTTGGCGGTGTTGCCGGTCCTGAAAAAAAATTTTAAGGACGAC‑AAAACGGGCTTAAgcg                                                                                        >  2:920236/1‑150 (MQ=255)
                                                                   gTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTTCCTGGGCGGTGTTGCCGATCCTGGAAAAGATTTTTAAGGACGGC‑AAACCGGGCTTAAACgcgg                                                                                     >  2:3268268/1‑150 (MQ=255)
                                                                    tAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAAAATTATAAGGACGAC‑AAACCGGGCTTAAACgcggc                                                                                    >  2:719046/1‑150 (MQ=255)
                                                                       aaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGAt                                                                                 <  2:4330645/150‑1 (MQ=255)
                                                                       aaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCCGGGCGGTGGTGCCCGTCCTGAAAAAAAATTTTAAGGAGGAC‑AAACCCGGCTTAAGGgcggcgga                                                                                 >  2:3573235/1‑148 (MQ=255)
                                                                        aaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GG‑GGGTTTACCTTGGCGGTGTTGGCGGTCCTGTAAAAAATTTTTTAGGACGAA‑AAACCGGGCCTTAGCGGGGCGaag                                                                                >  2:3081030/1‑148 (MQ=255)
                                                                        aaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGGTCCTGCAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGGTg                                                                                >  2:820223/1‑150 (MQ=255)
                                                                        aaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATg                                                                                <  1:1083170/150‑1 (MQ=255)
                                                                         aaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGa                                                                               <  1:1557188/150‑1 (MQ=255)
                                                                         aaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAACGCGGCGGTGa                                                                               >  2:1361992/1‑150 (MQ=255)
                                                                              atCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGAt                                                                          <  2:279656/150‑1 (MQ=255)
                                                                              atCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGAt                                                                          <  1:3831123/150‑1 (MQ=255)
                                                                              atCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGAt                                                                          >  1:1419406/1‑150 (MQ=255)
                                                                              atCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAAATTTATAAGGGCGAC‑AAACCGGGCTTAAGCGCGGGGATGAAAGAt                                                                          >  2:536059/1‑150 (MQ=255)
                                                                                cACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATaa                                                                        <  2:1526285/150‑1 (MQ=255)
                                                                                 aCCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTTCCTGGGGGGTGGTGCCGGTCCCGAAAAAGAATTTTAAAGACGACAAAACCGGGTTTAA‑CGCGGGGGTGGAAGaaaaa                                                                       >  2:2035338/1‑149 (MQ=255)
                                                                                  cccGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAAcc                                                                      >  1:2645006/1‑150 (MQ=255)
                                                                                         gggCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑Att                                                              <  1:146611/150‑1 (MQ=255)
                                                                                            cTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTAACTGGGGGATGTTGGCGGTCCTGAAAAAGATTTATAAGGACGGC‑AAAACGGGGTTAAGGGCGGGGGTGAAAAATAACCTCGA‑ATttat                                                           >  2:3837816/1‑150 (MQ=255)
                                                                                                 gttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATa                                                      <  1:2531344/150‑1 (MQ=255)
                                                                                                  ttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATac                                                     <  1:3389617/150‑1 (MQ=255)
                                                                                                  ttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAAGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCCCGA‑ATTTATTAATac                                                     >  2:2586298/1‑150 (MQ=255)
                                                                                                   tcgtcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATacc                                                    <  2:3657095/150‑1 (MQ=255)
                                                                                                            ctgcAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTTCCTGGGGGGGGTTGCCGGTCCTGAAAAAGATTTATAAGGGCGGC‑AAAACGGGCTTAAGCGCGGGGGTGAAAGATAACCCCGA‑ATTTATTAATACCCCCCCGaac                                           >  2:3293561/1‑149 (MQ=255)
                                                                                                            ctgcAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTTCCCGGGCGATGTTGCCGATCCTGAAAAAAAATTTTAAAGAGGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGAAAAACTCCA‑ATTTATTAATAACCACCCCAAt                                           >  2:1901523/1‑150 (MQ=255)
                                                                                                               cAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGGTGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTg                                        >  1:3954803/1‑150 (MQ=255)
                                                                                                               cAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTg                                        <  1:4162648/150‑1 (MQ=255)
                                                                                                               cAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTg                                        >  1:3613314/1‑150 (MQ=255)
                                                                                                                aaGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTTCCGATCCTGAAAAAGATTTTTAAGGAAGAC‑AAACCGGGCTTAAGCGGGGCGATGAAAGAAAACCTCGA‑ATTTTTTAATACCCACCCGAATCTgg                                       >  2:3311440/1‑150 (MQ=255)
                                                                                                                aaGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTgg                                       >  1:4130230/1‑150 (MQ=255)
                                                                                                                aaGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTACCTGGGCGATGTTGCCCATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTgg                                       >  1:2382889/1‑150 (MQ=255)
                                                                                                                 aGCGAGCTTTAACTCCGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGt                                      <  1:58260/150‑1 (MQ=255)
                                                                                                                   cGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCg                                    <  1:403832/150‑1 (MQ=255)
                                                                                                                     aGCTTTAACTACGAAAGGATGCAGGCGG‑GG‑GGGTTTACCTGGGGGGTGTTGCCGAACCTGAAAAAGATTTATTAAGAAGCC‑AAACCGGGGTTAAGAACGGCGGAGAAAGATAAACTCGAGATTTTTTATTA‑CCACCCGCATCTGGGCggg                                  >  2:1803050/1‑149 (MQ=255)
                                                                                                                          tAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCt                             <  1:1605742/150‑1 (MQ=255)
                                                                                                                          tAACTACGAAAGGATGCAGGCGG‑GC‑GGGTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCt                             >  1:1311198/1‑150 (MQ=255)
                                                                                                                            aCTACGAAAGGATGCAGGCGG‑GC‑GGGTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTTCTGa                           >  1:2727128/1‑150 (MQ=255)
                                                                                                                             cTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGAt                          <  2:59601/150‑1 (MQ=255)
                                                                                                                              tACGAAAGGAAGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATg                         >  1:3905720/1‑150 (MQ=255)
                                                                                                                                cGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTCGGCGCTTGTGCCGGTCCTGGAAAAGATTTATAAGGACGAC‑CAACCCGGCTTAAACGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCCAAACTGGGCGGGTTCCTTATtgg                       >  2:108860/1‑150 (MQ=255)
                                                                                                                                     ggATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGttat                  <  2:2465301/150‑1 (MQ=255)
                                                                                                                                     ggATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGttat                  >  1:2548840/1‑150 (MQ=255)
                                                                                                                                       aTGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGttatta                <  1:2262506/150‑1 (MQ=255)
                                                                                                                                       aTGCAGGCGG‑GC‑GGGTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCCGGCTTAAGCGCGGCGGTGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGgtatta                >  2:2171830/1‑150 (MQ=255)
                                                                                                                                            ggcgg‑gc‑gGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTc           >  2:1656291/1‑150 (MQ=255)
                                                                                                                                               gg‑gc‑ggTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGAt        >  2:3152945/1‑150 (MQ=255)
                                                                                                                                               gg‑gc‑ggTTTTACCTCGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGAt        <  1:683054/150‑1 (MQ=255)
                                                                                                                                                  gc‑ggTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATgg      >  2:3899800/1‑150 (MQ=255)
                                                                                                                                                      gTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaag   >  2:1418385/1‑150 (MQ=255)
                                                                                                                                                      gTTTTACCTGGGCGATGTTACCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaag   <  1:4291401/150‑1 (MQ=255)
                                                                                                                                                       ttttACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaaga  >  1:1626912/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                       
AAGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGG‑GC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGA‑ATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGA  >  AP012306/2758760‑2759058

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