Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,745,327 A→G 8.2% intergenic (‑182/+65) fbp ← / ← glgP class 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,745,3270AG 8.2% 557.8 / 25.9 183intergenic (‑182/+65)fbp/glgPclass 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (84/84);  new base G (6/9);  total (90/93)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

Only 100 of 185 total aligned reads displayed.

ACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAAT  >  NC_002947/5745178‑5745476
                                                                                                                                                     |                                                                                                                                                     
aCCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTa                                                                                                                                                       <  2:1804889/150‑1 (MQ=255)
  ccACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACa                                                                                                                                                     <  1:603900/150‑1 (MQ=255)
  ccACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACa                                                                                                                                                     <  2:1762056/150‑1 (MQ=255)
          tttACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACATCGCAAGc                                                                                                                                             <  2:1607604/150‑1 (MQ=255)
               aGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGa                                                                                                                                        >  1:3158202/1‑150 (MQ=255)
                cgcaaggcggtattgaacacggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTACAGCGCAAGCCTTGaa                                                                                                                                       >  1:1430556‑M1/121‑150 (MQ=255)
                 ggAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAg                                                                                                                                      >  1:1305167/1‑150 (MQ=255)
                   agagTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCt                                                                                                                                    <  2:2530159/150‑1 (MQ=255)
                     agTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTg                                                                                                                                  <  1:4851335/150‑1 (MQ=255)
                     agTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTg                                                                                                                                  <  2:1649587/150‑1 (MQ=255)
                     agTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTg                                                                                                                                  <  2:3385666/150‑1 (MQ=255)
                     agTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTg                                                                                                                                  <  2:3081493/150‑1 (MQ=255)
                      gTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGc                                                                                                                                 <  2:3949938/150‑1 (MQ=255)
                      gTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCAAAGGTACAGCGCAAGCCTTGAAGCTTGc                                                                                                                                 <  1:5404961/150‑1 (MQ=255)
                        ggtattgaacacggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTGTCTCGCAAGCCTTGAAGCTTGCga                                                                                                                               >  1:391887‑M1/113‑150 (MQ=255)
                        ggtattgaacacggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTGTCTCGCAAGCCTTGAAGCTTGCga                                                                                                                               >  1:4387059‑M1/113‑150 (MQ=255)
                        gACGGTCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCTGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                               <  2:5370574/150‑1 (MQ=255)
                        gACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                               >  2:1264408/1‑150 (MQ=255)
                        gACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                               <  2:5370587/150‑1 (MQ=255)
                        gACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                               <  2:658501/150‑1 (MQ=255)
                          tattgaacacggcccggatgggctggttctcgtcggaccgcaccattcgcgaagatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTGTCGCGCAAGCCTTGAAGCTTGCgaga                                                                                                                             >  2:3952138‑M1/111‑150 (MQ=255)
                            ttgaacacggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGaca                                                                                                                           <  1:2089539‑M1/42‑1 (MQ=255)
                             tgaacacggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgctCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGacac                                                                                                                          >  2:5191916‑M1/108‑150 (MQ=255)
                             gCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCgctgtac                                                                                                                          <  2:823692/150‑7 (MQ=37)
                             gCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCgctgtac                                                                                                                          <  2:860325/150‑7 (MQ=37)
                               tAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCtt                                                                                                                                          >  1:5142653/1‑132 (MQ=255)
                                 gCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgt                                                                                                                      <  1:2193552/150‑1 (MQ=255)
                                 gCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgt                                                                                                                      <  1:2614612/150‑1 (MQ=255)
                                  acggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgtg                                                                                                                     >  2:5113509‑M1/103‑150 (MQ=255)
                                   cggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgtgtccattcgcgggcgcgctCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTgg                                                                                                                    >  1:4149953‑M1/102‑150 (MQ=255)
                                       acggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgctCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGacac                                                                                                                          >  1:967922‑M1/98‑140 (MQ=255)
                                       acggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgctCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGacac                                                                                                                          <  2:967922‑M1/43‑1 (MQ=255)
                                        atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCg                                                                                                               >  2:118756/1‑150 (MQ=255)
                                           atgggcgggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACCGGTACAGCGCACGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGc                                                                                                            >  2:291702‑M1/94‑150 (MQ=255)
                                             aCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGt                                                                                                          <  2:36921/150‑1 (MQ=255)
                                             aCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGt                                                                                                          <  2:3158202/150‑1 (MQ=255)
                                                   cACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcgc                                                                                                    >  2:2125467/1‑150 (MQ=255)
                                                     cGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGa                                                                                                  >  2:1215639/1‑150 (MQ=255)
                                                      ggCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                 <  2:2851189/150‑1 (MQ=255)
                                                      ggCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCACGCCcgctcc                                                                                                 <  2:2382487/150‑4 (MQ=255)
                                                         tctcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACga                                                                                              >  2:3816054/1‑149 (MQ=255)
                                                          ctcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGc                                                                                             >  1:2823668/1‑150 (MQ=255)
                                                           tcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGTCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcc                                                                                            >  1:1741440/1‑150 (MQ=255)
                                                           tcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCTTCAAGGCTTGCgctgt                                                                                                                            >  1:1857165/1‑114 (MQ=38)
                                                           tcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcc                                                                                            >  1:3233842/1‑150 (MQ=255)
                                                             aGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcctc                                                                                          >  2:1916909/1‑150 (MQ=255)
                                                              gACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                               >  2:1354813/1‑112 (MQ=255)
                                                              gACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                               <  1:1354814/112‑1 (MQ=255)
                                                                 ggCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTGCCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                                    >  1:3808390/1‑120 (MQ=18)
                                                                 gcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcgcgctCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                 <  1:3625728‑M1/68‑1 (MQ=255)
                                                                    acaATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCaa                                                                                   >  1:1976683/1‑150 (MQ=255)
                                                                      aaTAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                               >  1:1175456/1‑104 (MQ=255)
                                                                      aaTAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATc                                                                                 >  2:3444804/1‑150 (MQ=255)
                                                                      aaTAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATc                                                                                 >  2:3444568/1‑150 (MQ=255)
                                                                            tCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTg                                                                                                                                         >  2:1248435/1‑88 (MQ=255)
                                                                            tCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTg                                                                                                                                         <  1:1248435/88‑1 (MQ=255)
                                                                              ggTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCt                                                                         >  1:2418230/1‑150 (MQ=255)
                                                                                      aGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtc                                                                 >  1:312349/1‑150 (MQ=255)
                                                                                        ccTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCTTCAAGGCTTGCgc                                                                                                                               >  1:4552783/1‑85 (MQ=37)
                                                                                        ccTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCTTCAAGGCTTGCgc                                                                                                                               <  2:4552783/86‑2 (MQ=37)
                                                                                        ccTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCGGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtctc                                                               >  1:1830919/1‑150 (MQ=255)
                                                                                         ctgcGCCTGGCGGTTGTTTTGTTTGTTTGGCCATTCGGGGGCGCGCTCGAACCAACAGGTAAAGCGAACGCATTCATGCTTGCTAGACACGTGGGGGAGGGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCg                                                              <  1:3451485/147‑1 (MQ=255)
                                                                                          ttCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAAt                                                                                  <  2:3570810/129‑1 (MQ=255)
                                                                                          ttCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATc                                                                                 <  2:4889136/130‑1 (MQ=255)
                                                                                          ttCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCgg                                                             >  2:5145049/1‑150 (MQ=255)
                                                                                           tCGCCTGGCGGGTGGTTTGTTTGTGTGTCCATTCGGGGGCGCGATCGCCCCCACTGGTACAGCGAATGCCTTGACGCTGGGGACACACCTGTGGGAGCGGGGGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGt                                                            <  1:5257845/150‑1 (MQ=255)
                                                                                            cGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTg                                                           <  2:3873298/150‑1 (MQ=255)
                                                                                              tggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgt                                                                                                                      <  1:1182409‑M1/47‑1 (MQ=255)
                                                                                              tggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgt                                                                                                                      >  2:1182409‑M1/43‑89 (MQ=255)
                                                                                              cgcggCGGGTGTTGTGGTTGTGTGCCCATGCGGGGGCGCGCGCGCCCCCCCCGGTACCGCGCACGCCTTGAAGCTTGCGACACACGTGTGGGCGCGGGCGTCCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGc                                                         <  1:1124074/147‑1 (MQ=255)
                                                                                              ccAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGc                                                         >  2:3284027/1‑150 (MQ=255)
                                                                                               cAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGc                                                                                                                                 <  1:2552923/77‑1 (MQ=255)
                                                                                               aaGGCTGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCa                                                        >  2:560332/2‑150 (MQ=255)
                                                                                                      gtggtCTGTTTGTGTGTCCTTTCGGGGGCGCGCTCGCTCCCACTGGTACAGCGAACGCCTTGAAGCTTGCGAGACACGTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcgc                                                 <  1:2125467/150‑1 (MQ=255)
                                                                                                                tgtgtgTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCg                                       <  1:3135913/150‑1 (MQ=255)
                                                                                                                tgtgtgGCCTTTGGGGGGGGCGCGCGCTCCCACTGGTACCGCGCACGCATTCACGGTGGCGCCCCACGGGGGGGCGGGGGGGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCg                                       <  1:2374927/150‑1 (MQ=255)
                                                                                                                  tgtgtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGGGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGt                                     <  2:4362573/150‑1 (MQ=255)
                                                                                                                   gtgtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTc                                    >  2:2014965/1‑150 (MQ=255)
                                                                                                                    tgtCCAGTCGCGGGCGCGCTCGCTCCCACTTGTACAGCGCACGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTcc                                   <  1:4972881/150‑1 (MQ=255)
                                                                                                                       ccATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGACGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                                <  2:4092505/150‑1 (MQ=255)
                                                                                                                                   gcgcTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGccc                    >  2:2561713/1‑150 (MQ=255)
                                                                                                                                     gccCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTg                                                           >  1:5038326‑M1/4‑109 (MQ=255)
                                                                                                                                      ctcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATc                 <  2:3963782/150‑1 (MQ=255)
                                                                                                                                             CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGGCTCTGCGCGAATGGTGCGGTCCGACGAGACCCAGCCCATccgggccg          <  2:4598881‑M1/150‑1 (MQ=255)
                                                                                                                                             CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccg          <  1:2300652‑M1/150‑1 (MQ=255)
                                                                                                                                             CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccg          <  2:4540529‑M1/150‑1 (MQ=255)
                                                                                                                                             CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccg          <  1:4906939‑M1/150‑1 (MQ=255)
                                                                                                                                             CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccg          <  1:5203666‑M1/150‑1 (MQ=255)
                                                                                                                                             CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGAGCGCGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtct                                                                >  2:151200‑M1/1‑96 (MQ=255)
                                                                                                                                             CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGAGCGCGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtct                                                                <  1:151200‑M1/96‑1 (MQ=255)
                                                                                                                                              CACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGc                                        <  2:1706413‑M1/119‑1 (MQ=255)
                                                                                                                                              CACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGc                                        >  1:1706413‑M1/1‑119 (MQ=255)
                                                                                                                                               ACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcgc                                                 >  2:5071176‑M1/1‑109 (MQ=255)
                                                                                                                                                 AGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGc                                                                                             <  1:3367204‑M1/63‑1 (MQ=255)
                                                                                                                                                 AGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGt                                          >  2:2155273‑M1/1‑114 (MQ=255)
                                                                                                                                                 AGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGtt      >  2:2277305‑M1/1‑150 (MQ=255)
                                                                                                                                                 AGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGtt      <  1:5369757‑M1/150‑1 (MQ=255)
                                                                                                                                                  GGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGt                                                            <  2:3344882‑M1/95‑1 (MQ=255)
                                                                                                                                                  GGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGt                                                            >  1:3344882‑M1/1‑95 (MQ=255)
                                                                                                                                                     aCAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAAt  >  2:4228557/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCCACGCAGTTTACAGGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAAT  >  NC_002947/5745178‑5745476

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