Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,322,454 G→C 5.3% A251G (GCC→GGC)  PP_2041 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,322,4540GC 5.3% 629.6 / 18.3 188A251G (GCC→GGC) PP_2041LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (87/91);  new base C (4/6);  total (91/97)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 188 total aligned reads displayed.

TTTCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCTCGATGC  >  NC_002947/2322306‑2322601
                                                                                                                                                    |                                                                                                                                                    
tttCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGc                                                                                                                                                     >  2:5257836/1‑150 (MQ=255)
tttCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGc                                                                                                                                                     >  2:5023175/1‑150 (MQ=255)
tttCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGc                                                                                                                                                     >  1:3407826/1‑150 (MQ=255)
tttCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGc                                                                                                                                                     >  1:5280245/1‑150 (MQ=255)
tttCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGc                                                                                                                                                     <  2:1997005/150‑1 (MQ=255)
      gcgcTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGgcgc                                                                                                                                              <  2:1823079/150‑1 (MQ=255)
           cgctGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCa                                                                                                                                         <  1:5172752/150‑1 (MQ=255)
              tGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGAt                                                                                                                                      >  1:4602278/1‑150 (MQ=255)
               gcgcACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATc                                                                                                                                     <  2:3992483/150‑1 (MQ=255)
                  cACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCcag                                                                                                                                  >  2:4176719/1‑150 (MQ=255)
                   aCCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCcagc                                                                                                                                 <  2:3215888/150‑1 (MQ=255)
                    ccGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCcagca                                                                                                                                >  1:4681374/1‑150 (MQ=255)
                     cGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCcagcag                                                                                                                               >  1:2061143/1‑150 (MQ=255)
                      gCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCagcagc                                                                                                                              <  2:3502790/150‑1 (MQ=255)
                      gCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCagcagc                                                                                                                              >  2:1598311/1‑150 (MQ=255)
                            accaGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCggcggc                                                                                                                        <  1:3874108/150‑1 (MQ=255)
                             ccaGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGcc                                                                                                                       >  1:4453/1‑150 (MQ=255)
                              cagcagTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCa                                                                                                                      >  2:1118589/1‑150 (MQ=255)
                                gcagTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGt                                                                                                                    >  2:4645865/1‑150 (MQ=255)
                                 cagTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGtt                                                                                                                   <  2:2619492/150‑1 (MQ=255)
                                   gTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGc                                                                                                                 >  2:1704847/1‑150 (MQ=255)
                                      cgcgcACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGa                                                                                                              <  2:2090833/150‑1 (MQ=255)
                                        cgcACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGAtg                                                                                                            <  2:2175285/150‑1 (MQ=255)
                                        cgcACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCAGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGAtg                                                                                                            >  2:1646347/1‑150 (MQ=255)
                                         gcACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGAtgg                                                                                                           >  1:4701797/1‑150 (MQ=255)
                                          cACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTGGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGAtggt                                                                                                          >  2:2954435/1‑150 (MQ=255)
                                             cAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGgc                                                                                                       >  1:1152603/1‑150 (MQ=255)
                                             cAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGgc                                                                                                       <  2:1245058/150‑1 (MQ=255)
                                             cAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGgc                                                                                                       <  2:1245057/150‑1 (MQ=255)
                                             cAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGgc                                                                                                       <  2:4784226/150‑1 (MQ=255)
                                                  cgcgCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCgcca                                                                                                  <  2:3406099/150‑1 (MQ=255)
                                                   gcgcATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCgccag                                                                                                 >  2:2468697/1‑150 (MQ=255)
                                                     gcATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCCGCCGCTGCTGGATCTGGCCGCGCCTGCCAgc                                                                                                                    >  1:129566/1‑128 (MQ=37)
                                                      cATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCgccaggcc                                                                                              >  2:3686728/1‑150 (MQ=255)
                                                      cATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCgccaggcc                                                                                              >  2:4957026/1‑150 (MQ=255)
                                                          cAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCa                                                                                          <  1:4309516/150‑1 (MQ=255)
                                                          cAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCa                                                                                          >  2:919615/1‑150 (MQ=255)
                                                              cggcggcgAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGa                                                                                      <  1:5257836/150‑1 (MQ=255)
                                                               ggcggcgAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTTCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGAc                                                                                     <  1:4432746/150‑1 (MQ=255)
                                                                 cggcgAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCCGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCa                                                                                   <  2:4955589/150‑1 (MQ=255)
                                                                   gcgAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAgg                                                                                 <  1:2529368/150‑1 (MQ=255)
                                                                     gAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCCTGCCAGCGCGCGGctga                                                                                                                                        <  1:3221635/93‑4 (MQ=38)
                                                                     gAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCCTGCCAGCGCGCGGctga                                                                                                                                        >  2:3221535/1‑90 (MQ=38)
                                                                      aaCGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTg                                                                              <  1:3197022/150‑1 (MQ=255)
                                                                      aaCGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTg                                                                              >  2:2984881/1‑150 (MQ=255)
                                                                        cGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGt                                                                            <  2:3892832/150‑1 (MQ=255)
                                                                           cAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGAt                                                                         >  2:730272/1‑150 (MQ=255)
                                                                             gCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATAc                                                                       >  2:965665/1‑150 (MQ=255)
                                                                               ggCTTGTGCTGCGAGCGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCg                                                                     <  1:3060018/150‑1 (MQ=255)
                                                                               ggCTTGTGCTGCGAGCGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCg                                                                     <  1:329854/150‑1 (MQ=255)
                                                                                   tgtgCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCggggg                                                                 >  1:2624010/1‑150 (MQ=255)
                                                                                    gtgCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGc                                                                <  2:3309596/150‑1 (MQ=255)
                                                                                     tgctgcGAGTGGGCCAAGGCGATGGTGCGGGTGAGGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGcc                                                               <  2:4409440/150‑1 (MQ=255)
                                                                                       ctgcGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGa                                                             >  1:3945251/1‑150 (MQ=255)
                                                                                        tgcGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGAc                                                            >  2:119180/1‑150 (MQ=255)
                                                                                        tgcGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGAc                                                            >  2:4607794/1‑150 (MQ=255)
                                                                                        tgcGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGAc                                                            >  2:4292382/1‑150 (MQ=255)
                                                                                        tgcGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCCTGCCAGCGCGCGGc                                                                                                                                           <  2:3075311/71‑1 (MQ=37)
                                                                                          cGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGActg                                                          <  2:383227/150‑1 (MQ=255)
                                                                                           gAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGActgc                                                         <  1:4368196/150‑1 (MQ=255)
                                                                                            aGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGActgct                                                        <  1:1704847/150‑1 (MQ=255)
                                                                                              tGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGt                                                      <  2:3638245/150‑1 (MQ=255)
                                                                                              tGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGt                                                      >  1:1780807/1‑150 (MQ=255)
                                                                                                ggCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTc                                                    <  1:4059460/150‑1 (MQ=255)
                                                                                                 gCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCg                                                   >  2:4846281/1‑150 (MQ=255)
                                                                                                  ccAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGt                                                  >  2:641642/1‑150 (MQ=255)
                                                                                                    aaGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAg                                                >  1:2843983/1‑150 (MQ=255)
                                                                                                     aGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGc                                               <  1:1446698/150‑1 (MQ=255)
                                                                                                       gCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAg                                                                                                                                                   <  2:2337777/49‑1 (MQ=255)
                                                                                                       gCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCa                                                                                          <  1:2282769/105‑1 (MQ=255)
                                                                                                       gCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCggg                                                                   <  2:266372/128‑1 (MQ=255)
                                                                                                       gCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCggg                                                                   >  1:266372/1‑128 (MQ=255)
                                                                                                       gCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTg                                             <  1:2343460/150‑1 (MQ=255)
                                                                                                         gATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGa                                           <  1:740670/150‑1 (MQ=255)
                                                                                                            ggTGCGGGTGATGCGCGGTTCCCCTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTc                                        <  2:3224044/150‑1 (MQ=255)
                                                                                                              tGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCa                                      <  1:1231557/150‑1 (MQ=255)
                                                                                                                cGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAgt                                    <  2:51333/150‑1 (MQ=255)
                                                                                                                cGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAgt                                    <  1:1242223/150‑1 (MQ=255)
                                                                                                                    tGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCca                                                                                                                                   <  1:2992364/51‑1 (MQ=255)
                                                                                                                    tGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCca                                                                                                                                   >  2:2992364/1‑51 (MQ=255)
                                                                                                                              ttCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTc                      <  2:3865253/150‑1 (MQ=255)
                                                                                                                                ccACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCa                    <  2:1415240/150‑1 (MQ=255)
                                                                                                                                 cACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTc                                        >  2:216707/1‑129 (MQ=255)
                                                                                                                                 cACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCAt                   >  1:5412618/1‑150 (MQ=255)
                                                                                                                                 cACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCAt                   >  1:4506528/1‑150 (MQ=255)
                                                                                                                                 cACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCAt                   <  1:3734842/150‑1 (MQ=255)
                                                                                                                                  aCTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATg                  <  2:4099197/150‑1 (MQ=255)
                                                                                                                                    tATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGt                                                      >  2:3431695/1‑112 (MQ=255)
                                                                                                                                    tATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGt                                                      <  1:3431695/112‑1 (MQ=255)
                                                                                                                                     aTCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCa               <  1:2149098/150‑1 (MQ=255)
                                                                                                                                      tCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAg              >  1:5199123/1‑150 (MQ=255)
                                                                                                                                       tggCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTGTCCATGGCAGc             >  1:4952443/3‑150 (MQ=255)
                                                                                                                                       cAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGc             >  1:1129803/1‑150 (MQ=255)
                                                                                                                                         gCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCt           >  1:1303329/1‑150 (MQ=255)
                                                                                                                                         gCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGGGTCACACCGGTTTCCATGGCAGCCt           >  1:1848095/1‑150 (MQ=255)
                                                                                                                                             cgcgcTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCATCCTGCTc       <  1:5023175/150‑1 (MQ=255)
                                                                                                                                               cgcTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCTCGa     <  1:641642/150‑1 (MQ=255)
                                                                                                                                                gcTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCTCGAt    <  1:1934054/150‑1 (MQ=255)
                                                                                                                                                  tGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCTCGATgc  <  1:3109559/150‑1 (MQ=255)
                                                                                                                                                  tGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCTCGATgc  >  2:2268504/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTCCAGCGCTCGCTGCGCACCGCATCCACCAGCAGTGCGCGCACCAGGTCGCGCATGCAGGCGGCGGCGAACGACAGCGGCTTGTGCTGCGAGTGGGCCAAGGCGATGGTGCGGGTGATGCGCGGTTCCACTATCAGCCGCGCGCTGGCAGGCGCGGCCAGATCCAGCAGCGGCGGCCAGTTGCTGATGGTGGCCGCCAGGCCAGCACGGACCAGGCTGGTGATACCGGGGGCCGACTGCTGTTCGTAGCTGGAGTCCAGTGTCACACCGGTTTCCATGGCAGCCTGCTCGATGC  >  NC_002947/2322306‑2322601

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: