Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,576,195 G→C 8.3% G301R (GGC→CGC)  glgX → glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,576,1950GC 8.3% 374.0 / 17.9 118G301R (GGC→CGC) glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base G (54/54);  new base C (6/4);  total (60/58)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

CGCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGCGAAAT  >  NC_002947/4576046‑4576343
                                                                                                                                                     |                                                                                                                                                    
cGCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGcgcg                                                                                                                                                      <  2:174071/150‑1 (MQ=255)
 gCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCgg                                                                                                                                                     >  2:3242604/1‑150 (MQ=255)
     ccGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATc                                                                                                                                                 <  1:2275830/150‑1 (MQ=255)
       gcgcTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGa                                                                                                                                               <  1:1917849/150‑1 (MQ=255)
       gcgcTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGa                                                                                                                                               <  1:3957387/150‑1 (MQ=255)
         gcTACCTGGCCAGCGGCAAGATCGCCGAGTTAAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTATACAACCACACCGCCGATGGAAACGAGCGCGGCCCGACCCTGTCTATGCCGCGCATCGACa                                                                                                                                             <  1:4647376/150‑1 (MQ=255)
          cTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACaa                                                                                                                                            >  1:1589621/1‑150 (MQ=255)
          cTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCCGCGCATCGACag                                                                                                                                            <  2:2105073/150‑2 (MQ=255)
           tACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAAc                                                                                                                                           >  2:2813297/1‑150 (MQ=255)
           tACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAAc                                                                                                                                           >  1:1859476/1‑150 (MQ=255)
              cTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAATGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACTACGcc                                                                                                                                        <  1:5195413/150‑1 (MQ=255)
                  ccAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGt                                                                                                                                    >  2:846527/1‑150 (MQ=255)
                   cAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGAc                                                                                                                                              <  1:3641947/139‑1 (MQ=255)
                   cAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGAc                                                                                                                                              >  2:3641947/1‑139 (MQ=255)
                        gCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTAcc                                                                                                                              >  1:2904098/1‑150 (MQ=255)
                         cAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCg                                                                                                                             >  2:802886/1‑150 (MQ=255)
                          aaGATCGCCGAGTTCATGGAGTTGGTCGCGCACCGGCACGACGCGGGGCTTGCGGTGATCCTGGAGGTGGTCTACAACCACACCGCCTACGGCAACGCGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGc                                                                                                                            <  1:5115937/150‑1 (MQ=255)
                           aGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGcc                                                                                                                           >  1:5377249/1‑150 (MQ=255)
                                 ccGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATc                                                                                                                                                 <  2:2717869/122‑1 (MQ=255)
                                 ccGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATc                                                                                                                                                 >  1:2717869/1‑122 (MQ=255)
                                     gTTCAAGGAGATGGTCGCGCACCTGCACGCCGCGGGGCTTGAGGTGATCCTGGAGGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGga                                                                                                                 <  2:1469719/150‑1 (MQ=255)
                                     gTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGga                                                                                                                 >  2:1316119/1‑150 (MQ=255)
                                     gTTAAAGGAGATGGTCGCGCACCTGCACGACGCGGGGATTGTGGTGATCCTGGACGTGGTCTACAACCACACCGCCGACGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGga                                                                                                                 <  1:2578805/150‑1 (MQ=255)
                                         aaGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGTAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGgacgac                                                                                                             >  1:3416505/1‑150 (MQ=255)
                                         aaGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGgacgac                                                                                                             >  1:1722592/1‑150 (MQ=255)
                                             agaTGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAgc                                                                                                         <  2:2752794/150‑1 (MQ=255)
                                             agaTGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAgc                                                                                                         >  2:4244688/1‑150 (MQ=255)
                                              gaTGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAgcg                                                                                                        >  2:558936/1‑150 (MQ=255)
                                                  gTCGCGCACCTGAACGACGGGGGGTTTGAGGTGATCCTGGACGTGTTCTACAACAACACCGCCGAAGGAAACGAGCGCGGCCCTACCCTGTCGATGCGCGGCATCTACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGctac                                                                                                    <  1:4759184/150‑1 (MQ=255)
                                                    cgcgCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGctacta                                                                                                  <  2:4973949/150‑1 (MQ=255)
                                                     gcgcACCTGCACGACGCGGGGCTTGTGGTGATCCTGGACGTGTTCTACAACCACACCGCCGAAGGAAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCtactac                                                                                                 <  2:2382652/150‑1 (MQ=255)
                                                       gcACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACACGCGCTACTACAt                                                                                               >  1:4333351/1‑150 (MQ=255)
                                                        cACCTGCAGGACGGGGGGGTTGTGTTGAGCCTGTACGTGTTCTACAACCACACCGCGGATGGCAACGAGCGGGGCCCGACCCTGTCTATGCGGGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATc                                                                                              <  1:232167/150‑1 (MQ=255)
                                                        cACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATc                                                                                              <  2:1476682/150‑1 (MQ=255)
                                                           cTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAAc                                                                                           >  1:897625/1‑150 (MQ=255)
                                                              cACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACCAGCGCTACTACATCAACGAc                                                                                        >  1:2180718/1‑150 (MQ=255)
                                                                cgccgcgGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCTAAGGCAACGTGCGCGGCCCGACCCTGTCTATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTc                                                                                      <  2:4800605/147‑1 (MQ=255)
                                                                   cgcgGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTccgg                                                                                   >  2:4973429/1‑150 (MQ=255)
                                                                    gcgGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACCCGCGCTACTACATCAACGACTccggc                                                                                  >  1:967703/1‑150 (MQ=255)
                                                                            tGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCaa                                                                          <  2:1098643/150‑1 (MQ=255)
                                                                                gttaTCGTGGACGTGGTCTACACCCACACCGCCGAAGGAAAGGCGCGCGGCCCGCCCCTGTCGATGCGCGGCATCGAAAACGCCTCTTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACAcc                                                                      <  1:2577051/147‑1 (MQ=255)
                                                                                gTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACAcc                                                                      >  2:2065563/1‑150 (MQ=255)
                                                                                   aTCCTGTACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGAAAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTg                                                                   <  1:3153999/150‑1 (MQ=255)
                                                                                   aTCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGctac                                                                                                    <  2:3111460/117‑1 (MQ=255)
                                                                                    tCCTGGACGTGGTCTACAACCACACCGCCGACGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTgg                                                                  <  1:4464391/150‑1 (MQ=255)
                                                                                     ccTGGACGTGTTCTACAACCACACCGCCTAAGGAAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGa                                                                 <  2:2904098/150‑1 (MQ=255)
                                                                                      cTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGCGCGCGGCCCGACCGTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGAc                                                                <  2:2028990/150‑1 (MQ=255)
                                                                                       tGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGAcc                                                               <  2:3764520/150‑1 (MQ=255)
                                                                                        ggACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCt                                                              >  2:2050289/1‑150 (MQ=255)
                                                                                        ggACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCt                                                              >  2:2830561/1‑150 (MQ=255)
                                                                                          tggTGGTCTACAACCACCCCGCCGATGGCAACGCGGGCGGCCCGACCGTGTCTAGGCGGGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGa                                                            <  1:3219979/148‑1 (MQ=255)
                                                                                          aCGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCTACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGa                                                            <  2:2388669/150‑1 (MQ=255)
                                                                                          aCGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGa                                                            <  2:3505246/150‑1 (MQ=255)
                                                                                           cGTGTTCTACAACCCCACCGCCTATGGAAACGAGCGCGGCCCGACCCTGTGTATGCGCGGCATCGAAAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAg                                                           <  1:3383082/150‑1 (MQ=255)
                                                                                           cGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAg                                                           >  2:1248773/1‑150 (MQ=255)
                                                                                                 cTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAccc                                                     >  1:748587/1‑150 (MQ=255)
                                                                                                 atACAACCAAACCCCCTATGCAAACGCGCGCGCCCCGACCCTGTCGAGGCGCGCCAACCAAAACCCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAccc                                                     <  1:1597592/149‑1 (MQ=255)
                                                                                                 aaaaaaCCACACCGCCTAAGCAAACGAGCGCGGCCCCACCCTGTCTATGCGCGCAAACAACAACGCATCGTACTACCGCCTTATGCCGGACGACAAGCGCTACTACAACAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAccc                                                     <  1:1597757/146‑1 (MQ=255)
                                                                                                  tACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAcccc                                                    >  2:2922547/1‑150 (MQ=255)
                                                                                                  tACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCAcccc                                                    >  2:1063894/1‑150 (MQ=255)
                                                                                                      aCCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCg                                                >  2:4471510/1‑150 (MQ=255)
                                                                                                       ccACACCGCCTAAGGAAACGAGCGGGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACACGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGt                                               <  1:3767395/150‑1 (MQ=255)
                                                                                                       ccACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACCAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGt                                               >  2:4859745/1‑150 (MQ=255)
                                                                                                        cacaCCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGtg                                              >  2:2586739/1‑150 (MQ=255)
                                                                                                         acacCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGCCAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGtgc                                             <  1:2677927/150‑1 (MQ=255)
                                                                                                         acacCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCCGCGCATCGAc                                                                                                                                              >  2:2601151/1‑53 (MQ=37)
                                                                                                           acCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGtgctg                                           <  2:4028837/150‑1 (MQ=255)
                                                                                                           acCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGtgctg                                           <  2:4017229/150‑1 (MQ=255)
                                                                                                              gccgAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAg                                        >  1:4419457/1‑150 (MQ=255)
                                                                                                               ccgAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGt                                       <  2:3551627/150‑1 (MQ=255)
                                                                                                                 taagcaaaCCAGCGCGGCCCCACCCTGTCGATGCGCGCAAACCAAAACGCCTCGTACTACCGCCTCATGCCGGACCACAAGCGCTACTACAACAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTg                                     <  1:1165528/144‑1 (MQ=255)
                                                                                                                 gacggCCACGCGCGCGGCCCGCCCCTGGCGAGGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTCCATCCACGACTCCGGCACCGGCAACACCCGGGACCTGAGCCACCCCTGCGTGCTGCAGTTg                                     <  1:4054593/147‑1 (MQ=255)
                                                                                                                 gAAGGCAACGCGCGCGGCCCGACCCTGTCGAGGCGGGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTg                                     <  2:1589621/150‑1 (MQ=255)
                                                                                                                 gAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTg                                     <  1:4054461/150‑1 (MQ=255)
                                                                                                                 caaGGAAACGAGCGCGGCCCGACCCTGTCCATGCGCGCAATCCACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCCCAACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTg                                     <  2:3743170/149‑1 (MQ=255)
                                                                                                                  aaGGCAACGCGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTgg                                    <  1:377532/150‑1 (MQ=255)
                                                                                                                   cggCAACGCGGGCGGCCCGCCCCTGTCGATGCGCGGAATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGt                                   <  1:2547685/149‑1 (MQ=255)
                                                                                                                    ggCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACACGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTc                                  >  1:5078474/1‑150 (MQ=255)
                                                                                                                      cAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCAc                                >  1:2519875/1‑150 (MQ=255)
                                                                                                                       aaCGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCAcc                               >  2:136706/1‑150 (MQ=255)
                                                                                                                          gAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGAc                            >  2:916780/1‑150 (MQ=255)
                                                                                                                                  cccGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgcttcg                    >  2:4463312/1‑150 (MQ=255)
                                                                                                                                   ccGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGAAGCCGGACGACAATCGCTACTACATCAACGACTCCGGCACCGGCAACACCCAGGACCTGAGCCACCCCTGCTTGCTGCAGTTGGTAACCGACtcgcttcgc                   >  1:2197075/1‑150 (MQ=255)
                                                                                                                                    cGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACtcgcttcgct                  >  1:756724/1‑150 (MQ=255)
                                                                                                                                     gACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACCAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTa                 >  2:5115696/1‑150 (MQ=255)
                                                                                                                                        gttGTCGATGCCGCGCATCGACAGCGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTg              >  1:3230943/3‑150 (MQ=255)
                                                                                                                                        gttGTCGATGCCGCGCATCGACAGCGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGCCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTg              >  1:1201447/3‑150 (MQ=255)
                                                                                                                                        ccTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTg              >  2:2033669/1‑150 (MQ=255)
                                                                                                                                        ccTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTg              <  2:5377249/150‑1 (MQ=255)
                                                                                                                                         ttGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTgg             >  2:461978/2‑150 (MQ=255)
                                                                                                                                         ttGTCGATGCCGCGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGtg                                              <  2:3892486/116‑1 (MQ=255)
                                                                                                                                         ttGTCGATGCCGCGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTgg             >  2:3647203/2‑150 (MQ=255)
                                                                                                                                         cTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTgg             <  2:2266362/150‑1 (MQ=255)
                                                                                                                                            tCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGcc          <  1:4463312/150‑1 (MQ=255)
                                                                                                                                             cGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTg                                     <  2:2062835/122‑1 (MQ=255)
                                                                                                                                             cGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTg                                     >  1:2062835/1‑122 (MQ=255)
                                                                                                                                              gATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCgg        >  1:2712066/1‑150 (MQ=255)
                                                                                                                                               ttGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTCCATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGc       <  1:3621560/149‑1 (MQ=255)
                                                                                                                                                   gcgGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACACGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACACCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGCGaaa   >  1:2295384/1‑150 (MQ=255)
                                                                                                                                                    cgGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGCGAAAt  <  2:3416505/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCACCCGCGCTACCTGGCCAGCGGCAAGATCGCCGAGTTCAAGGAGATGGTCGCGCACCTGCACGACGCGGGGCTTGAGGTGATCCTGGACGTGGTCTACAACCACACCGCCGAAGGCAACGAGCGCGGCCCGACCCTGTCGATGCGCGGCATCGACAACGCCTCGTACTACCGCCTGATGCCGGACGACAAGCGCTACTACATCAACGACTCCGGCACCGGCAACACCCTGGACCTGAGCCACCCCTGCGTGCTGCAGTTGGTCACCGACTCGCTTCGCTACTGGGCCGGCGAAAT  >  NC_002947/4576046‑4576343

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: