Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,603,823 T→C 7.7% intergenic (‑71/+60) dctB ← / ← bglX C4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,603,8230TC 7.7% 831.0 / 39.0 272intergenic (‑71/+60)dctB/bglXC4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase
Reads supporting (aligned to +/- strand):  ref base T (100/151);  new base C (10/11);  total (110/162)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.70e-01

Only 100 of 287 total aligned reads displayed.

AGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGCAA  >  NC_002947/1603674‑1603972
                                                                                                                                                     |                                                                                                                                                      
aGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTActct                                                                                                                                                        <  2:7475659/150‑1 (MQ=255)
 gCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTg                                                                                                                                                       <  2:4590451/150‑1 (MQ=255)
  ccGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTgc                                                                                                                                                      >  1:5884276/1‑150 (MQ=255)
   cGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTgca                                                                                                                                                     <  1:1552238/150‑1 (MQ=255)
   cGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTgca                                                                                                                                                     <  2:1553212/150‑1 (MQ=255)
   cGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTgca                                                                                                                                                     <  1:4590990/150‑1 (MQ=255)
   cGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTgca                                                                                                                                                     <  2:4404231/150‑1 (MQ=255)
   cGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTgca                                                                                                                                                     <  1:3216532/150‑1 (MQ=255)
       gcTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCAc                                                                                                                                                 <  1:7932268/150‑1 (MQ=255)
       gcTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCAc                                                                                                                                                 <  2:6113557/150‑1 (MQ=255)
       gcTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCAc                                                                                                                                                 <  2:4370811/150‑1 (MQ=255)
       gcTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCAc                                                                                                                                                 <  1:3124885/150‑1 (MQ=255)
         tGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCt                                                                                                                                               <  1:2646298/150‑1 (MQ=255)
             tGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAg                                                                                                                                           <  1:6814728/150‑1 (MQ=255)
             tGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAg                                                                                                                                           <  1:2825451/150‑1 (MQ=255)
             tGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAg                                                                                                                                           <  1:5580653/150‑1 (MQ=255)
             tGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCGGCCTTTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAg                                                                                                                                           <  2:7020854/150‑1 (MQ=37)
               cGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTGCTGCAGAGTACCTGTAGgc                                                                                                                                         <  1:6482940/150‑1 (MQ=38)
                 gTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGgcgc                                                                                                                                       <  2:6741764/150‑1 (MQ=255)
                 gTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGgcgc                                                                                                                                       >  2:4785382/1‑150 (MQ=255)
                   aCCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTTTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCgg                                                                                                                                    <  1:938938/150‑1 (MQ=37)
                    cccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGTAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGt                                                                                                                                   <  2:2427741/150‑1 (MQ=255)
                    cccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCAGCTGTAGGCGCGGt                                                                                                                                   <  2:6617580/150‑1 (MQ=255)
                    cccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGt                                                                                                                                   <  1:739729/150‑1 (MQ=255)
                    cccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGt                                                                                                                                   <  2:7172637/150‑1 (MQ=255)
                     ccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGtt                                                                                                                                  <  2:1170232/150‑1 (MQ=255)
                       aGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGtttt                                                                                                                                <  2:3786756/150‑1 (MQ=255)
                        gCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTa                                                                                                                               >  1:7758020/1‑150 (MQ=255)
                        gCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTa                                                                                                                               <  1:4334959/150‑1 (MQ=255)
                            gATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCcg                                                                                                                           >  2:6358285/1‑150 (MQ=255)
                               cagcagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGa                                                                                                                        >  1:8109835/1‑150 (MQ=255)
                                 gcagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGaaa                                                                                                                      >  2:3446228/1‑150 (MQ=255)
                                  cagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGaaaa                                                                                                                     <  2:2811399/150‑1 (MQ=255)
                                   agGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAg                                                                                                                    >  1:8551944/1‑150 (MQ=255)
                                     ggTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAATGc                                                                                                                  >  2:7696751/1‑150 (MQ=255)
                                        gATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCgg                                                                                                               >  1:3547657/1‑150 (MQ=255)
                                           cAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAAc                                                                                                            >  2:7246921/1‑150 (MQ=255)
                                           cAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCGGCCTTTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCAggccc                                                                                                            <  2:1212312/150‑4 (MQ=17)
                                             gCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCgcctgtt                                                                                                          <  2:2407265/150‑7 (MQ=255)
                                                ccAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcg                                                                                                       >  2:2962323/1‑150 (MQ=255)
                                                 cAGACGCAGGGCGCGAATGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcgg                                                                                                      <  1:293768/150‑1 (MQ=255)
                                                  aGACGCAGGGCGCGAAAGGAAAATGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCAGCCTCTTCGCGGGTAAAACCGCGCCAACAGGTACACTGCAGCACCAGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcggc                                                                                                     >  2:5092566/1‑150 (MQ=39)
                                                      gCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCttt                                                                                                 >  2:7364548/1‑150 (MQ=255)
                                                        aGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTATCCGCGAAAAGGCCGGAACAGGCGGCCttttt                                                                                               <  1:5074120/150‑1 (MQ=255)
                                                        aGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCttttt                                                                                               <  1:1810241/150‑1 (MQ=255)
                                                          ggCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                             <  2:5884276/150‑1 (MQ=255)
                                                          ggCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                             <  2:2319945/150‑1 (MQ=255)
                                                          ggCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                             <  1:3844843/150‑1 (MQ=255)
                                                           gcgcGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAg                                                                                            >  2:250305/1‑150 (MQ=255)
                                                               gAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTc                                                                                              >  2:863210/1‑144 (MQ=255)
                                                                    gAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAAt                                                                                  <  1:7581783/150‑1 (MQ=255)
                                                                         gggCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCaa                                                                              >  1:8609872/1‑150 (MQ=255)
                                                                            cATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGc                                                                           >  2:3668258/1‑150 (MQ=255)
                                                                                gggggCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTctg                                                                       <  2:2156557/150‑1 (MQ=255)
                                                                                gggggCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTctg                                                                       <  1:790313/150‑1 (MQ=255)
                                                                                 ggggCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGGACAGGCGGCCTTTTTCAGTCAAAGCCATTCAAAGCTctgc                                                                      >  1:5404085/1‑150 (MQ=255)
                                                                                 ggggCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTctgc                                                                      <  1:4144005/150‑1 (MQ=255)
                                                                                    gCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTctgctgc                                                                   >  2:8124201/1‑150 (MQ=255)
                                                                                           cGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCg                                                            >  2:1806441/1‑150 (MQ=255)
                                                                                           cGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCg                                                            >  2:6229487/1‑150 (MQ=255)
                                                                                            ggtggtGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGc                                                           <  2:5348915/150‑1 (MQ=255)
                                                                                              tggtggGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCt                                                         <  1:6455831/150‑1 (MQ=255)
                                                                                                gtggGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTgg                                                       <  1:3642526/150‑1 (MQ=255)
                                                                                                  gggATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGa                                                     <  1:1453369/150‑1 (MQ=255)
                                                                                                   ggATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGaa                                                    <  1:3069673/150‑1 (MQ=255)
                                                                                                   ggATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGaa                                                    <  1:3430567/150‑1 (MQ=255)
                                                                                                       gTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCa                                                >  1:3331018/1‑150 (MQ=255)
                                                                                                        tCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAg                                               <  1:1246441/150‑1 (MQ=255)
                                                                                                         cGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGc                                              <  1:6329927/150‑1 (MQ=255)
                                                                                                         cGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGc                                              <  1:7032944/150‑1 (MQ=255)
                                                                                                          tgttccggcctTTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGcc                                             >  1:2302506/8‑150 (MQ=17)
                                                                                                               ctgcctCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTATGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTTCTGCTGCAACGCCTGGGAATCCAGCCCTAcc                                        >  1:6907615/1‑150 (MQ=255)
                                                                                                                   ctctTCGCGGGTAAACCCGCGCCTACAGGTAATCTGCAGCACCGGGAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCATTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCAGGGAATCCAGCCCTACCAGGa                                    >  1:3994042/1‑150 (MQ=38)
                                                                                                                    ttttCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGAc                                   >  1:572433/3‑150 (MQ=21)
                                                                                                                      ttCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACAt                                 >  2:2089916/1‑150 (MQ=255)
                                                                                                                       tCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACAtt                                >  2:4166896/1‑150 (MQ=255)
                                                                                                                       tCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTAAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACAtt                                <  1:5288082/150‑1 (MQ=255)
                                                                                                                         gcgGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGa                              <  1:1239317/150‑1 (MQ=255)
                                                                                                                          cgGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGaa                             <  1:2043825/150‑1 (MQ=255)
                                                                                                                           gggTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAAc                            <  1:1934263/150‑1 (MQ=255)
                                                                                                                              tAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTcc                         >  2:2291938/1‑150 (MQ=255)
                                                                                                                               aaaCCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTccc                        <  2:7613296/150‑1 (MQ=255)
                                                                                                                                   ccGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTgg                    >  1:4482735/1‑150 (MQ=255)
                                                                                                                                    cgcgCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCAGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGt                   >  1:8748360/1‑150 (MQ=25)
                                                                                                                                     gcgcCTACCGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCCAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACTTTGAACTCCCCTGGtt                  <  2:1034895/150‑1 (MQ=255)
                                                                                                                                     gcgcCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGtt                  >  1:4757316/1‑150 (MQ=38)
                                                                                                                                     gcgcCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGtt                  >  1:2361535/1‑150 (MQ=255)
                                                                                                                                     gcgcCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACAGTGAACTCCCCTGGtt                  >  2:7915530/1‑150 (MQ=255)
                                                                                                                                       gcCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTaaa                                                                                        >  1:8158930/1‑79 (MQ=14)
                                                                                                                                       gcCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTaaa                                                                                        <  2:8158930/79‑1 (MQ=14)
                                                                                                                                        ccTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCACGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTACACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCgg               >  2:3066576/1‑150 (MQ=37)
                                                                                                                                         cTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGc              >  1:6262091/1‑150 (MQ=255)
                                                                                                                                         cTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGc              >  2:1903553/1‑150 (MQ=255)
                                                                                                                                            cAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAAc           <  1:4962185/150‑1 (MQ=255)
                                                                                                                                                tACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGc       >  1:7657530/1‑150 (MQ=255)
                                                                                                                                                tACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGc       >  2:3096290/1‑150 (MQ=255)
                                                                                                                                                   tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGc    >  2:689940/1‑150 (MQ=255)
                                                                                                                                                   tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGc    >  2:1390557/1‑150 (MQ=255)
                                                                                                                                                   tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGc    >  2:2262801/1‑150 (MQ=255)
                                                                                                                                                     tGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGCaa  >  2:3967320/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGCAA  >  NC_002947/1603674‑1603972

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