Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,056,245 C→G 5.9% intergenic (+16/+27) oprB‑III → / ← PP_3571 carbohydrate‑selective porin/acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,056,2450CG 5.9% 1254.2 / 48.0 372intergenic (+16/+27)oprB‑III/PP_3571carbohydrate‑selective porin/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (191/159);  new base G (11/11);  total (202/170)
Fisher's exact test for biased strand distribution p-value = 8.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 373 total aligned reads displayed.

CTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCACACCAC  >  NC_002947/4056097‑4056394
                                                                                                                                                    |                                                                                                                                                     
cTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACt                                                                                                                                                      >  2:2375937/1‑150 (MQ=255)
 tACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTg                                                                                                                                                     >  1:1724329/1‑150 (MQ=255)
 tACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTg                                                                                                                                                     >  1:3663637/1‑150 (MQ=255)
    ggCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCgg                                                                                                                                                  >  2:776069/1‑150 (MQ=255)
      cTACCAGGCCACGCCGTGGCTGAACCTGCGCCCGGACGTGCAATACATCCTCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGcc                                                                                                                                                <  2:2184679/150‑1 (MQ=255)
      cTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGcc                                                                                                                                                <  2:8616273/150‑1 (MQ=255)
      cTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGcc                                                                                                                                                <  1:5465774/150‑1 (MQ=255)
       tACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGccc                                                                                                                                               >  2:4831539/1‑150 (MQ=255)
       tACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGccc                                                                                                                                               >  2:3995299/1‑150 (MQ=255)
        aCCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGcccc                                                                                                                                              >  2:8085867/1‑150 (MQ=255)
          cAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCAGTGCGGccccgg                                                                                                                                            <  2:7123275/150‑3 (MQ=255)
          cAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCAGTGCGGccccgg                                                                                                                                            <  2:3535721/150‑3 (MQ=255)
          cAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCAGTGCGGccccgg                                                                                                                                            <  1:4315368/150‑3 (MQ=255)
          cAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCAt                                                                                                                                            >  1:2420760/1‑150 (MQ=255)
           aGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATc                                                                                                                                           <  2:3793094/150‑1 (MQ=255)
            ggCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCTGGGGCCGCAGTGCGGccccggcg                                                                                                                                          <  1:7974498/150‑5 (MQ=255)
            ggCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCagcg                                                                                                                                          <  1:1870637/150‑4 (MQ=255)
              ccACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATccgc                                                                                                                                        <  2:8665917/150‑1 (MQ=255)
                aCGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATccgccg                                                                                                                                      >  2:5489434/1‑150 (MQ=255)
                  gCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGgc                                                                                                                                    >  2:7942889/1‑150 (MQ=255)
                  gCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCACCCGCCGgc                                                                                                                                    <  2:6189980/150‑1 (MQ=255)
                        gCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAAcc                                                                                                                              >  2:3716728/1‑150 (MQ=255)
                        gCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAAcc                                                                                                                              >  1:7590375/1‑150 (MQ=255)
                             aCCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCaa                                                                                                                         <  1:5983821/150‑1 (MQ=255)
                                tGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGcc                                                                                                                      >  2:6843385/1‑150 (MQ=255)
                                tGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGcc                                                                                                                      >  2:1665710/1‑150 (MQ=255)
                                  cgcCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTg                                                                                                                    <  2:4103271/150‑1 (MQ=255)
                                    cccGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGtt                                                                                                                  <  1:8244519/150‑1 (MQ=255)
                                       ggACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGa                                                                                                               <  2:732194/150‑1 (MQ=255)
                                         aCGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAgg                                                                                                             <  1:7039629/150‑1 (MQ=255)
                                         aCGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAgg                                                                                                             <  2:4339683/150‑1 (MQ=255)
                                          cGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGt                                                                                                            <  2:8453769/150‑1 (MQ=255)
                                           gTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTa                                                                                                           <  2:594244/150‑1 (MQ=255)
                                              cAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTc                                                                                                        >  2:3398120/1‑150 (MQ=255)
                                               aaTACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTcc                                                                                                       >  1:7321950/1‑150 (MQ=255)
                                               aaTACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGATTGAGGTAGTcc                                                                                                       <  2:5543443/150‑1 (MQ=255)
                                                aTACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTcca                                                                                                      >  1:2311537/1‑150 (MQ=255)
                                                aTACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTcca                                                                                                      >  2:8243975/1‑150 (MQ=255)
                                                aTACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTcca                                                                                                      >  1:7893088/1‑150 (MQ=255)
                                                 tACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTccac                                                                                                     <  1:4478652/150‑1 (MQ=255)
                                                   catcatCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTccacca                                                                                                   <  2:6084582/150‑1 (MQ=255)
                                                    atcatcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAg                                                                                                  >  1:4191557/1‑150 (MQ=255)
                                                    atcatcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAg                                                                                                  >  2:168686/1‑150 (MQ=255)
                                                    atcatcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAg                                                                                                  >  2:2938559/1‑150 (MQ=255)
                                                       atcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCAGTGCGGccccgg                                                                                                                                            <  1:8737128/105‑3 (MQ=255)
                                                       atcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCt                                                                                               >  2:4904738/1‑150 (MQ=255)
                                                       atcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCt                                                                                               >  1:1971296/1‑150 (MQ=255)
                                                        tcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTc                                                                                              <  1:2893014/150‑1 (MQ=255)
                                                            gCCAGGGGCGTTCTCGGGCAAAGACGTCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTcagca                                                                                          <  2:3882548/150‑1 (MQ=255)
                                                               aGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCCCCGGCGCAACCCACAAGCCGGGTTGAGGTAGTCCACCAGCCACAgcagca                                                                                       >  1:2224494/1‑150 (MQ=255)
                                                                  ggCGTCCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATgg                                                                                    >  2:970611/1‑150 (MQ=255)
                                                                    cGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCa                                                                                  <  1:6072039/150‑1 (MQ=255)
                                                                     gTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCAt                                                                                 <  1:4243665/150‑1 (MQ=255)
                                                                       tctcGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATca                                                                               >  2:6709567/1‑150 (MQ=255)
                                                                        ctcGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATcac                                                                              >  2:1967630/1‑150 (MQ=255)
                                                                         tcGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCGGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATcaca                                                                             >  2:450127/1‑150 (MQ=255)
                                                                          cGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAg                                                                            <  2:229873/150‑1 (MQ=255)
                                                                           gggCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGc                                                                           <  2:2676113/150‑1 (MQ=255)
                                                                               aaaGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCcgcg                                                                       >  1:6907349/1‑150 (MQ=255)
                                                                                 aGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCGCAAGCCTGTTTGAGGTAGTCCACCCGCCTCAGCAGCATGGCCTCACAGCCGCGCa                                                                     >  1:1665892/1‑150 (MQ=255)
                                                                                     aTCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTg                                                                 >  2:6478843/1‑150 (MQ=255)
                                                                                      tCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGc                                                                <  2:1448376/150‑1 (MQ=255)
                                                                                       cGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACAGCt                                                               >  1:5322330/1‑150 (MQ=255)
                                                                                        gACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTc                                                              >  2:7165001/1‑150 (MQ=255)
                                                                                        gACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTc                                                              <  2:2551164/150‑1 (MQ=255)
                                                                                          cAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGa                                                            <  1:2636247/150‑1 (MQ=255)
                                                                                           aaCGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGac                                                           <  1:6444467/150‑1 (MQ=255)
                                                                                           aaCGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGac                                                           >  2:5438725/1‑150 (MQ=255)
                                                                                            aCGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGaca                                                          <  1:1607527/150‑1 (MQ=255)
                                                                                            aCGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGaca                                                          <  1:6372691/150‑1 (MQ=255)
                                                                                             cGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGacac                                                         >  2:1309959/1‑150 (MQ=255)
                                                                                              gCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACt                                                        <  1:1327671/150‑1 (MQ=255)
                                                                                              gCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACt                                                        >  2:3591189/1‑150 (MQ=255)
                                                                                              gCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCCCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATTACAGCCGCGCAACTGCTCGACACt                                                        >  2:4023561/1‑150 (MQ=255)
                                                                                                aCTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGa                                                      <  1:6220053/150‑1 (MQ=255)
                                                                                                aCTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGa                                                      <  1:5002596/150‑1 (MQ=255)
                                                                                                  tCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACa                                                    >  2:2763100/1‑150 (MQ=255)
                                                                                                  tCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACa                                                    >  2:6275323/1‑150 (MQ=255)
                                                                                                          ggCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCACGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATc                                            <  1:5109271/150‑1 (MQ=255)
                                                                                                                cAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCtt                                      <  2:5135011/150‑1 (MQ=255)
                                                                                                                 aGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTtg                                     <  1:3541150/150‑1 (MQ=255)
                                                                                                                      aaGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGcc                                >  2:6959794/1‑150 (MQ=255)
                                                                                                                       aGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGccc                               <  1:5689027/150‑1 (MQ=255)
                                                                                                                         gCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTg                             <  2:1586236/150‑1 (MQ=255)
                                                                                                                          cGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTgg                            <  2:6370688/150‑1 (MQ=255)
                                                                                                                           gACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGt                           <  1:5179321/150‑1 (MQ=255)
                                                                                                                             ccTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTcc                         <  2:1189010/150‑1 (MQ=255)
                                                                                                                              cTTCTGAACCGCCGGGGCCGCACTGCGGCCCCCTCCGCCGGCGCCACCCTCCAGCCTGGTTGCGGTAGTCCCCCAGCCTCAGCAGCATGGCCTCACAGCCGCGCAACTGCTCGACACTGACAAACTCCTCCGGCTTGTGGCGCTGGTCCa                        >  2:4231392/1‑150 (MQ=255)
                                                                                                                               ttCTGAACCGCCGGGGCCGCACGGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCAt                       >  1:4005580/1‑150 (MQ=255)
                                                                                                                                       cgccgGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAg               >  2:4723329/1‑150 (MQ=255)
                                                                                                                                         ccgGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAggg             >  2:3110194/1‑150 (MQ=255)
                                                                                                                                         atggggCCGCAGTGCGGCCCCGTCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAggg             >  1:113045/3‑150 (MQ=255)
                                                                                                                                         atggggCCGCAGTGCGGCCCCGTCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAggg             >  1:7204167/3‑150 (MQ=255)
                                                                                                                                           ggggCCGCAGTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGcc           >  1:4762230/1‑150 (MQ=255)
                                                                                                                                               ccGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGcac       >  1:4374031/1‑150 (MQ=255)
                                                                                                                                                 gCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCacac     >  2:85494/1‑150 (MQ=255)
                                                                                                                                                 gCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCacac     >  1:90810/1‑150 (MQ=255)
                                                                                                                                                   aCTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCa                <  1:7712064/137‑1 (MQ=255)
                                                                                                                                                   aCTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCAcacca   >  2:2100507/1‑150 (MQ=255)
                                                                                                                                                    cTGCGGTCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCAcaccac  <  1:3380821/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCACACCAC  >  NC_002947/4056097‑4056394

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: