Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,049,008 A→T 13.3% L231Q (CTG→CAG)  glnL ← sensory histidine kinase in two‑component regulatory system with GlnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,049,0080AT13.3% 1001.1 / 123.0 331L231Q (CTG→CAG) glnLsensory histidine kinase in two‑component regulatory system with GlnG
Reads supporting (aligned to +/- strand):  ref base A (133/152);  new base T (21/23);  total (156/175)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

Only 100 of 336 total aligned reads displayed.

CCGCCAGCCGGTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCGGCCCCAACAGACGGTCGA  >  CP009273/4048859‑4049155
                                                                                                                                                     |                                                                                                                                                   
ccgccAGCCGGTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCa                                                                                                                                                     <  2:3669975/150‑1 (MQ=255)
   ccAGCCGGTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGtt                                                                                                                                                  >  1:1347189/1‑150 (MQ=255)
   ccAGCCGGTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGtt                                                                                                                                                  <  1:192912/150‑1 (MQ=255)
     aGCCGGTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTcc                                                                                                                                                <  1:4728256/150‑1 (MQ=255)
     aGCCGGTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTcc                                                                                                                                                <  1:1385192/150‑1 (MQ=255)
         ggTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCAGTGCGCACAATATTCAGCAAGACCGGTTCGAGTTGATCCGGGTCGGGCGCCCGTGCCGGtc                                                                                                                                            >  2:4884168/1‑149 (MQ=255)
          gTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAg                                                                                                                                           <  1:1624204/150‑1 (MQ=255)
               acTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTgg                                                                                                                                      <  2:4950106/149‑1 (MQ=255)
                cTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTggg                                                                                                                                     >  2:4000763/1‑150 (MQ=255)
                 tCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGa                                                                                                                                    <  2:1996177/150‑1 (MQ=255)
                  cGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGAt                                                                                                                                   <  2:2209071/150‑1 (MQ=255)
                  cGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGAt                                                                                                                                   <  2:242301/150‑1 (MQ=255)
                     cGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCTGTTCCGGTAGGCTGGGATCGt                                                                                                                                <  2:419486/150‑1 (MQ=255)
                     cGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGt                                                                                                                                >  2:1052851/1‑150 (MQ=255)
                        gtAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAAt                                                                                                                             >  2:4543066/1‑150 (MQ=255)
                        gtAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAAt                                                                                                                             <  2:4106781/150‑1 (MQ=255)
                              tCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGaa                                                                                                                       <  2:628105/150‑1 (MQ=255)
                               cAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAAt                                                                                                                      <  2:1217093/150‑1 (MQ=255)
                                 gTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCa                                                                                                                    <  2:5022795/150‑1 (MQ=255)
                                  ttGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCaa                                                                                                                   >  1:4225662/1‑150 (MQ=255)
                                  ttGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCaa                                                                                                                   <  1:4080643/150‑1 (MQ=255)
                                  ttGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCaa                                                                                                                   >  1:1218761/1‑150 (MQ=255)
                                   tGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAAc                                                                                                                  >  1:3860074/1‑150 (MQ=255)
                                       aCGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCa                                                                                                              >  2:4708733/1‑150 (MQ=255)
                                           ggTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGtt                                                                                                          <  2:3393000/150‑1 (MQ=255)
                                            gTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTg                                                                                                         >  1:2265993/1‑150 (MQ=255)
                                               cGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCTGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTcc                                                                                                      <  2:4086433/150‑1 (MQ=255)
                                               cGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTcc                                                                                                      <  1:3600643/150‑1 (MQ=255)
                                                gggTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCTGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCg                                                                                                     >  1:734493/1‑150 (MQ=255)
                                                gggTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCg                                                                                                     <  2:4228396/150‑1 (MQ=255)
                                                gggTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCg                                                                                                     >  2:1624941/1‑150 (MQ=255)
                                                   tACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCa                                                                                                  >  2:3758811/1‑150 (MQ=255)
                                                     cGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGt                                                                                                <  1:3106827/150‑1 (MQ=255)
                                                      gCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGtt                                                                                               >  2:277348/1‑150 (MQ=255)
                                                      gCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGtt                                                                                               >  2:3963059/1‑150 (MQ=255)
                                                          aaTGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCAt                                                                                           <  1:3514802/150‑1 (MQ=255)
                                                          aaTGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCAt                                                                                           <  2:3108828/150‑1 (MQ=255)
                                                           aTGATTTCACCGCCTTCCGTCCCAATAGCCTGTAGCTCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTTATCAGTGTCGTGCGCCAGTTCCGGTAGTCTGGGATCGTAATCACTAATCAACCTAACGTTGTCCGGCAGTTCAATc                                                                                          <  2:803401/150‑1 (MQ=255)
                                                             gATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATAGTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGa                                                                                        <  2:1273371/150‑1 (MQ=255)
                                                                ttCACCGCCTTCCGGCCCCAGCGCCTGGAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGacac                                                                                     <  1:92146/150‑1 (MQ=255)
                                                                    ccGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGc                                                                                 >  1:3219882/1‑150 (MQ=255)
                                                                      gCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGt                                                                               <  1:2564833/150‑1 (MQ=255)
                                                                       ccTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGtt                                                                              >  1:2897215/1‑150 (MQ=255)
                                                                        cTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTa                                                                             >  2:3972081/1‑150 (MQ=255)
                                                                         ttCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTAc                                                                            <  2:1422737/150‑1 (MQ=255)
                                                                         ttCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTAc                                                                            >  1:164772/1‑150 (MQ=255)
                                                                         tgccGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTAc                                                                            <  2:4254623/148‑1 (MQ=255)
                                                                            cagCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCAc                                                                         <  1:2804159/148‑1 (MQ=255)
                                                                            cGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCAc                                                                         <  1:1732745/150‑1 (MQ=255)
                                                                            cGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCAc                                                                         >  1:2504156/1‑150 (MQ=255)
                                                                                cccAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGt                                                                     <  1:2060739/150‑1 (MQ=255)
                                                                                  cAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTc                                                                   <  2:3570611/150‑1 (MQ=255)
                                                                                  cAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTc                                                                   <  1:1052851/150‑1 (MQ=255)
                                                                                  cAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTc                                                                   >  1:3300383/1‑150 (MQ=255)
                                                                                       ccTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCa                                                              >  1:5069354/1‑150 (MQ=255)
                                                                                         tGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACt                                                            <  2:1789833/150‑1 (MQ=255)
                                                                                         tGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACt                                                            <  2:4398661/150‑1 (MQ=255)
                                                                                          gTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACtt                                                           <  1:4313649/150‑1 (MQ=255)
                                                                                          gTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACtt                                                           <  2:4312458/150‑1 (MQ=255)
                                                                                          gTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACtt                                                           <  2:2163252/150‑1 (MQ=255)
                                                                                          gTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACtt                                                           <  1:154519/150‑1 (MQ=255)
                                                                                               gcATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGa                                                      <  1:3730043/150‑1 (MQ=255)
                                                                                               gcATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGa                                                      <  1:2693498/150‑1 (MQ=255)
                                                                                                 aTTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAAt                                                    <  2:2460075/150‑1 (MQ=255)
                                                                                                   tGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATAc                                                  >  2:3625053/1‑150 (MQ=255)
                                                                                                    gcgcACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACt                                                 <  1:5178999/150‑1 (MQ=255)
                                                                                                    gcgcACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACt                                                 <  1:2829109/150‑1 (MQ=255)
                                                                                                       cacaATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTc                                              >  2:5105720/1‑150 (MQ=255)
                                                                                                          aaTATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGt                                           >  1:844142/1‑150 (MQ=255)
                                                                                                           atatTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTa                                          <  2:2931329/150‑1 (MQ=255)
                                                                                                           atatTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTa                                          <  1:1867069/150‑1 (MQ=255)
                                                                                                              ttCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAAcg                                       <  2:2259347/150‑1 (MQ=255)
                                                                                                                cAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAAcgcg                                     >  2:1374976/1‑150 (MQ=255)
                                                                                                                 aGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAAcgcgc                                    >  1:1130612/1‑150 (MQ=255)
                                                                                                                  gCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGGTACCACGCGTTCAGCCACTGTGTGACTACTTTCGGTCAcgcgcg                                   >  2:2985217/1‑150 (MQ=255)
                                                                                                                     aGACCTGTTCAATTTGATCCGGGTCGTGCGCCTGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTAc                                >  1:3045172/1‑150 (MQ=255)
                                                                                                                      gACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTAcc                               <  1:1427543/150‑1 (MQ=255)
                                                                                                                       aCCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCg                              <  1:5224803/150‑1 (MQ=255)
                                                                                                                         cTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCggg                            >  1:3123513/1‑150 (MQ=255)
                                                                                                                         cTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCggg                            >  1:4851802/1‑150 (MQ=255)
                                                                                                                            ttCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAg                         <  1:298583/150‑1 (MQ=255)
                                                                                                                              cAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCt                       >  2:1041379/1‑150 (MQ=255)
                                                                                                                              cAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCt                       >  1:3878597/1‑150 (MQ=255)
                                                                                                                              cAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCt                       <  1:2123770/150‑1 (MQ=255)
                                                                                                                                aTTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGc                     >  1:3422167/1‑150 (MQ=255)
                                                                                                                                aTTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGc                     >  1:3188183/1‑150 (MQ=255)
                                                                                                                                aTTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTGCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTga                     >  2:4928494/1‑149 (MQ=255)
                                                                                                                                 tttGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCg                    <  1:66675/150‑1 (MQ=255)
                                                                                                                                  ttGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCgg                   <  2:561448/150‑1 (MQ=255)
                                                                                                                                     aTCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCGGccc                >  2:3216776/1‑150 (MQ=255)
                                                                                                                                       ccGGGTCGTGCGCCTGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCGGCCCCa              >  2:2875293/1‑150 (MQ=255)
                                                                                                                                         gggTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCGGCCCCAAc            <  1:4254188/150‑1 (MQ=255)
                                                                                                                                         gggTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCGGCCCCAAc            >  2:1162023/1‑150 (MQ=255)
                                                                                                                                             cGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCGGCCCCAACAGAc        >  1:2309597/1‑150 (MQ=255)
                                                                                                                                              gTGCGCCTGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCGGCCCCAACAGACg       >  2:4740232/1‑150 (MQ=255)
                                                                                                                                              gTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCGGCCCCAACAGACg       <  1:4209934/150‑1 (MQ=255)
                                                                                                                                                gcgcCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCGGCCCCAACAGACGGt     <  1:1952708/150‑1 (MQ=255)
                                                                                                                                                  gcCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCGGCCCCAACAGACGGTCg   >  2:1720408/1‑150 (MQ=255)
                                                                                                                                                  gcCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCGGCCCCAACAGACGGTCg   <  1:3996596/150‑1 (MQ=255)
                                                                                                                                                   ccAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCGGCCCCAACAGACGGTCGa  >  2:4768534/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGCCAGCCGGTAGCGCTCGCCGTGTAAGGTCAGTTGAAACGCGGTGCGGGTACGCAGAATGATTTCACCGCCTTCCGGCCCCAGCGCCTGTAGCGCATTGCGCACAATATTCAGCAAGACCTGTTCAATTTGATCCGGGTCGTGCGCCAGTTCCGGTAGGCTGGGATCGTAATCACGAATCAACCGCACGTTGTCCGGCAGTTCCATCGACACCAGCGTTACCACGCGTTCAGCCACTTTGTGAATACTTTCGGTAACGCGCGTACCGGGCAGCTGCGGCCCCAACAGACGGTCGA  >  CP009273/4048859‑4049155

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: