Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,430,349 G→C 5.5% S189S (TCG→TCC PP_3045 → ATP‑dependent protease ClpP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,430,3490GC 5.5% 236.1 / 3.0 73S189S (TCG→TCCPP_3045ATP‑dependent protease ClpP
Reads supporting (aligned to +/- strand):  ref base G (29/40);  new base C (2/2);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

CGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATGGGG  >  NC_002947/3430202‑3430491
                                                                                                                                                   |                                                                                                                                              
cGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCg                                                                                                                                                >  5:145990/1‑148 (MQ=255)
 ggTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGa                                                                                                                                              <  4:38107/149‑1 (MQ=255)
   tGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACTTCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAgg                                                                                                                                            <  1:196534/149‑1 (MQ=255)
       gCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCCACGGCCTCCGAGGCCGt                                                                                                                                        >  8:173684/1‑148 (MQ=255)
        ccGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTTTAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTg                                                                                                                                       <  8:225457/149‑1 (MQ=255)
        ccGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGACGAAAGTTGGCTCACGGCCTCGGAGGCCGTg                                                                                                                                       >  6:188657/1‑149 (MQ=255)
          gAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCCGAGGCCGt                                                                                                                                        >  3:94144/1‑146 (MQ=255)
          gAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCCGAGGCCGTga                                                                                                                                      <  4:94144/148‑2 (MQ=255)
            ggACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGgcgc                                                                                                                                   <  4:297312/149‑1 (MQ=255)
            ggACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCCACGGCCTCCGAGGCCGTgagcca                                                                                                                                  <  7:173684/149‑6 (MQ=255)
            ggACCTTCGTCGGGTGGCGGACGTGCGTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGgcgc                                                                                                                                   <  6:64955/149‑1 (MQ=255)
               ccTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTgg                                                                                                                                <  4:159379/149‑1 (MQ=255)
                cTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTggg                                                                                                                               >  7:48855/1‑149 (MQ=255)
                 ttCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTgggg                                                                                                                              >  3:158517/1‑149 (MQ=255)
                  tcgtcgGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGc                                                                                                                             >  2:57836/1‑149 (MQ=255)
                  tcgtcgGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGc                                                                                                                             >  2:60281/1‑149 (MQ=255)
                     tcgGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTgg                                                                                                                          <  6:155009/149‑1 (MQ=255)
                            gCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGAt                                                                                                                     >  1:91855/1‑147 (MQ=255)
                            gCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGa                                                                                                                   <  1:292633/149‑1 (MQ=255)
                              ggACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAgg                                                                                                                 <  2:91855/149‑1 (MQ=255)
                               gACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGt                                                                                                                <  7:310199/149‑1 (MQ=255)
                                  gTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGGCgg                                                                                                             <  3:44520/149‑1 (MQ=255)
                                   tGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGCTGAGGTCgg                                                                                                             >  5:57166/1‑148 (MQ=255)
                                        gATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAAt                                                                                                         <  5:356024/147‑1 (MQ=255)
                                         aTCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGgt                                                                                                      <  8:48855/149‑1 (MQ=255)
                                              gCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCGGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTc                                                                                                   >  5:331238/1‑147 (MQ=255)
                                              gCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCaa                                                                                                 <  1:60281/149‑1 (MQ=255)
                                                cTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAgg                                                                                               >  3:88354/1‑149 (MQ=255)
                                                cTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAgg                                                                                               >  5:381517/1‑149 (MQ=255)
                                                 ttCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGt                                                                                              >  7:170121/1‑149 (MQ=255)
                                                    gAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTc                                                                                            <  1:57836/148‑1 (MQ=255)
                                                     aGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAg                                                                                          <  7:122488/149‑1 (MQ=255)
                                                         aaTGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATg                                                                                                                    <  6:310258/119‑1 (MQ=255)
                                                         aaTGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATg                                                                                                                    >  5:310258/1‑119 (MQ=255)
                                                         aaTGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGgtg                                                                                                     >  5:116541/1‑134 (MQ=255)
                                                         aaTGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGgtg                                                                                                     <  6:116541/134‑1 (MQ=255)
                                                         aaTGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCaa                                                                                                 <  8:295833/138‑1 (MQ=255)
                                                         aaTGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCaa                                                                                                 >  7:295833/1‑138 (MQ=255)
                                                         aaTGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGACGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGt                                                                                      <  5:188657/149‑1 (MQ=255)
                                                            gATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGcc                                                                                   >  2:222870/1‑149 (MQ=255)
                                                                gCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGCTGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTgg                                                                               <  6:57166/149‑1 (MQ=255)
                                                                gCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTgg                                                                               <  8:13251/149‑1 (MQ=255)
                                                                gCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTgg                                                                               <  3:227204/149‑1 (MQ=255)
                                                                gCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTg                                                                                <  6:168308/148‑1 (MQ=255)
                                                                    ccTATAAGGCCAAGGCTCCTAACAGCGATGACGCGGAGCGGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCGGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCaa                                                                           <  8:278112/149‑1 (MQ=255)
                                                                    ccTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTTGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCaa                                                                           <  6:144349/149‑1 (MQ=255)
                                                                     cTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAg                                                                          <  6:145990/149‑1 (MQ=255)
                                                                       ataAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGgc                                                                        >  7:123378/1‑149 (MQ=255)
                                                                           ggCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCgcgg                                                                    <  6:381517/149‑1 (MQ=255)
                                                                                aGGCTCCGAACATCGATGACGCGGAGTTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATg                                                               >  7:162212/1‑149 (MQ=255)
                                                                                     ccGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGcgcg                                                          >  7:36024/1‑149 (MQ=255)
                                                                                      cGAACATCGATGACGCGGAGTTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGgcgc                                                           <  8:162212/147‑1 (MQ=255)
                                                                                           aTCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCg                                                    <  8:123378/149‑1 (MQ=255)
                                                                                                    gcgGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATAcccc                                           >  5:94392/1‑149 (MQ=255)
                                                                                                    gcgGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATAcccc                                           >  8:309543/1‑149 (MQ=255)
                                                                                                      ggAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAg                                                                                                <  3:310280/94‑1 (MQ=255)
                                                                                                      ggAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAg                                                                                                >  4:310280/1‑94 (MQ=255)
                                                                                                             agaCGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGccc                                   >  3:294628/1‑148 (MQ=255)
                                                                                                               aCGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTAt                                >  5:287648/1‑149 (MQ=255)
                                                                                                                cGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTAtt                               <  4:158517/149‑1 (MQ=255)
                                                                                                                cGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTAtt                               <  6:287648/149‑1 (MQ=255)
                                                                                                                   cTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGa                            >  8:136246/1‑149 (MQ=255)
                                                                                                                      gTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGa                            <  4:294628/146‑1 (MQ=255)
                                                                                                                      gTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGa                            <  7:309543/146‑1 (MQ=255)
                                                                                                                       tCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTAGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAg                        >  1:42027/1‑149 (MQ=255)
                                                                                                                              cGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAgg                 >  1:5916/1‑149 (MQ=255)
                                                                                                                               gAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGc                >  6:256335/1‑149 (MQ=255)
                                                                                                                               gAAAGTTGGCTCACGGACTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGc                <  4:88354/149‑1 (MQ=255)
                                                                                                                                   gTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAg             >  1:363671/1‑148 (MQ=255)
                                                                                                                                       gCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAgg         <  7:136246/148‑1 (MQ=255)
                                                                                                                                         tCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTa           <  4:345762/144‑1 (MQ=255)
                                                                                                                                         tCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTa           >  3:345762/1‑144 (MQ=255)
                                                                                                                                            cGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATggg   <  8:36024/149‑1 (MQ=255)
                                                                                                                                             ggCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTGAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATgggg  <  1:112777/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
CGGTGATGCCGAGGACCTTCGTCGGGTGGCGGACGTGCTTGATCAGGCCTTCGAGGCAATGATCGCAGCCTATAAGGCCAAGGCTCCGAACATCGATGACGCGGAGCTGAGACGTCTGGTCAATGCCGAAAGTTGGCTCACGGCCTCGGAGGCCGTGGCGCTGGGGCTGGCCGATGAGGTCGGTAATGGTGTCAAGGTTCAGGCGTGCCTTGGGCAAGGCGCGGCGATGGCGCGTTACCGCAATACCCCGCAAGCCCTATTGGATCAGCTCAAGGCCAGTAGGGATGGGG  >  NC_002947/3430202‑3430491

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