Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 372,660 C→G 6.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 6.0% 786.6 / 29.9 232intergenic (+77/‑102)dgcA/dgcBdimethylglycine dehydrogenase subunit/dimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base C (98/119);  new base G (7/7);  total (106/126)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 235 total aligned reads displayed.

CGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCT  >  NC_002947/372511‑372808
                                                                                                                                                     |                                                                                                                                                    
cGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAacac                                                                                                                                                      <  1:2879325/150‑1 (MQ=255)
 ggCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACt                                                                                                                                                     >  1:2995647/1‑150 (MQ=255)
 ggCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACt                                                                                                                                                     <  2:309607/150‑1 (MQ=255)
  gCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTg                                                                                                                                                    >  1:2631103/1‑150 (MQ=255)
  gCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTg                                                                                                                                                    >  1:821940/1‑150 (MQ=255)
  gCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTg                                                                                                                                                    >  1:3091251/1‑150 (MQ=255)
    gACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGtt                                                                                                                                                  <  2:3804823/150‑1 (MQ=255)
    gACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGtt                                                                                                                                                  <  1:1970875/150‑1 (MQ=255)
       tGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGcc                                                                                                                                               <  2:3003281/150‑1 (MQ=255)
        gCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCt                                                                                                                                              <  1:3134047/150‑1 (MQ=255)
        gCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCt                                                                                                                                              <  1:1538908/150‑1 (MQ=255)
        gCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCt                                                                                                                                              <  1:1166862/150‑1 (MQ=255)
          gTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGa                                                                                                                                            >  1:3213524/1‑150 (MQ=255)
          gTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGa                                                                                                                                            >  2:3273334/1‑150 (MQ=255)
           tGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAc                                                                                                                                           <  1:1919298/150‑1 (MQ=255)
            ggCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAcc                                                                                                                                          <  2:4213552/150‑1 (MQ=255)
               ccAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGa                                                                                                                                       >  2:524522/1‑150 (MQ=255)
                cAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                                      <  2:962469/150‑1 (MQ=255)
                cAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                                      <  1:2120344/150‑1 (MQ=255)
                cAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                                      <  2:882734/150‑1 (MQ=255)
                     cAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCAtt                                                                                                                                 <  1:1663651/150‑1 (MQ=255)
                      aaCGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTc                                                                                                                                >  2:3426426/1‑150 (MQ=255)
                        cGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAAGAGCGGGTTCACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTcgc                                                                                                                              <  1:1699493/150‑1 (MQ=18)
                          tGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCgg                                                                                                                            <  1:4470320/150‑1 (MQ=21)
                          tGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg                                                                                                                            <  1:2279544/150‑1 (MQ=255)
                          tGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg                                                                                                                            <  2:2102641/150‑1 (MQ=255)
                          tGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACGTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGTATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg                                                                                                                            <  2:4166504/150‑1 (MQ=255)
                                   cgATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACGACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCg                                                                                                                   <  2:1333727/150‑1 (MQ=255)
                                      tCTACGATGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGctc                                                                                                                <  2:2817941/150‑1 (MQ=255)
                                        tACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCtt                                                                                                              <  2:2972340/150‑1 (MQ=255)
                                          cgacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTAc                                                                                                            <  1:3542654/150‑1 (MQ=255)
                                           gacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACa                                                                                                           <  2:4522941/150‑1 (MQ=255)
                                           gacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACa                                                                                                           >  2:2195317/1‑150 (MQ=255)
                                           gacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACa                                                                                                           <  1:1792586/150‑1 (MQ=255)
                                           gacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACa                                                                                                           <  2:1843588/150‑1 (MQ=255)
                                           gacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCAGTCGCGGGTGAACCCGCGCTTACa                                                                                                           >  1:3726157/1‑150 (MQ=255)
                                            acgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACAt                                                                                                          <  2:2666873/150‑1 (MQ=255)
                                                cTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATgggg                                                                                                      >  1:1830595/1‑150 (MQ=255)
                                                cTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATgggg                                                                                                      >  2:2416606/1‑150 (MQ=255)
                                                       cTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGtct                                                                                               >  2:303673/1‑150 (MQ=255)
                                                       cTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGtct                                                                                               >  1:3243179/1‑150 (MQ=255)
                                                       cTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGtct                                                                                               <  2:2118233/150‑1 (MQ=255)
                                                       cTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACGGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACCTGGGGCCAGtct                                                                                               >  2:3407740/1‑150 (MQ=255)
                                                           gCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGcc                                                                                           <  1:4408217/150‑1 (MQ=255)
                                                            cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGccc                                                                                          >  1:1625645/1‑150 (MQ=255)
                                                             aaGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCct                                                                                         >  1:2519890/1‑150 (MQ=255)
                                                              aGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCctc                                                                                        >  2:2554596/1‑150 (MQ=255)
                                                               ggACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCctct                                                                                       >  2:2715091/1‑150 (MQ=255)
                                                                gACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTg                                                                                      <  2:1436444/150‑1 (MQ=255)
                                                                  cTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGt                                                                                    <  2:3943058/150‑1 (MQ=255)
                                                                  cTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGt                                                                                    <  1:131636/150‑1 (MQ=255)
                                                                   ttCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTc                                                                                   <  1:3662433/150‑1 (MQ=255)
                                                                   ttCGGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTc                                                                                   >  1:464215/1‑150 (MQ=255)
                                                                    tCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGt                                                                                    >  1:2880617/1‑148 (MQ=255)
                                                                    tCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCg                                                                                  <  1:940141/150‑1 (MQ=255)
                                                                     cTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCgg                                                                                 <  1:2702489/150‑1 (MQ=255)
                                                                       gATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGgcg                                                                               <  1:3945146/150‑1 (MQ=255)
                                                                          cGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCaa                                                                            <  1:4122885/150‑1 (MQ=255)
                                                                           gTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAg                                                                           >  2:3682796/1‑150 (MQ=255)
                                                                              aaCAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAAt                                                                        >  2:3608322/1‑150 (MQ=255)
                                                                                    cgccgcCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCt                                                                  <  2:4001025/150‑1 (MQ=255)
                                                                                          cATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgt                                                            <  2:1909159/150‑1 (MQ=255)
                                                                                          cATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgt                                                            <  2:1750260/150‑1 (MQ=255)
                                                                                          cATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgt                                                            <  1:617862/150‑1 (MQ=255)
                                                                                          cATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTAGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgt                                                            >  1:541071/1‑150 (MQ=255)
                                                                                           aTCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgtg                                                           >  1:975547/1‑150 (MQ=255)
                                                                                           aTCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgtg                                                           <  1:1578852/150‑1 (MQ=255)
                                                                                           aTCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgtg                                                           <  1:3760162/150‑1 (MQ=255)
                                                                                              gggCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCACGAATCCAGCTGTTTGTGGGa                                                        >  1:1108418/1‑150 (MQ=255)
                                                                                              gggCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGa                                                        <  2:1589537/150‑1 (MQ=255)
                                                                                              gggCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGa                                                        <  1:4289394/150‑1 (MQ=255)
                                                                                              gggCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGAATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGa                                                        <  2:1515749/150‑1 (MQ=255)
                                                                                                gcgcGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCAGTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCGGTTTGTGGGAGc                                                      >  1:3862118/1‑150 (MQ=255)
                                                                                                   cgATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTc                                                   <  2:1108418/150‑1 (MQ=255)
                                                                                                   cgATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTc                                                   <  2:310322/150‑1 (MQ=255)
                                                                                                     aTCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCACGACTCCCGCTGGGTGTGGGAGCCTccc                                                 >  1:3390547/1‑150 (MQ=39)
                                                                                                     aTCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTccc                                                 >  2:2611572/1‑150 (MQ=255)
                                                                                                        ccTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATg                                              >  1:4454397/1‑150 (MQ=255)
                                                                                                          tGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCGAGCTGTTTGTGGGAGCCTCCCATGtt                                            >  2:2501196/1‑150 (MQ=255)
                                                                                                             aGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGaa                                         <  2:4377854/150‑1 (MQ=255)
                                                                                                              ggAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACTTGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAac                                        >  2:3767234/1‑150 (MQ=255)
                                                                                                               gAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAaca                                       >  2:2786501/1‑150 (MQ=255)
                                                                                                                    ggTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCtt                                  >  2:3122355/1‑150 (MQ=255)
                                                                                                                      ttcaCCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCGAGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGACTCCCGCGGTGTGTGGGAGCCTCCCCTGTTGAACACCCTTCt                                >  2:4416878/4‑150 (MQ=17)
                                                                                                                       ttACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTa                               <  1:3562416/150‑1 (MQ=255)
                                                                                                                          cccGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTATGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTACTAccc                            >  1:2424069/1‑150 (MQ=21)
                                                                                                                          cccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGTATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTAccc                            >  1:971060/1‑150 (MQ=255)
                                                                                                                           ccGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCGGGTAAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCa                           >  2:4045365/1‑150 (MQ=21)
                                                                                                                           ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCa                           >  1:3527544/1‑150 (MQ=255)
                                                                                                                                     gTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGtt                 >  2:399505/1‑150 (MQ=255)
                                                                                                                                         ccACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCt             <  2:2338538/149‑1 (MQ=255)
                                                                                                                                            cTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGcc          >  1:2235663/1‑150 (MQ=255)
                                                                                                                                             tGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTgg                                                          <  1:4318347/101‑1 (MQ=255)
                                                                                                                                             tGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTgg                                                          >  2:4318347/1‑101 (MQ=255)
                                                                                                                                               acaacaCTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCtt       >  2:2319888/1‑150 (MQ=255)
                                                                                                                                                 aacaCTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTgg     >  1:3898747/1‑150 (MQ=255)
                                                                                                                                                 aacaCTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTgg     >  1:4102378/1‑150 (MQ=255)
                                                                                                                                                 aacaCTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTgg     >  2:643651/1‑150 (MQ=255)
                                                                                                                                                 aacaCTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTgg     <  1:1385651/150‑1 (MQ=255)
                                                                                                                                                 aacaCTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTgg     >  1:1030990/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGGCCT  >  NC_002947/372511‑372808

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: