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,1430AT98.2% 971.5 / 4.0 282I124F (ATT→TTT) fimBtyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base T (137/139);  total (141/141)
Fisher's exact test for biased strand distribution p-value = 6.84e-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 290 total aligned reads displayed.

TTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATAC  >  CP009273/4530996‑4531291
                                                                                                                                                   |                                                                                                                                                    
ttttCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCAtttt                                                                                                                                                    <  2:2903761/150‑1 (MQ=255)
 tttCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCAtttta                                                                                                                                                   <  2:656629/150‑1 (MQ=255)
     aacaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtctc                                                                                                                                               >  2:916972/1‑150 (MQ=255)
     aacaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtctc                                                                                                                                               <  1:3667194/150‑1 (MQ=255)
     aacaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtctc                                                                                                                                               <  1:965351/150‑1 (MQ=255)
     aacaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtctc                                                                                                                                               >  2:393492/1‑150 (MQ=255)
      acaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCg                                                                                                                                              <  2:1840175/150‑1 (MQ=255)
       caacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGa                                                                                                                                             >  1:1647349/1‑150 (MQ=255)
        aacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGAc                                                                                                                                            >  1:1360029/1‑150 (MQ=255)
         acGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACt                                                                                                                                           >  1:754829/1‑150 (MQ=255)
           gCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTc                                                                                                                                         >  2:313662/1‑150 (MQ=255)
             aCCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCg                                                                                                                                       <  1:3796223/150‑1 (MQ=255)
               ccGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCggt                                                                                                                                     >  2:780184/1‑150 (MQ=255)
                cGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCTTACTTCGTACCCGCAAGCAGAGAGCGAGAGGGTATTTTTATCACCTAAGGGGAATACGCTTACACGGCAACAGTTATACCAAATTATCTCTACTTCAggtg                                                                                                                                    >  2:2659539/1‑150 (MQ=255)
                cGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCggtg                                                                                                                                    >  2:1606367/1‑150 (MQ=255)
                   tATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTa                                                                                                                                 <  1:1972538/150‑1 (MQ=255)
                    aTTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTaa                                                                                                                                <  1:4563047/150‑1 (MQ=255)
                    aTTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTaa                                                                                                                                <  2:1491263/150‑1 (MQ=255)
                     ttGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAAt                                                                                                                               >  1:2317141/1‑150 (MQ=255)
                       gAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGc                                                                                                                             <  1:2547172/150‑1 (MQ=255)
                        aaTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGcc                                                                                                                            <  2:2426402/150‑1 (MQ=255)
                             agaagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGtt                                                                                                                       <  1:3014170/150‑1 (MQ=255)
                              gaagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTg                                                                                                                      >  2:4259553/1‑150 (MQ=255)
                               aagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGt                                                                                                                     >  1:3711825/1‑150 (MQ=255)
                                 gTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCa                                                                                                                   >  1:72063/1‑150 (MQ=255)
                                   tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACt                                                                                                                 >  2:1287275/1‑150 (MQ=255)
                                   tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACt                                                                                                                 >  2:1910835/1‑150 (MQ=255)
                                   tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACt                                                                                                                 <  1:2707736/150‑1 (MQ=255)
                                   tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACt                                                                                                                 <  1:4033421/150‑1 (MQ=255)
                                    cAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTg                                                                                                                <  2:1282612/150‑1 (MQ=255)
                                     aGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTgg                                                                                                               >  2:4538357/1‑150 (MQ=255)
                                         tttAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGAt                                                                                                           >  1:3543068/1‑150 (MQ=255)
                                         tttAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGAt                                                                                                           >  2:3550695/1‑150 (MQ=255)
                                          ttAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTCTCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGAtt                                                                                                          <  2:1229090/150‑1 (MQ=255)
                                               aaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATcc                                                                                                     >  1:4119933/1‑150 (MQ=255)
                                               aaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATcc                                                                                                     <  1:2767938/150‑1 (MQ=255)
                                                aaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCg                                                                                                    >  2:3609567/1‑150 (MQ=255)
                                                   tGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGcac                                                                                                 <  1:3981675/150‑1 (MQ=255)
                                                       tGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGt                                                                                             <  1:2735580/150‑1 (MQ=255)
                                                         aGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCTCTGGAGATTCATCCGCACATGTTa                                                                                           <  2:3818762/150‑1 (MQ=255)
                                                         aGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTa                                                                                           >  1:2596134/1‑150 (MQ=255)
                                                         aGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTa                                                                                           >  1:2596354/1‑150 (MQ=255)
                                                         aGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTa                                                                                           >  1:1928959/1‑150 (MQ=255)
                                                            aTCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGc                                                                                        >  2:4404819/1‑150 (MQ=255)
                                                              ccGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCa                                                                                      <  2:3161153/150‑1 (MQ=255)
                                                                   cTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgt                                                                                 >  1:2296663/1‑150 (MQ=255)
                                                                    ttCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgtg                                                                                <  1:1846525/150‑1 (MQ=255)
                                                                       gTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTgg                                                                             >  2:3260245/1‑150 (MQ=255)
                                                                        tACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGt                                                                            <  1:4488135/150‑1 (MQ=255)
                                                                         aCCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGtt                                                                           >  1:2416380/1‑150 (MQ=255)
                                                                           ccGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTTGAGATTCATCCGCACATGTTACGCCATTCGTGTGGtttt                                                                         <  2:1377570/150‑1 (MQ=255)
                                                                           ccGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGtttt                                                                         >  2:1097272/1‑150 (MQ=255)
                                                                           ccGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGtttt                                                                         >  2:3428060/1‑150 (MQ=255)
                                                                              cATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCt                                                                      <  2:1201040/150‑1 (MQ=255)
                                                                              cATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCt                                                                      <  1:1306381/150‑1 (MQ=255)
                                                                               aTGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCtt                                                                     <  1:483404/150‑1 (MQ=255)
                                                                                tGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCttt                                                                    >  2:4586817/1‑150 (MQ=255)
                                                                                tGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCttt                                                                    <  1:3143162/150‑1 (MQ=255)
                                                                                 gCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTg                                                                   <  2:3353777/150‑1 (MQ=255)
                                                                                   tGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGc                                                                 >  2:2550557/1‑150 (MQ=255)
                                                                                    gagaGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCTTGTGGTTTTGCTTTGGCg                                                                >  2:306262/1‑150 (MQ=255)
                                                                                     agagCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGa                                                               <  1:2864881/150‑1 (MQ=255)
                                                                                       agCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAat                                                             <  2:2736314/150‑1 (MQ=255)
                                                                                        gCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAata                                                            <  2:1339743/150‑1 (MQ=255)
                                                                                        gCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAata                                                            <  2:407793/150‑1 (MQ=255)
                                                                                          gAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATg                                                          >  2:2193836/1‑150 (MQ=255)
                                                                                           aGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATgg                                                         <  1:1030376/150‑1 (MQ=255)
                                                                                            gTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATATCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATggg                                                        <  1:2154571/150‑1 (MQ=255)
                                                                                             tGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGa                                                       <  2:2001256/150‑1 (MQ=255)
                                                                                             tGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGa                                                       <  1:1108254/150‑1 (MQ=255)
                                                                                             tGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGa                                                       <  1:1159548/150‑1 (MQ=255)
                                                                                                    ttttATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATa                                                <  2:2959489/150‑1 (MQ=255)
                                                                                                      ttATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATAcg                                              <  2:3483082/150‑1 (MQ=255)
                                                                                                          cACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGAc                                          >  1:640485/1‑150 (MQ=255)
                                                                                                           aCGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACt                                         <  2:2287092/150‑1 (MQ=255)
                                                                                                             gTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTa                                       <  2:1156886/150‑1 (MQ=255)
                                                                                                             gTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTa                                       <  1:2926027/150‑1 (MQ=255)
                                                                                                                aGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATcc                                    >  1:100330/1‑150 (MQ=255)
                                                                                                                  gggAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAg                                  <  2:3316463/150‑1 (MQ=255)
                                                                                                                      aTCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGAtt                              <  2:1676444/150‑1 (MQ=255)
                                                                                                                      aTCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGAtt                              <  1:2085419/150‑1 (MQ=255)
                                                                                                                        ccGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTAt                            >  1:1546039/1‑150 (MQ=255)
                                                                                                                        ccGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTAt                            >  2:1561572/1‑150 (MQ=255)
                                                                                                                         cGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATc                           >  2:3379787/1‑150 (MQ=255)
                                                                                                                          gCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCt                          >  2:2052105/1‑150 (MQ=255)
                                                                                                                             ttCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTgg                       >  2:4182511/1‑150 (MQ=255)
                                                                                                                             ttCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTgg                       >  1:1196972/1‑150 (MQ=255)
                                                                                                                             ttCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTgg                       <  1:4264097/150‑1 (MQ=255)
                                                                                                                               ctcGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGc                     >  2:3981741/1‑150 (MQ=255)
                                                                                                                               ctcGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGc                     >  2:279634/1‑150 (MQ=255)
                                                                                                                                tcGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCa                    >  1:787122/1‑150 (MQ=255)
                                                                                                                                tcGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCa                    >  1:2050444/1‑150 (MQ=255)
                                                                                                                                 cGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCAt                   >  2:1248691/1‑150 (MQ=255)
                                                                                                                                  ggCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATc                  >  2:3768305/1‑150 (MQ=255)
                                                                                                                                     aaCAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCa               <  2:470537/150‑1 (MQ=255)
                                                                                                                                        aGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAata            >  1:3373729/1‑150 (MQ=255)
                                                                                                                                        aGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAata            >  1:3804787/1‑150 (MQ=255)
                                                                                                                                            ttACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCg        <  2:1936280/150‑1 (MQ=255)
                                                                                                                                              aCCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTc      >  1:1470194/1‑150 (MQ=255)
                                                                                                                                               ccATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCa     >  1:3821411/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATAC  >  CP009273/4530996‑4531291

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: