Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 565,299 T→G 9.7% N314K (AAT→AAG ybcK → DLP12 prophage; putative recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273565,2990TG 9.7% 1009.6 / 93.0 357N314K (AAT→AAGybcKDLP12 prophage; putative recombinase
Reads supporting (aligned to +/- strand):  ref base T (157/165);  new base G (22/13);  total (179/178)
Fisher's exact test for biased strand distribution p-value = 1.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 362 total aligned reads displayed.

GCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAAT  >  CP009273/565150‑565448
                                                                                                                                                     |                                                                                                                                                     
  tGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAAGAc                                                                                                                                                     >  2:2974158/1‑150 (MQ=255)
     cTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCAt                                                                                                                                                  >  2:2881012/1‑150 (MQ=255)
      tattatCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATg                                                                                                                                                 >  2:1348106/1‑150 (MQ=255)
      tattatCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATg                                                                                                                                                 <  1:262339/150‑1 (MQ=255)
       attatCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGa                                                                                                                                                >  2:54493/1‑150 (MQ=255)
         tatCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGAtt                                                                                                                                              >  1:4217869/1‑150 (MQ=255)
         tatCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGAtt                                                                                                                                              <  1:4148765/150‑1 (MQ=255)
          atCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTg                                                                                                                                             >  2:851615/1‑150 (MQ=255)
              cAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCa                                                                                                                                         >  1:1007214/1‑150 (MQ=255)
              cAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTGATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCa                                                                                                                                         <  2:4409963/150‑1 (MQ=255)
               agagTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCAt                                                                                                                                        <  1:2852231/150‑1 (MQ=255)
                 agTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGc                                                                                                                                      >  1:3457193/1‑150 (MQ=255)
                 agTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGc                                                                                                                                      <  2:202298/150‑1 (MQ=255)
                 agTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGc                                                                                                                                      <  2:3159542/150‑1 (MQ=255)
                  gTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCg                                                                                                                                     <  2:3310555/150‑1 (MQ=255)
                   tCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCgg                                                                                                                                    <  2:1803688/150‑1 (MQ=255)
                    cATATCAGATGATTTGTTTTACGCTGTACAGGAAAGTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGt                                                                                                                                   <  1:1827181/150‑1 (MQ=255)
                     atatCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTa                                                                                                                                  <  1:4376912/150‑1 (MQ=255)
                     atatCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTa                                                                                                                                  <  1:325066/150‑1 (MQ=255)
                         cAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTg                                                                                                                              >  2:2099299/1‑150 (MQ=255)
                          aGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACCGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTgg                                                                                                                             <  2:1522418/150‑1 (MQ=255)
                            atgatTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGaa                                                                                                                           >  2:4923848/1‑150 (MQ=255)
                              gatTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGt                                                                                                                         <  1:538777/150‑1 (MQ=255)
                              gatTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGt                                                                                                                         >  1:451750/1‑150 (MQ=255)
                                  tGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGACGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGc                                                                                                                     <  1:4678964/150‑1 (MQ=255)
                                         cGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCta                                                                                                              <  1:5040871/150‑1 (MQ=255)
                                         cGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCta                                                                                                              <  1:1795888/150‑1 (MQ=255)
                                         cGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCta                                                                                                              >  1:2993658/1‑150 (MQ=255)
                                           cTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCtatt                                                                                                            <  1:4742940/150‑1 (MQ=255)
                                            tGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCtatta                                                                                                           >  1:1234883/1‑150 (MQ=255)
                                             gTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTCGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCtattat                                                                                                          <  1:354359/150‑1 (MQ=255)
                                              tACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATg                                                                                                         <  2:5187804/150‑1 (MQ=255)
                                              tACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATg                                                                                                         <  1:1446427/150‑1 (MQ=255)
                                                cAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGtt                                                                                                       <  1:5221782/150‑1 (MQ=255)
                                                 aGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGttt                                                                                                      >  1:3745209/1‑150 (MQ=255)
                                                    aaaTTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTc                                                                                                   <  2:4446722/150‑1 (MQ=255)
                                                    aaaTTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTc                                                                                                   >  2:2216706/1‑150 (MQ=255)
                                                    aaaTTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTc                                                                                                   <  1:1145649/150‑1 (MQ=255)
                                                    aaaTTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTc                                                                                                   >  1:2508219/1‑150 (MQ=255)
                                                     aaTTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTcc                                                                                                  >  2:5540035/1‑150 (MQ=255)
                                                      aTTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCg                                                                                                 <  1:5305348/150‑1 (MQ=255)
                                                        tCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGAt                                                                                               <  1:2352031/150‑1 (MQ=255)
                                                              ggCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGagaag                                                                                         <  1:4553676/150‑1 (MQ=255)
                                                                      ttGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATc                                                                                 >  2:4079741/1‑150 (MQ=255)
                                                                         gTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGAt                                                                              >  1:495763/1‑150 (MQ=255)
                                                                         gTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGAt                                                                              >  1:1261716/1‑150 (MQ=255)
                                                                         gTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGAt                                                                              >  1:5538867/1‑150 (MQ=255)
                                                                         gTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAAGACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGAt                                                                              <  2:3255567/150‑1 (MQ=255)
                                                                         ctATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAAGACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGAt                                                                              <  2:5292160/149‑1 (MQ=255)
                                                                           aTTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGAtgt                                                                            <  2:1882900/150‑1 (MQ=255)
                                                                             tAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGa                                                                          <  2:4905142/150‑1 (MQ=255)
                                                                              aGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGAc                                                                         <  1:5353109/150‑1 (MQ=255)
                                                                                  aTAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAAGACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGGTGAGAAGATTACATCGATGTGACAGGc                                                                     >  2:731644/1‑150 (MQ=255)
                                                                                    agtagCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCa                                                                   >  1:60823/1‑150 (MQ=255)
                                                                                    agtagCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCa                                                                   >  2:2378665/1‑150 (MQ=255)
                                                                                         cAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAAt                                                              <  2:2262261/150‑1 (MQ=255)
                                                                                           aGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATaa                                                            <  1:2266036/150‑1 (MQ=255)
                                                                                           aGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAAGACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATaa                                                            <  1:5471285/150‑1 (MQ=255)
                                                                                            gAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATaaa                                                           <  1:4889711/150‑1 (MQ=255)
                                                                                            gAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATaaa                                                           >  2:2516362/1‑150 (MQ=255)
                                                                                              aTCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATaacaa                                                         >  1:4049905/1‑150 (MQ=255)
                                                                                                ccTATTCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAaga                                                       <  2:3947484/150‑1 (MQ=255)
                                                                                                 cTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAagag                                                      <  1:2659554/150‑1 (MQ=255)
                                                                                                 cTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAAGACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAagag                                                      >  2:2305574/1‑150 (MQ=255)
                                                                                                     gCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGAttt                                                  <  1:825342/150‑1 (MQ=255)
                                                                                                     gCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGAttt                                                  <  2:375456/150‑1 (MQ=255)
                                                                                                        aataaATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGt                                               >  2:1719392/1‑150 (MQ=255)
                                                                                                            aaTCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGAt                                           <  1:2889893/150‑1 (MQ=255)
                                                                                                            aaTCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGAt                                           <  1:4618514/150‑1 (MQ=255)
                                                                                                            aaTCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGAt                                           <  2:603608/150‑1 (MQ=255)
                                                                                                            aaTCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGAt                                           <  2:4217869/150‑1 (MQ=255)
                                                                                                             aTCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGAtt                                          >  2:1963577/1‑150 (MQ=255)
                                                                                                             aTCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGAtt                                          <  1:1535/150‑1 (MQ=255)
                                                                                                             aTCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGAtt                                          <  1:1594778/150‑1 (MQ=255)
                                                                                                                   ttCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTataata                                    >  2:2626587/1‑150 (MQ=255)
                                                                                                                    tCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAatat                                   <  1:5190228/150‑1 (MQ=255)
                                                                                                                     cGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATc                                  >  2:579494/1‑150 (MQ=255)
                                                                                                                       aaCAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCAt                                >  2:1789815/1‑150 (MQ=255)
                                                                                                                       aaCAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCAt                                >  2:1125621/1‑150 (MQ=255)
                                                                                                                       aaCAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCAt                                >  2:2864664/1‑150 (MQ=255)
                                                                                                                          aGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTaa                             <  2:1007214/150‑1 (MQ=255)
                                                                                                                           gTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAAt                            >  2:2663052/1‑150 (MQ=255)
                                                                                                                            ttATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATg                           <  1:216428/150‑1 (MQ=255)
                                                                                                                             tATGAAGTGTGGGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGa                          >  2:1668949/1‑150 (MQ=255)
                                                                                                                              aTGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGaa                         >  1:20872/1‑150 (MQ=255)
                                                                                                                              aTGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGaa                         >  1:2001290/1‑150 (MQ=255)
                                                                                                                                  aGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGc                     >  2:29789/1‑150 (MQ=255)
                                                                                                                                    tgtgAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCtt                   <  1:3952469/150‑1 (MQ=255)
                                                                                                                                     gtgAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCttt                  <  2:3343438/150‑1 (MQ=255)
                                                                                                                                       gAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGTAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCttttt                <  2:4682988/150‑1 (MQ=255)
                                                                                                                                       gAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCttttt                <  2:1268518/150‑1 (MQ=255)
                                                                                                                                       gAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCttttt                <  2:4602854/150‑1 (MQ=255)
                                                                                                                                        aGGCTTGTGGTAAGACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTa               >  1:952011/1‑150 (MQ=255)
                                                                                                                                          gCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAAt             <  2:5061410/150‑1 (MQ=255)
                                                                                                                                             tgtgGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGt          >  1:1907937/1‑150 (MQ=255)
                                                                                                                                             tgtgGTAAGACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGt          <  2:2766070/150‑1 (MQ=255)
                                                                                                                                               tgGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAg        >  2:3084743/1‑150 (MQ=255)
                                                                                                                                                  tAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCaa     <  1:3927758/150‑1 (MQ=255)
                                                                                                                                                    aTACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCaaaa   >  1:3327227/1‑150 (MQ=255)
                                                                                                                                                     tACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAAt  >  2:4827068/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTGGCTATTATCCCAGAGTCATATCAGATGATTTGTTTTACGCTGTACAGGAAATTCGGTTGGCACCTTTTGGTATTAGCAATAGTAGCAAGAATCCTATGCTAATAAATCTACTTCGAACAGTTATGAAGTGTGAGGCTTGTGGTAATACCATGATTGTTCATGCGGTATCTGGAAGTTTGCATGGCTATTATGTTTGTCCGATGAGAAGATTACATCGATGTGACAGGCCATCAATAAAAAGAGATTTGGTTGATTATAATATCATTAATGAATTGCTTTTTAATTGTAGCAAAAT  >  CP009273/565150‑565448

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