Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123062,758,9090TG23.7% 576.0 / 16.0 211G33G (GGT→GGG) agaDN‑acetylgalactosamine‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (30/131);  new base G (50/0);  total (80/131)
Fisher's exact test for biased strand distribution p-value = 9.17e-28
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 286 total aligned reads displayed.

GCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGA  >  AP012306/2758762‑2759058
                                                                                                                                                    |                                                                                                                                                     
gcgcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGttt                                                                                                                                                     >  2:3465481/1‑150 (MQ=255)
 cgcAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGgttt                                                                                                                                                    >  1:4206404/1‑150 (MQ=255)
    aaCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTAcc                                                                                                                                                 <  1:2236710/150‑1 (MQ=255)
    aaCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTAcc                                                                                                                                                 >  2:1301620/1‑150 (MQ=255)
     aCCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCt                                                                                                                                                >  2:3791917/1‑150 (MQ=255)
       cGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCCgg                                                                                                                                              >  2:1415666/1‑150 (MQ=255)
        gCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTggg                                                                                                                                             <  2:144444/150‑1 (MQ=255)
            cccGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGAt                                                                                                                                         <  2:329128/150‑1 (MQ=255)
             ccGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATg                                                                                                                                        <  1:1387065/150‑1 (MQ=255)
             ccGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGGGATg                                                                                                                                        >  2:2806261/1‑150 (MQ=255)
              cGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTTCCTGGGCGATGt                                                                                                                                       >  1:2969595/1‑150 (MQ=255)
                ttGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTACTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTg                                                                                                                                     <  2:3088983/150‑1 (MQ=255)
                 tGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGc                                                                                                                                    <  1:572930/150‑1 (MQ=255)
                 tGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGGTGTTGc                                                                                                                                    >  1:3072577/1‑150 (MQ=255)
                  gCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGcc                                                                                                                                   <  1:4313623/150‑1 (MQ=255)
                  gCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGGTGtttgc                                                                                                                                   >  1:969167/1‑147 (MQ=255)
                     aGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGGTGTTGCCCAt                                                                                                                                >  1:1482063/1‑150 (MQ=255)
                     aGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGAt                                                                                                                                <  1:1448315/150‑1 (MQ=255)
                      gTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTTCCGATc                                                                                                                               >  2:1762639/1‑150 (MQ=255)
                           aaTGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGATGTTGCCGATCCTGaa                                                                                                                          >  2:2211252/1‑150 (MQ=255)
                               aagaagAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAg                                                                                                                      >  1:4146191/1‑150 (MQ=255)
                                 gaagaGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGG‑GGGTTTACCTGGGCGATGTTGCCGAGCCTGAAAAAGAt                                                                                                                    >  1:299626/1‑150 (MQ=255)
                                   agagGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTTCCGATCCTGAAAAAGAatt                                                                                                                  >  1:1034741/1‑150 (MQ=255)
                                      ggACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTtata                                                                                                               <  1:2127349/150‑1 (MQ=255)
                                      ggACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTtata                                                                                                               >  1:2820346/1‑150 (MQ=255)
                                      ggACTACAGCAAAGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGGTGTTGCCCGACCTGTAAAAGATTtata                                                                                                               >  2:3184089/1‑150 (MQ=255)
                                        aCTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGGCGATCCTGAAAAAGATTTATAAg                                                                                                             >  2:3688899/1‑150 (MQ=255)
                                        aCTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTTGGCGATGTTTTCCGTCCTGGAAAAGATTTATaaa                                                                                                             >  2:438488/1‑149 (MQ=255)
                                           aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGgac                                                                                                          <  1:2778092/150‑1 (MQ=255)
                                           aCAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGgac                                                                                                          <  1:3723858/150‑1 (MQ=255)
                                             aGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGgacga                                                                                                        <  2:2997785/150‑1 (MQ=255)
                                             aGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGgacga                                                                                                        >  1:1279936/1‑150 (MQ=255)
                                             aGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCCGGGCGATGTTGCCGATCCTGAAAAAAATTTATAAAgacga                                                                                                        >  2:3028874/1‑150 (MQ=255)
                                              gCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGgacgac                                                                                                       <  1:1162949/150‑1 (MQ=255)
                                               cAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCCGGGCGATGTTGCCGATCCTGCAAAAAATTTATAAAGACGACa                                                                                                      >  1:325844/1‑150 (MQ=255)
                                                aaTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAAATTTTTAAGGACGACaa                                                                                                     >  2:3385783/1‑150 (MQ=255)
                                                 aTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACaaa                                                                                                    <  2:1804770/150‑1 (MQ=255)
                                                 aTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCCATCCTGAAAAAGATTTATAAAGACGGCaaa                                                                                                    >  1:2923637/1‑150 (MQ=255)
                                                 aTGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGGTGTTGCCCATCCTGAAAAAAATTTTTAAGGAAGACaaa                                                                                                    >  2:2370204/1‑150 (MQ=255)
                                                    ggATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTTCCTGGGCGATGTTGCCGATCCTTAAAAAGATTTATTAAGACGACAAACCg                                                                                                 >  1:4134913/1‑150 (MQ=255)
                                                    ggATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTTACTGGGCGGTGTTTGCGGTACTTAAAAAAAATTTTAAGGACGCCaaaacc                                                                                                 >  2:1856722/1‑147 (MQ=255)
                                                    ggATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTGCCGGTCCTGAAAAAGATTTTTTAGGACGACAAACCg                                                                                                 >  1:1566181/1‑150 (MQ=255)
                                                     gATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGG‑GGGTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAAGACGACAAAcccg                                                                                                >  2:910612/1‑148 (MQ=255)
                                                      aTCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGG‑GGGTTTAACTGGGGGAAGTTTCCGATACTGCAAAAGATTTATAAAGACGACAAACCcgg                                                                                               >  2:815835/1‑150 (MQ=255)
                                                       tCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGAGGCAGGCGGGC‑GGGGTTGCCAGGGCGGTGAGGTCGGGCTTGAAGAAGAACTATAAAGAGGATGAACCgggg                                                                                              >  2:879688/1‑149 (MQ=255)
                                                        cTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGATGTTTCCGGTCCTGAAAAAGAATTATAAGGAAGACAAAACCGGCt                                                                                             >  2:1378686/1‑150 (MQ=255)
                                                         tGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCtt                                                                                            <  1:3757759/150‑1 (MQ=255)
                                                           aaaTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCTCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCTGGCTGGCGGGTTTTATC‑AGGCGATTTTATCGTTGCTGAAAAAGAATTTTAAGGGGGAAAAACCtgggtcta                                                                                          >  2:4270261/1‑145 (MQ=255)
                                                            aaTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAAGACGACAAACCGGGCTTAAg                                                                                         >  1:2442268/1‑150 (MQ=255)
                                                             aTCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAgc                                                                                        <  1:2648271/150‑1 (MQ=255)
                                                              tCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAgcg                                                                                       <  1:4412212/150‑1 (MQ=255)
                                                                  tAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCgcggc                                                                                   <  1:2031352/150‑1 (MQ=255)
                                                                    aaaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGG‑GGGTTTACCTGGGGGGTGTTGCCGACCCCGGAAAAGATTTTTAAGGGCGAAAAAACGGGGTTTAGCgcggcgg                                                                                 >  2:554950/1‑149 (MQ=255)
                                                                     aaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGAt                                                                                <  1:1663538/150‑1 (MQ=255)
                                                                     aaaaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGAt                                                                                <  1:367172/150‑1 (MQ=255)
                                                                       aaGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTCCCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGa                                                                              <  2:2131209/150‑1 (MQ=255)
                                                                        aGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGaa                                                                             >  1:550233/1‑150 (MQ=255)
                                                                         gATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGaaa                                                                            <  1:1016085/150‑1 (MQ=255)
                                                                           tatCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGa                                                                          <  1:3880537/150‑1 (MQ=255)
                                                                            atCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTCGGCGATGTTGCCGGTCCTGAAAAAAATTTATAAAGACGACAAACCGGGCTTAAACGCGGCGATGGAAGaa                                                                         >  2:1297222/1‑149 (MQ=255)
                                                                               aCCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTTGGCGATGTTGGCGATCCTGAAAAAGATTTATAAAGACGCCAAAACGGGCCTTAACGCGGCGGTGAAAGATAAc                                                                      >  2:1924834/1‑150 (MQ=255)
                                                                                     cTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTAACTGGGCGGGGTTTCCGATCCTGAAAAAGATTTATTAAGACGACAAAACGGGCCTTAACGCGGCGATGAAAGATAAAATCGaa                                                                >  2:1970678/1‑150 (MQ=255)
                                                                                         gCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAAATTTATAAGGACGACAAACCCGGCTTAAGCGCGGCGATGAAAGATAAACTCGAATtta                                                            >  1:3947030/1‑150 (MQ=255)
                                                                                           tttCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAAGACGACAAACCGGGCTTAAGCGCGGCGATGAAAAATAACCTCGAATttatt                                                          >  2:4034475/1‑150 (MQ=255)
                                                                                            ttcgttcgTCGCTGCTGCAAGCTAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGGAAGATAACCTCCAATttatta                                                         >  1:1640504/1‑150 (MQ=255)
                                                                                            ttcgttcgTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATttatta                                                         <  2:1095933/150‑1 (MQ=255)
                                                                                              cgttcgtCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAAt                                                       <  1:4092674/150‑1 (MQ=255)
                                                                                                   gtcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGAGCGGCGGTGAAAGATAACCTCGAATTTATTAATaccca                                                  >  2:2581745/1‑150 (MQ=255)
                                                                                                   gtcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGGTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGGTGAAAGATAACCTCGAATTTATTAATaccca                                                  >  2:3874288/1‑150 (MQ=255)
                                                                                                    tcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATacccac                                                 <  1:2366596/150‑1 (MQ=255)
                                                                                                      gctgctGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATacccaccc                                               <  2:2442561/150‑1 (MQ=255)
                                                                                                          ctgcAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAAt                                           <  1:2681320/150‑1 (MQ=255)
                                                                                                              aaGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGGTGTTGCCGATCCTGAAAAAGATTTTTAAGGACGACAAACCGGGCTTAAGCGCGGGGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTgg                                       >  2:1586208/1‑150 (MQ=255)
                                                                                                              aaGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGGTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAGATTTTTAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTgg                                       >  2:2257463/1‑150 (MQ=255)
                                                                                                              aaGCGAGCTTTAACTACGAAAGGATGCAGGC‑GGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGt                                      >  1:3227219/1‑150 (MQ=255)
                                                                                                                   aGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGa                                  <  1:68367/150‑1 (MQ=255)
                                                                                                                    gCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGAt                                 <  1:257160/150‑1 (MQ=255)
                                                                                                                     cTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGAtt                                <  1:1716547/150‑1 (MQ=255)
                                                                                                                        tAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCt                             <  1:536982/150‑1 (MQ=255)
                                                                                                                          aCTACGAAAGGATGCAGGCGGGG‑GGGTTTTCCTAGGCGGTGGTGCCAAACCTGAAAAAGATTTATAAGGACGCCAAACCGGGATTAAGCGCGGGGGTGAAAGATAAACTCGGGTTTTTTAATACCCACACGGATTTGGTCGGGTTCCTgt                           >  2:4098426/1‑149 (MQ=255)
                                                                                                                           cTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGAt                          <  1:658644/150‑1 (MQ=255)
                                                                                                                            tACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATg                         <  1:3404120/150‑1 (MQ=255)
                                                                                                                             aCGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATgg                        <  1:1100222/150‑1 (MQ=255)
                                                                                                                             aCGAAAGGATGCAGGCGGGC‑GGGTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGGAGACAAACCGGGCTTAAACGCGGGGATGGAAGATAACCCCGAATTTTTTTATTACCCCCCGAATCTGGTTGGATTCCTGATgg                        >  2:683379/1‑150 (MQ=255)
                                                                                                                              cGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATggg                       <  1:3775161/150‑1 (MQ=255)
                                                                                                                                 aaGGATGCAGGCGGGC‑GGTTTTTCCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGCACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATTCCCCCCCGAATCTGGTCGGATTCCTGATGGGGtt                    <  2:1232550/150‑1 (MQ=255)
                                                                                                                                 aaGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGtt                    <  1:1928441/150‑1 (MQ=255)
                                                                                                                                 aaGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGtt                    <  1:1678551/150‑1 (MQ=255)
                                                                                                                                     aTGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGttatta                <  2:1004346/150‑1 (MQ=255)
                                                                                                                                       gCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAAt              <  1:1364096/150‑1 (MQ=255)
                                                                                                                                             gggc‑ggTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGAt        <  1:3354954/150‑1 (MQ=255)
                                                                                                                                              ggc‑ggTTTTACCTGGGCGATGTTGCAGATCCTGAAAAAGATTTATAAGGACGCCAAACCGGGCTTAAGCGCGGCGATGAAAGATACCCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATg       <  2:1019157/150‑1 (MQ=255)
                                                                                                                                               gc‑ggTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATgg      <  2:292942/150‑1 (MQ=255)
                                                                                                                                                c‑ggTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGga     <  1:1876752/150‑1 (MQ=255)
                                                                                                                                                  ggTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaa    <  2:75572/150‑1 (MQ=255)
                                                                                                                                                  ggTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaa    >  1:2735934/1‑150 (MQ=255)
                                                                                                                                                   gTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaag   >  2:238730/1‑150 (MQ=255)
                                                                                                                                                    ttttACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaaga  <  2:461832/150‑1 (MQ=255)
                                                                                                                                                    ttttACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaaga  >  2:3842658/1‑150 (MQ=255)
                                                                                                                                                    ttttACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGgaaga  >  2:2656911/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGCAACCGCAGCCCGTTGCCAGTAAAAATGAAGAAGAGGACTACAGCAATGGGATCTGAAATCAGTAAAAAAGATATCACCCGTCTGGGCTTTCGTTCGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGC‑GGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGA  >  AP012306/2758762‑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: ‑