Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,531,451 A→G 15.9% intergenic (+75/‑403) fimB → / → fimE tyrosine recombinase/inversion of on/off regulator of fimA/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,4510AG15.9% 792.1 / 129.8 282intergenic (+75/‑403)fimB/fimEtyrosine recombinase/inversion of on/off regulator of fimA/tyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base A (103/134);  new base G (29/16);  total (132/150)
Fisher's exact test for biased strand distribution p-value = 1.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

Only 100 of 302 total aligned reads displayed.

ACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTC  >  CP009273/4531302‑4531600
                                                                                                                                                     |                                                                                                                                                     
 ccGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATttat                                                                                                                                                      >  1:4367464/1‑150 (MQ=255)
  cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATttatt                                                                                                                                                     >  2:737208/1‑150 (MQ=255)
  cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATttatt                                                                                                                                                     <  1:2008531/150‑1 (MQ=255)
  cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTGTATAGAGATAGTTGAGATACCAGGGATGGTGTTCAAAACTACAGTCTTCTTACGATATATCAAAAATGATTAACCATttgtt                                                                                                                                                     >  2:3648579/1‑150 (MQ=255)
   gCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATttatta                                                                                                                                                    <  2:1302017/150‑1 (MQ=255)
    ccAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTaa                                                                                                                                                   <  2:714194/150‑1 (MQ=255)
     cAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAt                                                                                                                                                  <  1:2503259/150‑1 (MQ=255)
      aGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAtt                                                                                                                                                 <  2:5032182/150‑1 (MQ=255)
       gCAATGCAGTGCGTTTTTACGGCATCTGTGATAGAGCCAGAGGACGACAGCGTCACGCTTTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                                                                                                                <  1:4462728/150‑1 (MQ=255)
       gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTg                                                                                                                                                >  1:268242/1‑150 (MQ=255)
       gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                                                                                                                >  1:583788/1‑150 (MQ=255)
          aTGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATg                                                                                                                                             >  1:3586294/1‑150 (MQ=255)
                ggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTg                                                                                                                                       >  2:549640/1‑150 (MQ=255)
                ggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTg                                                                                                                                       >  2:5023642/1‑150 (MQ=255)
                     ttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAgt                                                                                                                                  <  2:3248678/150‑1 (MQ=255)
                         aCGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGa                                                                                                                              >  1:3915788/1‑150 (MQ=255)
                          cGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGAt                                                                                                                             >  2:1840911/1‑150 (MQ=255)
                          cGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGAt                                                                                                                             <  1:656537/150‑1 (MQ=255)
                          cGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGAt                                                                                                                             <  1:4617197/150‑1 (MQ=255)
                           ggCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGAtt                                                                                                                            >  2:4134114/1‑150 (MQ=255)
                             cATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCt                                                                                                                          >  1:4745752/1‑150 (MQ=255)
                              aTCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCtt                                                                                                                         >  1:1000642/1‑150 (MQ=255)
                              aTCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCtt                                                                                                                         >  1:1116523/1‑150 (MQ=255)
                                cTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTaa                                                                                                                       >  2:2514145/1‑150 (MQ=255)
                                cTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTaa                                                                                                                       <  1:1146092/150‑1 (MQ=255)
                                 tGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAAt                                                                                                                      <  1:1080480/150‑1 (MQ=255)
                                  gggATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAAtt                                                                                                                     <  2:3157246/150‑1 (MQ=255)
                                           ccAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCataata                                                                                                            >  1:3617228/1‑150 (MQ=255)
                                            cAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCAtaataa                                                                                                           >  2:4932082/1‑150 (MQ=255)
                                              gagGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAAc                                                                                                         >  2:882552/1‑150 (MQ=255)
                                                ggACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACAt                                                                                                       <  2:3063003/150‑1 (MQ=255)
                                                cgacgacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACAt                                                                                                       >  2:784434/2‑150 (MQ=255)
                                                 gacgacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACAtt                                                                                                      <  1:4169973/150‑1 (MQ=255)
                                                 gacgacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACAtt                                                                                                      >  1:225845/1‑150 (MQ=255)
                                                   cgacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTaa                                                                                                    <  2:3658374/150‑1 (MQ=255)
                                                     acAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGt                                                                                                  <  1:5106289/150‑1 (MQ=255)
                                                       aGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTa                                                                                                >  2:3508331/1‑150 (MQ=255)
                                                         cGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAc                                                                                              <  2:1928861/150‑1 (MQ=255)
                                                          gTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAcc                                                                                             >  1:2830604/1‑150 (MQ=255)
                                                            cACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCat                                                                                           <  1:31578/150‑1 (MQ=255)
                                                            cACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCat                                                                                           <  1:1492100/150‑1 (MQ=255)
                                                             aCGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCata                                                                                          <  2:4816919/150‑1 (MQ=255)
                                                              cGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCatat                                                                                         <  1:1741945/150‑1 (MQ=255)
                                                              cGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCatat                                                                                         <  1:2301016/150‑1 (MQ=255)
                                                               gCTGTTTTATAGAGATAGTTGAGATACCAGGGATGCTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATc                                                                                        >  1:4295929/1‑150 (MQ=255)
                                                                  gTTTTATAGAGATAGTTGAGATACCAGGGATTGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCAt                                                                                     <  2:3986700/150‑1 (MQ=255)
                                                                    tttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATAc                                                                                   <  1:2770034/150‑1 (MQ=255)
                                                                     ttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACa                                                                                  >  1:504474/1‑150 (MQ=255)
                                                                     ttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACa                                                                                  <  2:1631242/150‑1 (MQ=255)
                                                                      tataGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATAAAGTTAACCATATCCATACaa                                                                                 >  2:3833256/1‑150 (MQ=255)
                                                                         agagATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAAt                                                                              <  1:3787260/150‑1 (MQ=255)
                                                                          gagaTAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATa                                                                             <  1:2996513/150‑1 (MQ=255)
                                                                          gagaTAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATa                                                                             <  1:5070634/150‑1 (MQ=255)
                                                                               aGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATg                                                                        <  1:759498/150‑1 (MQ=255)
                                                                                   gagaTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGttt                                                                    >  2:2132529/1‑150 (MQ=255)
                                                                                         ccAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttc                                                              <  1:3280875/150‑1 (MQ=255)
                                                                                          cAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttct                                                             >  2:667464/1‑150 (MQ=255)
                                                                                           atggATGGTGTTAAAAAATACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTATTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttctt                                                            <  2:4830236/148‑1 (MQ=255)
                                                                                           aGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttctt                                                            >  1:2910336/1‑150 (MQ=255)
                                                                                           aGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttctt                                                            >  2:3181626/1‑150 (MQ=255)
                                                                                            gggATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTACTTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttcttc                                                           <  1:4963014/150‑1 (MQ=255)
                                                                                            gggATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttcttc                                                           <  1:1918295/150‑1 (MQ=255)
                                                                                             ggATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCa                                                          <  1:2376665/150‑1 (MQ=255)
                                                                                              gATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCaa                                                         >  2:4015007/1‑150 (MQ=255)
                                                                                              gATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCaa                                                         <  2:3223465/150‑1 (MQ=255)
                                                                                                    gtTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaataaa                                                   >  2:537433/1‑150 (MQ=255)
                                                                                                      tAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAaataaata                                                 <  2:3759328/150‑1 (MQ=255)
                                                                                                      tAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAaataaata                                                 >  1:1765322/1‑150 (MQ=255)
                                                                                                            cTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAATTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACaaa                                           >  1:3134528/1‑150 (MQ=255)
                                                                                                             tACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACaaaa                                          >  2:432651/1‑150 (MQ=255)
                                                                                                              aCATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAt                                         >  1:2146607/1‑150 (MQ=255)
                                                                                                                attattCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtca                                       >  2:4537904/1‑150 (MQ=255)
                                                                                                                attattCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtca                                       <  2:3576375/150‑1 (MQ=255)
                                                                                                                 ttattCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcat                                      <  1:1436682/150‑1 (MQ=255)
                                                                                                                   attCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcattc                                    <  1:2607043/150‑1 (MQ=255)
                                                                                                                   attCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcattc                                    <  1:2606958/150‑1 (MQ=255)
                                                                                                                    ttattACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcattca                                   <  2:2669761/147‑1 (MQ=255)
                                                                                                                    ttCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcattca                                   <  2:232260/150‑1 (MQ=255)
                                                                                                                    ttCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcattca                                   <  2:2829751/150‑1 (MQ=255)
                                                                                                                      cTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATa                                 <  2:2248274/150‑1 (MQ=255)
                                                                                                                         aCGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAt                              <  1:2724361/150‑1 (MQ=255)
                                                                                                                           gATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAttt                            <  1:524480/150‑1 (MQ=255)
                                                                                                                              atatCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTaca                         <  1:3540626/150‑1 (MQ=255)
                                                                                                                                  cAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATc                     >  2:3708917/1‑150 (MQ=255)
                                                                                                                                  cAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATc                     >  2:1136735/1‑150 (MQ=255)
                                                                                                                                  cAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATc                     >  1:1593034/1‑150 (MQ=255)
                                                                                                                                     aaaTGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACt                  >  1:4001042/1‑150 (MQ=255)
                                                                                                                                      aaTGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACtt                 >  2:2281414/1‑150 (MQ=255)
                                                                                                                                      aaTGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACtt                 >  1:1614169/1‑150 (MQ=255)
                                                                                                                                       aTGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTa                <  2:4976012/150‑1 (MQ=255)
                                                                                                                                            tAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAtc           <  2:2100636/150‑1 (MQ=255)
                                                                                                                                              aCCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAtctc         >  2:2753663/1‑150 (MQ=255)
                                                                                                                                               ccATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTcc        >  2:4680779/1‑150 (MQ=255)
                                                                                                                                                cATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCt       <  2:1566658/150‑1 (MQ=255)
                                                                                                                                                  tttATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGt     >  2:3617377/1‑150 (MQ=255)
                                                                                                                                                   ttattaATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGtt    <  1:1265379/150‑1 (MQ=255)
                                                                                                                                                    tattaATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGttt   >  2:2888227/1‑150 (MQ=255)
                                                                                                                                                    tattaATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGttt   >  2:724299/1‑150 (MQ=255)
                                                                                                                                                     attaATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTc  <  2:2993535/150‑1 (MQ=255)
                                                                                                                                                     attaATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAATTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTc  <  2:3420048/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTC  >  CP009273/4531302‑4531600

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: