Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123062,758,9090TG26.9% 656.1 / 25.4 257G33G (GGT→GGG) agaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (46/142);  new base G (69/0);  total (115/142)
Fisher's exact test for biased strand distribution p-value = 6.50e-32
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 353 total aligned reads displayed.

AAGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGA  >  AP012306/2758760‑2759058
                                                                                                                                                      |                                                                                                                                                         
aaGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGt                                                                                                                                                           >  1:4984237/1‑150 (MQ=255)
    gcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTa                                                                                                                                                       >  1:1374176/1‑150 (MQ=255)
    gcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTa                                                                                                                                                       >  2:4356697/1‑150 (MQ=255)
      aaCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTAcc                                                                                                                                                     >  1:1078253/1‑150 (MQ=255)
        ccGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCtt                                                                                                                                                   >  1:93098/1‑149 (MQ=255)
         cGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTgg                                                                                                                                                  <  2:2230057/150‑1 (MQ=255)
             gCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGa                                                                                                                                              >  1:2991458/1‑150 (MQ=255)
              cccGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGAt                                                                                                                                             <  1:4328516/150‑1 (MQ=255)
               ccGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGATg                                                                                                                                            >  2:3124834/1‑150 (MQ=255)
               ccGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGGTg                                                                                                                                            >  1:2609633/1‑150 (MQ=255)
                cGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑t                                                                                                                                          <  2:305215/150‑1 (MQ=255)
                cGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑t                                                                                                                                          <  2:3138496/150‑1 (MQ=255)
                   tGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTTCCTCGCCGATT‑TTGc                                                                                                                                       >  1:1856527/1‑150 (MQ=255)
                    gCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGGTG‑TTGcc                                                                                                                                      >  2:3450518/1‑150 (MQ=255)
                     ccAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGG‑GGGTTTGACTGGGCGCTG‑tatccg                                                                                                                                     >  2:2858151/1‑145 (MQ=255)
                     ccAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCg                                                                                                                                     <  1:3223228/150‑1 (MQ=255)
                        gTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATc                                                                                                                                  <  2:1185215/150‑1 (MQ=255)
                         tAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATG‑TTTCCCATcc                                                                                                                                 >  2:4018334/1‑150 (MQ=255)
                           aaaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATG‑TTTCCGGTCCTg                                                                                                                               >  2:65939/1‑150 (MQ=255)
                            aaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGGTG‑TTTCCGGTCCTGa                                                                                                                              >  2:4799568/1‑150 (MQ=255)
                              aTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATT‑GTGCCGATCCTGgaa                                                                                                                            >  2:1977268/1‑150 (MQ=255)
                               tGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTATCTGGGCGGTG‑TTGCCGATCCTGaaaa                                                                                                                           >  2:2697884/1‑150 (MQ=255)
                               tGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGaaaa                                                                                                                           <  1:316512/150‑1 (MQ=255)
                               tGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGATG‑TTGCCGATCCTGaaaa                                                                                                                           >  2:4195494/1‑150 (MQ=255)
                               tGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGATG‑TTGCCGATCCTGaaaa                                                                                                                           >  2:265123/1‑150 (MQ=255)
                                gaagaaGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGaaaaa                                                                                                                          >  1:3203943/1‑150 (MQ=255)
                                gaagaaGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGaaaaa                                                                                                                          <  1:4314694/150‑1 (MQ=255)
                                gaagaaGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGaaaaa                                                                                                                          <  2:4032119/150‑1 (MQ=255)
                                 aagaagAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAg                                                                                                                         <  1:2863253/150‑1 (MQ=255)
                                 aagaagAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAg                                                                                                                         >  1:63007/1‑150 (MQ=255)
                                  agaagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGATG‑TTGCCGGTCCTGAAAAAGa                                                                                                                        >  2:1810451/1‑150 (MQ=255)
                                   gaagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGAt                                                                                                                       <  2:4957627/150‑1 (MQ=255)
                                   gaagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTAGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGAt                                                                                                                       <  1:1301475/150‑1 (MQ=255)
                                       agGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGGGGGG‑TTGGCGATCCTGGAAAAGAttttt                                                                                                                   >  2:2092406/1‑148 (MQ=255)
                                        ggACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTTCCTGGGCGATG‑TTGCCGGCCCTGAAAAAGAATtata                                                                                                                  >  2:240793/1‑150 (MQ=255)
                                        ggACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTAACTGGGCGATG‑GTGCCGATCCTGAAAAAGATTtata                                                                                                                  >  1:370801/1‑150 (MQ=255)
                                        ggACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGCTTTACCTTGGCGATG‑TTTTCGCACCTGGAAAAGAATttta                                                                                                                  >  1:3835574/1‑150 (MQ=255)
                                           cTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAgg                                                                                                               <  1:3394525/150‑1 (MQ=255)
                                           cTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATG‑TTTCCGATCCTGAAAAAAAATTTTAAgg                                                                                                               >  2:169588/1‑150 (MQ=255)
                                             aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGgac                                                                                                             <  1:2390665/150‑1 (MQ=255)
                                             aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGATG‑TTTTCGATCCTGAAAAAAAATTTTAAAgac                                                                                                             >  1:5223112/1‑150 (MQ=255)
                                                 cAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGG‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑a                                                                                                        >  2:306085/1‑150 (MQ=255)
                                                  aaTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATG‑TTGCCGATCCTGTAAAAGATTTATAAGGACGAC‑aa                                                                                                       >  1:1662166/1‑150 (MQ=255)
                                                   aTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑aaa                                                                                                      <  1:3666692/150‑1 (MQ=255)
                                                   aTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑aaa                                                                                                      <  1:570835/150‑1 (MQ=255)
                                                     gtgATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGCTTTCCTTCGCGCTG‑GTGGCGATCCTCAAAAAAAATTATAAGGAAGAa‑aaaac                                                                                                    >  2:364232/3‑148 (MQ=255)
                                                             aaaTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTaa                                                                                           <  2:2350932/150‑1 (MQ=255)
                                                             aaaTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAAGACGAC‑AAACCGGG‑CTTaa                                                                                           >  1:4453321/1‑150 (MQ=255)
                                                                 cAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTTCTTGGGGGTGG‑TTG‑CGGTCCTTAAAAAAAATTTTAAAGAAGACAAAACCGGG‑CTTaacgc                                                                                        >  2:5248354/1‑147 (MQ=255)
                                                                 cAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAgcgc                                                                                       <  1:3599846/150‑1 (MQ=255)
                                                                 cAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑CGGTTTAACTGGGCGATT‑TTGCCGATCCTGGAAAAGATTTTAAAAGGCGAC‑AAAACCGG‑GTTAAgcgc                                                                                       >  2:4664720/1‑150 (MQ=255)
                                                                   gTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCgcgg                                                                                     <  1:1017358/150‑1 (MQ=255)
                                                                    tAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCgcggc                                                                                    <  2:4134831/150‑1 (MQ=255)
                                                                    tAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCgcggc                                                                                    >  1:414136/1‑150 (MQ=255)
                                                                       aaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGGCGAC‑AAACCGGG‑CTTAAGCGCGGCGAt                                                                                 >  2:4880389/1‑150 (MQ=255)
                                                                         aaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑Ga                                                                              <  1:1909646/150‑1 (MQ=255)
                                                                           gATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTACTGGGCGATGTTTGCCGA‑ACTGGAAAAGGTTTTTAAAGAAGAC‑AAAACGGG‑CTTAAACGCGGCGGT‑Ggaa                                                                            >  2:5054494/1‑150 (MQ=255)
                                                                            atatCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGGT‑GAAAg                                                                           >  1:3030480/1‑150 (MQ=255)
                                                                             tatCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGa                                                                          <  1:202600/150‑1 (MQ=255)
                                                                             taccccccGTCTGGGCTTTCGTTCGGCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGa                                                                          <  2:2545075/145‑1 (MQ=255)
                                                                              atCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGAt                                                                         <  1:2645060/150‑1 (MQ=255)
                                                                               tCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATa                                                                        <  1:2836431/150‑1 (MQ=255)
                                                                               tCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATG‑TTGACGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATa                                                                        >  2:3171618/1‑150 (MQ=255)
                                                                                  cccGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGTTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTTAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAAcc                                                                     <  2:25925/150‑1 (MQ=255)
                                                                                   ccGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCt                                                                    <  1:4773736/150‑1 (MQ=255)
                                                                                   ccGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCt                                                                    >  1:4690437/1‑150 (MQ=255)
                                                                                    cGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCCGGGGGGTG‑GTGCCCATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CCTAAGCGCGGCGAG‑TTAAAATAACCTc                                                                   >  2:1290054/1‑150 (MQ=255)
                                                                                     gTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGGTG‑TTGCCGGTCCTGAAAAAGATTTTTAAGGAAGGC‑AAACCGGG‑CTTAAGCGCGGCGGT‑GAAAGATAAACTCg                                                                  >  2:3825440/1‑150 (MQ=255)
                                                                                      tCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGa                                                                 >  1:1503918/1‑150 (MQ=255)
                                                                                        tGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTTCCTGGGCGGTG‑TTGCCGGTCCTGAAAAAAATTTTTAAAGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAAACTCGaaa                                                               >  2:2152114/1‑149 (MQ=255)
                                                                                            cTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCCGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATttat                                                           >  1:4379601/1‑150 (MQ=255)
                                                                                              ttcgttcgTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATG‑TTGCCGATCCCGAAAAAGATTTATAAAGGCGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGGTAACCTCGAATttttta                                                         >  2:1430278/1‑150 (MQ=255)
                                                                                               tcgttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTaa                                                        <  1:828114/150‑1 (MQ=255)
                                                                                                cgttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGG‑GGGTTTTTCTGGGCGGTGTTTCCCGTCCCGGAAAAAAATTTTTAA‑GACGACAAAACCGGG‑TTTAA‑CGCGGCGGTGGAAAAAAAACCCCAAATTTTTTaa                                                        >  2:3873211/1‑150 (MQ=255)
                                                                                                    cgtcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATaccc                                                   <  1:1616541/150‑1 (MQ=255)
                                                                                                     gtcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCCGGGCGATG‑TTGCCGGTCCTGAAAAAGATTTATAAGGACGAC‑AAAACCGG‑CCTAAACGCGGCGGT‑GGAAGATAACCTCGAAATTATTAATaacca                                                  >  2:624876/1‑150 (MQ=255)
                                                                                                       cgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATacccacc                                                >  1:259668/1‑150 (MQ=255)
                                                                                                       cgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATacccacc                                                <  1:4859425/150‑1 (MQ=255)
                                                                                                          tgctgcAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTTGGGGATG‑TTTTCGATCCTGAAAAAGATTTTTAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAAACTCCAATTTATTAATAACCCCCCGa                                             >  2:1852520/1‑150 (MQ=255)
                                                                                                            ctgcAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAAt                                           <  1:2068116/150‑1 (MQ=255)
                                                                                                                aaGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTgg                                       <  1:2508971/150‑1 (MQ=255)
                                                                                                                   cGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGGGATG‑TTGCCCATCCTGCAAAAGAATTATAAAGAAGAC‑AAAACCGG‑CTTAAGCGCCGCGGT‑GAAAGAAAAACTCGAAATTATTTATAACCACCCCAATCTGGTCg                                    >  2:5180189/1‑150 (MQ=255)
                                                                                                                   cGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCg                                    >  1:1720594/1‑150 (MQ=255)
                                                                                                                     aGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGCTG‑TTGCCGATTCCGAAAAAGATTTATAAGGACGAC‑AAACCGGGCCTTAAG‑GCGGGGGT‑GAAAGATAACCTCGAATTTTTTAATAACCACCCGAATCTgggtggg                                  >  2:4530736/1‑145 (MQ=255)
                                                                                                                            aCTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCCAATCTGGTCGGATTCCTGa                           >  1:2798575/1‑150 (MQ=255)
                                                                                                                            aCTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCCGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGGCGGATTCCTGa                           >  1:397572/1‑150 (MQ=255)
                                                                                                                             cTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGAt                          <  1:1244373/150‑1 (MQ=255)
                                                                                                                             cTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGGGGTG‑GTGCCGGTCCTGAAAAAGATTTATAAGGACGAC‑AAAACGGG‑GTTAAGCGGGGCGGA‑GAAAGAAAAACTCGAAATTATTAATAACCACCCGAAACTGGTCGGAATTCTGAt                          >  2:759657/1‑150 (MQ=255)
                                                                                                                                      gATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCCCCCGAATCTGGTCGGATTCCTGATGGGGttatt                 <  2:3586421/150‑1 (MQ=255)
                                                                                                                                       aTGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGttatta                <  2:1024968/150‑1 (MQ=255)
                                                                                                                                         gCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAAt              <  1:2291201/150‑1 (MQ=255)
                                                                                                                                          cAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAAtt             <  1:3881962/150‑1 (MQ=255)
                                                                                                                                          cAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAAtt             >  1:882588/1‑150 (MQ=255)
                                                                                                                                            ggcgggc‑gGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTc           <  1:3130806/150‑1 (MQ=255)
                                                                                                                                              cgggc‑ggTTTTACCTGGGCGATG‑TTGCCGTTCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGa         >  2:3290523/1‑150 (MQ=255)
                                                                                                                                              cgggc‑ggTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGa         >  2:3557066/1‑150 (MQ=255)
                                                                                                                                               gggc‑ggTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGAt        >  1:1430098/1‑150 (MQ=255)
                                                                                                                                                ggc‑ggTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATg       >  1:5109601/1‑150 (MQ=255)
                                                                                                                                                     gTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaag   >  2:16101/1‑150 (MQ=255)
                                                                                                                                                      ttttACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaaga  >  1:3445962/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                         
AAGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATG‑TTGCCGATCCTGAAAAAGATTTATAAGGACGAC‑AAACCGGG‑CTTAAGCGCGGCGAT‑GAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGA  >  AP012306/2758760‑2759058

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