Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,577,914:1 +A 5.8% intergenic (+79/+30) PP_4057 → / ← glgB membrane protein/1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,577,9141.A 5.8% 544.8 / 24.3 137intergenic (+79/+30)PP_4057/glgBmembrane protein/1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base . (60/69);  new base A (5/3);  total (65/72)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 141 total aligned reads displayed.

CAGTGATGTGGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCC  >  NC_002947/4577769‑4578063
                                                                                                                                                  |                                                                                                                                                     
cagTGATGTGGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCcg                                                                                                                                                   >  1:26760/1‑150 (MQ=255)
          gATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggcc                                                                                                                                         <  2:403175‑M2/150‑6 (MQ=255)
           ggccctgaacctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggctgtactggccccttcgcgggtaaACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCa                                                                                                                                        >  1:896945‑M1/82‑150 (MQ=255)
           aTGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCa                                                                                                                                        >  1:1886934/1‑150 (MQ=255)
             gACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGt                                                                                                                                      <  1:1034324/150‑1 (MQ=255)
             gACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGt                                                                                                                                      >  1:1455258/1‑150 (MQ=255)
             gACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGt                                                                                                                                      >  1:1358933/1‑150 (MQ=255)
               cAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgg                                                                                                                                            <  1:6874‑M2/142‑3 (MQ=255)
               cAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgg                                                                                                                                            >  2:6874‑M2/1‑140 (MQ=255)
                     ctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggctgtaccggccccttcgcgggtaaACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCaa                                                                                                                              >  2:300221‑M1/72‑150 (MQ=255)
                     ccTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCaa                                                                                                                              <  1:2403812/150‑1 (MQ=255)
                      cTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAAc                                                                                                                             <  1:1711396/150‑1 (MQ=255)
                           cGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATc                                                                                                                        <  2:901527/150‑1 (MQ=255)
                            gggTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGAACACCAAACCgctac                                                                                                                       >  2:2352496/1‑146 (MQ=255)
                            gggTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTAGAGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCa                                                                                                                       >  1:1013572/1‑150 (MQ=255)
                              gTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCAGCTCCTACACGGGTCACCGCCGGGGGCGGGGCCTTGTTAGGAGCGGCACTCCCGCGAAGTGTCCGGTTCAGCCAACCGagcggg                                                                                                                     >  2:2343278/1‑145 (MQ=255)
                               ttGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggc                                                                                                                                          <  2:215604‑M2/128‑5 (MQ=255)
                               ttGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggc                                                                                                                                          >  1:215604‑M2/1‑124 (MQ=255)
                                      gCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAg                                                                                                                      <  2:216130/141‑1 (MQ=255)
                                      gCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGc                                                                                                                    >  2:588677/1‑143 (MQ=255)
                                      gCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGc                                                                                                                    >  2:588676/1‑143 (MQ=255)
                                      gCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGNCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGc                                                                                                                    <  1:588676/143‑1 (MQ=255)
                                          ccTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGa                                                                                                         <  2:1455258/150‑1 (MQ=255)
                                            ttGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATc                                                                                                       <  2:1197359/150‑1 (MQ=255)
                                               gggTGCGGATGCGGTGGTATGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAac                                                                                                    <  2:182535/150‑1 (MQ=255)
                                               gggTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAac                                                                                                    <  1:403175/150‑1 (MQ=255)
                                                 gTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGtct                                                                                                                                                        <  2:832516/97‑1 (MQ=255)
                                                 gTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGtct                                                                                                                                                        >  1:832516/1‑97 (MQ=255)
                                                 gTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAacac                                                                                                  >  2:353842/1‑150 (MQ=255)
                                                  tGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACAcc                                                                                                 <  1:353842/150‑1 (MQ=255)
                                                           ggtggtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcggc                                                                                        >  2:341899/1‑150 (MQ=255)
                                                            gtggtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCa                                                                                                            <  1:1769666/129‑1 (MQ=255)
                                                             tggtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcggcgg                                                                                      >  2:1607767/1‑150 (MQ=255)
                                                               gtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCa                                                                                                                       >  2:1986518/1‑115 (MQ=255)
                                                               gtAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCa                                                                                                                       <  1:1986518/115‑1 (MQ=255)
                                                                    cAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTc                                                                               >  1:1532798/1‑150 (MQ=255)
                                                                      gggCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAg                                                                             <  2:1358933/150‑1 (MQ=255)
                                                                       ggccggccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcggc                                                                                        <  1:1312414/138‑1 (MQ=255)
                                                                       ggccggccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcggc                                                                                        >  2:1312414/1‑138 (MQ=255)
                                                                       ggccggccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAgg                                                                            <  2:26760/150‑1 (MQ=255)
                                                                       ggccggccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAgg                                                                            <  1:2210722/150‑1 (MQ=255)
                                                                       ggccggccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGAATGACAGGAAACGGGTCAGCGGCCAATTGCGGACGCCGGACACCCTACTAGTCGGGCTTGATGATGAACACCCCCAACGGCGGCAGGTTCAgg                                                                            <  2:1465163/150‑1 (MQ=255)
                                                                          cggccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggccggccctgcctaccaccgcatccgcaccccaaggctgccaatcaacccgttcaggtcattgtcat                                                                         <  2:2395323‑M2/150‑70 (MQ=255)
                                                                           tggccctttcgcgggagacccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtTTACCCGCGAAGGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGcc                                                                         >  1:896945‑M2/70‑150 (MQ=255)
                                                                            gccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGCTGAACACACCCAACGTCGGCAGGGTCAGGGCCAg                                                                       >  1:748970/1‑150 (MQ=255)
                                                                            gccCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAg                                                                       >  1:489317/1‑150 (MQ=255)
                                                                             ccccTCCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggccagccctgcctaccaccgcatccgcaccccaaggctgccaatcaacccgttcaggtcattgtcatcca                                                                      <  2:221586‑M2/150‑73 (MQ=255)
                                                                             caccTTCGCGGGGAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGc                                                                      >  2:958510/3‑150 (MQ=255)
                                                                               ccTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCa                                                                                                                                        >  2:278373/1‑82 (MQ=255)
                                                                               ccTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCa                                                                                                                                        <  1:278373/82‑1 (MQ=255)
                                                                                cTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGAt                                                                   <  2:489317/150‑1 (MQ=255)
                                                                                 tttcgcgggagacccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtTTACCCGCGAAGGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGAt                                                                   >  2:2350860‑M2/64‑150 (MQ=255)
                                                                                   cgcgGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAAt                                                                <  1:2081391/150‑1 (MQ=255)
                                                                                      gggTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGa                                                                                                         >  1:1736228/1‑106 (MQ=255)
                                                                                      gggTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGa                                                                                                         <  2:1736228/106‑1 (MQ=255)
                                                                                       ggTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCt                                                            <  1:1076141/150‑1 (MQ=255)
                                                                                        gTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTg                                                           >  2:1964073/1‑150 (MQ=255)
                                                                                            ACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCg                                                       <  1:1443417‑M1/150‑1 (MQ=255)
                                                                                             cccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggaacgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGAAGGTTCAGGGCCAGCGATAATGGCTGCCCGt                                                      >  1:1078394‑M2/52‑150 (MQ=255)
                                                                                              ccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGtg                                                     >  1:1100230‑M2/51‑150 (MQ=255)
                                                                                              ccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGtg                                                     >  1:1278386‑M2/51‑150 (MQ=255)
                                                                                              ccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGtg                                                     >  1:1277454‑M2/51‑150 (MQ=255)
                                                                                              ccGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggccggccctgcctaccaccgcatccgcaccccaaggctgccaatcaac                                                                                            >  1:823925‑M2/1‑61 (MQ=255)
                                                                                              ccGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAgggccggccctgcctaccaccgcatccgcaccccaaggctgccaatcaac                                                                                            <  2:823925‑M2/111‑51 (MQ=255)
                                                                                              CCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTg                                                     <  1:2189234‑M1/150‑1 (MQ=255)
                                                                                               cgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGtgg                                                    <  2:1100230‑M2/101‑1 (MQ=255)
                                                                                               CGCTCCTACGCGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACATCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTgg                                                    >  1:749421‑M1/1‑150 (MQ=255)
                                                                                               CGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAGCGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTgg                                                    >  2:1013572‑M1/1‑150 (MQ=255)
                                                                                                GCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGc                                                   <  1:2384904‑M1/150‑1 (MQ=255)
                                                                                                    ctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtTTACCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGtggctgg                                                >  2:1864300‑M2/45‑150 (MQ=255)
                                                                                                    cTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAg                                                                             >  2:1073798/1‑120 (MQ=255)
                                                                                                     tacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGcc                                              <  2:1277454‑M2/107‑1 (MQ=255)
                                                                                                     tacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGcc                                              <  2:1272607‑M2/107‑1 (MQ=255)
                                                                                                     tacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGcc                                              <  2:1278386‑M2/107‑1 (MQ=255)
                                                                                                     tACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGcc                                              >  2:1064188/1‑150 (MQ=255)
                                                                                                          gggTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAAcgg                                                                                         >  2:255380/1‑102 (MQ=255)
                                                                                                          gggTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATgg                                                              <  2:1268656/129‑1 (MQ=255)
                                                                                                          gggTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATgg                                                              >  1:1268656/1‑129 (MQ=255)
                                                                                                          aggcaccgcacccggcggtgacccgcgtaggaacgggtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGa                                                                                                         >  2:1349054‑M2/39‑86 (MQ=255)
                                                                                                             caccgcacccggcggtgacccgtgtaggagcgggtTTACCCGCGAAGGGGCCGGCCCTGCCTACCGATCAGGCGGGCTTCATGATCAACACTCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGt                                       <  1:611615‑M2/115‑1 (MQ=255)
                                                                                                                cgcacccggcggtgacccgtgtaggagcgggtTTACCCGCGAAGGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAgg                                                                            >  1:2284462‑M2/33‑110 (MQ=255)
                                                                                                                cgcacccggcggtgacccgtgtaggagcgggtTTACCCGCGAAGGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAgg                                                                            <  2:2284462‑M2/78‑1 (MQ=255)
                                                                                                                ccgccgGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACt                                   >  1:833199/1‑150 (MQ=255)
                                                                                                                 cgccgGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTg                                  <  2:1597553/150‑1 (MQ=255)
                                                                                                                  cacccggcggtgacccgtgtaggagcgggtTTACCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACCCTCCCCACGGCGACAATTTCATCGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTg                                  <  1:1864300‑M2/120‑1 (MQ=255)
                                                                                                                   ccgGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTCCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTc                                         >  2:1311574/1‑141 (MQ=255)
                                                                                                                   ccgGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTCCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTc                                         <  1:1311574/141‑1 (MQ=255)
                                                                                                                       gtgcggtgcCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCAcc                            >  1:1136665/1‑150 (MQ=255)
                                                                                                                                cTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCgaga                   >  2:2077821/1‑150 (MQ=255)
                                                                                                                                cTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCgaga                   >  2:2269025/1‑150 (MQ=255)
                                                                                                                                  gCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGAtt                 >  2:1358915/1‑150 (MQ=255)
                                                                                                                                    gTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGc               <  1:518828/150‑1 (MQ=255)
                                                                                                                                        gAACGGGTCTCCCGCGAAAgggccggccctgcctaccaccgcatccgcaccccaaggctgccaatcaacccgttcaggtcattgtcatccacatcactgctgtaatcggcgc                                                 >  2:2336547‑M2/1‑19 (MQ=255)
                                                                                                                                        gAACGGGTCTCCCGCGAAAgggccggccctgcctaccaccgcatccgcaccccaaggctgccaatcaacccgttcaggtcattgtcatccacatcactgctgtaatcggcgc                                                 <  1:2336547‑M2/112‑94 (MQ=255)
                                                                                                                                          aCGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAg                                                                             >  2:333207/1‑82 (MQ=255)
                                                                                                                                          aCGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAg                                                                             <  1:333207/82‑1 (MQ=255)
                                                                                                                                              gtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGa     <  1:1064188‑M2/148‑1 (MQ=255)
                                                                                                                                              gtCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGAAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGa     <  2:1078394‑M2/148‑1 (MQ=255)
                                                                                                                                                ctcCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGAcc   <  1:2269025/150‑1 (MQ=255)
                                                                                                                                                 tcCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGCTCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGAccc  <  1:1330143/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CAGTGATGTGGATGACAATGACCTGAACGGGTTGATTGGCAGCCTTGGGGTGCGGATGCGGTGGTAGGCAGGGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCC  >  NC_002947/4577769‑4578063

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