Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,051,296 C→G 5.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 5.7% 274.0 / 6.6 87intergenic (+16/+27)oprB‑III/PP_3571carbohydrate‑selective porin/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base C (37/45);  new base G (2/3);  total (39/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01

GCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCACAC  >  NC_002947/4051153‑4051442
                                                                                                                                               |                                                                                                                                                  
gCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGc                                                                                                                                              <  2:1176762/150‑1 (MQ=255)
gCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGACCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACtgctgc                                                                                                                                              <  2:668646/150‑1 (MQ=255)
gCTACCAGGCCACGCCATGGCTGAACCTGCGACCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGTCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGACGCACTGcttc                                                                                                                                              <  2:1843872/150‑4 (MQ=255)
     cAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCAGTGCGGccccgg                                                                                                                                         <  2:2096648/150‑3 (MQ=255)
     cAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCAt                                                                                                                                         >  1:9207/1‑150 (MQ=255)
       ggCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATcc                                                                                                                                       <  2:814631/150‑1 (MQ=255)
               cATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGgcgc                                                                                                                               >  2:1378151/1‑150 (MQ=255)
                     tGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCt                                                                                                                         >  1:811229/1‑150 (MQ=255)
                        aCCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCaa                                                                                                                      <  1:900796/150‑1 (MQ=255)
                        aCCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTCCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTTAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCaa                                                                                                                      <  1:1574610/150‑1 (MQ=255)
                         ccTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAg                                                                                                                     >  2:1059865/1‑150 (MQ=255)
                          cTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCTGGGGCCGCAGTGCGGccccggcg                                                                                                                                       <  2:2207674/131‑5 (MQ=255)
                          cTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCTGGGGCCGCAGTGCGGccccggcg                                                                                                                                       >  1:2207674/1‑127 (MQ=255)
                           tGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAAcc                                                                                                                           >  1:335217/1‑142 (MQ=255)
                           tGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAAcc                                                                                                                           <  2:335217/142‑1 (MQ=255)
                                   gACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAg                                                                                                           <  1:2057638/150‑1 (MQ=255)
                                   gACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAg                                                                                                           <  2:2296075/150‑1 (MQ=255)
                                   gACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAg                                                                                                           <  2:2308744/150‑1 (MQ=255)
                                   gACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGGCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAg                                                                                                           >  2:56673/1‑150 (MQ=255)
                                    gatctcAATACATCATCGAGCCAGGGGCGTTCTCGGAAACATACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAgg                                                                                                          <  1:613572/145‑1 (MQ=255)
                                      gTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTa                                                                                                        <  2:9207/150‑1 (MQ=255)
                                          aaTACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGGGGCGCGAACCACCGGCTCAACCCTCAAGCAGGATTGAGGTaaatc                                                                                                    >  2:613572/1‑146 (MQ=255)
                                           aTACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTcca                                                                                                   <  2:148522/150‑1 (MQ=255)
                                            tACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTccac                                                                                                  <  1:1059865/150‑1 (MQ=255)
                                              caacaccGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACACCGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTccacca                                                                                                <  1:799433/144‑1 (MQ=255)
                                               atcatcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTcca                                                                                                   >  1:20513/1‑146 (MQ=255)
                                               atcatcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTcca                                                                                                   <  2:20513/146‑1 (MQ=255)
                                                tcatcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCTGGGGCCGCAGTGCGGccccggcg                                                                                                                                       <  2:739876/109‑5 (MQ=255)
                                                tcatcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCTGGGGCCGCAGTGCGGccccggcg                                                                                                                                       >  1:739876/1‑105 (MQ=255)
                                                tcatcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCaa                                                                                                                             >  1:349694/1‑119 (MQ=255)
                                                tcatcGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCaa                                                                                                                             <  2:349694/119‑1 (MQ=255)
                                                tcaacGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCGTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGc                                                                                              <  2:495476/150‑1 (MQ=255)
                                                     gAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTcag                                                                                         <  1:1725702/150‑1 (MQ=255)
                                                     gAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTcag                                                                                         <  2:455562/150‑1 (MQ=255)
                                                     gAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTcag                                                                                         <  1:720260/150‑1 (MQ=255)
                                                       gCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTcagca                                                                                       <  2:1990793/150‑1 (MQ=255)
                                                         cAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCCACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCagcagc                                                                                     >  1:2135923/1‑150 (MQ=255)
                                                           ggggCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCAt                                                                                   <  1:1378151/150‑1 (MQ=255)
                                                            gggCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATg                                                                                  >  2:1431638/1‑150 (MQ=255)
                                                             ggCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATgg                                                                                 >  2:615548/1‑150 (MQ=255)
                                                               cGTTCTCGGGCAAAGACATCGACAACGCACTCGCGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCa                                                                               <  2:182548/150‑1 (MQ=255)
                                                                   ctcGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCTGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATcac                                                                           <  1:231384/150‑1 (MQ=255)
                                                                      gggCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAg                                                                                                           <  2:1736165/115‑1 (MQ=255)
                                                                      gggCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAg                                                                                                           >  1:1736165/1‑115 (MQ=255)
                                                                         cAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAATCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCcgc                                                                     <  1:615548/150‑1 (MQ=255)
                                                                         cAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCcgc                                                                     >  2:1921129/1‑150 (MQ=255)
                                                                                aTCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGt                                                                                                      >  2:1966494/1‑110 (MQ=255)
                                                                                aTCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGt                                                                                                      <  1:1966494/110‑1 (MQ=255)
                                                                                aTCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTg                                                              >  1:960073/1‑150 (MQ=255)
                                                                                     cAACGCACGCGTGGTCGGCCTGCAGGTCACAGCGTCCTTCTGTTCCGCTGGGTTCGCGCTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGa                                                         <  2:2135923/150‑1 (MQ=255)
                                                                                      aaCGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGac                                                        >  2:318286/1‑150 (MQ=255)
                                                                                       aCGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGaca                                                       <  1:874949/150‑1 (MQ=255)
                                                                                          cACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTg                                                    <  1:842016/150‑1 (MQ=255)
                                                                                               gTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGCGCCGCACTGCGGCCCCATCCGCCGGCCCAACCCTCAAGCCGGTTTGACGTAGAACACCAGGCCAAGCAGCAGAACATCGCATCCGCGCAACTGCTCGACACTGACaaa                                               >  1:730873/1‑150 (MQ=255)
                                                                                                tGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAAACGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAAc                                              <  2:895756/150‑1 (MQ=255)
                                                                                                   tCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCa                                                                               <  2:2189847/114‑1 (MQ=255)
                                                                                                   tCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCa                                                                               >  1:2189847/1‑114 (MQ=255)
                                                                                                      gCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATcc                                        >  2:1856315/1‑150 (MQ=255)
                                                                                                              gTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTtgtg                                <  1:1431638/150‑1 (MQ=255)
                                                                                                               tCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTgg                               >  1:1268762/1‑150 (MQ=255)
                                                                                                                cAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTtgtg                                >  2:438187/1‑148 (MQ=255)
                                                                                                                cAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTtgtg                                <  1:438187/148‑1 (MQ=255)
                                                                                                                 aaGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAg                                                                                               >  2:1212978/1‑84 (MQ=255)
                                                                                                                 aaGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGcc                             <  2:1428379/150‑1 (MQ=255)
                                                                                                                 aaGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGACCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAg                                                                                               <  1:1212978/84‑1 (MQ=255)
                                                                                                                   ggCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATc                                                                             <  2:185719/100‑1 (MQ=255)
                                                                                                                   ggCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATc                                                                             >  1:185719/1‑100 (MQ=255)
                                                                                                                    gCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTg                          >  1:1889484/1‑150 (MQ=255)
                                                                                                                      gACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATg                                                                                  <  2:2235319/92‑1 (MQ=255)
                                                                                                                      gACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATg                                                                                  >  1:2235319/1‑92 (MQ=255)
                                                                                                                      gACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTc                                                           <  1:719583/115‑1 (MQ=255)
                                                                                                                      gACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTc                                                           >  2:719583/1‑115 (MQ=255)
                                                                                                                          ttCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCAt                    <  2:811229/150‑1 (MQ=255)
                                                                                                                            cTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCAtgc                  >  1:838554/1‑150 (MQ=255)
                                                                                                                               aaCCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGa                                                         <  2:930656/108‑1 (MQ=255)
                                                                                                                               aaCCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGa                                                         >  1:930656/1‑108 (MQ=255)
                                                                                                                                aCCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGcc              >  2:1853283/1‑150 (MQ=255)
                                                                                                                                 ccgccgGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCa             <  1:318286/150‑1 (MQ=255)
                                                                                                                                  cgccgGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAg            >  1:2019192/1‑150 (MQ=255)
                                                                                                                                  cgccgGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAg            >  1:1123118/1‑150 (MQ=255)
                                                                                                                                   gacgGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAgg           >  1:2215020/3‑150 (MQ=255)
                                                                                                                                      ggggCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTTAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGcc        >  1:2300055/1‑150 (MQ=255)
                                                                                                                                         gCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGca     <  1:1853283/150‑1 (MQ=255)
                                                                                                                                           cGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATca                                                                            >  2:2334121/1‑77 (MQ=255)
                                                                                                                                           cGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATca                                                                            <  1:2334121/77‑1 (MQ=255)
                                                                                                                                           cGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGcaca   >  1:968775/1‑150 (MQ=255)
                                                                                                                                            gcccTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCacac  <  1:56673/147‑1 (MQ=255)
                                                                                                                                            gCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCacac  >  1:1940247/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GCTACCAGGCCACGCCATGGCTGAACCTGCGCCCGGACGTGCAATACATCATCGAGCCAGGGGCGTTCTCGGGCAAAGACATCGACAACGCACTCGTGGTCGGCCTGCAGGTCAAGGCGACCTTCTGAACCGCCGGGGCCGCACTGCGGCCCCATCCGCCGGCGCAACCCTCAAGCCTGTTTGAGGTAGTCCACCAGCCTCAGCAGCATGGCATCACAGCCGCGCAACTGCTCGACACTGACAAACTCATCCGGCTTGTGGCCCTGGTCCATGCTGCCAGGGCCGCACAC  >  NC_002947/4051153‑4051442

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