Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 518,820 C→A 100% A35E (GCG→GAG)  rhsD → Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273518,8200CA98.1% 678.4 / 2.6 211A35E (GCG→GAG) rhsDRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (132/75);  total (134/77)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 231 total aligned reads displayed.

TTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGA  >  CP009273/518673‑518966
                                                                                                                                                   |                                                                                                                                                  
ttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtggcgt                                                                                                                                                  >  1:977531/1‑150 (MQ=255)
ttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtggagt                                                                                                                                                  <  2:335712/150‑1 (MQ=255)
ttgCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAACGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCggcgtggagt                                                                                                                                                  <  2:4760524/150‑1 (MQ=255)
    tAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTc                                                                                                                                              >  2:3049714/1‑150 (MQ=255)
      acacAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCgg                                                                                                                                            >  2:3232351/1‑150 (MQ=255)
        acaGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtg                                                                                                                                          >  2:2731847/1‑150 (MQ=255)
         caGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtgt                                                                                                                                         <  2:2249702/150‑1 (MQ=255)
          aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtgtg                                                                                                                                        >  1:1108310/1‑150 (MQ=255)
          aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGgtgtg                                                                                                                                        >  1:174682/1‑150 (MQ=255)
               tatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCcgg                                                                                                                                   >  2:3622736/1‑150 (MQ=255)
               tatGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCcgg                                                                                                                                   >  1:508583/1‑150 (MQ=255)
                atGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACCCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCcggg                                                                                                                                  >  1:2256725/1‑150 (MQ=255)
                       aTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATg                                                                                                                           >  2:145581/1‑150 (MQ=255)
                        tAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGa                                                                                                                          >  2:2140868/1‑150 (MQ=255)
                         aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGAc                                                                                                                         >  2:924474/1‑150 (MQ=255)
                         aaCTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGAc                                                                                                                         >  2:3641330/1‑150 (MQ=255)
                             aaaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCg                                                                                                                     <  1:466034/150‑1 (MQ=255)
                              aaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCgg                                                                                                                    <  1:3878716/150‑1 (MQ=255)
                              aaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCgg                                                                                                                    >  2:4549944/1‑150 (MQ=255)
                              aaaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCgg                                                                                                                    >  1:26858/1‑150 (MQ=255)
                               aaGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCggg                                                                                                                   >  2:3006666/1‑150 (MQ=255)
                                aGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGc                                                                                                                  >  2:4478408/1‑150 (MQ=255)
                                 gggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCa                                                                                                                 >  2:3154452/1‑150 (MQ=255)
                                 gggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCa                                                                                                                 >  2:2828771/1‑150 (MQ=255)
                                 gggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCa                                                                                                                 <  1:4551590/150‑1 (MQ=255)
                                  ggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCaa                                                                                                                >  2:475089/1‑150 (MQ=255)
                                      cTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCg                                                                                                            <  2:3971586/150‑1 (MQ=255)
                                        ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCTGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGt                                                                                                          >  1:4020259/1‑150 (MQ=255)
                                        ttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGt                                                                                                          >  2:420255/1‑150 (MQ=255)
                                          atatGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGTCTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTaa                                                                                                        >  1:2582629/1‑150 (MQ=255)
                                          atatGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTaa                                                                                                        >  1:903186/1‑150 (MQ=255)
                                              gAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATcc                                                                                                    <  2:1204646/150‑1 (MQ=255)
                                               aGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCg                                                                                                   >  2:2006295/1‑150 (MQ=255)
                                                 cGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgct                                                                                                 <  1:3655406/150‑1 (MQ=255)
                                                  ggAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCgctg                                                                                                >  2:850616/1‑150 (MQ=255)
                                                     aaaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGctgctg                                                                                             >  2:1335674/1‑150 (MQ=255)
                                                      aaCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTgg                                                                                            <  1:2344572/150‑1 (MQ=255)
                                                       aCCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTggg                                                                                           >  1:3256992/1‑150 (MQ=255)
                                                        ccAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTgggg                                                                                          >  1:1690040/1‑150 (MQ=255)
                                                        ccAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTgggg                                                                                          >  1:4074267/1‑150 (MQ=255)
                                                        ccAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTgggg                                                                                          <  2:423934/150‑1 (MQ=255)
                                                          aGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGc                                                                                        >  1:3324729/1‑150 (MQ=255)
                                                            cggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGa                                                                                      >  1:1563363/1‑150 (MQ=255)
                                                              gcgcGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAg                                                                                    >  1:3154861/1‑150 (MQ=255)
                                                              gcgcGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAg                                                                                    >  2:1677230/1‑150 (MQ=255)
                                                                 cgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTCAATCCGCTGCTGGGGGCGAAGGtg                                                                                 >  2:4307703/1‑150 (MQ=255)
                                                                 cgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtg                                                                                 >  1:166456/1‑150 (MQ=255)
                                                                   tCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTGAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgct                                                                               <  1:2928832/150‑1 (MQ=255)
                                                                   tCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgct                                                                               <  1:4181017/150‑1 (MQ=255)
                                                                   tCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgct                                                                               <  2:1555721/150‑1 (MQ=255)
                                                                    cAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgctg                                                                              >  2:3133509/1‑150 (MQ=255)
                                                                     aGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGtgctgc                                                                             >  2:2796487/1‑150 (MQ=255)
                                                                        gagaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCg                                                                          <  2:1885990/150‑1 (MQ=255)
                                                                        gagaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCg                                                                          <  2:1112627/150‑1 (MQ=255)
                                                                         agaTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTCAATCCGCTGCTGGGGGCGAAGGTGCTGCCCgg                                                                         >  1:3439318/1‑150 (MQ=255)
                                                                                aCTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACg                                                                  >  1:2469343/1‑150 (MQ=255)
                                                                                aCTCAGTATGGCGGTCCCATGGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACg                                                                  <  2:1396642/150‑1 (MQ=255)
                                                                                   cAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGAc                                                               >  2:2605220/1‑150 (MQ=255)
                                                                                      tATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCtt                                                            >  2:2775144/1‑150 (MQ=255)
                                                                                      tATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCtt                                                            <  1:792600/150‑1 (MQ=255)
                                                                                        tGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgc                                                          >  2:3779586/1‑150 (MQ=255)
                                                                                        tGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgc                                                          >  1:4482649/1‑150 (MQ=255)
                                                                                         ggCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcg                                                         <  2:4154951/150‑1 (MQ=255)
                                                                                         ggCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTgcg                                                         <  1:4658134/150‑1 (MQ=255)
                                                                                            ggTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTg                                                      >  1:3645297/1‑150 (MQ=255)
                                                                                            ggTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTg                                                      >  1:2403303/1‑150 (MQ=255)
                                                                                               cccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgccc                                                   <  2:445521/150‑1 (MQ=255)
                                                                                               cccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgccc                                                   >  2:1977958/1‑150 (MQ=255)
                                                                                                ccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccg                                                  >  2:621967/1‑150 (MQ=255)
                                                                                                ccATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccg                                                  >  2:776333/1‑150 (MQ=255)
                                                                                                 cATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccgg                                                 >  2:2509579/1‑150 (MQ=255)
                                                                                                 cATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccgg                                                 >  2:2681741/1‑150 (MQ=255)
                                                                                                     gTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcccg                                             >  1:1048679/1‑150 (MQ=255)
                                                                                                     gTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTgcccggcccg                                             <  2:1646455/150‑1 (MQ=255)
                                                                                                       ccAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCt                                           >  1:2323944/1‑150 (MQ=255)
                                                                                                         aGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCAGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGc                                         >  1:1253570/1‑150 (MQ=255)
                                                                                                          gggTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGcc                                        <  1:4689975/150‑1 (MQ=255)
                                                                                                           gggtCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGGGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCg                                       <  2:833866/147‑1 (MQ=255)
                                                                                                                  cAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtc                                >  1:4395833/1‑150 (MQ=255)
                                                                                                                   aGGTGTAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtct                               >  2:1266218/1‑150 (MQ=255)
                                                                                                                    ggTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctc                              <  1:1342062/150‑1 (MQ=255)
                                                                                                                      tgtAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATtctctc                            <  2:4843409/150‑1 (MQ=255)
                                                                                                                        tAAGAATTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTccc                          >  1:4940094/1‑150 (MQ=255)
                                                                                                                             aTTGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCAcc                     >  1:4025767/1‑150 (MQ=255)
                                                                                                                              ttGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCt                    >  1:4296686/1‑150 (MQ=255)
                                                                                                                              ttGGCGCGCCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCt                    <  1:3726277/150‑1 (MQ=255)
                                                                                                                                 gcgcgcCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAca                 <  2:3322489/150‑1 (MQ=255)
                                                                                                                                  cgcgcCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCCTTCTCTCCCGCACCTAcag                >  2:4676898/1‑150 (MQ=255)
                                                                                                                                  cgcgcCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcag                >  2:4792099/1‑150 (MQ=255)
                                                                                                                                  cgcgcCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGGAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcag                >  2:3077458/1‑150 (MQ=255)
                                                                                                                                   gcgcCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagc               >  1:2141024/1‑150 (MQ=255)
                                                                                                                                    cgcCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagca              <  1:2585979/150‑1 (MQ=255)
                                                                                                                                    cgcCCACCGGCGTGGAGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagca              <  2:4074267/150‑1 (MQ=255)
                                                                                                                                            ggcgtggcgtGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCgg      <  1:183854/150‑1 (MQ=255)
                                                                                                                                            ggcgtggagtGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCgg      >  2:3392993/1‑150 (MQ=255)
                                                                                                                                             gcgtggagtgCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGa     <  1:829071/150‑1 (MQ=255)
                                                                                                                                              cgtggagtgCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGAc    >  1:587269/1‑150 (MQ=255)
                                                                                                                                               gtggagtgCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACg   <  2:2200932/150‑1 (MQ=255)
                                                                                                                                               gtggagtgCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACg   <  2:2256725/150‑1 (MQ=255)
                                                                                                                                                tggagtgCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGa  >  1:4238289/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTGCTCGGTGTGTCCGGGCGGGATGACTTCGGGCAACCCGGTAAATCCGCTGCTGGGGGCGAAGGTGCTGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGA  >  CP009273/518673‑518966

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: ‑