Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123061,393,3790AC23.7% 752.9 / 40.4 281T92P (ACC→CCC) nthDNA glycosylase and apyrimidinic (AP) lyase
Reads supporting (aligned to +/- strand):  ref base A (180/33);  new base C (0/66);  total (180/101)
Fisher's exact test for biased strand distribution p-value = 1.81e-39
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 318 total aligned reads displayed.

CGACCGATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGCACATT  >  AP012306/1393230‑1393528
                                                                                                                                                       |                                                                                                                                                     
 gACCGATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAc                                                                                                                                                      >  1:2539701/1‑150 (MQ=255)
  aCCGATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAAcc                                                                                                                                                     <  1:2104098/150‑1 (MQ=255)
     gATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGc                                                                                                                                                  >  1:3079126/1‑150 (MQ=255)
      aTGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAGAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGcc                                                                                                                                                 <  2:2842406/150‑1 (MQ=255)
      aTGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGcc                                                                                                                                                 <  1:229851/150‑1 (MQ=255)
          cAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt                                                                                                                                             >  1:3401050/1‑150 (MQ=255)
          cAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTAt                                                                                                                                             <  2:3455604/150‑1 (MQ=255)
           aGTTTTAATAAGGCGACGGCGAAACTCCACCCGGGGCCAAACACTCCCGCAGGGATTTTTAACGTGGGCGTGGAAGGGGTGAAAACTTTAAAAAAAACATTGGGTTTTTAAAACAGCAAAGCAGAAAATATCAACAAACCCTGCCGTATc                                                                                                                                            <  2:1288173/150‑1 (MQ=255)
           aGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAAAACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGCGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAACGCAGAAAATATCATCAAAACCTGCCGTATc                                                                                                                                            >  2:3966456/1‑150 (MQ=255)
               ttAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGc                                                                                                                                        >  1:1069754/1‑150 (MQ=255)
                     aGGCGACGGCGAACCTCTCCCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGagc                                                                                                                                  <  1:1751638/150‑1 (MQ=255)
                     aGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGagc                                                                                                                                  >  2:1689089/1‑150 (MQ=255)
                     aGGCGACGGCGAAACTCTACCCGGGGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTGGAAGGGGTGAAAACCTATATCAAAACGATGGGGTTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGagc                                                                                                                                  <  2:1393289/150‑1 (MQ=255)
                       gcggcggcgAATCTCCCCCAGGGGCGGAACACGCCTGCAGCGTTGCTTGACCTGGGCGTTAAGGGGTAAAAACCCGAAAACAAACCGATTGGCCTTTAAAACAGCAAGGGAGAAAATATCAAAAAACCCTGCCGTATCTTGCTGGagcag                                                                                                                                <  2:828953/150‑1 (MQ=255)
                       gCGAGGGCGAACCTCTCCCCGGTGGCAAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTAAAACCTTATATAAAAACGTTTGGTTTTTATAACAGCAAAGAAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGagcag                                                                                                                                <  2:2984474/150‑1 (MQ=255)
                          aCGGGGAAACTATCCCCGGGGGCGAATACGCCTGAAGGGATGCTTGAACTGGGGGTTGAAGGGGTGAAAACCTATAAAAAAACGATTGGGTTTTATAAAAGAAAAGAAGAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCAt                                                                                                                             <  2:1080823/150‑1 (MQ=255)
                            ggCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATaa                                                                                                                           >  2:1828106/1‑150 (MQ=255)
                            ggCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATaa                                                                                                                           >  2:3999027/1‑150 (MQ=255)
                             gCGAACCTCTCCCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAACCCTATATAAAAACGATTGGGCTTTATAACAGAAAAGCAAAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAAt                                                                                                                          <  1:933063/150‑1 (MQ=255)
                               gAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATgg                                                                                                                        >  2:1989681/1‑150 (MQ=255)
                                 aaCTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCg                                                                                                                      >  2:2498402/1‑150 (MQ=255)
                                 aaCTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCg                                                                                                                      >  1:3789704/1‑150 (MQ=255)
                                 aaCTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCg                                                                                                                      >  2:1906561/1‑150 (MQ=255)
                                  aCTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGa                                                                                                                     >  2:3305771/1‑150 (MQ=255)
                                  aCTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGa                                                                                                                     >  1:3528596/1‑150 (MQ=255)
                                    tctACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAgg                                                                                                                   >  1:3732945/1‑150 (MQ=255)
                                      tACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGtt                                                                                                                 >  2:1492786/1‑150 (MQ=255)
                                             tGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAg                                                                                                          <  1:1164811/150‑1 (MQ=255)
                                              ggCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGa                                                                                                         >  1:1238942/1‑150 (MQ=255)
                                                  aaTCCGCCTGAAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAACCCTATAACAAAACGATTGGGTTTTATAACAGCAAAGAAGAAAATATCAAAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGt                                                                                                     <  2:1370613/150‑1 (MQ=255)
                                                     aCGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGtgct                                                                                                  >  1:2913006/1‑150 (MQ=255)
                                                                gATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGc                                                                                       >  2:3753087/1‑150 (MQ=255)
                                                                 aTGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGcc                                                                                      <  1:1338347/150‑1 (MQ=255)
                                                                  tGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGccc                                                                                     >  1:3115089/1‑150 (MQ=255)
                                                                   gCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCt                                                                                    >  1:1430345/1‑150 (MQ=255)
                                                                    cTTGAACTGGGCGTTGAAGGGGTGAAACCTTATATAAAAACGTTGGGGCTTTAAAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTg                                                                                   <  2:2373036/150‑1 (MQ=255)
                                                                     ttGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGc                                                                                  >  1:2621309/1‑150 (MQ=255)
                                                                      tGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGcc                                                                                 >  2:1642114/1‑150 (MQ=255)
                                                                                gtggAAGGGGTAAAAACTTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGt                                                                       <  1:4020739/147‑1 (MQ=255)
                                                                                 ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc                                                                      >  1:1303214/1‑150 (MQ=255)
                                                                                 ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc                                                                      >  1:250236/1‑150 (MQ=255)
                                                                                 ttGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTc                                                                      >  2:2416009/1‑150 (MQ=255)
                                                                                  tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCg                                                                     >  1:3604829/1‑150 (MQ=255)
                                                                                  tGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCg                                                                     >  2:3511073/1‑150 (MQ=255)
                                                                                   gAAGGGGTAAAAACCTATATCAAAACGATTGGGCTTTAAAACAGCAAGGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGt                                                                    <  1:3334108/150‑1 (MQ=255)
                                                                                      ggggTGAAACCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTaaa                                                                 <  1:3013380/150‑1 (MQ=255)
                                                                                       gggTGAAACCCTATATCAAAACGATTGGGCTTTATAACAGAAAAGAAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTaaaa                                                                <  1:2778260/150‑1 (MQ=255)
                                                                                       gggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTaaaa                                                                >  1:2281544/1‑150 (MQ=255)
                                                                                       gggTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTaaaa                                                                >  2:2256347/1‑150 (MQ=255)
                                                                                         gTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACa                                                              >  1:992021/1‑150 (MQ=255)
                                                                                          tGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAg                                                             >  2:872669/1‑150 (MQ=255)
                                                                                           gAAACCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGc                                                            <  1:3841093/150‑1 (MQ=255)
                                                                                              aaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCaa                                                         >  2:906249/1‑150 (MQ=255)
                                                                                              aaCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCaa                                                         >  1:295334/1‑150 (MQ=255)
                                                                                               aCCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAc                                                        >  2:115601/1‑150 (MQ=255)
                                                                                                cccTATATCAAACCGATTGGGTTTTATAACAGAAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAc                                                        <  1:463545/149‑1 (MQ=255)
                                                                                                  cTATATCAAACCGATTGGGCTTTATAACAGAAAAGCAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcgt                                                      <  1:3441109/150‑1 (MQ=255)
                                                                                                   tatataAAACCGATTGGGTTTTATAACACAAAAGAAAAAAATACCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcgtc                                                     <  2:2799660/150‑1 (MQ=255)
                                                                                                   tatataAAAACGATTGGGTTTTATAACAGAAAAGAAAAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcgtc                                                     <  2:1195251/150‑1 (MQ=255)
                                                                                                   tataTCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcgtc                                                     >  1:1891086/1‑150 (MQ=255)
                                                                                                   aatatAAAAACAATTGGTTTTTATAACAGCAAAGAAAAAAATATCACAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcgtc                                                     <  2:1596743/149‑1 (MQ=255)
                                                                                                    atatCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcgtcg                                                    >  2:2162180/1‑150 (MQ=255)
                                                                                                    aaaaaaaaaCGTTTGGTTTTTAAAACAGAAAAGAAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAAcgtcg                                                    <  2:3055367/145‑1 (MQ=255)
                                                                                                       tCAAAACGATTGGGCTTTATAACAGCAAAGCAAAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTAt                                                 <  1:2154975/150‑1 (MQ=255)
                                                                                                         aaaaCGATTGGTCTTTAAAAAAGCAAAGAAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa                                               <  2:1107559/150‑1 (MQ=255)
                                                                                                         aaaaCGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTa                                               >  1:1199318/1‑150 (MQ=255)
                                                                                                            aCGTTTGGGTTTTATAACAGAAAAGAAGAAAATATCAACAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAac                                            <  2:1728864/150‑1 (MQ=255)
                                                                                                               aTTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACt                                         >  2:641633/1‑150 (MQ=255)
                                                                                                               aTTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACt                                         >  1:751722/1‑150 (MQ=255)
                                                                                                               aTTGGCTTTTAAAACAGCAAAGCAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACt                                         <  2:759727/150‑1 (MQ=255)
                                                                                                                ttGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTg                                        >  1:645850/1‑150 (MQ=255)
                                                                                                                 tgggtcTTTAAAACAGAAAAGAAAAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTg                                        <  2:3772006/145‑1 (MQ=255)
                                                                                                                 tGGGTTTTATAACAGCAAAGCAAAAAATATCAACAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGc                                       <  2:4077023/150‑1 (MQ=255)
                                                                                                                    gCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCAtt                                    >  1:1052611/1‑150 (MQ=255)
                                                                                                                      cTTTATAACAGCAAAGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTc                                   <  1:1568388/150‑1 (MQ=255)
                                                                                                                          ataACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCt                               >  1:875731/1‑150 (MQ=255)
                                                                                                                          ataACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCt                               >  2:2624472/1‑150 (MQ=255)
                                                                                                                          aaaaCAGAAAAGAAAAAAATATCATAAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCt                               <  1:1701775/148‑1 (MQ=255)
                                                                                                                           aaaCAGCAAAGCAGAAAATATCATAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTg                              <  2:2457150/149‑1 (MQ=255)
                                                                                                                            aaCAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTgg                             >  1:5127/1‑150 (MQ=255)
                                                                                                                              cAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGcc                           >  1:1415132/1‑150 (MQ=255)
                                                                                                                              cAGAAAAGAAAAAAATATCAACAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGcc                           <  2:2285110/150‑1 (MQ=255)
                                                                                                                               aGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCg                          >  1:2371531/1‑150 (MQ=255)
                                                                                                                               aGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCg                          >  1:2987584/1‑150 (MQ=255)
                                                                                                                                 cAAGGCAGAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGAc                        <  1:1492507/150‑1 (MQ=255)
                                                                                                                                    gaaaGAAAATATCATAAAAGCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTAt                     <  2:2469090/147‑1 (MQ=255)
                                                                                                                                        gAAAATATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCt                 <  1:556061/150‑1 (MQ=255)
                                                                                                                                         aaaaTATAACAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTg                <  2:3032531/150‑1 (MQ=255)
                                                                                                                                           aaTATCATCAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTc              <  1:719731/150‑1 (MQ=255)
                                                                                                                                           aaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTc              >  1:2601963/1‑150 (MQ=255)
                                                                                                                                           aaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTc              >  1:3560456/1‑150 (MQ=255)
                                                                                                                                           aaTATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTc              >  1:1975822/1‑150 (MQ=255)
                                                                                                                                           aaTATCACAAAACCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGTAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTc              <  1:2074371/150‑1 (MQ=255)
                                                                                                                                                 atcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACg        <  1:1643420/150‑1 (MQ=255)
                                                                                                                                                 atcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACg        >  1:1279014/1‑150 (MQ=255)
                                                                                                                                                  tcaAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGc       >  1:473016/1‑150 (MQ=255)
                                                                                                                                                    aaacCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGcac     <  2:3650203/150‑1 (MQ=255)
                                                                                                                                                    aaaaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGcac     >  1:3268646/1‑150 (MQ=255)
                                                                                                                                                      aaCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGCACAt   >  2:3243803/1‑150 (MQ=255)
                                                                                                                                                       aCCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGCACAtt  >  2:534322/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CGACCGATGTCAGTGTTAATAAGGCGACGGCGAAACTCTACCCGGTGGCGAATACGCCTGCAGCGATGCTTGAACTGGGCGTTGAAGGGGTGAAAACCTATATCAAAACGATTGGGCTTTATAACAGCAAAGCAGAAAATATCATCAAAACCTGCCGTATCTTGCTGGAGCAGCATAATGGCGAGGTTCCGGAAGATCGTGCTGCGCTTGAAGCCCTGCCCGGCGTAGGTCGTAAAACAGCCAACGTCGTATTAAACACTGCATTCGGCTGGCCGACTATTGCTGTCGACACGCACATT  >  AP012306/1393230‑1393528

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