Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,393,3790AC23.9% 823.8 / 61.8 319T92P (ACC→CCC) nthDNA glycosylase and apyrimidinic (AP) lyase
Reads supporting (aligned to +/- strand):  ref base A (214/28);  new base C (0/76);  total (214/105)
Fisher's exact test for biased strand distribution p-value = 0.00e+00
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 368 total aligned reads displayed.

CGACCGATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGCACATT  >  AP012306/1393230‑1393528
                                                                                                                                                     |                                                                                                                                                     
 gACCGATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAc                                                                                                                                                      >  1:2422484/1‑150 (MQ=255)
 gACCGATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAc                                                                                                                                                      >  1:273757/1‑150 (MQ=255)
 gACCGATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAc                                                                                                                                                      >  1:1204923/1‑150 (MQ=255)
 gACCGATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAc                                                                                                                                                      >  2:3207356/1‑150 (MQ=255)
   ccGATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt                                                                                                                                                    >  1:1109219/1‑150 (MQ=255)
   ccGATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt                                                                                                                                                    <  1:171619/150‑1 (MQ=255)
   ccGATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCt                                                                                                                                                    >  1:1413465/1‑150 (MQ=255)
    cGATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTg                                                                                                                                                   <  1:1972129/150‑1 (MQ=255)
     gATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGc                                                                                                                                                  >  2:104265/1‑150 (MQ=255)
     gATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGc                                                                                                                                                  <  2:4020603/150‑1 (MQ=255)
      aTGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGcc                                                                                                                                                 <  1:1547929/150‑1 (MQ=255)
        gTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGt                                                                                                                                               >  2:3056622/1‑150 (MQ=255)
         tCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTa                                                                                                                                              >  2:1095010/1‑150 (MQ=255)
         tCAGTGTTAATAAGGCGACGGCAAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTAAAAACCTATATCAAAACGTTGGGGCTTTATAACAGAAAAGCAAAAAATATCATAAAACCCTGCCGTa                                                                                                                                              <  2:2167343/150‑1 (MQ=255)
          cAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt                                                                                                                                             >  2:2755418/1‑150 (MQ=255)
            gtgtTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCt                                                                                                                                           >  1:3076342/1‑150 (MQ=255)
            gtgtTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCt                                                                                                                                           >  2:258546/1‑150 (MQ=255)
             tgtTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCtt                                                                                                                                          >  1:4426547/1‑150 (MQ=255)
              tttAAAAAGGCGCCGGCAAAACTCTCCCCGGGGGCGAATACGCCTGCAGCGATGTTTGAACTGGGCGTTGAAGGGGTAAAAACCTATATAAAAACGATTGGGTTTTATAACAGCAAAGAAAAAAATATCAAAAAACCCTGCCGTATCTTg                                                                                                                                         <  2:3252153/149‑1 (MQ=255)
               ttAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGc                                                                                                                                        >  2:750609/1‑150 (MQ=255)
                         gACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATTTCATCAAAACCTGCCGTATCTTGCTGGAgcagca                                                                                                                              >  2:1328582/1‑150 (MQ=255)
                            ggCGAACCCCTCCCGGGTGGCAAATCCGCCTGCAGCGATGCTTGAACTGGGCGTTAAAGGGGTAAAAACCTAAATAAAAACGTTGGGGCTTTAAAACAGAAAAAAAAAAAAAATCAACAAACCCTGCCGTATCTTGCTGGAGCAGCATaa                                                                                                                           <  2:3628199/150‑1 (MQ=255)
                             gCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt                                                                                                                          >  2:3409522/1‑150 (MQ=255)
                             gCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTAAAACCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt                                                                                                                          <  1:3708691/150‑1 (MQ=255)
                              cGAAGCTCTCCCCGGTGGCGAATCCGCCTGCGGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAACCCTATATCAAAACGATTGGGCTTAATAACAGCAAAGCAGAAAATATCATCAAAGCCTGCCGTATCTTGCTGGAGCAGCATAATg                                                                                                                         <  1:3407784/150‑1 (MQ=255)
                              cGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATg                                                                                                                         >  2:1941999/1‑150 (MQ=255)
                               gAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATgg                                                                                                                        >  1:3129012/1‑150 (MQ=255)
                                    tctACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg                                                                                                                   >  1:1062146/1‑150 (MQ=255)
                                     ctACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGt                                                                                                                  >  2:3879323/1‑150 (MQ=255)
                                           ggTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGa                                                                                                            >  1:1653833/1‑150 (MQ=255)
                                               gcaaaTCCGCCTGCAGCGATGCTTGAACTGGGCGTTGAGGGGGTAAAACCCTATATAAAAACGATTGGGTTTTATAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGAt                                                                                                        <  2:288980/147‑1 (MQ=255)
                                                cGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATc                                                                                                       >  1:1690984/1‑150 (MQ=255)
                                                       gCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgctgc                                                                                                >  1:3160854/1‑150 (MQ=255)
                                                         cTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTgcgc                                                                                              >  2:4315695/1‑150 (MQ=255)
                                                                    cTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg                                                                                   >  2:210036/1‑150 (MQ=255)
                                                                     ttGAACGGGGCGTTGAAGGGGTAAAACCCTATACAAAACCGTTTGGGTTTTATAAAACAAAAGAAAAAAATATCATAAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGc                                                                                  <  2:4461686/150‑1 (MQ=255)
                                                                           ggggCGTTGAGGGGGTAAAACCCTATATCAAAACGATTGGGCTTTATAACAGAAAAGCAGAAAATATCAAAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCg                                                                            <  2:4516720/149‑1 (MQ=255)
                                                                             ggCGTTGAAGGGGTAAAAACCTATATAAAAACGATTGGGTTTTATAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTa                                                                          <  2:1803004/150‑1 (MQ=255)
                                                                               cGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAgg                                                                        >  2:2733337/1‑150 (MQ=255)
                                                                                gTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt                                                                       >  1:155599/1‑150 (MQ=255)
                                                                                 ttGAGGGGGTGAAACCCTATATCAAAACGATTGGGTTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc                                                                      <  2:1866297/150‑1 (MQ=255)
                                                                                 ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc                                                                      >  2:1620075/1‑150 (MQ=255)
                                                                                 ttGAAGGGGTAAAAACCTATATAAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc                                                                      <  1:4286875/150‑1 (MQ=255)
                                                                                  tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCg                                                                     >  1:2015033/1‑150 (MQ=255)
                                                                                    aaGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTa                                                                   >  1:3220871/1‑150 (MQ=255)
                                                                                        ggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAAc                                                               >  2:1508886/1‑150 (MQ=255)
                                                                                            aaaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGcc                                                           >  1:2527255/1‑150 (MQ=255)
                                                                                             aaaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCa                                                          >  1:871451/1‑150 (MQ=255)
                                                                                              aaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCaa                                                         >  1:528574/1‑150 (MQ=255)
                                                                                               aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAc                                                        >  1:4235259/1‑150 (MQ=255)
                                                                                                ccTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTCCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcg                                                       >  2:61525/1‑150 (MQ=255)
                                                                                                 cTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcgt                                                      <  1:1141497/150‑1 (MQ=255)
                                                                                                 cTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcgt                                                      >  2:3483497/1‑150 (MQ=255)
                                                                                                 cTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcgt                                                      >  1:2951056/1‑150 (MQ=255)
                                                                                                  tatataAAAACGATTGGGTTTTATAACAGCAAAGCAGAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcgtc                                                     <  1:3380254/150‑1 (MQ=255)
                                                                                                  tataTCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcgtc                                                     >  2:698836/1‑150 (MQ=255)
                                                                                                        aaaaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa                                               >  1:4111151/1‑150 (MQ=255)
                                                                                                          aaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTaaa                                             >  1:2137952/1‑150 (MQ=255)
                                                                                                          aaCGATGGGTTTTTAAAACAGAAAACGAAAAAATATCACAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTaaa                                             <  2:3324160/150‑1 (MQ=255)
                                                                                                           aCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAac                                            >  1:1860008/1‑150 (MQ=255)
                                                                                                             gATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAacac                                          >  2:1771152/1‑150 (MQ=255)
                                                                                                              aTTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACt                                         >  2:2263733/1‑150 (MQ=255)
                                                                                                                gggtttttAAAACAAAAAACCAAAAAATTTCACCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGc                                       <  2:3372842/145‑1 (MQ=255)
                                                                                                                    cTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTc                                   >  1:279532/1‑150 (MQ=255)
                                                                                                                     tttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCg                                  >  2:3705241/1‑150 (MQ=255)
                                                                                                                      ttATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCgg                                 >  2:2617739/1‑150 (MQ=255)
                                                                                                                      ttAAAACAGAAAAGAAAAAAATATCATAAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGATGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCgg                                 <  2:4147038/150‑1 (MQ=255)
                                                                                                                       tataACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGc                                >  1:3488683/1‑150 (MQ=255)
                                                                                                                       tataACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGc                                >  2:2027297/1‑150 (MQ=255)
                                                                                                                          aaCAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTgg                             >  1:4013290/1‑150 (MQ=255)
                                                                                                                          aaCAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTgg                             >  2:3571610/1‑150 (MQ=255)
                                                                                                                           aCAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGc                            >  2:312411/1‑150 (MQ=255)
                                                                                                                           aCAGCAAAGAAAAAAATATCACAAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGc                            <  2:1207368/150‑1 (MQ=255)
                                                                                                                            cAGAAAAGAAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGcc                           <  2:2943089/150‑1 (MQ=255)
                                                                                                                             aGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCg                          >  2:1663589/1‑150 (MQ=255)
                                                                                                                                 aagaaaaaaaTATCAACAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTa                      <  2:2313370/144‑1 (MQ=255)
                                                                                                                                  aGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTAt                     >  1:3588007/1‑150 (MQ=255)
                                                                                                                                    cAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTg                   >  2:4389364/1‑150 (MQ=255)
                                                                                                                                    cAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTg                   >  2:1427022/1‑150 (MQ=255)
                                                                                                                                     aGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGc                  >  2:856357/1‑150 (MQ=255)
                                                                                                                                     aGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGc                  >  1:2640674/1‑150 (MQ=255)
                                                                                                                                       aaaaTATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTg                <  1:285464/150‑1 (MQ=255)
                                                                                                                                       aaaaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTg                >  2:3782357/1‑150 (MQ=255)
                                                                                                                                        aaaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGt               >  2:2814456/1‑150 (MQ=255)
                                                                                                                                        aaaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGt               >  2:2122300/1‑150 (MQ=255)
                                                                                                                                        aaaTATCATCAAAACCTGCCGAATCTTGCTGGAGCAGCATAATGGCGAGGTTTCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGt               >  1:4192448/1‑150 (MQ=255)
                                                                                                                                           tatCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGa            >  2:2472118/1‑150 (MQ=255)
                                                                                                                                           tatCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGa            >  2:3304/1‑150 (MQ=255)
                                                                                                                                           tatCATCAAAACATGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGa            >  2:4592318/1‑150 (MQ=255)
                                                                                                                                            atcatcAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGac           >  1:3907062/1‑150 (MQ=255)
                                                                                                                                              catcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGacac         <  1:4542258/150‑1 (MQ=255)
                                                                                                                                               atcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACg        >  1:799864/1‑150 (MQ=255)
                                                                                                                                               atcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGCAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACg        >  1:1314751/1‑150 (MQ=255)
                                                                                                                                                tcaAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGc       <  2:4341276/150‑1 (MQ=255)
                                                                                                                                                tcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGc       >  1:3815902/1‑150 (MQ=255)
                                                                                                                                                 caAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGca      >  2:4083270/1‑150 (MQ=255)
                                                                                                                                                  aaaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGcac     >  1:2583092/1‑150 (MQ=255)
                                                                                                                                                   aaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGcaca    >  1:739509/1‑150 (MQ=255)
                                                                                                                                                   aaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGcaca    >  1:2289235/1‑150 (MQ=255)
                                                                                                                                                     aCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGCACAtt  >  2:2405557/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGACCGATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGCACATT  >  AP012306/1393230‑1393528

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