Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,531,451 A→G 9.6% 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,4510AG 9.6% 445.0 / 32.5 145intergenic (+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 (46/85);  new base G (6/8);  total (52/93)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 159 total aligned reads displayed.

ACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTT  >  CP009273/4531302‑4531599
                                                                                                                                                     |                                                                                                                                                    
aCCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATtta                                                                                                                                                      <  2:843457/150‑1 (MQ=255)
aCCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATtta                                                                                                                                                      <  2:818428/150‑1 (MQ=255)
  cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATttatt                                                                                                                                                    <  2:2738908/150‑1 (MQ=255)
  cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATttatt                                                                                                                                                    <  2:124790/150‑1 (MQ=255)
  cGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATttatt                                                                                                                                                    >  1:4786850/1‑150 (MQ=255)
   gCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATttatta                                                                                                                                                   <  2:1045091/150‑1 (MQ=255)
     cAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAt                                                                                                                                                 <  2:1186628/150‑1 (MQ=255)
      aGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAtt                                                                                                                                                <  1:3035260/150‑1 (MQ=255)
      aGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAAtt                                                                                                                                                <  2:3793334/150‑1 (MQ=255)
        cAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCCCGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACAGTATTCTTACGATATATCAAAAATGATTCACCATGTATTAATTGa                                                                                                                                              >  1:139342/1‑150 (MQ=255)
        cAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGa                                                                                                                                              >  1:1875872/1‑150 (MQ=255)
           tGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGt                                                                                                                                           <  1:3369720/150‑1 (MQ=255)
            gCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGtt                                                                                                                                          >  2:3156607/1‑150 (MQ=255)
             cAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTa                                                                                                                                         <  1:4463970/150‑1 (MQ=255)
             cAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTa                                                                                                                                         <  1:996986/150‑1 (MQ=255)
              aGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTAt                                                                                                                                        >  2:1965181/1‑150 (MQ=255)
              aGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTAt                                                                                                                                        >  1:1573754/1‑150 (MQ=255)
                   gTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTaa                                                                                                                                   <  1:2951056/150‑1 (MQ=255)
                    tttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAg                                                                                                                                  >  1:1286710/1‑150 (MQ=255)
                    tttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAg                                                                                                                                  >  1:1286663/1‑150 (MQ=255)
                     ttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAgt                                                                                                                                 <  1:4088815/150‑1 (MQ=255)
                       ttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAgtgt                                                                                                                               <  1:1826237/150‑1 (MQ=255)
                       ttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTAGAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAgtgt                                                                                                                               <  2:948202/150‑1 (MQ=255)
                          cGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGAt                                                                                                                            <  1:2464330/150‑1 (MQ=255)
                           ggCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGAtt                                                                                                                           <  2:695885/150‑1 (MQ=255)
                           ggCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGAtt                                                                                                                           >  1:2957101/1‑150 (MQ=255)
                             cATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCt                                                                                                                         <  1:1779881/150‑1 (MQ=255)
                                 tGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAAt                                                                                                                     <  1:2371131/150‑1 (MQ=255)
                                 tGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAAt                                                                                                                     <  1:2447618/150‑1 (MQ=255)
                                   ggATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTa                                                                                                                   <  1:2258337/150‑1 (MQ=255)
                                     aTAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATc                                                                                                                 >  2:4542635/1‑150 (MQ=255)
                                       agagCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCat                                                                                                               <  2:3320661/150‑1 (MQ=255)
                                        gagCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCata                                                                                                              >  2:38018/1‑150 (MQ=255)
                                         agCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCataa                                                                                                             <  2:2221764/150‑1 (MQ=255)
                                         agCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCataa                                                                                                             <  1:3799932/150‑1 (MQ=255)
                                          gCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCataat                                                                                                            <  1:476937/150‑1 (MQ=255)
                                          gCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCataat                                                                                                            <  1:3726658/150‑1 (MQ=255)
                                           ccAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCataata                                                                                                           <  1:4522724/150‑1 (MQ=255)
                                             agagGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGTGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATaaa                                                                                                         <  1:4049272/150‑1 (MQ=255)
                                             agagGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATaaa                                                                                                         <  2:4348385/150‑1 (MQ=255)
                                             agagGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATaaa                                                                                                         <  2:4348541/150‑1 (MQ=255)
                                              gagGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAAc                                                                                                        <  1:1124382/150‑1 (MQ=255)
                                               agGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACa                                                                                                       >  1:1279616/1‑150 (MQ=255)
                                                ggACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACAt                                                                                                      <  2:2729460/150‑1 (MQ=255)
                                                    gacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAg                                                                                                  <  1:3531819/150‑1 (MQ=255)
                                                    gacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAg                                                                                                  <  1:3531671/150‑1 (MQ=255)
                                                    gacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAg                                                                                                  <  2:1680044/150‑1 (MQ=255)
                                                              cGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCatat                                                                                        <  1:3845601/150‑1 (MQ=255)
                                                              cGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCatat                                                                                        <  2:1859028/150‑1 (MQ=255)
                                                                 tGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCa                                                                                     <  2:1927509/150‑1 (MQ=255)
                                                                 tGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCa                                                                                     <  2:682062/150‑1 (MQ=255)
                                                                  gTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCAt                                                                                    <  1:742661/150‑1 (MQ=255)
                                                                    tttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATAc                                                                                  >  2:327809/1‑150 (MQ=255)
                                                                    tttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAAAATTAAGTTAACAATATCCATAc                                                                                  >  2:3261576/1‑150 (MQ=255)
                                                                     ttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACa                                                                                 >  1:4298803/1‑150 (MQ=255)
                                                                     ttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACa                                                                                 >  1:4298815/1‑150 (MQ=255)
                                                                     ttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACa                                                                                 <  2:1786611/150‑1 (MQ=255)
                                                                     ttATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACa                                                                                 >  1:4301423/1‑150 (MQ=255)
                                                                      tataGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCAATACaa                                                                                >  2:3251626/1‑150 (MQ=255)
                                                                      tataGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAAAATTAAGTTAACCATATCCATACaa                                                                                >  2:3250125/1‑150 (MQ=255)
                                                                             aTAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACaa                                                                         <  1:4300325/150‑1 (MQ=255)
                                                                             aTAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACaa                                                                         <  2:187200/150‑1 (MQ=255)
                                                                                   gagaTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGttt                                                                   <  2:3393724/150‑1 (MQ=255)
                                                                                     gaTACCATGGATTGTGTTAAAAACTACATTATTCTTAAGATATATCAAATATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAACATATCCATACAAAATACAATGGTTTAt                                                                 <  1:1892804/150‑1 (MQ=255)
                                                                                     gaTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTAt                                                                 <  1:1048582/150‑1 (MQ=255)
                                                                                     gaTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTAt                                                                 <  1:3120613/150‑1 (MQ=255)
                                                                                     gaTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTAt                                                                 <  1:3388020/150‑1 (MQ=255)
                                                                                      aTACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATg                                                                <  1:846397/150‑1 (MQ=255)
                                                                                           aGGGATTGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttctt                                                           <  2:2926662/150‑1 (MQ=255)
                                                                                           aGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGttctt                                                           >  2:3763138/1‑150 (MQ=255)
                                                                                             ggATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCa                                                         >  2:3171696/1‑150 (MQ=255)
                                                                                                 ggTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaat                                                     <  2:1322961/150‑1 (MQ=255)
                                                                                                 ggTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaat                                                     <  2:1322959/150‑1 (MQ=255)
                                                                                                  gtgtTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAaaata                                                    >  1:4316424/1‑150 (MQ=255)
                                                                                                       aaaaaCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAataaataa                                               >  2:1090292/1‑150 (MQ=255)
                                                                                                        aaaaCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAAtaaataaa                                              <  1:981963/150‑1 (MQ=255)
                                                                                                          aaCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACa                                            >  1:1531380/1‑150 (MQ=255)
                                                                                                          aaCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACa                                            >  1:4170943/1‑150 (MQ=255)
                                                                                                           aCTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACaa                                           <  1:3496423/150‑1 (MQ=255)
                                                                                                             tACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACaaaa                                         >  2:3030157/1‑150 (MQ=255)
                                                                                                              aCATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAt                                        <  1:4745574/150‑1 (MQ=255)
                                                                                                                attattCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtca                                      <  2:3307144/150‑1 (MQ=255)
                                                                                                                 ttattCTTACGATATATCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcat                                     <  2:688139/150‑1 (MQ=255)
                                                                                                                 ttattCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAAtcat                                     <  1:105668/150‑1 (MQ=255)
                                                                                                                       ttACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATaa                               <  2:760538/150‑1 (MQ=255)
                                                                                                                        tACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATaaa                              >  2:764484/1‑150 (MQ=255)
                                                                                                                         aCGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAt                             <  2:3449402/150‑1 (MQ=255)
                                                                                                                          cGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAtt                            >  2:468940/1‑150 (MQ=255)
                                                                                                                          cGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAAtt                            >  2:470065/1‑150 (MQ=255)
                                                                                                                               tatCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTacac                       <  1:3302963/150‑1 (MQ=255)
                                                                                                                                atCAAAAATGATTAACCATTTGTTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTAcaca                      >  2:131888/1‑150 (MQ=255)
                                                                                                                                atCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTAcaca                      <  1:4542635/150‑1 (MQ=255)
                                                                                                                                  cAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATc                    >  2:3617574/1‑150 (MQ=255)
                                                                                                                                     aaaTGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACt                 >  2:3393507/1‑150 (MQ=255)
                                                                                                                                       aTGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTa               <  2:1718643/150‑1 (MQ=255)
                                                                                                                                        tGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTaa              <  1:3680715/150‑1 (MQ=255)
                                                                                                                                           ttAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAt           <  1:4310320/150‑1 (MQ=255)
                                                                                                                                             aaCCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAAtct         >  1:1776615/1‑150 (MQ=255)
                                                                                                                                                  tttATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGt    >  1:2145904/1‑150 (MQ=255)
                                                                                                                                                    tattaATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGttt  <  1:4693813/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCATATCCATACAAAATACAATGGTTTATGTTCTTCAAAATAAATAAACAAAATCATTCATAAATTTACACATCACTTAAAATCTCCTGTTT  >  CP009273/4531302‑4531599

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: