Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,464,043 A→T 5.9% intergenic (+102/‑251) gtrS → / → tfaS serotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage/pseudogene, CPS‑53 (KpLE1) prophage; tail fiber assembly protein fragment;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,464,0430AT 5.9% 983.0 / 37.3 287intergenic (+102/‑251)gtrS/tfaSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage/pseudogene, CPS‑53 (KpLE1) prophage; tail fiber assembly protein fragment;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (152/118);  new base T (9/8);  total (161/126)
Fisher's exact test for biased strand distribution p-value = 8.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 290 total aligned reads displayed.

TATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACATATACTAACTGATT  >  CP009273/2463894‑2464191
                                                                                                                                                     |                                                                                                                                                    
    aaaTAATTCAAAGAAATCAATACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt                                                                                                                                                  >  2:3169094/1‑150 (MQ=255)
    aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt                                                                                                                                                  >  1:29231/1‑150 (MQ=255)
    aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt                                                                                                                                                  >  1:2529721/1‑150 (MQ=255)
    aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt                                                                                                                                                  >  1:180713/1‑150 (MQ=255)
    aaaTAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTAt                                                                                                                                                  >  1:525850/1‑150 (MQ=255)
     aataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGAGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATc                                                                                                                                                 >  1:2626290/1‑150 (MQ=255)
      ataatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCa                                                                                                                                                <  1:2745634/150‑1 (MQ=255)
       taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaa                                                                                                                                               >  1:2982329/1‑150 (MQ=255)
       taatTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCaa                                                                                                                                               >  1:4262318/1‑150 (MQ=255)
         atTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAAt                                                                                                                                             >  1:1645071/1‑150 (MQ=255)
          ttCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATg                                                                                                                                            <  1:4784225/150‑1 (MQ=255)
               aGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTg                                                                                                                                       <  2:2283241/150‑1 (MQ=255)
                gAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGa                                                                                                                                      >  2:642566/1‑150 (MQ=255)
                 aaaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAt                                                                                                                                     >  1:4455777/1‑150 (MQ=255)
                 aaaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAt                                                                                                                                     <  2:199503/150‑1 (MQ=255)
                 aaaTCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGAt                                                                                                                                     >  1:3591032/1‑150 (MQ=255)
                     cAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGTCTATCAAATGCCCTGATTCgt                                                                                                                                 >  2:4004442/1‑150 (MQ=255)
                     cAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgt                                                                                                                                 >  2:1376085/1‑150 (MQ=255)
                     cAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgt                                                                                                                                 >  2:457104/1‑150 (MQ=255)
                       aGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCgttg                                                                                                                               >  1:846424/1‑150 (MQ=255)
                        gACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGTCTATCAAATGCCCTGATTCgttgt                                                                                                                              >  2:3094309/1‑150 (MQ=255)
                          cAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTAGTGACTATCAAATGCCCTGATGCGTGGttt                                                                                                                            >  2:32810/1‑150 (MQ=255)
                             tAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCACTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGTCTATCAAATGCCCTGATTCGTTGTTTAtt                                                                                                                         >  1:2828667/1‑150 (MQ=255)
                                cGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTa                                                                                                                      >  1:2158869/1‑150 (MQ=255)
                                  aTGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAat                                                                                                                    <  2:2900368/150‑1 (MQ=255)
                                  aTGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAat                                                                                                                    <  2:1565487/150‑1 (MQ=255)
                                   tGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAata                                                                                                                   <  1:41715/150‑1 (MQ=255)
                                         gCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATctcct                                                                                                             >  1:1716678/1‑150 (MQ=255)
                                          cACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATctcctc                                                                                                            <  2:3409984/150‑1 (MQ=255)
                                          cACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATctcctc                                                                                                            <  2:4985039/150‑1 (MQ=255)
                                           aCTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCctct                                                                                                           <  1:4604377/150‑1 (MQ=255)
                                            cTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTa                                                                                                          >  1:4832667/1‑150 (MQ=255)
                                             tAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTAt                                                                                                         >  2:231981/1‑150 (MQ=255)
                                             tAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTAt                                                                                                         >  2:5134110/1‑150 (MQ=255)
                                              aaGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATc                                                                                                        <  2:509636/150‑1 (MQ=255)
                                               aGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCt                                                                                                       <  1:1543301/150‑1 (MQ=255)
                                                gTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTg                                                                                                      <  1:5014804/150‑1 (MQ=255)
                                                   tAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAg                                                                                                   >  1:2689092/1‑150 (MQ=255)
                                                     aaTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGTCTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGAc                                                                                                 <  1:2225112/150‑1 (MQ=255)
                                                     aaTTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGAc                                                                                                 <  2:4884858/150‑1 (MQ=255)
                                                       ttGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGa                                                                                               >  1:4723420/1‑150 (MQ=255)
                                                        tGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGAt                                                                                              >  2:3684165/1‑150 (MQ=255)
                                                         gcgcgcCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGTCTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATa                                                                                             <  1:868097/150‑1 (MQ=255)
                                                           gcgcCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAAc                                                                                           <  1:1263679/150‑1 (MQ=255)
                                                              ccAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTa                                                                                        >  2:4306004/1‑150 (MQ=255)
                                                               cAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTaa                                                                                       <  2:4292502/150‑1 (MQ=255)
                                                                aaTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTaaa                                                                                      >  1:3950160/1‑150 (MQ=255)
                                                                 aTCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAAt                                                                                     <  2:2697232/150‑1 (MQ=255)
                                                                   cgtGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGc                                                                                   >  2:3571195/3‑150 (MQ=255)
                                                                     tGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGTCTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGcat                                                                                 >  2:2623374/1‑150 (MQ=255)
                                                                     tGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGcat                                                                                 >  2:2476689/1‑150 (MQ=255)
                                                                        cgcgcACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGcatcat                                                                              >  2:2041768/1‑150 (MQ=255)
                                                                            cacaAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGcc                                                                          <  1:2753755/150‑1 (MQ=255)
                                                                              caAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCca                                                                        <  1:1051449/150‑1 (MQ=255)
                                                                                aGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCcaca                                                                      >  2:1860282/1‑150 (MQ=255)
                                                                                      aataCCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAATTGCATCATTGCCCACAACAAAc                                                                >  2:2971453/1‑150 (MQ=255)
                                                                                        taCCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAAccc                                                              >  1:2321081/1‑150 (MQ=255)
                                                                                          ccAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAAccccc                                                            >  1:362919/1‑150 (MQ=255)
                                                                                           cAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCt                                                           >  1:1361479/1‑150 (MQ=255)
                                                                                            aaCCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTc                                                          >  2:1332954/1‑150 (MQ=255)
                                                                                              ccTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCa                                                        <  1:698711/150‑1 (MQ=255)
                                                                                               cTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGTCTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAg                                                       <  2:1743565/150‑1 (MQ=255)
                                                                                                tAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGTCTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGa                                                      >  2:1674750/1‑150 (MQ=255)
                                                                                                tAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGa                                                      >  1:606679/1‑150 (MQ=255)
                                                                                                    ttCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGAGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCa                                                  >  1:4918017/1‑150 (MQ=255)
                                                                                                        cctcctcTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTg                                              >  2:4400215/1‑150 (MQ=255)
                                                                                                          tcctcTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCt                                            >  2:3887864/1‑150 (MQ=255)
                                                                                                           cctcTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTg                                           >  2:145117/1‑150 (MQ=255)
                                                                                                             tctTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGat                                         >  2:4148727/1‑150 (MQ=255)
                                                                                                               ttAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGatat                                       <  2:1146103/150‑1 (MQ=255)
                                                                                                                tAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGAtatt                                      <  2:2778538/150‑1 (MQ=255)
                                                                                                                tAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGAtatt                                      <  1:1132498/150‑1 (MQ=255)
                                                                                                                 agagTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGAtatta                                     <  2:763686/150‑1 (MQ=255)
                                                                                                                  gagTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGAtattat                                    <  2:3208292/150‑1 (MQ=255)
                                                                                                                    gTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTatca                                  >  1:4650042/1‑150 (MQ=255)
                                                                                                                     tGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTatcat                                 <  1:920558/150‑1 (MQ=255)
                                                                                                                         tataTCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAAc                             >  1:3429133/1‑150 (MQ=255)
                                                                                                                         tataTCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAAc                             <  1:3908319/150‑1 (MQ=255)
                                                                                                                          atatCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACa                            <  1:3399265/150‑1 (MQ=255)
                                                                                                                          atatCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACa                            >  1:2106935/1‑150 (MQ=255)
                                                                                                                          atatCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACa                            <  1:4101207/150‑1 (MQ=255)
                                                                                                                           tatCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAg                           <  2:3472946/150‑1 (MQ=255)
                                                                                                                            atCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGt                          <  2:1825394/150‑1 (MQ=255)
                                                                                                                               tcCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTgat                       <  1:3103771/150‑1 (MQ=255)
                                                                                                                               tcCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTgat                       <  1:4659290/150‑1 (MQ=255)
                                                                                                                                ccTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTgata                      >  1:3462247/1‑150 (MQ=255)
                                                                                                                                 cTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTgatag                     <  2:2594097/150‑1 (MQ=255)
                                                                                                                                      aGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACa                >  1:2474709/1‑150 (MQ=255)
                                                                                                                                       gAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACat               >  1:3030220/1‑150 (MQ=255)
                                                                                                                                         agtGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACatat             <  2:4866989/148‑1 (MQ=255)
                                                                                                                                         aTTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACatat             <  1:749698/150‑1 (MQ=255)
                                                                                                                                           tGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACATATAc           >  2:1902285/1‑150 (MQ=255)
                                                                                                                                           tGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACATATAc           <  2:1219514/150‑1 (MQ=255)
                                                                                                                                             cGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACATATACTa         >  1:1483818/1‑150 (MQ=255)
                                                                                                                                              ggTATTGACTATCAAATGCCCTGATTAGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACATATACTaa        <  1:3169094/150‑1 (MQ=255)
                                                                                                                                                tATTGTCTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACATATACTAACt      <  2:3746962/150‑1 (MQ=255)
                                                                                                                                                 aTTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACATATACTAACTg     <  1:1951570/150‑1 (MQ=255)
                                                                                                                                                  ttGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACATATACTAACTGa    >  2:462604/1‑150 (MQ=255)
                                                                                                                                                  ttGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACATATACTAACTGa    <  1:951234/150‑1 (MQ=255)
                                                                                                                                                    gACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACATATACTAACTGAtt  >  1:3507305/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATAATACCAACCTAATTTCTCCTCCTCTTAGAGTGACTATATCTCCTGATAGAATTGCGGTATTGACTATCAAATGCCCTGATTCGTTGTTTATTGTAATATCTCCTCTATCTGCAGACGATAACTTAAATGCATCATTGCCCACAACAAACCCCCTCCAGAACCAAGTGCTGATATTATCATCAACAGTGATAGATACATATACTAACTGATT  >  CP009273/2463894‑2464191

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: