Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,672,556 C→G 6.6% A252G (GCC→GGC)  PP_4979 → substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,672,5560CG 6.6% 273.3 / 7.4 86A252G (GCC→GGC) PP_4979substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (41/38);  new base G (3/3);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCAGTACCTGGGCCTGGTGCGCAAGCCTGGGCGCGAGGCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGATGTGC  >  NC_002947/5672411‑5672705
                                                                                                                                                 |                                                                                                                                                     
gcgcAGTACCTGGGCCTGGTGCGCAAGCCTGGGCGCGAGGCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGtt                                                                                                                                                   <  1:1078413/150‑1 (MQ=255)
gcgcAGTACCTGGGCCTGGTGCGCAAGCCTGGGCGCGAGGCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGtt                                                                                                                                                   <  1:352731/150‑1 (MQ=255)
         cTGGGCCTGGTGCGCAAGCCTGGGCGCGAGGCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGGCGTTGagaat                                                                                                                                             <  1:2035723/147‑4 (MQ=255)
         cTGGGCCTGGTGCGCAAGCCTGGGCGCGAGGCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGGCGTTGagaat                                                                                                                                             >  2:2035723/1‑144 (MQ=255)
         cTGGGCCTGGTGCGCAAGCCTGGGCGCGAGGCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGa                                                                                                                                                 <  2:1182229/143‑1 (MQ=255)
         cTGGGCCTGGTGCGCAAGCCTGGGCGCGAGGCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGa                                                                                                                                                 >  1:1182229/1‑143 (MQ=255)
                      gcAAGCCTGGGCGCGAGTCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAgcagca                                                                                                                                <  2:356706/147‑1 (MQ=255)
                      gcAAGCCTGGGCGCGAGTCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAgcagca                                                                                                                                >  1:356706/1‑147 (MQ=255)
                        aaGCCTGGGCGCGAGGCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAgcg                                                                                                                           >  2:2034376/1‑150 (MQ=255)
                          gCCTGGGCGCGAGGCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAgcagca                                                                                                                                           >  2:2306259/1‑132 (MQ=255)
                          gCCTGGGCGCGAGGCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAgcagca                                                                                                                                           <  1:2306259/132‑1 (MQ=255)
                           ccTGGGCGCGAGGCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGa                                                                                                                        <  1:2034376/150‑1 (MQ=255)
                                  gcgAGGCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGctct                                                                                                                 <  1:432311/150‑1 (MQ=255)
                                     aGGCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGTAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGtt                                                                                                              <  2:1048359/150‑1 (MQ=255)
                                        cACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTtata                                                                                                           <  1:1961613/150‑1 (MQ=255)
                                           tGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAagc                                                                                                                                    <  2:251878/122‑1 (MQ=255)
                                           tGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAagc                                                                                                                                    >  1:251878/1‑122 (MQ=255)
                                           tGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATCCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTc                                                                                                        >  2:640063/1‑150 (MQ=255)
                                            ggCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCg                                                                                                       <  2:2348119/150‑1 (MQ=255)
                                            ggCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCg                                                                                                       <  1:972666/150‑1 (MQ=255)
                                             gccagtgCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCgg                                                                                                      <  1:742075/144‑1 (MQ=255)
                                             gCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTAAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCTCGATAGCTCTGTTATACTCgg                                                                                                      <  1:2031009/150‑1 (MQ=255)
                                              ccGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCggg                                                                                                     >  2:551876/1‑150 (MQ=255)
                                              ccGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCggg                                                                                                     >  1:1025511/1‑150 (MQ=255)
                                                     cTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGGCGTTGagaat                                                                                                                                             >  1:2129900/1‑100 (MQ=255)
                                                     cTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGGCGTTGagaat                                                                                                                                             <  2:2129900/103‑4 (MQ=255)
                                                     cTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTcgc                                                                                              >  2:485975/1‑148 (MQ=255)
                                                      tttGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTcccg                                                                                             >  2:739/1‑150 (MQ=255)
                                                           cGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGACGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCgg                                                                                        >  2:1027086/1‑150 (MQ=255)
                                                            gAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGc                                                                                       >  2:1335836/1‑150 (MQ=255)
                                                            gAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGc                                                                                       <  1:1339156/150‑1 (MQ=255)
                                                                  cTGCATCGGTTCAAGCAGGAGCCGGCCCATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGc                                                                                 <  1:994326/150‑1 (MQ=255)
                                                                           ttCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATAc                                                                                                          >  1:261119/1‑116 (MQ=255)
                                                                           ttCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATAc                                                                                                          <  2:261119/116‑1 (MQ=255)
                                                                              aaGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCgg                                                                     >  2:416387/1‑150 (MQ=255)
                                                                               aGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGc                                                                    >  2:1802366/1‑150 (MQ=255)
                                                                                  aGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTc                                                                 <  1:243762/150‑1 (MQ=255)
                                                                                     aGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCa                                                              >  1:2385156/1‑150 (MQ=255)
                                                                                          gCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAgg                                                         >  2:1890464/1‑150 (MQ=255)
                                                                                          gCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAgg                                                         >  2:1913450/1‑150 (MQ=255)
                                                                                          gCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAgg                                                         >  1:1197060/1‑150 (MQ=255)
                                                                                            cTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTtata                                                                                                           >  2:286500/1‑98 (MQ=255)
                                                                                            cTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTtata                                                                                                           <  1:286500/98‑1 (MQ=255)
                                                                                              aTGCAGCGATCAGCAACCGTTATACCACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATc                                                     <  2:2109194/150‑1 (MQ=255)
                                                                                                gCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAGCGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGTTCGGCCTCGCAACAGGCATccc                                                   >  1:475781/1‑150 (MQ=255)
                                                                                                gCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATccc                                                   >  1:2285213/1‑150 (MQ=255)
                                                                                                   gCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAgcagca                                                                                                                                           <  2:670146/59‑1 (MQ=255)
                                                                                                   gCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAgcagca                                                                                                                                           >  1:670146/1‑59 (MQ=255)
                                                                                                        cAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTa                                           <  1:1802366/150‑1 (MQ=255)
                                                                                                         aGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTAc                                          <  1:1782336/150‑1 (MQ=255)
                                                                                                          gCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATa                                                                                                                      <  2:537803/73‑1 (MQ=255)
                                                                                                          gCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATa                                                                                                                      >  1:537803/1‑73 (MQ=255)
                                                                                                          gCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTAcc                                         >  2:1777919/1‑150 (MQ=255)
                                                                                                          gCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAACTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGa                                                                            >  1:869186/1‑115 (MQ=255)
                                                                                                          gCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAACTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGa                                                                            <  2:869186/115‑1 (MQ=255)
                                                                                                           cAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGa                                                                                                                                                 >  2:2325219/1‑45 (MQ=255)
                                                                                                           cAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGa                                                                                                                                                 <  1:2325219/45‑1 (MQ=255)
                                                                                                               cGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGc                                    >  1:621164/1‑150 (MQ=255)
                                                                                                               cGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAACATGCTTATACCCGGTGGTTCGCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGc                                    >  1:1977192/1‑150 (MQ=255)
                                                                                                                gTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGcc                                   >  2:1931397/1‑150 (MQ=255)
                                                                                                                gTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGcc                                   <  2:1025511/150‑1 (MQ=255)
                                                                                                                 ttATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGccc                                  >  1:2309655/1‑150 (MQ=255)
                                                                                                                       aaCCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCt                                                                  <  2:311484/112‑1 (MQ=255)
                                                                                                                       aaCCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCt                                                                  >  1:311484/1‑112 (MQ=255)
                                                                                                                       aaCCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTcgc                            >  2:200647/1‑150 (MQ=255)
                                                                                                                       aaCCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTcgc                            <  2:2028653/150‑1 (MQ=255)
                                                                                                                       aaCCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTcgc                            >  2:220985/1‑150 (MQ=255)
                                                                                                                            acacTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTcc                        >  2:703112/1‑149 (MQ=255)
                                                                                                                            acacTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTcc                        <  1:703112/149‑1 (MQ=255)
                                                                                                                              acTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCCCGCGCGGACGCTCGGCCTCGCAACAGGCATCCCGATc                                               <  1:1171326/124‑1 (MQ=255)
                                                                                                                              acTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATc                                               >  2:1171326/1‑124 (MQ=255)
                                                                                                                                tGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTc                                                                       >  1:1927807/1‑98 (MQ=255)
                                                                                                                                tGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTc                                                                       <  2:1927807/98‑1 (MQ=255)
                                                                                                                                tGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCt                   >  2:992806/1‑150 (MQ=255)
                                                                                                                                   aaaCATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGGCGCCCCTTCGCGTCCAGCCTCCg                <  1:1335836/150‑1 (MQ=255)
                                                                                                                                    aaCATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGt               <  2:621164/150‑1 (MQ=255)
                                                                                                                                     tgctgctcAACGGCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGtt              >  1:485596/6‑150 (MQ=255)
                                                                                                                                     aCATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGc                                                                                 <  1:723884/83‑1 (MQ=255)
                                                                                                                                     aCATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGc                                                                                 >  2:723884/1‑83 (MQ=255)
                                                                                                                                         tctcAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTccc          <  1:1777919/150‑1 (MQ=255)
                                                                                                                                           tcAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCgg        >  1:1517394/1‑150 (MQ=255)
                                                                                                                                            cAACGCCGTTGAGCAGCAGGAAAGCACGGGATCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGa       <  1:739/150‑1 (MQ=255)
                                                                                                                                             aaCGGCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGAt      <  2:485596/150‑1 (MQ=255)
                                                                                                                                             aaCGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGAt                                                >  1:1671791/1‑108 (MQ=255)
                                                                                                                                             aaCGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGAt                                                <  2:1671791/108‑1 (MQ=255)
                                                                                                                                               cGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGAtgt    >  1:1189993/1‑150 (MQ=255)
                                                                                                                                                 ccGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGATGTGc  <  1:1913450/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCGCAGTACCTGGGCCTGGTGCGCAAGCCTGGGCGCGAGGCACTGGCCGTGGCCTTTGCCGAAGAGCTGCAGCGGTTCAAGCAGGAGCCGGCCTATGCAGCGATCAGCAACCGTTATACAACCGACACTGGAAACATTCTCAACGCCGTTGAGCAGCAGGAAAGCAGCACAGCGCGATAGCTCTGTTATACTCGGGCCTTCCCGCCCGGCTCACGCCCGGACGCTCGGCCTCGCAACAGGCATCCCGATCGGTACCGACGCCCCTTCGCGTCCAGCCTCCGTTTCCCGGATGTGC  >  NC_002947/5672411‑5672705

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: