Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,432,366 A→G 8.3% intergenic (‑78/+83) PP_4769 ← / ← PP_4770 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,432,3660AG 8.3% 501.6 / 24.0 168intergenic (‑78/+83)PP_4769/PP_4770hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (68/86);  new base G (10/4);  total (78/90)
Fisher's exact test for biased strand distribution p-value = 9.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

Only 100 of 174 total aligned reads displayed.

CCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTAT  >  NC_002947/5432224‑5432515
                                                                                                                                              |                                                                                                                                                        
ccGTGGGCCTGGCGCTCGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGTTCTTTTTGACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGTGGGGCTGCTGCGCAGACCTTCGCGGTCACGCCCGCTCCCACAGGTATcgcgcg                                                                                                                                                  <  2:343011/150‑1 (MQ=255)
  gTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTa                                                                                                                                                <  1:4633836/150‑1 (MQ=255)
  gTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTa                                                                                                                                                <  1:2059543/150‑1 (MQ=255)
       ccTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACAccc                                                                                                                                           >  2:1721505/1‑150 (MQ=255)
        cTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCg                                                                                                                                          >  2:2899459/1‑150 (MQ=255)
         tGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGTGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGt                                                                                                                                         <  1:1866490/150‑1 (MQ=255)
          ggCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCGCAGGTATCGCGCGTACACCCGTg                                                                                                                                        >  2:1244640/1‑150 (MQ=255)
          ggCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTg                                                                                                                                        <  1:4052752/150‑1 (MQ=255)
          ggCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTg                                                                                                                                        <  1:1113385/150‑1 (MQ=255)
          ggCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTg                                                                                                                                        <  1:2165943/150‑1 (MQ=255)
           gcgcACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTgg                                                                                                                                       <  1:4838770/150‑1 (MQ=255)
           gcgcACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTgg                                                                                                                                       <  1:3906421/150‑1 (MQ=255)
             gcACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACGGGTGTACGCGCGATACCTGTGGGa                                                                                                                                     <  1:658425/150‑1 (MQ=38)
             gcACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGa                                                                                                                                     <  1:1150914/150‑1 (MQ=255)
               aCGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACGGGTGTACGCGCGATACCTGTGGGAGc                                                                                                                                   <  2:1168190/150‑1 (MQ=37)
                 ggTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACGGGTGTACGCGCGATACCTGTGGGAGCgg                                                                                                                                 <  2:4703501/150‑1 (MQ=37)
                 ggTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCgg                                                                                                                                 <  1:2887170/150‑1 (MQ=255)
                 ggTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCgg                                                                                                                                 <  1:357194/150‑1 (MQ=255)
                 ggTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCgg                                                                                                                                 <  1:1258062/150‑1 (MQ=255)
                 ggTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCgg                                                                                                                                 <  1:4043176/150‑1 (MQ=255)
                        tagAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCggg                                                                                                                                >  2:4001904/1‑144 (MQ=255)
                        tagAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGcc                                                                                                                          >  2:3995941/1‑150 (MQ=255)
                        tagAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGcc                                                                                                                          >  1:1346646/1‑150 (MQ=255)
                         agagGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGTACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGGGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGGGGGCGTGccc                                                                                                                         <  2:622613/150‑1 (MQ=255)
                          gagGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCCCAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCcg                                                                                                                        >  1:622121/1‑150 (MQ=255)
                          gagGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGcccc                                                                                                                        >  1:622189/1‑149 (MQ=255)
                              ttCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGaa                                                                                                                    <  1:2359731/150‑1 (MQ=255)
                                cTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAgg                                                                                                                  <  2:1831505/150‑1 (MQ=255)
                                  gCACATCGGCGACCAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAgg‑‑gg                                                                                                              >  1:1382961/1‑150 (MQ=255)
                                   cacaTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑Ggc                                                                                                             <  1:777834/150‑1 (MQ=255)
                                   cacaTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑Ggc                                                                                                             <  1:2209382/150‑1 (MQ=255)
                                      aTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGcaca                                                                                                          >  1:2084883/1‑150 (MQ=255)
                                      aTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGACGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGcaca                                                                                                          >  1:2084833/1‑150 (MQ=255)
                                        cGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCaaa                                                                                                        <  1:5011852/150‑1 (MQ=255)
                                        cGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCaaa                                                                                                        <  1:4838406/150‑1 (MQ=255)
                                          gcgcACAGCGCGATCTTTTTCACCTGGGAACCTCATGTTTGGCCGGGAGTGCCGGGGGGCTGCGGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCagcag                                                                                                      <  2:1119086/150‑5 (MQ=39)
                                          gCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTGTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCACAgc                                                                                                      >  2:5011205/1‑150 (MQ=255)
                                            gAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACCCCCGCTCCCACAGGTAACGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAgcgc                                                                                                    >  1:580253/1‑150 (MQ=255)
                                               cAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCGGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGcccc                                                                                                 <  1:976137/150‑1 (MQ=255)
                                               cAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGcccc                                                                                                 <  2:197457/150‑1 (MQ=255)
                                                 gcgcGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCg                                                                                               <  2:1640020/150‑1 (MQ=255)
                                                 gcgcGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCg                                                                                               <  2:531786/150‑1 (MQ=255)
                                                  cgcgATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGGCACGCCCGCTCCCACAGGTATCGTGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCgg                                                                                              >  2:3305477/1‑150 (MQ=255)
                                                  cgcgATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCgg                                                                                              >  2:899684/1‑150 (MQ=255)
                                                     gATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGaaa                                                                                           <  1:747703/150‑1 (MQ=255)
                                                     gATCTTTCTCACCTGGGAACCTCACTTTTGGCCGGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGaaa                                                                                           <  2:1192976/150‑1 (MQ=255)
                                                      aTCTTTTTCACCTGGGAACCTCATTTTGGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAAt                                                                                          <  1:446027/150‑1 (MQ=255)
                                                      aTCTTTTGCCCCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGGGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAAt                                                                                          <  1:446028/150‑1 (MQ=255)
                                                         tttttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGcc                                                                                       <  1:218627/150‑1 (MQ=255)
                                                         tttttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGcc                                                                                       >  1:3379393/1‑150 (MQ=255)
                                                          ttttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCca                                                                                      <  2:862265/150‑1 (MQ=255)
                                                          ttttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCca                                                                                      <  1:2072437/150‑1 (MQ=255)
                                                          ttttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCca                                                                                      <  2:4429395/150‑1 (MQ=255)
                                                                   gggAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGc                                                                             >  1:3236201/1‑150 (MQ=255)
                                                                     gAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCa                                                                           <  1:5116640/150‑1 (MQ=255)
                                                                     gAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCa                                                                           <  2:4558323/150‑1 (MQ=255)
                                                                      aaCCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAg                                                                          >  2:254989/1‑150 (MQ=255)
                                                                        ccTCATTCTGGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAc                                                                        <  1:451774/150‑1 (MQ=255)
                                                                         cTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAcc                                                                       >  2:4221104/1‑150 (MQ=255)
                                                                         cTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAcc                                                                       >  1:3238661/1‑150 (MQ=255)
                                                                         cTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAcc                                                                       >  1:2694855/1‑150 (MQ=255)
                                                                             tttttGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTc                                                                   <  2:1567160/150‑1 (MQ=255)
                                                                             tttttGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTc                                                                   <  2:1368898/150‑1 (MQ=255)
                                                                              ttttGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCa                                                                  >  2:2127834/1‑150 (MQ=255)
                                                                                 tGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGAc                                                               <  2:3836706/150‑1 (MQ=255)
                                                                                 tGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGAc                                                               >  1:383306/1‑150 (MQ=255)
                                                                                     cTGGTGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACCCCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTaa                                                           <  1:2609381/150‑1 (MQ=255)
                                                                                     cTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTaa                                                           <  2:4312982/150‑1 (MQ=255)
                                                                                     cTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTaa                                                           <  1:379638/150‑1 (MQ=255)
                                                                                          gTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCTTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtc                                                      <  1:3462158/150‑1 (MQ=255)
                                                                                            gCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACTGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCGCGGAAATGCCTGCAGTATGCGAGACCTGTCAGACGTAAGAAtcgt                                                    <  2:3963285/150‑1 (MQ=255)
                                                                                             ccGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTTCACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcgtc                                                   <  2:2737797/150‑1 (MQ=255)
                                                                                             ccGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcgtc                                                   <  2:276069/150‑1 (MQ=255)
                                                                                             ccGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcgtc                                                   <  2:3235904/150‑1 (MQ=255)
                                                                                             ccGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcgtc                                                   <  2:3236201/150‑1 (MQ=255)
                                                                                             ccGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcgtc                                                   <  2:44424/150‑1 (MQ=255)
                                                                                                gggggCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCgtcgtc                                                <  1:4339974/150‑1 (MQ=255)
                                                                                                 ggggCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTcc                                               >  2:2204405/1‑150 (MQ=255)
                                                                                                 ggggCTGCTGCGCAGACCTTCGCGGACACGCCCGCTACCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTcc                                               >  2:1698608/1‑150 (MQ=255)
                                                                                                              aGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCAc                                  <  1:4390224/150‑1 (MQ=255)
                                                                                                              aGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCAc                                  <  1:2389196/150‑1 (MQ=255)
                                                                                                                  cTTCGCGGACACGCCCGCTCCCACAGGTATAGCGCGTACACCCCTGGGATCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCATCAGTATGCCAGACCAGTCAGAAGTAAGAATAGACTTCCGAGAAGAACCCACCGTc                              >  1:2927311/1‑150 (MQ=255)
                                                                                                                     cgcgGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATc                           <  1:3796472/150‑1 (MQ=255)
                                                                                                                      gcgGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCg                          >  2:3737251/1‑150 (MQ=255)
                                                                                                                        ggACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCt                        >  1:1525189/1‑150 (MQ=255)
                                                                                                                        ggACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCt                        >  1:2416935/1‑150 (MQ=255)
                                                                                                                              gCCCGCTCCCACGGGTGTACGCGCGATACCTGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCaaa                  >  2:4052354/1‑150 (MQ=37)
                                                                                                                              gCCCGCTCCCACGGGTGTACGCGCGATACCTGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCaaa                  >  2:2812209/1‑150 (MQ=37)
                                                                                                                              gCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCaaa                  >  2:3709849/1‑150 (MQ=255)
                                                                                                                              gCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCaaa                  >  2:2400662/1‑150 (MQ=255)
                                                                                                                                ccGCTCCCACGGGTGTACGCGCGATACCTGTGGGAGCGGGCGTGTCCGCGAAGGTCTGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCaaa                  >  2:4527398/1‑150 (MQ=21)
                                                                                                                                  gCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAacca              >  1:3692574/1‑150 (MQ=255)
                                                                                                                                       cacaGGTATCGCGCGTACACCCATGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTc         >  2:4326827/1‑150 (MQ=255)
                                                                                                                                       cacaGGTATCGCGCGTACACCCATGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTc         >  2:1851070/1‑150 (MQ=255)
                                                                                                                                        acaGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCg        >  1:1642516/1‑150 (MQ=255)
                                                                                                                                        acaGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCg        >  1:1364102/1‑150 (MQ=255)
                                                                                                                                              aTCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTAt  >  2:2663200/1‑150 (MQ=255)
                                                                                                                                              aTCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTAt  >  2:149019/1‑150 (MQ=255)
                                                                                                                                              aTCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTAt  >  1:3559095/1‑150 (MQ=255)
                                                                                                                                              aTCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCACCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTTAGATTCGTCGTCCGAGAAGAACCCAGAGTCATCGCTGCGAAAACCACCGTCTGCGTAt  >  2:2662985/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                        
CCGTGGGCCTGGCGCACGGTGTAGTAGAGGTTCTGCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTATCGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGG‑‑GGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTAT  >  NC_002947/5432224‑5432515

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