Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,320,198 C→A 6.9% D33Y (GAT→TAT)  ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,320,1980CA 6.9% 714.7 / 6.1 217D33Y (GAT→TAT) ftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  ref base C (60/129);  new base A (26/2);  total (86/131)
Fisher's exact test for biased strand distribution p-value = 4.54e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 221 total aligned reads displayed.

CAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTGTGTTAAAAACAGCGTCAGGA  >  CP009273/3320049‑3320347
                                                                                                                                                     |                                                                                                                                                     
 aaCAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATcc                                                                                                                                                      <  1:2233262/150‑1 (MQ=255)
  aCAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCa                                                                                                                                                     <  1:4075210/150‑1 (MQ=255)
    aGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCAcc                                                                                                                                                   <  1:16603/150‑1 (MQ=255)
         aTCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGTCCTCTTGTTGGAAGGTAGAGTAATACACCATACg                                                                                                                                              >  1:2621507/1‑150 (MQ=255)
         aTCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTGATTGACCTCGTGTAGGAAGGTATATTAATACACCATAca                                                                                                                                              >  1:2621380/1‑149 (MQ=255)
          tCCAGTAATTTCGGATCCTGACCCGGAATGTAAGTGGTATAACGGGTACTATCTTTCGTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACgg                                                                                                                                             <  1:3007176/150‑1 (MQ=255)
          tCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTTAACTCTTGTAGGAAGGTATAGTAATACACCTTACgg                                                                                                                                             >  1:1306980/1‑150 (MQ=255)
             aGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGCAGGTAGAGTAATCCACCTTACGGCCa                                                                                                                                          >  2:2526418/1‑150 (MQ=255)
             aGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCa                                                                                                                                          <  1:4147016/150‑1 (MQ=255)
                aaTTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTa                                                                                                                                       <  2:143103/150‑1 (MQ=255)
                aaTTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTa                                                                                                                                       <  1:272618/150‑1 (MQ=255)
                 aTTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAg                                                                                                                                      <  2:2835145/150‑1 (MQ=255)
                  tttCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTTACCACTTGTAGGAAGGTAGATTAATCCACCTTACGGCCCATAGa                                                                                                                                     >  1:4342974/1‑150 (MQ=255)
                   ttCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGAc                                                                                                                                    <  1:3844083/150‑1 (MQ=255)
                    tCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACt                                                                                                                                   <  2:2050702/150‑1 (MQ=255)
                      ggATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCAATAGACTCg                                                                                                                                 >  2:1468151/1‑150 (MQ=255)
                       gATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCACTTGTAGGAAGGTAGAGTAATCCACCGTACGGCCAATAGACACGc                                                                                                                                >  1:2070593/1‑150 (MQ=255)
                        aTCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCt                                                                                                                               <  1:615206/150‑1 (MQ=255)
                        aTCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCt                                                                                                                               <  2:431435/150‑1 (MQ=255)
                        aTCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGAACTCTTGTAGGAAGGTATAGTAATCCACCATAAGGCCCTTATACTCGCt                                                                                                                               >  1:3553058/1‑150 (MQ=255)
                           cTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTGATTTAACACTTGTAGGAAGGTAGAGTAATCCACCATAAGGCCAATAGACTCGCAggg                                                                                                                            >  2:3273127/1‑150 (MQ=255)
                              aaCCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTTACCTCTTGTAGGAAGGTAGAGTAATACACCATACGGCCCTTAGACTCGCGGGGccc                                                                                                                         >  1:2892926/1‑150 (MQ=255)
                              aaCCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGccc                                                                                                                         <  2:2008400/150‑1 (MQ=255)
                                 cGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGTAGGTAGAGTAATCCACCATACGGCCAATAGACTCGCGGGGCCCaaa                                                                                                                      >  1:3629911/1‑150 (MQ=255)
                                     aTGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGctc                                                                                                                  >  2:270989/1‑150 (MQ=255)
                                      tGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGctct                                                                                                                 <  1:430438/150‑1 (MQ=255)
                                      tGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGctct                                                                                                                 <  1:321057/150‑1 (MQ=255)
                                        tAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGAACTCTTGTAGGAAGGTAGATTAATCCACCTTACGGCCAATAGACTCGCTGGGCCCAAAGCTCTgg                                                                                                               >  1:1411510/1‑150 (MQ=255)
                                          aGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTTAACTCTTTTAGGAAGGTAGAGTAATACACCATAAGGCCAATAGACTCGCTGGGCCCAAAGCTCTGGaa                                                                                                             >  1:1929344/1‑150 (MQ=255)
                                          aGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCACTTGTAGGAAGGTAGAGTAATCCACCGTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGaa                                                                                                             >  1:2400425/1‑150 (MQ=255)
                                                   aCGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACAt                                                                                                    <  1:3087760/150‑1 (MQ=255)
                                                    cGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATc                                                                                                   <  1:502793/150‑1 (MQ=255)
                                                    cGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATc                                                                                                   <  1:1339328/150‑1 (MQ=255)
                                                     ggTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCa                                                                                                  <  1:2913876/150‑1 (MQ=255)
                                                     ggTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCa                                                                                                  <  2:430893/150‑1 (MQ=255)
                                                      gTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAg                                                                                                 <  2:923199/150‑1 (MQ=255)
                                                      gTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAg                                                                                                 <  2:1464995/150‑1 (MQ=255)
                                                         aCTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGcac                                                                                              <  2:2862630/150‑1 (MQ=255)
                                                         aCTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGcac                                                                                              >  2:3935964/1‑150 (MQ=255)
                                                         aCTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTGTTTTAACGCTTGTAGGAAGGTAGATTAATACACCATACGGCCAATATACTCGCGGGGCCCCAAGCACTGGAAAAAAGACAACAGcac                                                                                              >  1:3816982/1‑150 (MQ=255)
                                                            aTCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGAACACATGTAGGAAGGTAGAGTAATCCACCATACGGCCCTTAGACACGCTGGGCCCCAAGCTCTGGAATACAGACATCAGCACAAc                                                                                           >  2:3480015/1‑150 (MQ=255)
                                                             tctttcttGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACg                                                                                          <  2:1720487/150‑1 (MQ=255)
                                                             tctttcttGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTAATGACCTCATGTAGGAAGGTAGAGTAATCCACCTTAAGGCCCATAGACACGCTGGGCCCCAAGCTCTGGAATACTGACCTCAGCACAACg                                                                                          >  2:841809/1‑150 (MQ=255)
                                                              ctttcttGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACgg                                                                                         <  1:2296156/150‑1 (MQ=255)
                                                                ttcttGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTTTTGTCCACTTTTTTGAATGTAGAGTAATACACCATACGGCCCATAGACACGCGGGGCCCAAATCACTGGGATAAAGACAACAGCACaaacgcc                                                                                       >  1:3318175/1‑145 (MQ=255)
                                                                  cttGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCACTTGTAGGAAGGTAGAGTAATACACCATAAGGCCATTATACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACCGCAAt                                                                                     >  2:3649702/1‑150 (MQ=255)
                                                                   ttGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTTTACTCATGTATGAAGGTATAGTAAAACACCATAAGGCCAATAGAATAGCTGGGGCCAAACCTCTGGTATAATGACATCAGCACCACGGccaat                                                                                    >  1:4160113/1‑146 (MQ=255)
                                                                       tAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATACACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGAcca                                                                                <  1:1201922/150‑1 (MQ=255)
                                                                       tAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATACACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGAcca                                                                                <  1:1709851/150‑1 (MQ=255)
                                                                        aaCGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGAccag                                                                               <  2:2789757/150‑1 (MQ=255)
                                                                           gTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGAccagcca                                                                            <  1:2804921/150‑1 (MQ=255)
                                                                           gTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGAccagcca                                                                            <  1:1071532/150‑1 (MQ=255)
                                                                           gTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGAccagcca                                                                            <  1:511309/150‑1 (MQ=255)
                                                                              gATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGATTAATCCACCATACGGCCAATAGACTCGCGGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCagag                                                                         >  1:1010309/1‑150 (MQ=255)
                                                                              gATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCagag                                                                         <  1:2258450/150‑1 (MQ=255)
                                                                                tttCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTTACCTCTTTTATGAAGGTAGAGTAATACAACTTACGGCCATTATACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAAACGCAATTACCAGCCAGAGTa                                                                       >  2:3700008/1‑150 (MQ=255)
                                                                                    aCGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAg                                                                   <  2:1988217/150‑1 (MQ=255)
                                                                                      gTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGt                                                                 <  2:2124316/150‑1 (MQ=255)
                                                                                        ccGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGttt                                                               <  2:3721226/150‑1 (MQ=255)
                                                                                         cGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGtttt                                                              <  2:4032320/150‑1 (MQ=255)
                                                                                          gTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGttttt                                                             <  1:1276908/150‑1 (MQ=255)
                                                                                           ttGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCACAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTc                                                            <  1:2510928/150‑1 (MQ=255)
                                                                                           ttGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTc                                                            <  1:2046706/150‑1 (MQ=255)
                                                                                           ttGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTc                                                            <  1:2941465/150‑1 (MQ=255)
                                                                                            tGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCg                                                           <  2:3557473/150‑1 (MQ=255)
                                                                                             gATACGCGCTTCACGAACCTGGTCGTTATTTTCCTCTTGTAGGAAGGTAGAGTAAACCACCTTACGGCCCTTAGACTCGCTGGGCCCAAAGCGCTGGAATACTGACATCAGCACAACGGCCATGACCAGCCAGAGTATTAGGTTTTTCGc                                                          >  2:3337711/1‑150 (MQ=255)
                                                                                                aCGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCAt                                                       <  2:624826/150‑1 (MQ=255)
                                                                                                aCGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCAt                                                       <  2:2091026/150‑1 (MQ=255)
                                                                                                 cgcgCTTCACGAACCTGGTCGTTATTTAACTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGAATCGCAGGGCCCAAAGCGCTGGAATACTTACATCAGCACAACGGCAATTAACAGCCAGAGTATTAGGTTTTTAGCCATg                                                      >  2:1017835/1‑150 (MQ=255)
                                                                                                    gcTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCa                                                   <  2:2621380/150‑1 (MQ=255)
                                                                                                    gcTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCa                                                   <  2:2621507/150‑1 (MQ=255)
                                                                                                     cTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCAc                                                  <  1:2564934/150‑1 (MQ=255)
                                                                                                     cTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCAc                                                  <  1:89797/150‑1 (MQ=255)
                                                                                                      ttCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACt                                                 <  1:2651873/150‑1 (MQ=255)
                                                                                                          cGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAg                                             <  2:2892926/150‑1 (MQ=255)
                                                                                                           gAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATACACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATAAAGACAAAAGCACAACGGCAAAGACCAGCCCGAGTATTATGTTTTTCGCCATGTAACTCaacg                                            >  2:1030560/1‑148 (MQ=255)
                                                                                                             aCCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGa                                          <  2:3811865/150‑1 (MQ=255)
                                                                                                                  gTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAAc                                     <  2:4339000/150‑1 (MQ=255)
                                                                                                                  gTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAAc                                     <  1:3861824/150‑1 (MQ=255)
                                                                                                                  gTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAAc                                     <  2:2166735/150‑1 (MQ=255)
                                                                                                                   tCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAAcc                                    <  2:3719836/150‑1 (MQ=255)
                                                                                                                    cGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACct                                   <  1:287254/150‑1 (MQ=255)
                                                                                                                       tATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCtt                                >  2:3124166/1‑150 (MQ=255)
                                                                                                                        aTTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCtta                               <  1:2754739/150‑1 (MQ=255)
                                                                                                                           gACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCttatta                            >  2:2941568/1‑150 (MQ=255)
                                                                                                                                ctTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACt                       >  2:4103814/1‑150 (MQ=255)
                                                                                                                                  tGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCACAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCTAGGGATTAACCTCTTATTACAACcgt                     <  2:1195447/150‑1 (MQ=255)
                                                                                                                                    tAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTgtgt                   >  1:2135593/1‑150 (MQ=255)
                                                                                                                                      ggAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTGTGTTa                 >  1:3167048/1‑150 (MQ=255)
                                                                                                                                       gAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTGTGTTaa                >  2:463051/1‑150 (MQ=255)
                                                                                                                                        aaGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTGTGTTaaa               <  1:3732199/150‑1 (MQ=255)
                                                                                                                                          ggTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTGTGTTaaaaa             >  2:3209019/1‑150 (MQ=255)
                                                                                                                                           gTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTGTGTTAAAAAc            >  1:2082602/1‑150 (MQ=255)
                                                                                                                                            tAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTGTGTTAAAAACa           <  1:4244446/150‑1 (MQ=255)
                                                                                                                                             agagTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTGTGTTAAAAACAg          >  2:113653/1‑150 (MQ=255)
                                                                                                                                                 tAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTGTGTTAAATACAGCGTc      <  1:1530975/150‑1 (MQ=255)
                                                                                                                                                   aTCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTGTGTTAAAAACAGCGTCAg    <  1:2126829/150‑1 (MQ=255)
                                                                                                                                                    tCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTGTGTTAAAAACAGCGTCAgg   <  1:47714/150‑1 (MQ=255)
                                                                                                                                                    tCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTGTGTTAAAAACAGCGTCAgg   <  2:417661/150‑1 (MQ=255)
                                                                                                                                                    tCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTGTGTTAAAAACAGCGTCAgg   >  2:618883/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTTATTACAACTGTGTTAAAAACAGCGTCAGGA  >  CP009273/3320049‑3320347

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: