Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,531,143 A→T 92.8% I124F (ATT→TTT)  fimB → tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,531,1430AT92.8% 967.5 / 46.2 286I124F (ATT→TTT) fimBtyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base A (10/11);  new base T (124/141);  total (134/152)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 300 total aligned reads displayed.

TTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACT  >  CP009273/4530996‑4531292
                                                                                                                                                   |                                                                                                                                                     
ttatCAACAACGCACCCGCTATTTAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCAtttt                                                                                                                                                     <  1:3280494/147‑1 (MQ=255)
  ttCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCAttttat                                                                                                                                                   <  2:1897682/150‑1 (MQ=255)
    caacaaCGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAATGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtct                                                                                                                                                 >  1:418114/1‑150 (MQ=255)
    caacaaCGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtct                                                                                                                                                 <  2:2823917/150‑1 (MQ=255)
    caacaaCGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtct                                                                                                                                                 <  1:3341952/150‑1 (MQ=255)
    caacaaCGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCAGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtct                                                                                                                                                 >  2:2246076/1‑150 (MQ=255)
     aacaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtctc                                                                                                                                                <  2:3562566/150‑1 (MQ=255)
      acaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCg                                                                                                                                               <  1:1022853/150‑1 (MQ=255)
           gCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTc                                                                                                                                          >  1:2282853/1‑150 (MQ=255)
            cACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTcc                                                                                                                                         <  1:352631/150‑1 (MQ=255)
            cACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTcc                                                                                                                                         <  1:2822909/150‑1 (MQ=255)
              cccGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCgg                                                                                                                                       >  1:2128661/1‑150 (MQ=255)
               ccGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCggt                                                                                                                                      <  1:1570188/150‑1 (MQ=255)
                cGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCggtg                                                                                                                                     >  1:3909281/1‑150 (MQ=255)
                 gCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCggtgg                                                                                                                                    <  2:3686283/150‑1 (MQ=255)
                 gCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCggtgg                                                                                                                                    <  2:857476/150‑1 (MQ=255)
                  cTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCggtggt                                                                                                                                   >  2:2740522/1‑150 (MQ=255)
                    aTTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCAttttat                                                                                                                                                   <  2:2323603/132‑1 (MQ=255)
                             agaagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGACTCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGtt                                                                                                                        >  2:2447827/1‑150 (MQ=255)
                               aagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGt                                                                                                                      >  2:2179168/1‑150 (MQ=255)
                                agTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGTGTTGTc                                                                                                                     <  1:185382/150‑1 (MQ=255)
                                agTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTc                                                                                                                     <  1:4018654/150‑1 (MQ=255)
                                 gTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCa                                                                                                                    >  1:4223701/1‑150 (MQ=255)
                                  ttCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCAc                                                                                                                   >  1:3842260/1‑150 (MQ=255)
                                   tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACt                                                                                                                  <  2:2335544/150‑1 (MQ=255)
                                   tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACt                                                                                                                  >  1:2995090/1‑150 (MQ=255)
                                   tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACt                                                                                                                  <  1:2347724/150‑1 (MQ=255)
                                       gCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGgag                                                                                                              >  2:2085199/1‑150 (MQ=255)
                                        cTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGgaga                                                                                                             >  1:3517620/1‑150 (MQ=255)
                                         tttAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGAt                                                                                                            >  1:2039071/1‑150 (MQ=255)
                                             aaaaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCAt                                                                                                        >  1:3407705/1‑150 (MQ=255)
                                              aaaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATc                                                                                                       <  2:4289688/150‑1 (MQ=255)
                                              aaaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATc                                                                                                       >  2:4076291/1‑150 (MQ=255)
                                               aaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATcc                                                                                                      >  2:1559051/1‑150 (MQ=255)
                                               aaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATcc                                                                                                      <  2:2362087/150‑1 (MQ=255)
                                                aaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCg                                                                                                     >  1:2567497/1‑150 (MQ=255)
                                                      ttGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATg                                                                                               <  2:1108262/150‑1 (MQ=255)
                                                      ttGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATg                                                                                               >  1:749718/1‑150 (MQ=255)
                                                       tGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATACGCACATGt                                                                                              >  1:8434/1‑150 (MQ=255)
                                                        gAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGtt                                                                                             >  1:17800/1‑150 (MQ=255)
                                                        gAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGtt                                                                                             >  1:16401/1‑150 (MQ=255)
                                                         aGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTa                                                                                            <  2:583013/150‑1 (MQ=255)
                                                            aTCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGc                                                                                         >  2:541095/1‑150 (MQ=255)
                                                               cGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCAt                                                                                      <  1:2507614/150‑1 (MQ=255)
                                                               cGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCAt                                                                                      >  2:634018/1‑150 (MQ=255)
                                                                gTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCAtt                                                                                     >  2:3253246/1‑150 (MQ=255)
                                                                gTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCAtt                                                                                     >  2:569441/1‑150 (MQ=255)
                                                                gTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCAtt                                                                                     <  2:485429/150‑1 (MQ=255)
                                                                 tACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTc                                                                                    >  2:1598791/1‑150 (MQ=255)
                                                                      cGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGtgtg                                                                               <  2:1850684/150‑1 (MQ=255)
                                                                       gTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTgg                                                                              <  2:1631919/150‑1 (MQ=255)
                                                                          cccGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGttt                                                                           <  2:2431942/150‑1 (MQ=255)
                                                                          cccGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGttt                                                                           <  1:1200322/150‑1 (MQ=255)
                                                                             gCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGc                                                                        >  2:1528594/1‑150 (MQ=255)
                                                                              cATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCt                                                                       <  2:293099/150‑1 (MQ=255)
                                                                                  cTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTgg                                                                   <  1:2322867/150‑1 (MQ=255)
                                                                                   tGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGc                                                                  <  2:2048458/150‑1 (MQ=255)
                                                                                    gagaGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCg                                                                 <  1:1987868/150‑1 (MQ=255)
                                                                                     agagCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGa                                                                >  1:1898874/1‑150 (MQ=255)
                                                                                     agagCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGATTTGGCGa                                                                >  2:3041943/1‑150 (MQ=255)
                                                                                        gCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAata                                                             >  2:1921450/1‑150 (MQ=255)
                                                                                             tGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGa                                                        >  2:3374282/1‑150 (MQ=255)
                                                                                              gggTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGaa                                                       <  2:3990498/150‑1 (MQ=255)
                                                                                              gggTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGaa                                                       >  1:2265671/1‑150 (MQ=255)
                                                                                              gggTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGaa                                                       >  1:4275600/1‑150 (MQ=255)
                                                                                              gggTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGaa                                                       <  2:1053645/150‑1 (MQ=255)
                                                                                               ggTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAAt                                                      >  2:2681823/1‑150 (MQ=255)
                                                                                                gTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATa                                                     <  1:444157/150‑1 (MQ=255)
                                                                                                  aTTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGa                                                   >  2:3503431/1‑150 (MQ=255)
                                                                                                  aTTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGa                                                   >  1:3606275/1‑150 (MQ=255)
                                                                                                    ttttATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATa                                                 <  1:3567923/150‑1 (MQ=255)
                                                                                                    ttttATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATa                                                 <  2:442758/150‑1 (MQ=255)
                                                                                                      ttATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATAcg                                               >  2:552981/1‑150 (MQ=255)
                                                                                                      ttATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATAcg                                               <  2:755317/150‑1 (MQ=255)
                                                                                                        aTCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATAcgcg                                             <  1:3311142/150‑1 (MQ=255)
                                                                                                             gTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTa                                        >  2:2440608/1‑150 (MQ=255)
                                                                                                             ctAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAACTCGACTTCCGGTGGTAATGCCGTGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCTTGTGGTTTTGCTTTGGCGAATATGGGAATATATACGCGACTTa                                        >  2:1367180/2‑150 (MQ=255)
                                                                                                              tAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTAt                                       <  2:3444508/150‑1 (MQ=255)
                                                                                                               aaGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATc                                      >  2:838700/1‑150 (MQ=255)
                                                                                                                aGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATcc                                     <  1:1190576/150‑1 (MQ=255)
                                                                                                                   ggAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAgg                                  >  2:2935422/1‑150 (MQ=255)
                                                                                                                     aaTCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGAt                                <  2:4223701/150‑1 (MQ=255)
                                                                                                                     aaTCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGAt                                <  2:1721273/150‑1 (MQ=255)
                                                                                                                     aaTCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCAGTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGAt                                >  1:1079745/1‑150 (MQ=255)
                                                                                                                        ccGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTAt                             >  2:3331009/1‑150 (MQ=255)
                                                                                                                              tctcGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTggg                       <  2:914371/150‑1 (MQ=255)
                                                                                                                               ctcGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGc                      <  2:2689271/150‑1 (MQ=255)
                                                                                                                                  ggCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATc                   <  2:2543940/150‑1 (MQ=255)
                                                                                                                                  ggCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATc                   >  1:497659/1‑150 (MQ=255)
                                                                                                                                        aGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAata             <  1:2399388/150‑1 (MQ=255)
                                                                                                                                        aGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAata             >  1:2205059/1‑150 (MQ=255)
                                                                                                                                          ttttACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATAtt           >  2:2298410/1‑150 (MQ=255)
                                                                                                                                          ttttACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATAtt           <  2:2131248/150‑1 (MQ=255)
                                                                                                                                           tttACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTc          <  1:1712557/150‑1 (MQ=255)
                                                                                                                                           tttACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTc          <  1:1297317/150‑1 (MQ=255)
                                                                                                                                            ttACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCg         <  1:1478768/150‑1 (MQ=255)
                                                                                                                                             tACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGt        >  1:3187718/1‑150 (MQ=255)
                                                                                                                                               ccATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCa      <  1:4251889/150‑1 (MQ=255)
                                                                                                                                               ccATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCa      <  1:2447827/150‑1 (MQ=255)
                                                                                                                                                   attatCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACt  <  1:541095/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACT  >  CP009273/4530996‑4531292

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: