Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,531,143 A→T 100% 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,1430AT96.0% 732.0 / 12.3 210I124F (ATT→TTT) fimBtyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base A (5/4);  new base T (95/106);  total (100/110)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 213 total aligned reads displayed.

TTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCAT  >  CP009273/4530998‑4531289
                                                                                                                                                 |                                                                                                                                                  
ttCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCAttttat                                                                                                                                                <  1:4867087/150‑1 (MQ=255)
 tCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtc                                                                                                                                               <  1:4669853/150‑1 (MQ=255)
 tCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtc                                                                                                                                               >  1:733726/1‑150 (MQ=255)
   aacaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtctc                                                                                                                                             >  1:3870576/1‑150 (MQ=255)
   aacaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtctc                                                                                                                                             >  1:1369977/1‑150 (MQ=255)
    accacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTTCTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTAtctct                                                                                                                                            <  1:949799/147‑2 (MQ=255)
    acaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCg                                                                                                                                            <  1:1404397/150‑1 (MQ=255)
    acaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCg                                                                                                                                            <  2:4842103/150‑1 (MQ=255)
     caacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCCTTTTATATCGa                                                                                                                                           >  1:2216359/1‑150 (MQ=255)
      aacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGAc                                                                                                                                          >  2:318086/1‑150 (MQ=255)
       acGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACt                                                                                                                                         >  2:539204/1‑150 (MQ=255)
         gCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTc                                                                                                                                       >  1:2219046/1‑150 (MQ=255)
           aCCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCg                                                                                                                                     <  1:570771/150‑1 (MQ=255)
               gCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCggtgg                                                                                                                                 <  2:2236063/150‑1 (MQ=255)
               gCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCggtgg                                                                                                                                 >  2:266187/1‑150 (MQ=255)
                 tATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTa                                                                                                                               <  1:865658/150‑1 (MQ=255)
                   ttGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAAt                                                                                                                             <  2:1626062/150‑1 (MQ=255)
                     gAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGc                                                                                                                           >  1:3742867/1‑150 (MQ=255)
                      aaTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGcc                                                                                                                          >  2:3981773/1‑150 (MQ=255)
                       aTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCg                                                                                                                         <  1:4898424/150‑1 (MQ=255)
                       aTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCg                                                                                                                         <  2:5076575/150‑1 (MQ=255)
                          aagaagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGt                                                                                                                      <  2:586891/150‑1 (MQ=255)
                          aagaagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGt                                                                                                                      >  1:2066321/1‑150 (MQ=255)
                           agaagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGtt                                                                                                                     <  1:933590/150‑1 (MQ=255)
                           agaagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGtt                                                                                                                     <  2:3405575/150‑1 (MQ=255)
                            gaagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTg                                                                                                                    >  2:3726154/1‑150 (MQ=255)
                                 tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACt                                                                                                               <  1:1703069/150‑1 (MQ=255)
                                  cAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTg                                                                                                              <  2:4302253/150‑1 (MQ=255)
                                     gCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGgag                                                                                                           <  2:3567890/150‑1 (MQ=255)
                                     gCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGgag                                                                                                           <  2:2273526/150‑1 (MQ=255)
                                     gCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGgag                                                                                                           <  2:1925050/150‑1 (MQ=255)
                                           aaaaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCAt                                                                                                     <  1:2450417/150‑1 (MQ=255)
                                           aaaaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCAt                                                                                                     >  1:4733274/1‑150 (MQ=255)
                                                cTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGca                                                                                                >  1:5073720/1‑150 (MQ=255)
                                                cTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGca                                                                                                >  1:3233166/1‑150 (MQ=255)
                                                 tGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGcac                                                                                               <  2:1574229/150‑1 (MQ=255)
                                                       aGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTa                                                                                         <  1:291928/150‑1 (MQ=255)
                                                        gTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTAc                                                                                        <  1:1001117/150‑1 (MQ=255)
                                                           tCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGcc                                                                                     >  1:1545541/1‑150 (MQ=255)
                                                            ccGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCa                                                                                    >  1:4813623/1‑150 (MQ=255)
                                                             cGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGccct                                                                                   >  1:1828456/1‑148 (MQ=255)
                                                              gTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCAtt                                                                                  <  2:734354/150‑1 (MQ=255)
                                                              gTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCAtt                                                                                  <  2:734093/150‑1 (MQ=255)
                                                               tACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTc                                                                                 >  1:2267569/1‑150 (MQ=255)
                                                                aCTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCg                                                                                >  1:3999886/1‑150 (MQ=255)
                                                                 cTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgt                                                                               >  2:327442/1‑150 (MQ=255)
                                                                  ttCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgtg                                                                              <  1:422893/150‑1 (MQ=255)
                                                                    cGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGtgtg                                                                            <  1:1470477/150‑1 (MQ=255)
                                                                      tACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGt                                                                          >  1:4260001/1‑150 (MQ=255)
                                                                       aCCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGtt                                                                         <  1:300874/150‑1 (MQ=255)
                                                                             aTGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGcgt                                                                   >  2:4470640/1‑148 (MQ=255)
                                                                               gCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTg                                                                 >  2:2039961/1‑150 (MQ=255)
                                                                                 tGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGc                                                               <  2:1213619/150‑1 (MQ=255)
                                                                                 tGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCAGTCGTGTGGTTTTGCTTTGGc                                                               >  1:1197163/1‑150 (MQ=255)
                                                                                  gagaGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCg                                                              >  2:3108667/1‑150 (MQ=255)
                                                                                   agagCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGa                                                             <  1:2739176/150‑1 (MQ=255)
                                                                                    gagCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGaa                                                            <  2:112159/150‑1 (MQ=255)
                                                                                     agCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAat                                                           <  2:3213254/150‑1 (MQ=255)
                                                                                       cGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAatat                                                         <  1:1970963/150‑1 (MQ=255)
                                                                                        gAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATg                                                        <  2:3987279/150‑1 (MQ=255)
                                                                                         aGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATgg                                                       <  2:3011222/150‑1 (MQ=255)
                                                                                          gTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATggg                                                      >  1:174800/1‑150 (MQ=255)
                                                                                          gTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATggg                                                      >  1:4608319/1‑150 (MQ=255)
                                                                                            gggTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGaa                                                    <  2:4640147/150‑1 (MQ=255)
                                                                                              gTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATa                                                  <  1:3419822/150‑1 (MQ=255)
                                                                                              gTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATa                                                  <  1:1153232/150‑1 (MQ=255)
                                                                                              gTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATa                                                  >  2:1131082/1‑150 (MQ=255)
                                                                                               tATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAg                                                 >  2:175799/1‑150 (MQ=255)
                                                                                                aTTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGa                                                >  1:643883/1‑150 (MQ=255)
                                                                                                aTTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGa                                                >  2:1880098/1‑150 (MQ=255)
                                                                                                aTTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGa                                                >  1:71829/1‑150 (MQ=255)
                                                                                                aTTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGa                                                >  1:2649779/1‑150 (MQ=255)
                                                                                                aTTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGa                                                >  1:3867452/1‑150 (MQ=255)
                                                                                                aTTCTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGa                                                >  1:2649725/1‑150 (MQ=255)
                                                                                                 tttttATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGAt                                               <  2:2357130/150‑1 (MQ=255)
                                                                                                   tttATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATAc                                             <  2:3742867/150‑1 (MQ=255)
                                                                                                      aTCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATAcgcg                                          >  1:4132800/1‑150 (MQ=255)
                                                                                                        cACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGAc                                        <  1:3638506/150‑1 (MQ=255)
                                                                                                           gTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTa                                     <  2:1997840/150‑1 (MQ=255)
                                                                                                           gTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTa                                     >  2:638686/1‑150 (MQ=255)
                                                                                                              aGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATcc                                  <  1:249151/150‑1 (MQ=255)
                                                                                                                      ccGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTAt                          <  2:4036192/150‑1 (MQ=255)
                                                                                                                       cGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATc                         >  1:3502610/1‑150 (MQ=255)
                                                                                                                        gCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCt                        <  1:539301/150‑1 (MQ=255)
                                                                                                                          tttCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTg                      <  1:1157708/150‑1 (MQ=255)
                                                                                                                          tttCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTg                      <  2:467455/150‑1 (MQ=255)
                                                                                                                          tttCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTg                      <  2:467097/150‑1 (MQ=255)
                                                                                                                              tcGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCa                  <  2:2174834/150‑1 (MQ=255)
                                                                                                                                ggCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGACTAGATACGCGACTTATCCAGGATTATCTTGGGCATc                >  2:630156/1‑150 (MQ=255)
                                                                                                                                ggCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATc                <  1:1043441/150‑1 (MQ=255)
                                                                                                                                ggCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATc                <  1:5061613/150‑1 (MQ=255)
                                                                                                                                 gCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCg               <  2:1668454/150‑1 (MQ=255)
                                                                                                                                   aaCAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGATTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCa             >  2:4638071/1‑150 (MQ=255)
                                                                                                                                   aaCAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCa             >  1:3851912/1‑150 (MQ=255)
                                                                                                                                      aGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAata          >  1:5164715/1‑150 (MQ=255)
                                                                                                                                        ttttACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATAtt        >  2:4135752/1‑150 (MQ=255)
                                                                                                                                           tACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGt     >  1:1405398/1‑150 (MQ=255)
                                                                                                                                           tACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGt     >  2:3900746/1‑150 (MQ=255)
                                                                                                                                            aCCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTc    <  1:5012658/150‑1 (MQ=255)
                                                                                                                                              cATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGTCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCAt  <  1:1314817/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCAT  >  CP009273/4530998‑4531289

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: