Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,287,332 C→G 5.3% G135A (GGC→GCC)  PP_3758 ← two component system response regulator/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,287,3320CG 5.3% 242.5 / 3.7 75G135A (GGC→GCC) PP_3758two component system response regulator/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (35/36);  new base G (2/2);  total (37/38)
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.29e-01

GCGTACGCACCTCGTGGGAGACGATCGACATGAAGTCGTCACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGAC  >  NC_002947/4287189‑4287480
                                                                                                                                               |                                                                                                                                                     
gcgTACGCACCTCGTGGGAGACGATCGACATGAAGTCGTCACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGc                                                                                                                                                  <  6:178371/149‑1 (MQ=255)
      gCACCTCGTGGGAGACGATCGACATGAAGTCGTCACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTg                                                                                                                                            >  3:230032/1‑149 (MQ=255)
       cACCTCGTGGGAGACGATCGACATGAAGTCGTCACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTgc                                                                                                                                           <  6:108260/149‑1 (MQ=255)
             gTGGGAGACGATCGACATGAAGTCGTCACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTgg                                                                                                                                      <  5:419/148‑1 (MQ=255)
              tGGGAGACGATCGACATGAAGTCGTCACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCg                                                                                                                                    <  4:230032/149‑1 (MQ=255)
               gggAGACGATCGACATGAAGTCGTCACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCgg                                                                                                                                   <  8:307074/149‑1 (MQ=255)
                ggAGACGATCGACATGAAGTCGTCACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCgg                                                                                                                                   >  6:240658/1‑148 (MQ=255)
                ggAGACGATCGACATGAAGTCGTCACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCgg                                                                                                                                   >  4:318074/1‑148 (MQ=255)
                 gagaCGATCGACATGAAGTCGTCACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGaa                                                                                                                                 <  5:245474/149‑1 (MQ=255)
                  agaCGATCGACATGAAGTCGTCACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAAc                                                                                                                                <  1:236143/149‑1 (MQ=255)
                                aaGTCGTCACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAacac                                                                                                                  >  4:144689/1‑149 (MQ=255)
                                aaGTCGTCACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAacac                                                                                                                  >  5:87557/1‑149 (MQ=255)
                                  gtcgtcACGCATGCGCCCGGCGTGCTCCAGCTCACCCCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACAcc                                                                                                                 >  3:229192/1‑148 (MQ=255)
                                  gtcgtcACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCg                                                                                                                <  3:144689/149‑1 (MQ=255)
                                   tcgtcACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGGCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACAcc                                                                                                                 >  6:278242/1‑147 (MQ=255)
                                   tcgtcACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGGCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACAcc                                                                                                                 >  6:273451/1‑147 (MQ=255)
                                   tcgtcACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGGCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCga                                                                                                               <  5:278242/149‑1 (MQ=255)
                                   tcgtcACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGGCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCga                                                                                                               <  5:273451/149‑1 (MQ=255)
                                     gtcACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCgag                                                                                                              <  3:318074/148‑1 (MQ=255)
                                      tcACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGAc                                                                                                            <  5:197784/149‑1 (MQ=255)
                                             tGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCtt                                                                                                     <  7:87935/149‑1 (MQ=255)
                                                        tGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCATCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAg                                                                                         <  6:34201/149‑1 (MQ=255)
                                                        tGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAg                                                                                         <  5:240658/149‑1 (MQ=255)
                                                          cTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGt                                                                                                                            >  5:151227/1‑113 (MQ=255)
                                                          cTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGt                                                                                                                            <  6:151227/113‑1 (MQ=255)
                                                           tCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGgt                                                                                      >  4:168309/1‑149 (MQ=255)
                                                           tCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGgt                                                                                      >  3:99405/1‑149 (MQ=255)
                                                            ccAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGgtg                                                                                     >  2:270321/1‑149 (MQ=255)
                                                            ccAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTCTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGCTCGACGAACACCGAGACCTTGAATTTTACCGCCAGGgtg                                                                                     >  2:174022/1‑149 (MQ=255)
                                                             cAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGc                                                                                                       <  6:256020/131‑1 (MQ=255)
                                                                     acCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGc                                                                            <  6:87557/149‑1 (MQ=255)
                                                                     acCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGc                                                                            >  8:168215/1‑149 (MQ=255)
                                                                       cGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAATCACCTTGCGCTGGCGGAACAGGTCGACGATCACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCtt                                                                          <  4:229192/149‑1 (MQ=255)
                                                                           ccACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTtgtg                                                                       >  7:107412/1‑148 (MQ=255)
                                                                            cgcatGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGCCGAGCGGCTTGTGc                                                                      <  7:299580/144‑1 (MQ=255)
                                                                            cACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCa                                                                     >  8:258701/1‑149 (MQ=255)
                                                                                  cAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGgaa                                                                 <  8:107412/147‑1 (MQ=255)
                                                                                  cAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGaaa                                                                 <  7:258701/147‑1 (MQ=255)
                                                                                   aGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTc                                                              >  4:124064/1‑149 (MQ=255)
                                                                                       ggctgAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACg                                                          >  1:187596/1‑149 (MQ=255)
                                                                                        gctgAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCAc                                                           <  3:168309/147‑1 (MQ=255)
                                                                                                   gcTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCg                                              >  5:59496/1‑149 (MQ=255)
                                                                                                     tCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGc                                                                      <  1:124512/123‑1 (MQ=255)
                                                                                                     tCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGc                                                                      >  2:124512/1‑123 (MQ=255)
                                                                                                     tCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCg                                              >  8:181230/1‑147 (MQ=255)
                                                                                                     tCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCg                                              >  7:17941/1‑147 (MQ=255)
                                                                                                      ccTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTtgtg                                                                       <  8:150987/121‑1 (MQ=255)
                                                                                                      ccTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTtgtg                                                                       >  7:150987/1‑121 (MQ=255)
                                                                                                      ccTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAg                                           >  2:31223/1‑149 (MQ=255)
                                                                                                          ttCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCt                                       <  3:124064/149‑1 (MQ=255)
                                                                                                                gcggcTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGGTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAg                                  <  1:174022/148‑1 (MQ=255)
                                                                                                                gcggcTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCACCGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAg                                  <  1:270321/149‑1 (MQ=255)
                                                                                                                          tCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCa                                                            <  2:337123/112‑1 (MQ=255)
                                                                                                                          tCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCa                                                            >  1:337123/1‑112 (MQ=255)
                                                                                                                          tCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACg                                                          >  1:191010/1‑114 (MQ=255)
                                                                                                                          tCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACg                                                          <  2:191010/114‑1 (MQ=255)
                                                                                                                          tCCAGGGCTTCCAGCTGCCGGCCAAGCACATTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCa                        >  5:45452/1‑148 (MQ=255)
                                                                                                                           ccAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCAtc                      >  4:290661/1‑149 (MQ=255)
                                                                                                                             aGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGaaa                                                                 <  7:179867/104‑1 (MQ=255)
                                                                                                                             aGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGaaa                                                                 >  8:179867/1‑104 (MQ=255)
                                                                                                                               ggCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTcgc                  >  5:244272/1‑149 (MQ=255)
                                                                                                                                gCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTcgcg                 >  4:201016/1‑149 (MQ=255)
                                                                                                                                gCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTcgcg                 <  2:187596/149‑1 (MQ=255)
                                                                                                                                 cTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCAtctacggc                >  3:33946/1‑144 (MQ=255)
                                                                                                                                       gCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCccgcc          >  4:293084/1‑149 (MQ=255)
                                                                                                                                          gCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGcc                                         >  2:133354/1‑115 (MQ=255)
                                                                                                                                          gCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGcc                                         <  1:133354/115‑1 (MQ=255)
                                                                                                                                           ccGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCgcgc      >  2:27123/1‑149 (MQ=255)
                                                                                                                                            cGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCcgcgtt     <  3:290661/149‑3 (MQ=255)
                                                                                                                                             ggCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCGGGCcgcgaga    <  4:99405/149‑4 (MQ=255)
                                                                                                                                             ggCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCGGGCcgcgaga    <  3:201016/149‑4 (MQ=255)
                                                                                                                                             ggCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTg    <  1:31223/149‑1 (MQ=255)
                                                                                                                                              gCCAAGCACCTTGCGCTGGGGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTcgcg                 >  3:135593/1‑135 (MQ=255)
                                                                                                                                              gCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTcgcg                 <  4:135593/135‑1 (MQ=255)
                                                                                                                                              gCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGa   >  1:269442/1‑149 (MQ=255)
                                                                                                                                               ccAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGAc  >  7:252290/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GCGTACGCACCTCGTGGGAGACGATCGACATGAAGTCGTCACGCATGCGCACGGCGTGCTCCAGCTCACACCGTGCCACCTGCAGCTGGCTGAGCAGCAGCTCCTGTTCCTGGCGGCTCTGCTCCAGGGCTTCCAGCTGCCGGCCAAGCACCTTGCGCTGGCGGAACAGGTCGACGAACACCGAGACCTTGCTTTTTACCGCCAGGGTGTCGAGCGGCTTGTGCAGAAAGTCCACGGCCCCGCTCTCGTAGCCCTTGAAGGCATAATTCATCTCGCGGCCCGCCGCGCTGAC  >  NC_002947/4287189‑4287480

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