Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,598,873 C→G 11.2% intergenic (‑70/+61) 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,598,8730CG11.2% 327.3 / 27.2 107intergenic (‑70/+61)dctB/bglXC4‑dicarboxylate transport sensor protein/beta‑D‑glucoside glucohydrolase
Reads supporting (aligned to +/- strand):  ref base C (37/58);  new base G (5/7);  total (42/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 113 total aligned reads displayed.

CAGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGC  >  NC_002947/1598724‑1599021
                                                                                                                                                     |                                                                                                                                                    
cAGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTActc                                                                                                                                                      <  2:5095171/150‑1 (MQ=255)
 aGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTActct                                                                                                                                                     <  2:3509744/150‑1 (MQ=255)
 aGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTActct                                                                                                                                                     <  2:2717601/150‑1 (MQ=255)
  gCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTg                                                                                                                                                    <  1:122857/150‑1 (MQ=255)
   ccGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTgc                                                                                                                                                   <  2:1703297/150‑1 (MQ=255)
   ccGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTgc                                                                                                                                                   <  1:3261387/150‑1 (MQ=255)
    cGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTgca                                                                                                                                                  <  1:472905/150‑1 (MQ=255)
        gcTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCAc                                                                                                                                              <  1:151569/150‑1 (MQ=255)
          tGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCt                                                                                                                                            >  1:2167724/1‑150 (MQ=255)
             gTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTa                                                                                                                                         <  1:4243774/150‑1 (MQ=255)
             gTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTa                                                                                                                                         <  2:626846/150‑1 (MQ=255)
             gTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTa                                                                                                                                         <  1:1125040/150‑1 (MQ=255)
              tGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTTTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAg                                                                                                                                        <  1:2742745/150‑1 (MQ=37)
              tGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAg                                                                                                                                        <  1:1190479/150‑1 (MQ=255)
              tGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGAGCCTGCCTCTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAg                                                                                                                                        <  1:1649745/150‑1 (MQ=37)
                    aCCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCgg                                                                                                                                  <  1:1483013/150‑1 (MQ=38)
                      ccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGTAAGGAATAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCGGCCTTTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCgggt                                                                                                                                <  2:3885164/150‑3 (MQ=21)
                      ccAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCGGCCTTTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCgggt                                                                                                                                <  2:3884704/150‑3 (MQ=25)
                       cAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGttt                                                                                                                               <  2:1701464/150‑1 (MQ=255)
                         gCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTa                                                                                                                             <  1:1802363/150‑1 (MQ=255)
                          cACGATCAGCAGGGTGCTCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTAc                                                                                                                            <  1:3082733/150‑1 (MQ=255)
                          cACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTAc                                                                                                                            <  2:4092991/150‑1 (MQ=25)
                                 agcagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTTTTCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGaa                                                                                                                     <  2:4629919/150‑1 (MQ=21)
                                   cagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGaaaa                                                                                                                   <  1:5001472/150‑1 (MQ=255)
                                   cagGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGaaaa                                                                                                                   <  1:4781581/150‑1 (MQ=255)
                                      ggTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCTGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGc                                                                                                                >  2:4655271/1‑150 (MQ=255)
                                      ggTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGc                                                                                                                <  2:3354050/150‑1 (MQ=255)
                                       gTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCCGTAGGCGCGGTTTTACCCGCGAAAAGGcc                                                                                                               <  2:3036636/150‑1 (MQ=255)
                                        tGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCg                                                                                                              <  2:2701949/150‑1 (MQ=255)
                                        tGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCg                                                                                                              <  2:4434732/150‑1 (MQ=255)
                                        tGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCg                                                                                                              <  2:3100853/150‑1 (MQ=255)
                                        tGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCg                                                                                                              <  2:1738154/150‑1 (MQ=255)
                                        tGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCg                                                                                                              >  1:4856438/1‑150 (MQ=255)
                                                cccAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggc                                                                                                      <  1:2083840/150‑1 (MQ=255)
                                                cccAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggc                                                                                                      <  2:803789/150‑1 (MQ=255)
                                                 ccAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcg                                                                                                     >  1:100756/1‑150 (MQ=255)
                                                   aGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAggcggc                                                                                                   <  2:2028636/150‑1 (MQ=255)
                                                    gACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGcc                                                                                                  >  1:2162602/1‑150 (MQ=255)
                                                        cAGGGCGCGAAAGGAAATGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCtttt                                                                                              <  1:1784231/150‑1 (MQ=255)
                                                         aGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCttttt                                                                                             <  2:4326218/150‑1 (MQ=255)
                                                         aGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCttttt                                                                                             >  2:1282819/1‑150 (MQ=255)
                                                          gggCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTc                                                                                            <  2:865840/150‑1 (MQ=255)
                                                          gggCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTc                                                                                            <  2:865858/150‑1 (MQ=255)
                                                           ggCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                           <  1:588979/150‑1 (MQ=255)
                                                           ggCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCa                                                                                           <  2:1163490/150‑1 (MQ=255)
                                                                gAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaa                                                                                      >  1:1124014/1‑150 (MQ=255)
                                                                 aaaGGCAACGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                     <  1:4829093/150‑1 (MQ=255)
                                                                 aaaGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                     <  2:3179376/150‑1 (MQ=255)
                                                                 aaaGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                     <  1:3117084/150‑1 (MQ=255)
                                                                 aaaGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                     >  2:5018764/1‑150 (MQ=255)
                                                                 aaaGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                     <  2:1322725/150‑1 (MQ=255)
                                                                 aaaGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                     <  2:1252197/150‑1 (MQ=255)
                                                                 aaaGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCaaa                                                                                     <  1:934592/150‑1 (MQ=255)
                                                                   aGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGc                                                                                   >  2:2456673/1‑150 (MQ=255)
                                                                    ggAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCa                                                                                  >  1:1634732/1‑150 (MQ=255)
                                                                     gAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCaa                                                                                 >  1:4336658/1‑150 (MQ=255)
                                                                      aaaaGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAAt                                                                                <  2:1538331/150‑1 (MQ=255)
                                                                      aaaaGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAAt                                                                                <  2:2334178/150‑1 (MQ=255)
                                                                          gggCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCaa                                                                            >  2:3201363/1‑150 (MQ=255)
                                                                           ggCATCGGGGGCGGTTCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTATACGCGCGCCTACAGGTACTCTGCAGCACCTGTATGCGCGTTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCaaa                                                                           <  2:1239375/150‑1 (MQ=255)
                                                                           ggCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCaaa                                                                           <  1:1282819/150‑1 (MQ=255)
                                                                               tCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAAGTCAAAGctc                                                                       <  2:3112361/150‑1 (MQ=255)
                                                                                     gCGGATCCGGTGGTGGGATGTCAGGGCCTGCCTCTACGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTctgctgc                                                                 <  1:3448036/150‑1 (MQ=255)
                                                                                         aTCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCa                                                             <  1:2164513/150‑1 (MQ=255)
                                                                                                   gggATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCCAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGa                                                   >  2:1671574/1‑150 (MQ=255)
                                                                                                   gggATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGa                                                   <  1:4965557/150‑1 (MQ=255)
                                                                                                   gggATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGa                                                   >  2:1671807/1‑150 (MQ=255)
                                                                                                    ggATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGaa                                                  >  2:4525927/1‑150 (MQ=255)
                                                                                                     gATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAAt                                                 <  2:3966047/150‑1 (MQ=255)
                                                                                                     gATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAAt                                                 <  1:4152066/150‑1 (MQ=255)
                                                                                                                  gcctCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTg                                    >  2:2126944/1‑150 (MQ=255)
                                                                                                                    ctctTCGCGGGTAAACCCGCGCCTACCGGTACTCTGCAGCACCTGTAGGCGCGGTATTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGa                                  <  1:4766717/150‑1 (MQ=255)
                                                                                                                    ctctTCGCGGGTAAACCCGCGCCTACCGGGACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGa                                  <  1:3255465/150‑1 (MQ=255)
                                                                                                                    ctctTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGa                                  <  1:4766903/150‑1 (MQ=255)
                                                                                                                     ttttCGCGGGTAAAACCGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGAc                                 >  2:4511842/3‑150 (MQ=21)
                                                                                                                        tCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACAtt                              <  1:4675276/150‑1 (MQ=255)
                                                                                                                          gcgGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGa                            >  2:1570365/1‑150 (MQ=255)
                                                                                                                            gggTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAAc                          <  1:530416/150‑1 (MQ=255)
                                                                                                                            gggTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAAc                          >  2:668537/1‑150 (MQ=255)
                                                                                                                                  aaccGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCt                    >  1:917812/3‑150 (MQ=25)
                                                                                                                                  aaccGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCt                    >  1:3021463/3‑150 (MQ=25)
                                                                                                                                  aaccGCGCCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCAGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCt                    >  2:4719767/3‑150 (MQ=25)
                                                                                                                                  aCCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCt                    >  2:4700923/1‑150 (MQ=255)
                                                                                                                                        gcCTACAGGTGCTGCAGAGTACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCg              >  1:4973629/1‑150 (MQ=38)
                                                                                                                                          cTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGc            >  1:727158/1‑150 (MQ=255)
                                                                                                                                          cTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGc            >  1:52654/1‑150 (MQ=255)
                                                                                                                                           tACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCa           >  1:595710/1‑150 (MQ=255)
                                                                                                                                           tACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCac          >  1:2488652/1‑149 (MQ=255)
                                                                                                                                            aCAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCaa          >  2:4748961/1‑150 (MQ=255)
                                                                                                                                              aGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAAcc        >  1:4619390/1‑150 (MQ=255)
                                                                                                                                              aGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAAcc        >  2:1816017/1‑150 (MQ=255)
                                                                                                                                                gTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCg      >  2:2819820/1‑150 (MQ=255)
                                                                                                                                                gTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCg      >  2:4005459/1‑150 (MQ=255)
                                                                                                                                                gTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACAGTGAACTCCCCTGGTTCGGCAACCCg      >  2:4005327/1‑150 (MQ=255)
                                                                                                                                                  aCTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCt    <  2:2621160/150‑1 (MQ=255)
                                                                                                                                                  aCTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCt    <  2:1441120/150‑1 (MQ=255)
                                                                                                                                                  aCTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCt    <  2:1441118/150‑1 (MQ=255)
                                                                                                                                                    tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGc  >  1:4465057/1‑150 (MQ=255)
                                                                                                                                                    tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGc  >  1:629737/1‑150 (MQ=255)
                                                                                                                                                    tctGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGc  >  2:634479/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGCCGGCGCTGAGTGCGGTACCCAGCACGATCAGCAGGGTGATCAGCCCCAGACGCAGGGCGCGAAAGGAAAAGGGCATCGGGGGCGGATCCGGTGGTGGGATGTCGGGGCCTGCCTCTTCGCGGGTAAACCCGCGCCTACAGGTACTCTGCAGCACCTGTAGGCGCGGTTTTACCCGCGAAAAGGCCGGAACAGGCGGCCTTTTTCAGTCAAAGCAATTCAAAGCTCTGCTGCTGCACCGCCTGGGAATCCAGCCCTACCTGGACATTGAACTCCCCTGGTTCGGCAACCCGCTGC  >  NC_002947/1598724‑1599021

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