Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,745,343 A→G 8.5% intergenic (‑198/+49) 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,3430AG 8.5% 527.8 / 26.0 180intergenic (‑198/+49)fbp/glgPclass 1 fructose‑1,6‑bisphosphatase/glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (88/72);  new base G (4/11);  total (92/88)
Fisher's exact test for biased strand distribution p-value = 5.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01

Only 100 of 183 total aligned reads displayed.

GGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCAC  >  NC_002947/5745194‑5745491
                                                                                                                                                     |                                                                                                                                                    
cgcaaggcggtattgaacacggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTACAGCGCAAGCCTTGaa                                                                                                                                                      >  1:1430556‑M1/121‑150 (MQ=255)
 ggAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAg                                                                                                                                                     >  1:1305167/1‑150 (MQ=255)
   agagTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCt                                                                                                                                                   <  2:2530159/150‑1 (MQ=255)
     agTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTg                                                                                                                                                 <  2:3385666/150‑1 (MQ=255)
     agTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTg                                                                                                                                                 <  2:3081547/150‑1 (MQ=255)
      gTGACGGGGTAGCGAGTGATGAGACTCGACACGGCTCTCAGAAGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCTTTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGc                                                                                                                                                <  2:1192275/150‑1 (MQ=255)
      gTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGc                                                                                                                                                <  2:1297117/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)
        gACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                                              <  1:4980543/150‑1 (MQ=255)
        gACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                                              <  2:5370587/150‑1 (MQ=255)
        gACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                                              <  1:2304514/150‑1 (MQ=255)
          tattgaacacggcccggatgggctggttctcgtcggaccgcaccattcgcgaagatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTGTCGCGCAAGCCTTGAAGCTTGCgaga                                                                                                                                            >  2:3952138‑M1/111‑150 (MQ=255)
            ttgaacacggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGaca                                                                                                                                          <  1:2089539‑M1/42‑1 (MQ=255)
             gCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCgctgtac                                                                                                                                         <  2:860325/150‑7 (MQ=37)
               tAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCt                                                                                                                                       >  1:2277305/1‑150 (MQ=255)
               tAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCt                                                                                                                                       >  1:4916295/1‑150 (MQ=255)
                  acggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgtg                                                                                                                                    >  2:2089539‑M1/103‑150 (MQ=255)
                  acggcacggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgtg                                                                                                                                    >  2:5113509‑M1/103‑150 (MQ=255)
                   gAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                                                   <  2:1685223/150‑1 (MQ=25)
                   gAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                                                   <  2:170643/150‑1 (MQ=25)
                       gatgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGc                                                                                                                               >  1:2382487/1‑150 (MQ=255)
                       acggatgggctggttctcgtcggaccgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgctCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGacac                                                                                                                                         <  2:967922‑M1/43‑1 (MQ=255)
                        atgaGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCg                                                                                                                              >  2:118756/1‑150 (MQ=255)
                          gagaCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCggg                                                                                                                            >  2:5257845/1‑150 (MQ=255)
                             aCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGt                                                                                                                         <  2:36921/150‑1 (MQ=255)
                             aCTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCACGCCCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTGGGagcgagcgc                                                                                                                         <  1:4771657/150‑2 (MQ=17)
                                  ggttctcgtctgaccgcaccattcgcgaagctgccaccgagatctggtaggcattggattgaggtgaggccgttcgcgggcgcgctCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcg                                                                                                                    >  2:1512995‑M1/87‑150 (MQ=255)
                                   cACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCcgc                                                                                                                   <  2:633523/150‑1 (MQ=255)
                                      ggCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                                <  2:2851189/150‑1 (MQ=255)
                                         tctcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTGTCTCGCAAGCTTCAAGGCTTGCgctgtac                                                                                                                                         >  1:3386572/1‑116 (MQ=21)
                                         tctcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACgg                                                                                                             >  2:2581869/1‑150 (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)
                                           tcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCTTCAAGGCTTGCgctgt                                                                                                                                           <  2:1857165/118‑5 (MQ=38)
                                           tcAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGcc                                                                                                           >  1:3233842/1‑150 (MQ=255)
                                              gACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                                              <  1:1354814/112‑1 (MQ=255)
                                              gACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                                              >  2:1354813/1‑112 (MQ=255)
                                                cgcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTgg                                                                                                                                   >  2:5408590‑M1/73‑121 (MQ=255)
                                                 ggCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTGCCTCGCAAGCTTCAAGGCTTGCGCTGTACCTGTgg                                                                                                                                   >  1:3808390/1‑120 (MQ=18)
                                                 gcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcgcgctCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                                <  1:3625728‑M1/68‑1 (MQ=255)
                                                 gcaccattcgcgaatatgccaccgagatctggaaggcattggattgaggtgaggccgttcgcgggcgcgctCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                                >  2:3625728‑M1/72‑139 (MQ=255)
                                                    acaATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCaa                                                                                                  >  1:1976683/1‑150 (MQ=255)
                                                      aaTAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                                              <  2:1175456/104‑1 (MQ=255)
                                                      aaTAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCga                                                                                                                                              >  1:1175456/1‑104 (MQ=255)
                                                      aaTAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATc                                                                                                >  2:3444568/1‑150 (MQ=255)
                                                           gTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATg                                                                                           >  2:3973329/1‑150 (MQ=255)
                                                              ggTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCt                                                                                        >  1:2418230/1‑150 (MQ=255)
                                                                tAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTc                                                                                      <  2:3241264/150‑1 (MQ=255)
                                                                      aGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtc                                                                                >  1:3155724/1‑150 (MQ=255)
                                                                      aGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtc                                                                                >  1:312349/1‑150 (MQ=255)
                                                                        ccTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCTTCAAGGCTTGCgc                                                                                                                                              <  2:4552783/86‑2 (MQ=37)
                                                                        ccTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGc                                                                                                                                                <  1:72781/84‑1 (MQ=255)
                                                                        ccTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGc                                                                                                                                                >  2:72781/1‑84 (MQ=255)
                                                                        ccTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCGGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtctc                                                                              >  1:1830919/1‑150 (MQ=255)
                                                                         ctgcGCCTGGCGGTTGTTTTGTTTGTTTGGCCATTCGGGGGCGCGCTCGAACCAACAGGTAAAGCGAACGCATTCATGCTTGCTAGACACGTGGGGGAGGGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCg                                                                             <  1:3451485/147‑1 (MQ=255)
                                                                          ttCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATc                                                                                                >  1:4889136/1‑130 (MQ=255)
                                                                          ttCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATc                                                                                                <  2:4889136/130‑1 (MQ=255)
                                                                           tCGCCTGGCGGGTGGTTTGTTTGTGTGTCCATTCGGGGGCGCGATCGCCCCCACTGGTACAGCGAATGCCTTGACGCTGGGGACACACCTGTGGGAGCGGGGGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGt                                                                           <  1:5257845/150‑1 (MQ=255)
                                                                            cGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTg                                                                          <  2:3873298/150‑1 (MQ=255)
                                                                              tggaaggcattggattgaggtgaggccgttcgcgggcacgccCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCtgt                                                                                                                                     <  1:1182409‑M1/47‑1 (MQ=255)
                                                                              cgcggCGGGTGTTGTGGTTGTGTGCCCATGCGGGGGCGCGCGCGCCCCCCCCGGTACCGCGCACGCCTTGAAGCTTGCGACACACGTGTGGGCGCGGGCGTCCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGc                                                                        <  1:1124074/147‑1 (MQ=255)
                                                                              cctggCGGGTGGTCTGGTTGTGTGTCCTTTCGCGGGCGCGCTCGCTCCCACAGGTACCGCGCACGCCTTGAAGCTTGCGAGACACCTGTGGGAGGGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGc                                                                        <  2:1976683/147‑1 (MQ=255)
                                                                              ccAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGc                                                                        >  2:3284027/1‑150 (MQ=255)
                                                                               cAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGc                                                                                                                                                <  1:2552923/77‑1 (MQ=255)
                                                                                          tCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCAGCTGCGCGAAt                                                            >  2:3524730/1‑150 (MQ=255)
                                                                                           cTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATg                                                           >  2:3023151/1‑150 (MQ=255)
                                                                                                tgtgtgTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCg                                                      <  1:3135913/150‑1 (MQ=255)
                                                                                                  tgtgtCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGGGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGt                                                    <  2:4362573/150‑1 (MQ=255)
                                                                                                       ccATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGACGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCcgac                                               <  2:4092505/150‑1 (MQ=255)
                                                                                                            cgggggCGCGCTCGCTCCCTCCGGTACCGCGCACGCCTTGAAGCTGGCGAGACACCGGGGTGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCTATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGaa                                          <  2:4866738/147‑1 (MQ=255)
                                                                                                                ggggcgcgcgcTCCACCCGGTACCGCGCACGCCTTGATGCTTGGGCGACCCGGGGGGGCGCGGGGGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAg                                      <  1:4148110/147‑1 (MQ=255)
                                                                                                                   gcgcTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGccc                                   >  2:2561713/1‑150 (MQ=255)
                                                                                                                     gctcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCAt                                 >  1:1981851/1‑150 (MQ=255)
                                                                                                                     gccCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTg                                                                          <  2:5038326‑M1/106‑1 (MQ=255)
                                                                                                                       tcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATcc                               >  2:1259850/1‑150 (MQ=255)
                                                                                                                       tcgctcCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATcc                               >  2:5174823/1‑150 (MQ=255)
                                                                                                                             CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccg                         <  1:5203666‑M1/150‑1 (MQ=255)
                                                                                                                             CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATccgtgccg                         <  2:4540529‑M1/150‑1 (MQ=255)
                                                                                                                             CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGAGCGCGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtct                                                                               >  2:151200‑M1/1‑96 (MQ=255)
                                                                                                                             CCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGAGCGCGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGAtct                                                                               <  1:151200‑M1/96‑1 (MQ=255)
                                                                                                                              CACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGc                                                       >  1:1706413‑M1/1‑119 (MQ=255)
                                                                                                                               ACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcgc                                                                <  1:5071176‑M1/109‑1 (MQ=255)
                                                                                                                                 AGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGt                                                         <  1:2155273‑M1/114‑1 (MQ=255)
                                                                                                                                 AGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGt                                                         >  2:2155273‑M1/1‑114 (MQ=255)
                                                                                                                                 AGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGtt                     <  1:5369757‑M1/150‑1 (MQ=255)
                                                                                                                                 AGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGtt                     >  2:2277305‑M1/1‑150 (MQ=255)
                                                                                                                                  GGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGt                                                                           >  1:3344882‑M1/1‑95 (MQ=255)
                                                                                                                                  GGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGt                                                                           <  2:3344882‑M1/95‑1 (MQ=255)
                                                                                                                                     aCAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAAt                 >  2:4228557/1‑150 (MQ=255)
                                                                                                                                      cAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATa                <  1:2581869/150‑1 (MQ=255)
                                                                                                                                         cgcAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCg             >  1:480220/1‑150 (MQ=255)
                                                                                                                                                cTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTgcg      >  1:3580805/1‑150 (MQ=255)
                                                                                                                                                 ttGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTcgc                                                                >  2:2075569/1‑91 (MQ=255)
                                                                                                                                                 ttGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATAGTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTgcgc     >  2:3055860/1‑150 (MQ=255)
                                                                                                                                                   gAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGcca   <  2:2418230/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGAGAGTGACGGGCTAGCGAGTGATGAGACTCGCCACGGCTCTCAGACGGCACAATAGGTCGGTAGTTCAGCCTTCGCCAGGCCGGTGGTCTGGTTGTGTGTCCATTCGCGGGCGCGCTCGCTCCCACAGGTACAGCGCAAGCCTTGAAGCTTGCGAGACACCTGTGGGAGCGGGCGTGCCCGCGAACGGCCTCACCTCAATCCAATGCCTTCCAGATCTCGGTGGCATATTCGCGAATGGTGCGGTCCGACGAGAACCAGCCCATCCGTGCCGTGTTCAATACCGCCTTGCGCCAC  >  NC_002947/5745194‑5745491

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: