Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,111,193 G→T 20.3% A231A (GCC→GCA)  wcaK ← colanic acid biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,111,1930GT20.3% 797.0 / 162.5 295A231A (GCC→GCA) wcaKcolanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (123/112);  new base T (24/36);  total (147/148)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 303 total aligned reads displayed.

GTGGAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAG  >  CP009273/2111044‑2111342
                                                                                                                                                     |                                                                                                                                                     
gTGGAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCg                                                                                                                                                       >  2:4237900/1‑150 (MQ=255)
  ggAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGc                                                                                                                                                     >  2:375137/1‑150 (MQ=255)
   gagaGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCa                                                                                                                                                    <  2:4328990/150‑1 (MQ=255)
   gagaGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCa                                                                                                                                                    <  1:4218331/150‑1 (MQ=255)
     gagCGCAATCACCTGAGACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAAc                                                                                                                                                  <  2:2850006/150‑1 (MQ=255)
      agCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCAAca                                                                                                                                                 >  2:2109206/1‑149 (MQ=255)
        cgcAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTGTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTc                                                                                                                                               >  2:169957/1‑150 (MQ=255)
             tCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTccagcc                                                                                                                                          >  1:4427957/1‑150 (MQ=255)
             tCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTccagcc                                                                                                                                          >  1:1352879/1‑150 (MQ=255)
                       cccTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGTGGCAACGTCCAGCCAGTGTTGAAc                                                                                                                                <  1:214559/150‑1 (MQ=255)
                       cccTTCATCGAGAAGGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAAc                                                                                                                                <  2:159199/150‑1 (MQ=255)
                       ccattCATCGAGAATGCGATTGACCACCCCGTCAAAGGCTTTTTCTTACGATTGTTGAGTGGTGCAGAGACGTTTGTCTAACGTTGCCAGTTCGCGCAGTGTTATGGCCACCGTTTTCTGTTTTGCGGCAACGTACAGCCTGTGTTGAAc                                                                                                                                <  2:1101481/147‑1 (MQ=255)
                          ttCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGc                                                                                                                             >  2:1810602/1‑150 (MQ=255)
                           tCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCa                                                                                                                            >  1:932034/1‑150 (MQ=255)
                             aTCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATa                                                                                                                          >  1:1848183/1‑150 (MQ=255)
                             aTCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATa                                                                                                                          >  2:1627671/1‑150 (MQ=255)
                             aTCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATa                                                                                                                          <  2:4185355/150‑1 (MQ=255)
                             aTCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATa                                                                                                                          >  1:1052362/1‑150 (MQ=255)
                               cGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGc                                                                                                                        >  1:468710/1‑150 (MQ=255)
                                gagaATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCt                                                                                                                       <  2:4522588/150‑1 (MQ=255)
                                 agaATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGTGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCACCGTCCAGCCAGTGTTGAACGGCATAGCTg                                                                                                                      >  2:3895920/1‑150 (MQ=255)
                                 agaATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTg                                                                                                                      <  2:1832777/150‑1 (MQ=255)
                                  gaATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTgg                                                                                                                     >  2:901408/1‑150 (MQ=255)
                                   aaTGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGc                                                                                                                    <  1:451464/150‑1 (MQ=255)
                                   aaTGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGc                                                                                                                    >  2:4273290/1‑150 (MQ=255)
                                    aggCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCg                                                                                                                   <  2:1349826/148‑1 (MQ=255)
                                     tGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCgg                                                                                                                  >  2:3117486/1‑150 (MQ=255)
                                        gATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGAGGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGa                                                                                                               >  2:2084543/1‑150 (MQ=255)
                                          ttGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAg                                                                                                             >  2:3979734/1‑150 (MQ=255)
                                             accaccCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGtct                                                                                                          >  2:1867665/1‑150 (MQ=255)
                                              ccaccCCGGCAACGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGtctt                                                                                                         <  1:2846577/150‑1 (MQ=255)
                                              ccaccCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGtctt                                                                                                         >  1:3493506/1‑150 (MQ=255)
                                              ccaccCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGtctt                                                                                                         >  1:1871573/1‑150 (MQ=255)
                                                accCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGtcttct                                                                                                       >  2:2723421/1‑150 (MQ=255)
                                                  cccGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCtgt                                                                                                     <  1:4307394/150‑1 (MQ=255)
                                                   ccGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCtgtg                                                                                                    <  2:2342200/150‑1 (MQ=255)
                                                    cGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCtgtgt                                                                                                   >  1:235537/1‑150 (MQ=255)
                                                     ggCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCtgtgtg                                                                                                  <  1:1500494/150‑1 (MQ=255)
                                                     ggCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCtgtgtg                                                                                                  >  2:566475/1‑150 (MQ=255)
                                                     ggCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCtgtgtg                                                                                                  >  2:3149399/1‑150 (MQ=255)
                                                      gCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTgtgg                                                                                                 <  1:3722495/150‑1 (MQ=255)
                                                      gCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTgtgg                                                                                                 >  1:3207374/1‑150 (MQ=255)
                                                       cAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTGTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCCGGGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTgtggt                                                                                                >  1:12727/1‑150 (MQ=255)
                                                        aaaGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAAAGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTgtggtg                                                                                               >  1:1545101/1‑150 (MQ=255)
                                                         aaGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGa                                                                                              >  2:1536456/1‑150 (MQ=255)
                                                         aaGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGa                                                                                              >  2:3814781/1‑150 (MQ=255)
                                                          aGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGAt                                                                                             <  1:725692/150‑1 (MQ=255)
                                                           ggCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTCTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATc                                                                                            >  1:1803147/1‑150 (MQ=255)
                                                             cTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGa                                                                                          >  2:526709/1‑150 (MQ=255)
                                                              tttttCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGAc                                                                                         >  1:405573/1‑150 (MQ=255)
                                                                        gCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCCGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCAcgcg                                                                               >  1:3160979/1‑150 (MQ=255)
                                                                        gCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCAcgcg                                                                               >  2:1711697/1‑150 (MQ=255)
                                                                         cTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCgg                                                                              >  1:2113967/1‑150 (MQ=255)
                                                                          ttgttgAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGt                                                                             <  1:2596173/150‑1 (MQ=255)
                                                                            gttgAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTAt                                                                           <  1:3048936/150‑1 (MQ=255)
                                                                            gttgAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTAt                                                                           >  2:2019721/1‑150 (MQ=255)
                                                                             ttgAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTTATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATc                                                                          <  2:457294/150‑1 (MQ=255)
                                                                              tgAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCg                                                                         <  1:2174530/150‑1 (MQ=255)
                                                                                   ggtgCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGcc                                                                    <  1:2493575/150‑1 (MQ=255)
                                                                                     tgCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCAt                                                                  <  1:402362/150‑1 (MQ=255)
                                                                                       ccGAGACGTTTGTCAAACGGTGCCAGTTAGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAAAGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGt                                                                >  1:150816/1‑150 (MQ=255)
                                                                                        cGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGtt                                                               >  1:717989/1‑150 (MQ=255)
                                                                                         gagaCGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTc                                                              >  1:4126911/1‑150 (MQ=255)
                                                                                                ttGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACtttt                                                       <  1:1043650/150‑1 (MQ=255)
                                                                                                 tGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTg                                                      >  1:3246635/1‑150 (MQ=255)
                                                                                                      aaCGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCggtg                                                 <  1:3585492/150‑1 (MQ=255)
                                                                                                       aCGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCggtgg                                                <  2:2927735/150‑1 (MQ=255)
                                                                                                       aCGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCggtgg                                                <  1:3447991/150‑1 (MQ=255)
                                                                                                        cGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCggtggt                                               <  2:2496643/150‑1 (MQ=255)
                                                                                                        cGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCggtggt                                               >  2:2513025/1‑150 (MQ=255)
                                                                                                         ggTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGgtggtg                                              <  2:1167346/150‑1 (MQ=255)
                                                                                                          gTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGa                                             >  2:1466810/1‑150 (MQ=255)
                                                                                                           tGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGat                                            <  2:2413045/150‑1 (MQ=255)
                                                                                                            gCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGata                                           >  1:985088/1‑150 (MQ=255)
                                                                                                            gCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGata                                           >  1:2810627/1‑150 (MQ=255)
                                                                                                              cAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATAtt                                         >  1:3382088/1‑150 (MQ=255)
                                                                                                               aGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATtg                                        >  1:656572/1‑150 (MQ=255)
                                                                                                                   cgcgCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATtgctgc                                    <  2:3727732/150‑1 (MQ=255)
                                                                                                                    gcgcAGGGTAATGGCCACCTTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGAACAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCg                                   <  2:3443120/150‑1 (MQ=255)
                                                                                                                        aGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTtc                               >  2:792581/1‑150 (MQ=255)
                                                                                                                          ggTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTtcat                             <  2:3941629/150‑1 (MQ=255)
                                                                                                                          ggTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTtcat                             >  2:184211/1‑150 (MQ=255)
                                                                                                                            tAATGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTtcatca                           >  2:4489695/1‑150 (MQ=255)
                                                                                                                               tGGCCACCGTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGAt                        <  2:1768147/150‑1 (MQ=255)
                                                                                                                                    aCCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAg                   >  1:1113170/1‑150 (MQ=255)
                                                                                                                                      cGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCt                 >  1:2020240/1‑150 (MQ=255)
                                                                                                                                       gTTTTCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTg                <  2:66023/150‑1 (MQ=255)
                                                                                                                                       gTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTg                >  2:2590267/1‑150 (MQ=255)
                                                                                                                                       gTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTg                >  1:364834/1‑150 (MQ=255)
                                                                                                                                         tttCTGTTGTGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGAc              >  2:3508221/1‑150 (MQ=255)
                                                                                                                                          ttCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGAcc             <  1:786654/150‑1 (MQ=255)
                                                                                                                                             tgttgtGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGAt          >  1:3003183/1‑150 (MQ=255)
                                                                                                                                               ttgtGCTGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTc        <  2:2674899/150‑1 (MQ=255)
                                                                                                                                               ttgtGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTc        <  1:1775611/150‑1 (MQ=255)
                                                                                                                                               ttgtGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTc        <  1:2688092/150‑1 (MQ=255)
                                                                                                                                                 gtgCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTcgc      >  1:1725795/1‑150 (MQ=255)
                                                                                                                                                 gtgCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTcgc      >  2:1008471/1‑150 (MQ=255)
                                                                                                                                                  tgCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTcgcg     <  1:566475/150‑1 (MQ=255)
                                                                                                                                                    cGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCa   >  1:1854409/1‑150 (MQ=255)
                                                                                                                                                    cGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCa   <  1:988102/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGGAGAGCGCAATCACCTGATACCCTTCATCGAGAATGCGATTGACCACCCCGGCAAAGGCTTTTTCATACGCTTGTTGAGTGGTGCCGAGACGTTTGTCAAACGGTGCCAGTTCGCGCAGGGTAATGGCCACCGTTTTCTGTTGTGCGGCAACGTCCAGCCAGTGTTGAACGGCATAGCTGGCGGTGAAGTCTTCTGTGTGGTGATCGACCAGCCACGCGGTATCGACGCCATGTTCCACTTTTGCGGTGGTGATATTGCTGCGTTTCATCAGATCAAAGCTGACCGATTCGCGCAG  >  CP009273/2111044‑2111342

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