Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,977,650 A→T 17.6% intergenic (+97/+62) aslB → / ← aslA putative AslA‑specific sulfatase‑maturating enzyme/putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,977,6500AT17.6% 856.0 / 167.0 323intergenic (+97/+62)aslB/aslAputative AslA‑specific sulfatase‑maturating enzyme/putative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base A (125/141);  new base T (24/33);  total (149/174)
Fisher's exact test for biased strand distribution p-value = 5.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 331 total aligned reads displayed.

GTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATC  >  CP009273/3977501‑3977799
                                                                                                                                                     |                                                                                                                                                     
gTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAACTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTa                                                                                                                                                       >  2:3135355/1‑150 (MQ=255)
gTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTa                                                                                                                                                       >  2:3284835/1‑150 (MQ=255)
gTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTa                                                                                                                                                       <  1:3156490/150‑1 (MQ=255)
    cccGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCt                                                                                                                                                   <  1:3487589/150‑1 (MQ=255)
      cGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTgg                                                                                                                                                 <  1:2423196/150‑1 (MQ=255)
      cGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTgg                                                                                                                                                 <  1:1772379/150‑1 (MQ=255)
        gCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGAt                                                                                                                                               <  1:4347341/150‑1 (MQ=255)
        gCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGAt                                                                                                                                               <  1:4511074/150‑1 (MQ=255)
        gCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGAt                                                                                                                                               <  1:2016904/150‑1 (MQ=255)
        gCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGAt                                                                                                                                               <  1:1439304/150‑1 (MQ=255)
         ccAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATa                                                                                                                                              >  2:3062108/1‑150 (MQ=255)
          cAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATaa                                                                                                                                             >  1:1540996/1‑150 (MQ=255)
          cAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATaa                                                                                                                                             >  2:445858/1‑150 (MQ=255)
            gCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGa                                                                                                                                           <  2:3641989/150‑1 (MQ=255)
            gCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATaata                                                                                                                                           <  2:325813/150‑3 (MQ=255)
            gCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGa                                                                                                                                           <  2:4117469/150‑1 (MQ=255)
              gACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATg                                                                                                                                         >  2:2948554/1‑150 (MQ=255)
                cATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCg                                                                                                                                       <  1:4197967/150‑1 (MQ=255)
                   tATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCa                                                                                                                                    >  1:4561321/1‑150 (MQ=255)
                      gcatgcGCATTTGCTGGTGGTGAGTCAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgca                                                                                                                                 >  1:4428666/1‑150 (MQ=255)
                         tgcGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcg                                                                                                                              <  2:1028985/150‑1 (MQ=255)
                         tgcGCAGTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcg                                                                                                                              >  2:740180/1‑150 (MQ=255)
                           cgcATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcgca                                                                                                                            >  1:159206/1‑150 (MQ=255)
                            gcATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAgcatcgcat                                                                                                                           >  1:4272582/1‑150 (MQ=255)
                                 tGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGc                                                                                                                      >  2:2699473/1‑150 (MQ=255)
                                   cTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCaa                                                                                                                    >  1:539351/1‑150 (MQ=255)
                                     ggtggtGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAg                                                                                                                  <  1:4320260/150‑1 (MQ=255)
                                        ggtgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCa                                                                                                               <  2:4493651/150‑1 (MQ=255)
                                        ggtgAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCa                                                                                                               <  2:3083613/150‑1 (MQ=255)
                                           gAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGAt                                                                                                            <  1:2175884/150‑1 (MQ=255)
                                               aagtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCa                                                                                                        >  2:1188845/1‑150 (MQ=255)
                                                agtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAg                                                                                                       <  2:3404598/150‑1 (MQ=255)
                                                 gtaGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGc                                                                                                      <  1:625757/150‑1 (MQ=255)
                                                     aaaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATa                                                                                                  <  2:3511060/150‑1 (MQ=255)
                                                     aaaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATa                                                                                                  <  1:2126318/150‑1 (MQ=255)
                                                     aaaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATa                                                                                                  >  2:3986282/1‑150 (MQ=255)
                                                     aaaTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATa                                                                                                  <  1:2487052/150‑1 (MQ=255)
                                                       aTCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAgc                                                                                                >  1:1375871/1‑150 (MQ=255)
                                                          ggcggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGcc                                                                                             >  2:3833640/1‑150 (MQ=255)
                                                           gcggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCg                                                                                            <  1:3267997/150‑1 (MQ=255)
                                                            cggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCATCGATAGCGCCgg                                                                                           <  2:1768719/150‑1 (MQ=255)
                                                             ggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGc                                                                                          >  2:1337473/1‑150 (MQ=255)
                                                             ggcCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTAAGCGATAGCGCCGGc                                                                                          >  1:664819/1‑150 (MQ=255)
                                                              gccgccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCt                                                                                         <  2:3936316/150‑1 (MQ=255)
                                                              gccgccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCt                                                                                         <  2:1103195/150‑1 (MQ=255)
                                                               ccgccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCtt                                                                                        <  2:2147955/150‑1 (MQ=255)
                                                               ccgccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCtt                                                                                        >  1:779675/1‑150 (MQ=255)
                                                                 gccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAg                                                                                      >  1:1636547/1‑150 (MQ=255)
                                                                  ccTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGt                                                                                     <  1:573588/150‑1 (MQ=255)
                                                                      cGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGa                                                                                 <  1:2863544/150‑1 (MQ=255)
                                                                      cGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGa                                                                                 <  2:2504113/150‑1 (MQ=255)
                                                                        gTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGAtt                                                                               >  1:3208306/1‑150 (MQ=255)
                                                                         ttGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGAttt                                                                              <  2:610676/150‑1 (MQ=255)
                                                                           gATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTaa                                                                            <  2:2268641/150‑1 (MQ=255)
                                                                            aTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAAt                                                                           >  1:1862495/1‑150 (MQ=255)
                                                                              tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCt                                                                         <  1:4086368/150‑1 (MQ=255)
                                                                              tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCt                                                                         <  1:2940793/150‑1 (MQ=255)
                                                                               gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTg                                                                        <  2:3569516/150‑1 (MQ=255)
                                                                                ccGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTgc                                                                       <  1:1936169/150‑1 (MQ=255)
                                                                                       cggcgTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCgtg                                                                >  1:1007707/1‑150 (MQ=255)
                                                                                       cggcgTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCgtg                                                                >  1:1004684/1‑150 (MQ=255)
                                                                                        ggcgTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCgtgg                                                               <  2:434148/150‑1 (MQ=255)
                                                                                          cgTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCgtggtg                                                             >  1:3921294/1‑150 (MQ=255)
                                                                                           gTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGAGTTAATCTGCGCGCGtggtgg                                                            >  2:3440071/1‑150 (MQ=255)
                                                                                              aaCGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCTCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGata                                                         <  1:1149544/150‑1 (MQ=255)
                                                                                              aaCGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGata                                                         <  1:3357528/150‑1 (MQ=255)
                                                                                               aCGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGatat                                                        <  2:1981621/150‑1 (MQ=255)
                                                                                               aCGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGatat                                                        <  1:4514098/150‑1 (MQ=255)
                                                                                                   cTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATAttttt                                                    >  2:2770299/1‑150 (MQ=255)
                                                                                                    ttATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATAtttttt                                                   <  2:847894/150‑1 (MQ=255)
                                                                                                    ttATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATAtttttt                                                   >  1:3587322/1‑150 (MQ=255)
                                                                                                      aTCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCa                                                 >  2:4402513/1‑150 (MQ=255)
                                                                                                           gCCTACCTGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGAtc                                            <  2:318666/150‑1 (MQ=255)
                                                                                                              tACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTcc                                         >  2:242564/1‑150 (MQ=255)
                                                                                                                 aTGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCata                                      >  2:315280/1‑150 (MQ=255)
                                                                                                                  tGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCatat                                     >  2:885163/1‑150 (MQ=255)
                                                                                                                     tCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATAcg                                  <  2:2355042/150‑1 (MQ=255)
                                                                                                                             aTTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCAtt                          <  2:4480950/150‑1 (MQ=255)
                                                                                                                             aTTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCAtt                          <  1:1860905/150‑1 (MQ=255)
                                                                                                                              ttCAATAAATTGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCAttt                         <  1:2430042/150‑1 (MQ=255)
                                                                                                                              ttCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCAttt                         <  1:2708210/150‑1 (MQ=255)
                                                                                                                               tCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCTGGATCTCCATATACGCGTGCATTTc                        <  2:1209953/150‑1 (MQ=255)
                                                                                                                               tCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTc                        <  1:373876/150‑1 (MQ=255)
                                                                                                                               tCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTc                        <  1:1617249/150‑1 (MQ=255)
                                                                                                                                cAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCg                       <  2:3417987/150‑1 (MQ=255)
                                                                                                                                cAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCg                       <  1:4031955/150‑1 (MQ=255)
                                                                                                                                   tAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTc                    <  1:4589307/150‑1 (MQ=255)
                                                                                                                                    aaaTTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCt                   >  1:3455809/1‑150 (MQ=255)
                                                                                                                                     aaTTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTg                  <  2:4272582/150‑1 (MQ=255)
                                                                                                                                      aTTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGt                 <  1:1188853/150‑1 (MQ=255)
                                                                                                                                       ttGCAGCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTa                <  2:1314910/150‑1 (MQ=255)
                                                                                                                                       ttGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCTGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTa                <  2:1471178/150‑1 (MQ=255)
                                                                                                                                        tGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAg               >  2:4197489/1‑150 (MQ=255)
                                                                                                                                        tGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAg               <  1:3062113/150‑1 (MQ=255)
                                                                                                                                            gCGTTCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTACAGGATCTCCATAAACGCGAGCAATTCGGTCTGTATCGGt           >  1:3361358/1‑150 (MQ=255)
                                                                                                                                              gTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTac         >  1:900813/1‑150 (MQ=255)
                                                                                                                                               ttCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTaca        >  1:2738903/1‑150 (MQ=255)
                                                                                                                                               ttCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTaca        <  2:4044449/150‑1 (MQ=255)
                                                                                                                                                tCTGTTGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTacac       <  2:3801602/150‑1 (MQ=255)
                                                                                                                                                 cTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACAcc      >  2:3701329/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCGCCCGGCCAGCGACATTATGCATGCGCATTTGCTGGTGGTGAGTAAGTAGAAATCGGCGGCCGCCTGCGGTTGATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACATGATCGTGCAAATTCAATAAATTGCAGCGTTCTGTAGGCTGGATAAGATGCGTCAGCATCGCATCCGGCAAAGGCAGATCTCAGCGATAGCGCCGGCTTAGTCAGATTTAATCTGCGCGCGTGGTGGATATTTTTTCAGGATCTCCATATACGCGTGCATTTCGGTCTGTAGCGGTACACCCATC  >  CP009273/3977501‑3977799

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: