Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 372,660 C→G 5.0% intergenic (+77/‑102) dgcA → / → dgcB dimethylglycine dehydrogenase subunit/dimethylglycine dehydrogenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947372,6600CG 5.0% 1009.4 / 30.4 297intergenic (+77/‑102)dgcA/dgcBdimethylglycine dehydrogenase subunit/dimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base C (131/150);  new base G (10/5);  total (142/155)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 308 total aligned reads displayed.

GGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCTG  >  NC_002947/372512‑372809
                                                                                                                                                    |                                                                                                                                                     
 gCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTg                                                                                                                                                     <  2:5850357/150‑1 (MQ=255)
   gACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGtt                                                                                                                                                   <  2:1351929/150‑1 (MQ=255)
   gACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGtt                                                                                                                                                   >  2:5435014/1‑150 (MQ=255)
    aCTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTg                                                                                                                                                  <  1:6039833/150‑1 (MQ=255)
     cTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGc                                                                                                                                                 <  1:3859237/150‑1 (MQ=255)
      tGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGcc                                                                                                                                                <  2:675634/150‑1 (MQ=255)
      tGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGcc                                                                                                                                                <  1:17337/150‑1 (MQ=255)
      tGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGcc                                                                                                                                                >  2:6433876/1‑150 (MQ=255)
               cAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                                       <  1:1564292/150‑1 (MQ=255)
                 gcgcAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGc                                                                                                                                     <  2:5322609/150‑1 (MQ=255)
                    cAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCAtt                                                                                                                                  <  2:2475947/150‑1 (MQ=255)
                     aaCGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTc                                                                                                                                 <  2:1567/150‑1 (MQ=255)
                       cGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTcgc                                                                                                                               <  2:6371761/150‑1 (MQ=255)
                       cGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTcgc                                                                                                                               <  2:2145852/150‑1 (MQ=255)
                        gTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTcgcg                                                                                                                              <  1:6497384/150‑1 (MQ=255)
                        gTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTcgcg                                                                                                                              >  1:2397472/1‑150 (MQ=255)
                        gTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTcgcg                                                                                                                              <  2:4359475/150‑1 (MQ=255)
                         tGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg                                                                                                                             <  1:6771087/150‑1 (MQ=255)
                         tGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg                                                                                                                             <  2:6553984/150‑1 (MQ=255)
                          gCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCggg                                                                                                                            <  1:747330/150‑1 (MQ=21)
                          gCATGCGGCGAACTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCggg                                                                                                                            >  1:219646/1‑150 (MQ=255)
                                ggcgATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAc                                                                                                                                            >  1:6693984/1‑128 (MQ=255)
                                     tCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGctc                                                                                                                 >  1:2944927/1‑150 (MQ=255)
                                      cTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGctct                                                                                                                <  1:4180693/150‑1 (MQ=255)
                                         cgacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTAc                                                                                                             >  1:3964540/1‑150 (MQ=255)
                                          gacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACa                                                                                                            >  2:2255479/1‑150 (MQ=255)
                                          gacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACa                                                                                                            >  1:2591714/1‑150 (MQ=255)
                                          gacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCAGTCGCGGGTGAACCCGCTCTTACa                                                                                                            >  2:3114897/1‑150 (MQ=255)
                                            cgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCggg                                                                                                                            <  2:6718616/132‑1 (MQ=255)
                                             gCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATgg                                                                                                         >  1:3280810/1‑150 (MQ=255)
                                               cTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATgggg                                                                                                       <  2:2643044/150‑1 (MQ=255)
                                                  gcgcCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCa                                                                                                    <  1:2255286/150‑1 (MQ=255)
                                                        gtgCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCg                                                                                              <  1:6200111/150‑1 (MQ=255)
                                                            aaGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGGCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCct                                                                                          >  1:1270963/1‑150 (MQ=255)
                                                             aGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCctc                                                                                         >  2:5162673/1‑150 (MQ=255)
                                                              ggACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCctct                                                                                        <  1:5279132/150‑1 (MQ=255)
                                                               gACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTg                                                                                       <  2:3371265/150‑1 (MQ=255)
                                                                  ttCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTc                                                                                    >  1:4732017/1‑150 (MQ=255)
                                                                       aTCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGgcgc                                                                               >  2:937369/1‑150 (MQ=255)
                                                                        tCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCa                                                                              >  2:4020579/1‑150 (MQ=255)
                                                                         cGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTggggtcaa                                                                             <  2:3863272/150‑6 (MQ=255)
                                                                         cGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCaa                                                                             <  2:818761/150‑1 (MQ=255)
                                                                             aaCAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAAt                                                                         <  1:3739202/150‑1 (MQ=255)
                                                                                aGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGATCCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCa                                                                      <  1:5418854/150‑1 (MQ=255)
                                                                                aGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGACTCCa                                                                      >  2:6613902/1‑150 (MQ=255)
                                                                                aGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCa                                                                      <  2:2796441/150‑1 (MQ=255)
                                                                                aGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCa                                                                      <  1:5211206/150‑1 (MQ=255)
                                                                                aGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCa                                                                      >  1:5083416/1‑150 (MQ=255)
                                                                                 gTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAg                                                                     <  1:5598489/150‑1 (MQ=255)
                                                                                  tCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGc                                                                    >  1:5893078/1‑150 (MQ=255)
                                                                                  tCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGc                                                                    >  1:1063084/1‑150 (MQ=255)
                                                                                   cgccgcCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCt                                                                   >  1:2752509/1‑150 (MQ=255)
                                                                                     ccgccATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGt                                                                 >  1:6257609/1‑150 (MQ=255)
                                                                                       gccATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACCCTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGttt                                                               <  1:1349376/150‑1 (MQ=255)
                                                                                       gccATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGttt                                                               <  1:227638/150‑1 (MQ=255)
                                                                                         cATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgt                                                             <  2:3630172/150‑1 (MQ=255)
                                                                                             gggCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGa                                                         <  2:2095275/150‑1 (MQ=255)
                                                                                             gggCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGa                                                         <  2:3026518/150‑1 (MQ=255)
                                                                                              ggCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAg                                                        >  1:4024038/1‑150 (MQ=255)
                                                                                               gcgcGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGc                                                       >  2:6234523/1‑150 (MQ=255)
                                                                                               gcgcGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGc                                                       <  2:1097930/150‑1 (MQ=255)
                                                                                                cgcgATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGcc                                                      <  2:4507520/150‑1 (MQ=255)
                                                                                                 gcgATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCt                                                     <  1:735043/150‑1 (MQ=255)
                                                                                                   gATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTcc                                                   >  1:369357/1‑150 (MQ=255)
                                                                                                   gATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTcc                                                   <  2:3567026/150‑1 (MQ=255)
                                                                                                   gATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTcc                                                   >  2:6065235/1‑150 (MQ=255)
                                                                                                   gATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTcc                                                   >  1:6536326/1‑150 (MQ=255)
                                                                                                         tGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGtt                                             <  1:4395325/150‑1 (MQ=255)
                                                                                                         tGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGtt                                             >  2:5746822/1‑150 (MQ=255)
                                                                                                          gTAGGATCCGGTATACCCGCGAAAGCGTCAGCACTCACAACAATGTACCCTCAACGGACGAATTCGAGGGTGAACCCGCTCTAAAATCGGGCCAGTATCGCCATCTGGTCGTCGCAAGAATCCAGCTGTTTGTGGTAGCCTCCCAGGTTg                                            >  2:2407474/1‑150 (MQ=255)
                                                                                                              gAGCGGGTTCACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCGGGTAAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAACCCAGCTGTTTGTGGGAGCCTCCCATGTTGAaca                                        >  2:3261616/1‑150 (MQ=17)
                                                                                                              gAGCGGGTTCACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTcgc                                                                                                                               >  2:5641140/1‑63 (MQ=255)
                                                                                                              gAGCGGGTTCACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTcgc                                                                                                                               <  1:5641140/63‑1 (MQ=255)
                                                                                                                gCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACAcc                                      <  2:3416116/150‑1 (MQ=255)
                                                                                                                    gTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGGTTGTGGGAGCCTCCCATGTTGAACACCCTTc                                  >  1:2437402/1‑150 (MQ=255)
                                                                                                                      tcaCCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTa                                >  2:2120135/3‑150 (MQ=21)
                                                                                                                        aCCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTAcc                              >  1:1003241/1‑150 (MQ=255)
                                                                                                                          ccGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCa                            >  1:1243014/1‑150 (MQ=21)
                                                                                                                          ccGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCGGGTAAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTGTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCa                            >  2:1991322/1‑150 (MQ=18)
                                                                                                                          ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCa                            >  1:1751907/1‑150 (MQ=255)
                                                                                                                          ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCa                            >  1:2112930/1‑150 (MQ=255)
                                                                                                                           cgcgAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTAGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCAt                           >  1:6405846/1‑150 (MQ=255)
                                                                                                                              gAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCct                        >  1:3028417/1‑150 (MQ=255)
                                                                                                                               actGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCctg                       <  2:4786079/148‑1 (MQ=255)
                                                                                                                                     tCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTc                 >  2:1515876/1‑150 (MQ=255)
                                                                                                                                        gCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCt              <  2:489349/150‑1 (MQ=255)
                                                                                                                                         gtcaggcaacaGTGTTGTCAGTGCTGACGCATTCGCGGGTAAACCCGCTCCTACA‑‑GGG‑‑ATCGCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCtt        >  2:1810058/7‑150 (MQ=17)
                                                                                                                                          aCTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGc            >  2:5513300/1‑150 (MQ=255)
                                                                                                                                          aCTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGc            >  2:5700061/1‑150 (MQ=255)
                                                                                                                                           cTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGcc           <  2:6839995/150‑1 (MQ=255)
                                                                                                                                            tGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCGGTGTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGccc          >  2:3495984/1‑150 (MQ=255)
                                                                                                                                             gACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCt         >  1:5469233/1‑150 (MQ=255)
                                                                                                                                               caacaCTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTg       <  2:50156/150‑1 (MQ=255)
                                                                                                                                                aacaCTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTgg      >  1:4976156/1‑150 (MQ=255)
                                                                                                                                                 acacTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTACTGTTCGCTGCCCTTGGc     >  2:374643/1‑150 (MQ=255)
                                                                                                                                                   acTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCt   >  1:1431436/1‑150 (MQ=255)
                                                                                                                                                   acTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCt   <  2:2944927/150‑1 (MQ=255)
                                                                                                                                                   acTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCt   <  1:989825/150‑1 (MQ=255)
                                                                                                                                                   acTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGGTTGTGGGAGCCTCCCATGTTGAACACCCGTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCt   >  2:6107026/1‑150 (MQ=255)
                                                                                                                                                    cTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCTg  <  1:5162673/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCTG  >  NC_002947/372512‑372809

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: