Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,051,296 C→G 7.7% 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,051,2960CG 7.7% 289.2 / 12.1 91intergenic (+16/+27)oprB‑III/PP_3571carbohydrate‑selective porin/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (43/41);  new base G (4/3);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

TACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCACACC  >  NC_002947/4051149‑4051443
                                                                                                                                                   |                                                                                                                                                   
tACGGCTACCAGGCCACGCCATGGCTGAACCTGCACCCGGACGTGCAATACATCATCTAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTg                                                                                                                                                   <  1:274669/150‑1 (MQ=255)
    gCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGTAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGc                                                                                                                                               <  2:2184874/150‑1 (MQ=255)
    gCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGc                                                                                                                                               >  2:1929446/1‑150 (MQ=255)
      tACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGccc                                                                                                                                             <  2:915257/150‑1 (MQ=255)
           ggCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCAGTGCGGccccggcg                                                                                                                                        <  2:847957/150‑5 (MQ=255)
                 gCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGgc                                                                                                                                  >  2:1762816/1‑150 (MQ=255)
                   cATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCCCCGGGGCCGCACTGCGGCCCCATCCGCCGgcgc                                                                                                                                >  2:207004/1‑150 (MQ=255)
                      ggCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAAc                                                                                                                             <  2:2009647/150‑1 (MQ=255)
                         tGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCt                                                                                                                          >  1:600185/1‑150 (MQ=255)
                             ccTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAg                                                                                                                      >  2:1665040/1‑150 (MQ=255)
                              cTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGc                                                                                                                     >  1:199134/1‑150 (MQ=255)
                                 cgcCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCCCCGGCGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTg                                                                                                                  <  1:2359759/150‑1 (MQ=255)
                                      ggACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGa                                                                                                             <  2:1400822/150‑1 (MQ=255)
                                       gACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAg                                                                                                            <  2:1496424/150‑1 (MQ=255)
                                       gACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAg                                                                                                            <  1:2254835/150‑1 (MQ=255)
                                          gTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTa                                                                                                         <  2:1182631/150‑1 (MQ=255)
                                            gCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGcc                                                                                                                                              >  1:915759/1‑111 (MQ=255)
                                            gCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGcc                                                                                                                                              <  2:915759/111‑1 (MQ=255)
                                            gCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGt                                                                                                       >  1:1063880/1‑150 (MQ=255)
                                              aaTACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTcc                                                                                                     >  1:1297346/1‑150 (MQ=255)
                                                  catcatCGAGCCAGGGGCGTTCTAACCGATAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTccacca                                                                                                 >  1:1269775/1‑150 (MQ=255)
                                                   atcatcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAg                                                                                                >  1:462641/1‑150 (MQ=255)
                                                   atcatcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAg                                                                                                >  1:1235975/1‑150 (MQ=255)
                                                    tcatcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCTACCCTCAAGCCTGTTTGCGGTAGTCCACCAGc                                                                                               <  1:1667016/150‑1 (MQ=255)
                                                     catcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGcc                                                                                              >  1:410937/1‑150 (MQ=255)
                                                      atcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCt                                                                                             >  2:1754873/1‑150 (MQ=255)
                                                      atcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCt                                                                                             >  1:1875875/1‑150 (MQ=255)
                                                      atcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCt                                                                                             >  1:249876/1‑150 (MQ=255)
                                                          aGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCt                                                                                                                          >  2:1438565/1‑117 (MQ=255)
                                                          aGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCt                                                                                                                          <  1:1438565/117‑1 (MQ=255)
                                                             cAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCagcagc                                                                                      >  2:493649/1‑150 (MQ=255)
                                                              aGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGc                                                                                                                                               <  2:378738/92‑1 (MQ=255)
                                                              aGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGc                                                                                                                                               >  1:378738/1‑92 (MQ=255)
                                                               ggggCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCAt                                                                                    <  1:1929446/150‑1 (MQ=255)
                                                               ggggCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCAt                                                                                    <  1:2331524/150‑1 (MQ=255)
                                                               ggggCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCAt                                                                                    <  2:68649/150‑1 (MQ=255)
                                                                gggCGTTCTCGGTCAAAGACATCGCCAACGCACGCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATg                                                                                   <  2:1508203/150‑1 (MQ=255)
                                                                gggCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATg                                                                                   <  1:1665040/150‑1 (MQ=255)
                                                                    gTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCACCCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTcagc                                                                                         <  1:1949164/140‑1 (MQ=255)
                                                                    gTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTcagc                                                                                         >  2:1949164/1‑140 (MQ=255)
                                                                      tctcGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATca                                                                             <  2:199134/150‑1 (MQ=255)
                                                                              aaaGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCcgcg                                                                     >  2:883297/1‑150 (MQ=255)
                                                                              aaaGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCGGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCcgcg                                                                     >  2:1643702/1‑150 (MQ=255)
                                                                                    aTCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAgg                                                                                                           <  2:854114/106‑1 (MQ=255)
                                                                                    aTCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAgg                                                                                                           >  1:854114/1‑106 (MQ=255)
                                                                                    aTCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTg                                                               >  2:1659954/1‑150 (MQ=255)
                                                                                      cGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGccc                                                                                                                                             >  1:230697/1‑70 (MQ=255)
                                                                                      cGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGccc                                                                                                                                             <  2:230697/70‑1 (MQ=255)
                                                                                          aaCGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGac                                                         >  2:1909777/1‑150 (MQ=255)
                                                                                           aCGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGaca                                                        <  1:1754873/150‑1 (MQ=255)
                                                                                           aCGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGaca                                                        <  1:1800448/150‑1 (MQ=255)
                                                                                                   gTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACaaa                                                <  2:1297346/150‑1 (MQ=255)
                                                                                                   gTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGCCCGCCCTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACaaa                                                <  2:2077566/150‑1 (MQ=255)
                                                                                                           ccTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCg                                        >  1:600111/1‑150 (MQ=255)
                                                                                                            cTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGttt                                                                                                               <  2:2137005/78‑1 (MQ=255)
                                                                                                            cTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGttt                                                                                                               >  1:2137005/1‑78 (MQ=255)
                                                                                                            cTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCgg                                       >  2:1220105/1‑150 (MQ=255)
                                                                                                                  gTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGGcg                                                                       <  1:1667524/112‑1 (MQ=255)
                                                                                                                  gTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCAGCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCcg                                                                       >  2:1667524/1‑112 (MQ=255)
                                                                                                                  ccttaGGCGACCTCCTGAACCGCCGGTCCCGTACTGCGGACCCATCCGCCGGCGCAACCCTCAAGCCTTTTTTAGGGAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTtgtg                                 <  2:419500/146‑1 (MQ=255)
                                                                                                                   tCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTgg                                >  2:782701/1‑150 (MQ=255)
                                                                                                                     aaGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCtt                                    >  1:2250586/1‑144 (MQ=255)
                                                                                                                     aaGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCtt                                    <  2:2250586/144‑1 (MQ=255)
                                                                                                                         cGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTgg                          <  1:1220105/150‑1 (MQ=255)
                                                                                                                          gACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGt                         <  1:883297/150‑1 (MQ=255)
                                                                                                                            ccTTCTGAACCGCCGGGGCCGCACTGCGGCGCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTcc                       >  1:2064462/1‑150 (MQ=255)
                                                                                                                            ccTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTcc                       <  2:410221/150‑1 (MQ=255)
                                                                                                                             cTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCa                      >  1:2188636/1‑150 (MQ=255)
                                                                                                                             cTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGAGAATGATTGATGTTCACTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCa                      <  2:1269775/150‑1 (MQ=255)
                                                                                                                              ttCTGAACCGCCGGGTCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTcag                                                                                          <  2:1736043/81‑1 (MQ=255)
                                                                                                                              ttCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTcag                                                                                          >  1:1736043/1‑81 (MQ=255)
                                                                                                                              ttCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCAt                     <  2:1063880/150‑1 (MQ=255)
                                                                                                                                     ccgccgGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCa              <  2:2188636/150‑1 (MQ=255)
                                                                                                                                     ccgccgGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCa              >  2:811284/1‑150 (MQ=255)
                                                                                                                                      ggatggggCCGCAGTGCGGCCCCGTCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATcc                                         >  2:324938/5‑122 (MQ=255)
                                                                                                                                      ggatggggCCGCAGTGCGGCCCCGTCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATcc                                         <  1:324938/118‑1 (MQ=255)
                                                                                                                                      ggatggggCCGCAGTGCGGCCCCGTCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAg             >  1:1743465/5‑150 (MQ=255)
                                                                                                                                      ggacgGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAg             >  1:685613/4‑150 (MQ=255)
                                                                                                                                      cgccgGGGCCGCACTGCGGCGCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAg             <  2:2064462/150‑1 (MQ=255)
                                                                                                                                       gacgGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTa                                                                                                         >  2:399393/3‑57 (MQ=255)
                                                                                                                                       gacgGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTa                                                                                                         <  1:399393/55‑1 (MQ=255)
                                                                                                                                        atggggCCGCAGTGCGGCCCCGTCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAggg           >  2:2374563/3‑150 (MQ=255)
                                                                                                                                            ggCCGCAGTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATgg                                                                                  >  1:762547/1‑75 (MQ=255)
                                                                                                                                            ggCCGCAGTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATgg                                                                                  <  2:762547/75‑1 (MQ=255)
                                                                                                                                               cGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTtgt                                  >  1:2389042/1‑120 (MQ=255)
                                                                                                                                               cGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTtgt                                  <  2:2389042/120‑1 (MQ=255)
                                                                                                                                               cGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACGTGCAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGcaca    >  1:2330653/1‑150 (MQ=255)
                                                                                                                                                gCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGTCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCacac   >  1:2365644/1‑150 (MQ=255)
                                                                                                                                                 cACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCgcgc                                                                    >  2:627801/1‑84 (MQ=255)
                                                                                                                                                 cACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCgcgc                                                                    <  1:627801/84‑1 (MQ=255)
                                                                                                                                                 cACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCAcacc  >  1:1165696/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TACGGCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCACACC  >  NC_002947/4051149‑4051443

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: