Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123062,758,9090TG24.6% 581.9 / 16.6 229G33G (GGT→GGGagaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (35/137);  new base G (56/0);  total (92/137)
Fisher's exact test for biased strand distribution p-value = 1.90e-29
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 305 total aligned reads displayed.

AAGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGA  >  AP012306/2758760‑2759058
                                                                                                                                                      |                                                                                                                                                     
aaGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGt                                                                                                                                                       <  2:2336016/150‑1 (MQ=255)
  gcgcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGttt                                                                                                                                                     <  2:1516293/150‑1 (MQ=255)
   cgcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGtttt                                                                                                                                                    <  1:4037766/150‑1 (MQ=255)
   cgcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGtttt                                                                                                                                                    <  2:447660/150‑1 (MQ=255)
   cgcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCAGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGgttt                                                                                                                                                    >  2:3621588/1‑150 (MQ=255)
    gcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTACAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGgtttt                                                                                                                                                   >  2:3177954/1‑149 (MQ=255)
    gcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTa                                                                                                                                                   <  2:210250/150‑1 (MQ=255)
       aCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCt                                                                                                                                                <  2:2011658/150‑1 (MQ=255)
           cAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGc                                                                                                                                            <  1:755476/150‑1 (MQ=255)
              cccGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCTGGGCGAt                                                                                                                                         >  2:1369142/1‑150 (MQ=255)
                 gTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTACTTGGGCggttt                                                                                                                                       >  2:1851962/1‑146 (MQ=255)
                 gTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGtt                                                                                                                                      <  2:919398/150‑1 (MQ=255)
                 gTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGtt                                                                                                                                      <  2:2934452/150‑1 (MQ=255)
                 gTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGtt                                                                                                                                      <  2:3563467/150‑1 (MQ=255)
                 gTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACAAAAGGATGCAGGCGGGCGGGTTTTCCTTGGCggagtt                                                                                                                                      >  1:945213/1‑145 (MQ=255)
                  ttGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTg                                                                                                                                     <  1:3303161/150‑1 (MQ=255)
                    gCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGcc                                                                                                                                   <  2:126052/150‑1 (MQ=255)
                    gCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGcc                                                                                                                                   <  1:2407252/150‑1 (MQ=255)
                     ccAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCg                                                                                                                                  <  1:1770587/150‑1 (MQ=255)
                         tAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATcc                                                                                                                              <  2:1454904/150‑1 (MQ=255)
                          aaaaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGTCGGGTTTTTCCGGGGGATGTTGTCGGTCCt                                                                                                                             >  2:2752380/1‑150 (MQ=255)
                           aaaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGCTTTTCCTGGGCGCTGTTGCCGATACTg                                                                                                                            >  1:1871736/1‑150 (MQ=255)
                            aaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGGGGGTTTTCCCGGGGGATGTTGGCGAACCTGa                                                                                                                           >  2:1923782/1‑150 (MQ=255)
                            aaaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGa                                                                                                                           <  1:2922150/150‑1 (MQ=255)
                             aaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGaa                                                                                                                          <  1:502059/150‑1 (MQ=255)
                             aaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGaa                                                                                                                          >  2:620861/1‑150 (MQ=255)
                             aaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGCCGGTTTTACCGGGGCGATGTTGCCGTTACTgga                                                                                                                          >  2:102299/1‑148 (MQ=255)
                              aTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGaaa                                                                                                                         <  2:1656363/150‑1 (MQ=255)
                                gaagaaGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGaaaaa                                                                                                                       <  2:2200772/150‑1 (MQ=255)
                                 aagaagAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGCTCCTGAAAAAg                                                                                                                      >  2:1781385/1‑150 (MQ=255)
                                    aagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTTCCTGGGCGATGTTTCCCGTCCTGCAAAAGAtt                                                                                                                   >  1:3409255/1‑150 (MQ=255)
                                    aagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTAACTGGGCGGTGTTTCCGATCCTGGAAAAAAtt                                                                                                                   >  2:2436190/1‑150 (MQ=255)
                                     agagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGAttt                                                                                                                  <  2:3202888/150‑1 (MQ=255)
                                     agagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCTGGGGGGTGTTTCCGGTCCTGaaaaaaaatt                                                                                                                  >  2:2534858/1‑145 (MQ=255)
                                      gagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGGTCCTGAAAAAGATTta                                                                                                                 >  1:3325866/1‑150 (MQ=255)
                                        ggACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAGATTtata                                                                                                               <  1:1027141/150‑1 (MQ=255)
                                        ggACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTtata                                                                                                               >  1:1609775/1‑150 (MQ=255)
                                         gACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATaa                                                                                                              <  2:1953150/150‑1 (MQ=255)
                                          aCTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTATGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAg                                                                                                             >  2:2460537/1‑150 (MQ=255)
                                           cTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGGAAAAGATTTATAAgg                                                                                                            >  1:3295900/1‑150 (MQ=255)
                                            tACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGGTGTTGCCGATCCTGAAAAAGATTTATAAGga                                                                                                           >  1:133924/1‑150 (MQ=255)
                                             aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGGGGGTTTTTCTGGGCGGGGTTGCCGAACCCTGAAAAGATTTTTTAGgag                                                                                                          >  2:3105353/1‑149 (MQ=255)
                                             aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTTCCTTGGGGATGTTGCCGATCCTGGAAAAAAATTATAAGgaa                                                                                                          >  2:3430042/1‑149 (MQ=255)
                                             aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGgac                                                                                                          <  1:1138023/150‑1 (MQ=255)
                                             aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTTGGCGATGTTGTCGATCCTGTAAAAGATTTATAAAgac                                                                                                          >  1:71742/1‑150 (MQ=255)
                                             aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGgac                                                                                                          >  1:3491443/1‑150 (MQ=255)
                                              cAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGTTCCTGGAAAAGATTTATAAGgacg                                                                                                         >  1:4079803/1‑150 (MQ=255)
                                                 cAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGGCGATCCTGGAAAAGATTTATAAGGACGACa                                                                                                      >  1:576740/1‑150 (MQ=255)
                                                  aaTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGCTTTACCTGGGCGATGTTGCCGATCCTGAAAAAAATTTATAAGGACGACaa                                                                                                     >  1:1656874/1‑150 (MQ=255)
                                                   aTGGGATCTGAAATCAGTAAAAAAGATATCGCCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACaaa                                                                                                    <  2:3131112/150‑1 (MQ=255)
                                                      ggATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCg                                                                                                 >  1:2289596/1‑150 (MQ=255)
                                                       gATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCgg                                                                                                >  1:1411898/1‑150 (MQ=255)
                                                       gATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCgg                                                                                                <  1:2170238/150‑1 (MQ=255)
                                                        aTCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTAACTGGGGGGTGTTGCCGGTACTGGAAAAAAATTATTAGGTCGAAAAACCcgg                                                                                               >  2:2831691/1‑150 (MQ=255)
                                                         tCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTTCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGc                                                                                              >  1:3723949/1‑150 (MQ=255)
                                                          cTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCt                                                                                             <  1:1512556/150‑1 (MQ=255)
                                                          cTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCt                                                                                             >  1:389837/1‑150 (MQ=255)
                                                           tGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCtt                                                                                            >  2:553872/1‑150 (MQ=255)
                                                             aaaTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTaa                                                                                          <  2:3781787/150‑1 (MQ=255)
                                                              aaTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAg                                                                                         >  1:2384899/1‑150 (MQ=255)
                                                                 cAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTTCCCGGGCGATGTTTCCGGTCCTTGAAAAGATTTTTAAAGACGACAAAACGGGCTTAAacgc                                                                                      >  2:777138/1‑150 (MQ=255)
                                                                 cAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAgcgc                                                                                      <  2:1138898/150‑1 (MQ=255)
                                                                        aaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATg                                                                               <  1:2863281/150‑1 (MQ=255)
                                                                        aaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGATTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATg                                                                               <  1:2965460/150‑1 (MQ=255)
                                                                          aGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCATGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGaa                                                                             <  1:359597/150‑1 (MQ=255)
                                                                          aGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCCGGCTTAAGCGCGGCGATGaa                                                                             >  1:87228/1‑150 (MQ=255)
                                                                          aGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGGGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGaa                                                                             >  2:2398982/1‑150 (MQ=255)
                                                                             tatCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTTCCCGGGCGATGTTGCCGATCCTGGAAAAGATTTATAAAGACGACAAAACGGGGTTAAACGCGGCGATGAAAGa                                                                          >  2:2565699/1‑150 (MQ=255)
                                                                             tatCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTTTAAGGACGACAAACCGGGCTTAAGCGCGGCGATGaaaaa                                                                          >  1:3171819/1‑148 (MQ=255)
                                                                              atCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGAt                                                                         <  1:1188683/150‑1 (MQ=255)
                                                                               tCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGGGGTGGTTCCGGACCTGGAAAAGATTTTTAAGGACGACAAACCCGGGTTAAACGCGGCGATGAAAGATa                                                                        >  2:1303803/1‑150 (MQ=255)
                                                                                cACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATaa                                                                       <  1:3456645/150‑1 (MQ=255)
                                                                                cACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATaa                                                                       <  1:1272777/150‑1 (MQ=255)
                                                                                  cccGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAAcc                                                                     <  1:2021827/150‑1 (MQ=255)
                                                                                  cccGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGGTGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAAcc                                                                     >  2:688322/1‑150 (MQ=255)
                                                                                   ccGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCt                                                                    <  1:452821/150‑1 (MQ=255)
                                                                                                cgttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAAt                                                       >  2:2780470/1‑150 (MQ=255)
                                                                                                cgttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCCAGCCTGAAAAAAATTTATAAGGACGACAAAACGGGCTTAAACGCGGGGGTGAAAGATAACCTCGAATTTTTTAAt                                                       >  2:1331867/1‑150 (MQ=255)
                                                                                                 gttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATa                                                      >  2:4061487/1‑150 (MQ=255)
                                                                                                     gtcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAAGACGACAAACCCGGCTTAAGCGCGGCGATGAAAGATAACCTCGAAATTATTAATaccca                                                  >  2:2601291/1‑150 (MQ=255)
                                                                                                        gctgctGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATacccaccc                                               >  1:3807218/1‑150 (MQ=255)
                                                                                                         cagctgcAAGCGAGCTTTAACTACGAAAGGATGCAGACGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATAAATACCCACCCg                                              <  2:3393977/148‑1 (MQ=255)
                                                                                                              gcAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCt                                         <  1:2529314/150‑1 (MQ=255)
                                                                                                              gcAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCt                                         <  1:1353546/150‑1 (MQ=255)
                                                                                                              gcAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCt                                         <  1:393640/150‑1 (MQ=255)
                                                                                                              gcAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTCCTGGGGGATGTTGCCGATCCTGAAAAAAATTTTTAAAGGCGACAAAACGGGCTTAAGCGCGGCGATGAAAGATAAACTCGAATTTTTTAATACCCACCCCAATCt                                         >  2:858019/1‑150 (MQ=255)
                                                                                                                 aGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTTTCTGGGCGGGGTTTTCGATCCCGAAAAAAATTTTTAAAGGCGACAAACCGGGGTTAAACGCGGGGGTGGAAAATAAACCCGGAATTATTTATACCCCCCCCaaaatggg                                      >  2:4148551/1‑144 (MQ=255)
                                                                                                                 aGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGt                                      >  1:3694268/1‑150 (MQ=255)
                                                                                                                        tttAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTc                               >  1:2706771/1‑150 (MQ=255)
                                                                                                                        tttAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAAGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCCCGAATTTATTAATACCCACCCGAACCTGGTCGGATTc                               <  2:3671191/150‑1 (MQ=255)
                                                                                                                            aCTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAATATAACCTCGAATTTATTAATACCCGCCCGAATCTGGTCGGATTCCTGa                           <  2:1487581/150‑1 (MQ=255)
                                                                                                                             cTACGAAAGGATGCAGGCGGGCGGGTTTACCTGGGGGGTGTTGCCGGTCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGGAAGATAACCTCGAATTTATTAATAACCACCCGAATCTGGTCGGATTCCTGAt                          >  2:3551172/1‑150 (MQ=255)
                                                                                                                              tACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATg                         <  2:2886884/150‑1 (MQ=255)
                                                                                                                                  aaaGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGt                     <  1:2593561/150‑1 (MQ=255)
                                                                                                                                           aGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAAttt            <  1:2215807/150‑1 (MQ=255)
                                                                                                                                           aGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAAttt            <  1:27780/150‑1 (MQ=255)
                                                                                                                                           aGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAAttt            >  2:908023/1‑150 (MQ=255)
                                                                                                                                            ggctggcgGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTc           <  2:3244159/150‑1 (MQ=255)
                                                                                                                                             gcgggcggTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCg          <  1:553872/150‑1 (MQ=255)
                                                                                                                                                     gTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaag   <  2:1699196/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AAGCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGA  >  AP012306/2758760‑2759058

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