Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,531,461 G→A 45.8% intergenic (+85/‑393) 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,4610GA45.8% 107.9 / inf 265intergenic (+85/‑393)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 G (75/68);  new base A (66/55);  total (141/124)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 268 total aligned reads displayed.

GCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTCCGCACTTT  >  CP009273/4531314‑4531608
                                                                                                                                                   |                                                                                                                                                   
gCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGtt                                                                                                                                                   <  1:4269016/150‑1 (MQ=255)
gCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATAtt                                                                                                                                                   <  1:1894903/150‑1 (MQ=255)
   gggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATattatt                                                                                                                                                <  1:4031815/150‑1 (MQ=255)
   gggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATattatt                                                                                                                                                >  2:1615485/1‑150 (MQ=255)
   gggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTAtt                                                                                                                                                >  2:823857/1‑150 (MQ=255)
    ggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTg                                                                                                                                               <  1:3840992/150‑1 (MQ=255)
    ggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTg                                                                                                                                               <  1:599465/150‑1 (MQ=255)
     gCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGt                                                                                                                                              <  2:2050444/150‑1 (MQ=255)
        tttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAg                                                                                                                                           >  2:3302174/1‑150 (MQ=255)
         ttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAgt                                                                                                                                          <  2:2416380/150‑1 (MQ=255)
         ttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAgt                                                                                                                                          >  2:3870579/1‑150 (MQ=255)
          tttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAgtg                                                                                                                                         <  1:3428060/150‑1 (MQ=255)
              cGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGAt                                                                                                                                     <  1:1561572/150‑1 (MQ=255)
                   tCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTa                                                                                                                                >  1:3329039/1‑150 (MQ=255)
                    cTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTaa                                                                                                                               <  2:1281586/150‑1 (MQ=255)
                     tGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAAt                                                                                                                              <  1:3585952/150‑1 (MQ=255)
                      gggATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAAtt                                                                                                                             <  1:458710/150‑1 (MQ=255)
                        gATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTAt                                                                                                                           <  1:2009416/150‑1 (MQ=255)
                        gATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTAt                                                                                                                           >  2:2942684/1‑150 (MQ=255)
                         aTAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATc                                                                                                                          >  1:2939178/1‑150 (MQ=255)
                           agagCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCat                                                                                                                        >  2:1587089/1‑150 (MQ=255)
                                   agGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACa                                                                                                                >  2:2994591/1‑150 (MQ=255)
                                   agGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACa                                                                                                                >  1:1307701/1‑150 (MQ=255)
                                       cgacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTaa                                                                                                            <  2:1772094/150‑1 (MQ=255)
                                        gacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAg                                                                                                           <  2:3491071/150‑1 (MQ=255)
                                         acAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGt                                                                                                          >  1:1763788/1‑150 (MQ=255)
                                         acAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGt                                                                                                          >  1:1516461/1‑150 (MQ=255)
                                           aGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTa                                                                                                        >  2:3030118/1‑150 (MQ=255)
                                           aGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTa                                                                                                        >  2:3030219/1‑150 (MQ=255)
                                             cGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAc                                                                                                      >  1:663308/1‑150 (MQ=255)
                                              gTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAcc                                                                                                     >  2:471122/1‑150 (MQ=255)
                                              gTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAcc                                                                                                     <  2:705603/150‑1 (MQ=255)
                                                 aCGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCata                                                                                                  <  1:1447813/150‑1 (MQ=255)
                                                  cGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCatat                                                                                                 <  2:1664874/150‑1 (MQ=255)
                                                  cGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCatat                                                                                                 >  2:3381682/1‑150 (MQ=255)
                                                    cTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATcc                                                                                               <  1:126216/150‑1 (MQ=255)
                                                        tttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATAc                                                                                           >  2:4124761/1‑150 (MQ=255)
                                                         ttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTGATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACa                                                                                          >  2:3641780/1‑150 (MQ=255)
                                                          tataGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACaa                                                                                         <  2:422446/150‑1 (MQ=255)
                                                           ataGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACaaa                                                                                        >  2:276434/1‑150 (MQ=255)
                                                            taGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACaaaa                                                                                       >  1:303281/1‑150 (MQ=255)
                                                             agagATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAAt                                                                                      <  1:4586817/150‑1 (MQ=255)
                                                                   aGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATg                                                                                <  2:3433823/150‑1 (MQ=255)
                                                                       gagaTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGttt                                                                            <  2:1760353/150‑1 (MQ=255)
                                                                       gagaTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGttt                                                                            <  1:4404819/150‑1 (MQ=255)
                                                                        agaTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTa                                                                           >  2:1049168/1‑150 (MQ=255)
                                                                         gaTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTAt                                                                          <  1:1058919/150‑1 (MQ=255)
                                                                            aCCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGtt                                                                       <  2:269476/150‑1 (MQ=255)
                                                                            aCCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGtt                                                                       <  1:4153647/150‑1 (MQ=255)
                                                                             ccAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttc                                                                      >  2:1904058/1‑150 (MQ=255)
                                                                              cAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttct                                                                     <  2:1051764/150‑1 (MQ=255)
                                                                                gggATGGTGTTAAAAACTACATTATTCTTACGATTTATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttcttc                                                                   <  2:653077/150‑1 (MQ=255)
                                                                                gggATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTTTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttcttc                                                                   >  2:2850772/1‑150 (MQ=255)
                                                                                  gATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCaa                                                                 >  2:469600/1‑150 (MQ=255)
                                                                                   aTGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCaaa                                                                >  1:1471071/1‑150 (MQ=255)
                                                                                    tGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCaaaa                                                               >  2:3921112/1‑150 (MQ=255)
                                                                                     ggTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaat                                                              <  1:4485337/150‑1 (MQ=255)
                                                                                     ggTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGCTTATGTTCTTCAaaat                                                              <  1:3379795/150‑1 (MQ=255)
                                                                                      gtgtTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaata                                                             >  1:443391/1‑150 (MQ=255)
                                                                                        gtTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaataaa                                                           >  2:1615741/1‑150 (MQ=255)
                                                                                           aaaaaCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGGTCTTCAAAataaataa                                                        >  2:4007183/1‑150 (MQ=255)
                                                                                           aaaaaCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAataaataa                                                        >  1:3473599/1‑150 (MQ=255)
                                                                                              aaCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACa                                                     <  2:756975/150‑1 (MQ=255)
                                                                                              aaCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACa                                                     >  1:3300100/1‑150 (MQ=255)
                                                                                              aaCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACa                                                     >  1:2607236/1‑150 (MQ=255)
                                                                                               aCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACaa                                                    >  1:3583330/1‑150 (MQ=255)
                                                                                                 tACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACaaaa                                                  >  1:3221387/1‑150 (MQ=255)
                                                                                                 tACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACaaaa                                                  >  1:2180036/1‑150 (MQ=255)
                                                                                                  aCATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAt                                                 <  2:3292381/150‑1 (MQ=255)
                                                                                                     ttattCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcat                                              <  1:3154775/150‑1 (MQ=255)
                                                                                                         tCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcattcat                                          <  2:20302/150‑1 (MQ=255)
                                                                                                            tACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATaaa                                       <  1:3359192/150‑1 (MQ=255)
                                                                                                            tACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATaaa                                       <  1:4296438/150‑1 (MQ=255)
                                                                                                             aCGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAt                                      >  1:4432134/1‑150 (MQ=255)
                                                                                                              cGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAtt                                     >  2:2442968/1‑150 (MQ=255)
                                                                                                                atatatCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTa                                   >  1:3016796/1‑150 (MQ=255)
                                                                                                                  atatCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTaca                                 >  2:4472186/1‑150 (MQ=255)
                                                                                                                  atatCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTaca                                 <  1:396861/150‑1 (MQ=255)
                                                                                                                   tatCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTAAATAATCATAATAATCATTAAGTTAACCATATCCATACAAAATACAAAGGATTATGATCTACAAAATAAATAAACAAAATAATTAATTAGTTTacac                                >  1:2524736/1‑150 (MQ=255)
                                                                                                                       aaaaaTGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCa                            <  1:3768400/150‑1 (MQ=255)
                                                                                                                        aaaaTGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCAc                           >  2:321066/1‑150 (MQ=255)
                                                                                                                        aaaaTGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCAc                           >  1:1289667/1‑150 (MQ=255)
                                                                                                                          aaTGATTAACAATTTATTAATTGATTTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAACATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACtt                         <  2:1510862/150‑1 (MQ=255)
                                                                                                                           aTGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTa                        <  2:107013/150‑1 (MQ=255)
                                                                                                                            tGATTAACCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTaa                       <  1:140744/150‑1 (MQ=255)
                                                                                                                                 aaCCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAtct                  >  1:11343/1‑150 (MQ=255)
                                                                                                                                  aCCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAtctc                 >  2:53290/1‑150 (MQ=255)
                                                                                                                                  aCCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAtctc                 <  1:3641226/150‑1 (MQ=255)
                                                                                                                                  aCCATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAtctc                 <  1:1796066/150‑1 (MQ=255)
                                                                                                                                   ccATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTcc                <  1:3740751/150‑1 (MQ=255)
                                                                                                                                   ccATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTcc                <  1:1647291/150‑1 (MQ=255)
                                                                                                                                    cATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCt               <  1:974960/150‑1 (MQ=255)
                                                                                                                                    cATTTATTAATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCt               <  1:276434/150‑1 (MQ=255)
                                                                                                                                       ttattaATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGtt            >  2:112744/1‑150 (MQ=255)
                                                                                                                                        tattaATTGATATTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGttt           >  1:748099/1‑150 (MQ=255)
                                                                                                                                         attaATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTc          <  2:4598823/150‑1 (MQ=255)
                                                                                                                                         attaATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTc          >  1:1279548/1‑150 (MQ=255)
                                                                                                                                          ttaATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTcc         >  2:3787436/1‑150 (MQ=255)
                                                                                                                                           taATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTCCg        >  2:4315108/1‑150 (MQ=255)
                                                                                                                                                gATATTATTGTAAGTGTGATTCTTAATTATAATTATAAACATTAAGTTAACCATATACATACAAAATACAATGGTTTATGTTCTTCAAAATAAATTAACAAAATCATTCATAAATTTACACATCCCTTAAAATCTACTGTTTCCGCACtt   <  2:2258193/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTTCCGCACTTT  >  CP009273/4531314‑4531608

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: