Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,769,806 A→G 5.5% D261G (GAC→GGC)  ctpA → carboxy‑terminal‑processing protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,769,8060AG 5.5% 1052.6 / 30.2 325D261G (GAC→GGC) ctpAcarboxy‑terminal‑processing protease
Reads supporting (aligned to +/- strand):  ref base A (158/149);  new base G (10/8);  total (168/157)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 329 total aligned reads displayed.

ATCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGAAATTG  >  NC_002947/5769658‑5769955
                                                                                                                                                    |                                                                                                                                                     
aTCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGAc                                                                                                                                                      <  2:408759/150‑1 (MQ=255)
     ccAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACtt                                                                                                                                                 >  1:8610646/1‑150 (MQ=255)
      cAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTc                                                                                                                                                <  2:4931866/150‑1 (MQ=255)
       aGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTcc                                                                                                                                               >  2:1639097/1‑150 (MQ=255)
        gTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGTCGGCCACTTcca                                                                                                                                              <  1:1980181/150‑2 (MQ=255)
        gTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCt                                                                                                                                              <  1:6948058/150‑1 (MQ=255)
         ttCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTg                                                                                                                                             >  2:233229/1‑150 (MQ=255)
          tCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGTCGGCCACTTccact                                                                                                                                            <  1:1122856/150‑4 (MQ=255)
          tCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGa                                                                                                                                            <  1:3416847/150‑1 (MQ=255)
           ccAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGAc                                                                                                                                           >  2:8273648/1‑150 (MQ=255)
            cAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGAcc                                                                                                                                          >  1:5007268/1‑150 (MQ=255)
             aGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCa                                                                                                                                         >  1:1332735/1‑150 (MQ=255)
               gTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCaaa                                                                                                                                       <  2:6387018/150‑1 (MQ=255)
               gTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCaaa                                                                                                                                       >  1:7561383/1‑150 (MQ=255)
                     aCCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTg                                                                                                                                 <  1:2505554/150‑1 (MQ=255)
                      ccGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGa                                                                                                                                <  1:1220888/150‑1 (MQ=255)
                        ggTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATc                                                                                                                              >  1:989964/1‑150 (MQ=255)
                        ggTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATc                                                                                                                              <  1:184623/150‑1 (MQ=255)
                         gTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATcc                                                                                                                             >  2:5908465/1‑149 (MQ=255)
                           gACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTc                                                                                                                           <  2:2800852/150‑1 (MQ=255)
                           gACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTc                                                                                                                           <  1:307511/150‑1 (MQ=255)
                             cGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTa                                                                                                                         <  1:7438112/150‑1 (MQ=255)
                               aaGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTaca                                                                                                                       <  1:5966672/150‑1 (MQ=255)
                                aGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTacac                                                                                                                      <  1:831860/150‑1 (MQ=255)
                                    ggCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCaaa                                                                                                                  <  2:8687765/150‑1 (MQ=255)
                                      cAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAgg                                                                                                                >  2:8442080/1‑150 (MQ=255)
                                      cAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAgg                                                                                                                >  2:5765344/1‑150 (MQ=255)
                                             cTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAACGTCTGATCGTCTACACCAAAGGCCGCATc                                                                                                         <  2:4201953/150‑1 (MQ=255)
                                               ggCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGc                                                                                                       <  1:4055539/150‑1 (MQ=255)
                                               ggCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGc                                                                                                       <  1:6700480/150‑1 (MQ=255)
                                                  cAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCaa                                                                                                    >  1:4610301/1‑150 (MQ=255)
                                                  cAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCaa                                                                                                    >  1:4271952/1‑150 (MQ=255)
                                                   aaGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAAc                                                                                                   >  2:2858278/1‑150 (MQ=255)
                                                      cTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTcc                                                                                                >  1:7012172/1‑150 (MQ=255)
                                                        gcgcAAGGACAACGGCAAGAAGCTGCGCGGTGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGa                                                                                              <  1:247097/150‑1 (MQ=255)
                                                                 cAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGttt                                                                                     <  1:4351379/150‑1 (MQ=255)
                                                                 cAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGttt                                                                                     <  1:3711082/150‑1 (MQ=255)
                                                                   aCGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTtct                                                                                   >  1:8184112/1‑150 (MQ=255)
                                                                      gCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCgg                                                                                <  1:7788151/150‑1 (MQ=255)
                                                                       cAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGc                                                                               <  1:2657726/150‑1 (MQ=255)
                                                                        aagaagCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGcc                                                                              >  1:7854373/1‑150 (MQ=255)
                                                                          gaagCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGa                                                                            >  1:450817/1‑150 (MQ=255)
                                                                           aagCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGAc                                                                           <  2:7731649/150‑1 (MQ=255)
                                                                               tGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCgg                                                                       >  1:8123056/1‑150 (MQ=255)
                                                                               tGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCgg                                                                       >  2:6793547/1‑150 (MQ=255)
                                                                                gcgcGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGc                                                                      <  2:361182/150‑1 (MQ=255)
                                                                                gcgcGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGc                                                                      >  1:3452273/1‑150 (MQ=255)
                                                                                gcgcGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGc                                                                      <  2:6591246/150‑1 (MQ=255)
                                                                                      ggtggtGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGc                                                                <  2:5247773/150‑1 (MQ=255)
                                                                                      ggtggtGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGc                                                                >  2:1274999/1‑150 (MQ=255)
                                                                                        tggtgCTGGACCTGCGCAACAACCCGGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCa                                                              <  2:8355524/150‑1 (MQ=255)
                                                                                            gCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAgcgc                                                          <  2:2922976/150‑1 (MQ=255)
                                                                                              tGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCgg                                                        <  1:2381592/150‑1 (MQ=255)
                                                                                              tGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCgg                                                        <  2:3261311/150‑1 (MQ=255)
                                                                                               ggACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCggg                                                       >  1:7230245/1‑150 (MQ=255)
                                                                                               ggACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCggg                                                       <  2:3559725/150‑1 (MQ=255)
                                                                                                gACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGc                                                      <  2:3506067/150‑1 (MQ=255)
                                                                                                 aCCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCg                                                     >  2:7414361/1‑150 (MQ=255)
                                                                                                    tGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGc                                                  <  2:4271952/150‑1 (MQ=255)
                                                                                                    tGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGc                                                  <  1:1773430/150‑1 (MQ=255)
                                                                                                      cgcAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCg                                                >  1:1569454/1‑150 (MQ=255)
                                                                                                      cgcAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCg                                                >  1:8456393/1‑150 (MQ=255)
                                                                                                           caacCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggt                                           >  2:6068728/1‑150 (MQ=255)
                                                                                                            aacCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggtg                                          <  2:3217986/150‑1 (MQ=255)
                                                                                                            aacCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggtg                                          <  1:7980984/150‑1 (MQ=255)
                                                                                                            aacCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggtg                                          <  2:7213127/150‑1 (MQ=255)
                                                                                                            aacCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggtg                                          <  2:1521720/150‑1 (MQ=255)
                                                                                                            aacCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtggtg                                          <  2:513502/150‑1 (MQ=255)
                                                                                                              cccAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTggtggt                                        >  1:3268170/1‑150 (MQ=255)
                                                                                                              cccAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTggtggt                                        >  2:5352944/1‑150 (MQ=255)
                                                                                                               ccAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTa                                       <  2:8371486/150‑1 (MQ=255)
                                                                                                               ccAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTa                                       >  2:8150393/1‑150 (MQ=255)
                                                                                                                  ggcggcGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAACGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTg                                    <  2:5827906/150‑1 (MQ=255)
                                                                                                                  ggcggcGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTg                                    <  1:4226334/150‑1 (MQ=255)
                                                                                                                    cggcgTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGAt                                  >  2:5045670/1‑150 (MQ=255)
                                                                                                                     ggcgTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCAACGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGTCAGCGCGGGCGTGCCGCTGGTGGTACTGATc                                 >  1:2980007/1‑150 (MQ=255)
                                                                                                                        gTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAAc                              <  1:5138304/150‑1 (MQ=255)
                                                                                                                        gTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAAc                              <  1:5869461/150‑1 (MQ=255)
                                                                                                                           ctgcAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGt                           <  2:6989541/150‑1 (MQ=255)
                                                                                                                            tgcAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTg                          >  1:5452631/1‑150 (MQ=255)
                                                                                                                             gcAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTgg                         <  2:8007988/150‑1 (MQ=255)
                                                                                                                              cAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGc                        <  2:4564240/150‑1 (MQ=255)
                                                                                                                               aaTCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCtg                                             >  2:2562174/1‑128 (MQ=255)
                                                                                                                                 tCGGCAGTGGATGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGTCCGCATCGCCAACTCCGAGCTTCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCTCTGGTGGTACTGATCAACTGTTGCAGt                     <  2:5540829/150‑1 (MQ=255)
                                                                                                                                 tCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGt                     <  1:5668101/150‑1 (MQ=255)
                                                                                                                                    gCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgca                  <  2:4429105/150‑1 (MQ=255)
                                                                                                                                      tcaggAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatc                >  1:7204767/4‑150 (MQ=255)
                                                                                                                                      tcaggAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatc                >  2:76544/4‑150 (MQ=255)
                                                                                                                                        aggAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatcgg              >  1:7333274/2‑150 (MQ=255)
                                                                                                                                          gAAGTGGTCGGCCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatcggca            >  1:4792221/1‑150 (MQ=255)
                                                                                                                                           aaGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTgcatcggcat           >  2:5490541/1‑150 (MQ=255)
                                                                                                                                              tGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGcatcggcatcgg        >  1:4690328/1‑150 (MQ=255)
                                                                                                                                               ggCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGa       >  1:6502160/1‑150 (MQ=255)
                                                                                                                                                gCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGaa      >  2:190981/1‑150 (MQ=255)
                                                                                                                                                 ccGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGaaa     <  2:6844574/150‑1 (MQ=255)
                                                                                                                                                 ccGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGaaa     <  2:1800813/150‑1 (MQ=255)
                                                                                                                                                 ccGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGaaa     <  2:8456393/150‑1 (MQ=255)
                                                                                                                                                  cGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGAAAt    >  1:6220653/1‑150 (MQ=255)
                                                                                                                                                  cGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGAAAt    >  2:4623054/1‑150 (MQ=255)
                                                                                                                                                   gACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGAAAtt   >  1:7640347/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATCACCCAGTTCCAGGTCAAGACCGGTGACGAAGTGGGCAAGGCCCTGGCCAAGCTGCGCAAGGACAACGGCAAGAAGCTGCGCGGGGTGGTGCTGGACCTGCGCAACAACCCAGGCGGCGTGCTGCAATCGGCAGTGGAAGTGGCCGACCACTTCCTGACCAAAGGTCTGATCGTCTACACCAAAGGCCGCATCGCCAACTCCGAGCTGCGTTTCTCGGCCGACCCGGCCGATGCCAGCGCGGGCGTGCCGCTGGTGGTACTGATCAACGGTGGCAGTGCATCGGCATCGGAAATTG  >  NC_002947/5769658‑5769955

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: