Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,531,143 A→T 87.9% 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,1430AT87.9% 949.3 / 102.2 296I124F (ATT→TTT) fimBtyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base A (22/14);  new base T (142/118);  total (164/132)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 306 total aligned reads displayed.

TTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATA  >  CP009273/4530997‑4531290
                                                                                                                                                  |                                                                                                                                                   
tttCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCAtttt                                                                                                                                                   <  2:948042/149‑1 (MQ=255)
 ttCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCAtatta                                                                                                                                                  <  2:4633857/149‑1 (MQ=255)
  tCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCAttttat                                                                                                                                                 <  1:526014/149‑1 (MQ=255)
   caacaaCGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtc                                                                                                                                                >  1:1750122/1‑149 (MQ=255)
   caacaaCGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtc                                                                                                                                                >  2:4128930/1‑149 (MQ=255)
    aacaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTAtct                                                                                                                                               >  2:532769/1‑149 (MQ=255)
        acGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGAc                                                                                                                                           <  2:2560399/149‑1 (MQ=255)
           cACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTc                                                                                                                                        >  1:3819346/1‑149 (MQ=255)
            aCCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTcc                                                                                                                                       <  2:3596436/149‑1 (MQ=255)
                gCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCggtg                                                                                                                                   >  1:369310/1‑149 (MQ=255)
                gCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCggtg                                                                                                                                   >  1:1835249/1‑149 (MQ=255)
                 cTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCggtgg                                                                                                                                  <  2:1279131/149‑1 (MQ=255)
                    ttGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTaa                                                                                                                               <  2:3950010/149‑1 (MQ=255)
                      gAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATg                                                                                                                             >  1:3222849/1‑149 (MQ=255)
                       aaTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGc                                                                                                                            <  1:845930/149‑1 (MQ=255)
                       aaTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGc                                                                                                                            >  1:3189399/1‑149 (MQ=255)
                       aaTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGc                                                                                                                            >  1:946043/1‑149 (MQ=255)
                       aaTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGc                                                                                                                            >  2:1379051/1‑149 (MQ=255)
                        aTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGcc                                                                                                                           <  2:2357277/149‑1 (MQ=255)
                        aTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGcc                                                                                                                           <  1:1958474/149‑1 (MQ=255)
                          aaaGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCgg                                                                                                                         >  1:1230262/1‑149 (MQ=255)
                           aagaagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCggg                                                                                                                        >  1:979582/1‑149 (MQ=255)
                            agaagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGt                                                                                                                       >  2:4073230/1‑149 (MQ=255)
                             gaagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGtt                                                                                                                      >  2:747732/1‑149 (MQ=255)
                              aagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTg                                                                                                                     >  1:2791762/1‑149 (MQ=255)
                              aagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGGGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTg                                                                                                                     >  1:2330878/1‑149 (MQ=255)
                                 ttCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCa                                                                                                                  <  2:1502374/149‑1 (MQ=255)
                                  tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCAc                                                                                                                 >  2:4324352/1‑149 (MQ=255)
                                   cAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACt                                                                                                                >  2:2614850/1‑149 (MQ=255)
                                    aGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTg                                                                                                               >  2:511230/1‑149 (MQ=255)
                                     ggCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTgg                                                                                                              >  1:1180563/1‑149 (MQ=255)
                                     ggCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTgg                                                                                                              >  2:196740/1‑149 (MQ=255)
                                     ggCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTgg                                                                                                              >  1:1883261/1‑149 (MQ=255)
                                      gCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGga                                                                                                             <  1:257385/149‑1 (MQ=255)
                                      gCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGga                                                                                                             <  1:569599/149‑1 (MQ=255)
                                       cTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGgag                                                                                                            >  1:1510405/1‑149 (MQ=255)
                                         ttAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGAt                                                                                                          <  2:2304922/149‑1 (MQ=255)
                                          tAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGAtt                                                                                                         <  2:791743/149‑1 (MQ=255)
                                          tAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGAtt                                                                                                         >  1:4644903/1‑149 (MQ=255)
                                            aaaaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCa                                                                                                       >  1:4282213/1‑149 (MQ=255)
                                            aaaaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGATTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCa                                                                                                       >  2:1129765/1‑149 (MQ=255)
                                              aaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATc                                                                                                     <  2:2099592/149‑1 (MQ=255)
                                               aaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATcc                                                                                                    >  1:569631/1‑149 (MQ=255)
                                               aaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATcc                                                                                                    >  1:2096597/1‑149 (MQ=255)
                                               aaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATcc                                                                                                    >  1:1644277/1‑149 (MQ=255)
                                                aCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGTAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCg                                                                                                   >  1:4007424/1‑149 (MQ=255)
                                                    gTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGcaca                                                                                               >  2:2900761/1‑149 (MQ=255)
                                                        aGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGtt                                                                                           <  1:132726/149‑1 (MQ=255)
                                                         gTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTa                                                                                          <  1:2693223/149‑1 (MQ=255)
                                                           aTCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACg                                                                                        >  1:457571/1‑149 (MQ=255)
                                                              cGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTCAGGGGAATCCGCTGTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCa                                                                                     <  1:2009018/149‑1 (MQ=255)
                                                                tACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCAtt                                                                                   >  1:1081991/1‑149 (MQ=255)
                                                                 aCTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTc                                                                                  >  1:3421274/1‑149 (MQ=255)
                                                                    tCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgtg                                                                               >  2:3453042/1‑149 (MQ=255)
                                                                      gTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGtgtg                                                                             >  1:2270052/1‑149 (MQ=255)
                                                                       tACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTgg                                                                            <  2:3271320/149‑1 (MQ=255)
                                                                       tACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTgg                                                                            <  2:489882/149‑1 (MQ=255)
                                                                       tACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTgg                                                                            >  2:3889624/1‑149 (MQ=255)
                                                                         cccGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGtt                                                                          >  2:2066333/1‑149 (MQ=255)
                                                                             cATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGc                                                                      <  1:729631/149‑1 (MQ=255)
                                                                                    agagCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCg                                                               >  2:3088341/1‑149 (MQ=255)
                                                                                     gagCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGa                                                              >  1:4242422/1‑149 (MQ=255)
                                                                                      agCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGaa                                                             <  1:419662/149‑1 (MQ=255)
                                                                                          aGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATg                                                         >  1:4416571/1‑149 (MQ=255)
                                                                                               gTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAAt                                                    <  2:705945/149‑1 (MQ=255)
                                                                                                     ttATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATAc                                              <  1:307430/149‑1 (MQ=255)
                                                                                                        tCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATAcgcg                                           >  2:180205/1‑149 (MQ=255)
                                                                                                           cGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACt                                        <  1:2136319/149‑1 (MQ=255)
                                                                                                            gTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACtt                                       >  2:980113/1‑149 (MQ=255)
                                                                                                             tAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTa                                      <  1:3488676/149‑1 (MQ=255)
                                                                                                             tAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTa                                      <  1:101826/149‑1 (MQ=255)
                                                                                                              aaGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTAt                                     <  2:3787545/149‑1 (MQ=255)
                                                                                                                ggggAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATcc                                   >  2:186938/1‑149 (MQ=255)
                                                                                                                ggggAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATcc                                   >  2:3390737/1‑149 (MQ=255)
                                                                                                                  ggAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAg                                 >  1:602347/1‑149 (MQ=255)
                                                                                                                  ggAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAg                                 >  2:2792611/1‑149 (MQ=255)
                                                                                                                     aTCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCGTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGAt                              >  2:153997/1‑149 (MQ=255)
                                                                                                                       ccGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTa                            <  2:4216760/149‑1 (MQ=255)
                                                                                                                       ccGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTa                            <  1:1440839/149‑1 (MQ=255)
                                                                                                                       ccGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTa                            <  1:4601586/149‑1 (MQ=255)
                                                                                                                         gCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATc                          <  1:2954542/149‑1 (MQ=255)
                                                                                                                         gCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATc                          <  1:1231580/149‑1 (MQ=255)
                                                                                                                          cTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCt                         >  1:3572031/1‑149 (MQ=255)
                                                                                                                           tttCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCtt                        <  1:181665/149‑1 (MQ=255)
                                                                                                                           tttCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCtt                        <  1:62644/149‑1 (MQ=255)
                                                                                                                            ttCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTg                       <  2:2270754/149‑1 (MQ=255)
                                                                                                                            ttCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTg                       <  2:215205/149‑1 (MQ=255)
                                                                                                                             tctcGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTgg                      <  2:3074067/149‑1 (MQ=255)
                                                                                                                             tctcGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTgg                      >  1:2885692/1‑149 (MQ=255)
                                                                                                                               tcGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGc                    <  2:1052367/149‑1 (MQ=255)
                                                                                                                                cGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCa                   <  2:3220143/149‑1 (MQ=255)
                                                                                                                                cGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCa                   <  2:1643352/149‑1 (MQ=255)
                                                                                                                                 ggCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCAt                  >  2:3660922/1‑149 (MQ=255)
                                                                                                                                    aaCAGTTTTACCATTTTATCTCGACTTCCTGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGAGAATATGGGAATAGATACGCGACTTATCAAGGATTATCTAGGGCATCGc               >  2:1185827/1‑149 (MQ=255)
                                                                                                                                        gTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAata           >  1:3894274/1‑149 (MQ=255)
                                                                                                                                           ttACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTc        >  1:3670674/1‑149 (MQ=255)
                                                                                                                                           ttACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTc        >  2:1349605/1‑149 (MQ=255)
                                                                                                                                            tACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCg       <  1:4035577/149‑1 (MQ=255)
                                                                                                                                                atatTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCAt   <  1:2774846/149‑1 (MQ=255)
                                                                                                                                                 tattatCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATa  >  2:2750498/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATA  >  CP009273/4530997‑4531290

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: