Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 372,647 A→C 5.3% intergenic (+64/‑115) 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,6470AC 5.3% 785.3 / 24.4 226intergenic (+64/‑115)dgcA/dgcBdimethylglycine dehydrogenase subunit/dimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base A (87/126);  new base C (7/5);  total (95/131)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 226 total aligned reads displayed.

GTACCGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCG  >  NC_002947/372502‑372795
                                                                                                                                                 |                                                                                                                                                    
gTACCGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCAc                                                                                                                                                  <  2:1809890/150‑1 (MQ=255)
   ccGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGa                                                                                                                                               <  1:373754/150‑1 (MQ=255)
   ccGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGa                                                                                                                                               <  2:4250758/150‑1 (MQ=255)
    cGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCACCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGac                                                                                                                                              <  2:4088348/150‑1 (MQ=255)
        tCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGacaaca                                                                                                                                          <  2:663192/150‑1 (MQ=255)
         cGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAacac                                                                                                                                         <  1:2879325/150‑1 (MQ=255)
           gCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTg                                                                                                                                       >  1:2631103/1‑150 (MQ=255)
             gACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGtt                                                                                                                                     <  1:1970875/150‑1 (MQ=255)
               cTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGGTCAGGCAACAGTGTTgt                                                                                                                                   <  2:3119323/150‑2 (MQ=39)
                tGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGcc                                                                                                                                  <  2:3100216/150‑1 (MQ=255)
                tGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGcc                                                                                                                                  <  2:3003267/150‑1 (MQ=255)
                 gCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCt                                                                                                                                 <  1:1538908/150‑1 (MQ=255)
                   gTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGa                                                                                                                               >  1:3213524/1‑150 (MQ=255)
                    tGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAc                                                                                                                              <  2:100009/150‑1 (MQ=255)
                     ggCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAcc                                                                                                                             <  2:4213552/150‑1 (MQ=255)
                       cccAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCgg                                                                                                                           >  2:397428/1‑150 (MQ=255)
                        ccAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGa                                                                                                                          <  1:4132089/150‑1 (MQ=255)
                        ccAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGa                                                                                                                          >  2:524522/1‑150 (MQ=255)
                         cAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                         >  2:3731835/1‑150 (MQ=255)
                         cAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                         <  2:962469/150‑1 (MQ=255)
                         cAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                         <  2:882734/150‑1 (MQ=255)
                         cAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                         <  1:2120344/150‑1 (MQ=255)
                                 cGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAAGAGCGGGTTCACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTcgc                                                                                                                 <  2:2479005/150‑1 (MQ=18)
                                   tGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg                                                                                                               <  1:2279544/150‑1 (MQ=255)
                                   tGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg                                                                                                               <  2:4486233/150‑1 (MQ=255)
                                   tGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg                                                                                                               <  1:1571340/150‑1 (MQ=255)
                                   tGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg                                                                                                               <  1:4411481/150‑1 (MQ=255)
                                   tGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACGTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGTATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg                                                                                                               <  2:4166504/150‑1 (MQ=255)
                                   tGCATGCGGCGATCTACGACGCCGTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCgg                                                                                                               <  2:3036417/150‑1 (MQ=21)
                                    gCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCggg                                                                                                              <  1:4473855/150‑1 (MQ=255)
                                            cgATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACGACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCg                                                                                                      <  2:1333727/150‑1 (MQ=255)
                                               tCTACGATGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGctc                                                                                                   <  2:2817941/150‑1 (MQ=255)
                                                 tACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATGGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCtt                                                                                                 <  1:2148175/150‑1 (MQ=255)
                                                 tACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCtt                                                                                                 <  2:2972340/150‑1 (MQ=255)
                                                   cgacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTAc                                                                                               <  1:411678/150‑1 (MQ=255)
                                                   cgacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTAc                                                                                               <  1:3542654/150‑1 (MQ=255)
                                                    gacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACa                                                                                              <  1:1792586/150‑1 (MQ=255)
                                                    gacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACa                                                                                              <  2:4522941/150‑1 (MQ=255)
                                                    gacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCAGTCGCGGGTGAACCCGCGCTTACa                                                                                              >  1:3726157/1‑150 (MQ=255)
                                                     acgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACAt                                                                                             >  2:3942773/1‑150 (MQ=255)
                                                     acgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACAt                                                                                             <  2:2666873/150‑1 (MQ=255)
                                                       gCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATgg                                                                                           <  2:3636789/150‑1 (MQ=255)
                                                       gCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATgg                                                                                           <  1:4327830/150‑1 (MQ=255)
                                                         cTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATgggg                                                                                         >  1:1830595/1‑150 (MQ=255)
                                                         cTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATgggg                                                                                         <  1:1599567/150‑1 (MQ=255)
                                                          ttGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGc                                                                                        <  2:3998445/150‑1 (MQ=255)
                                                             cgcCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAg                                                                                     <  1:1249077/150‑1 (MQ=255)
                                                                cTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGtct                                                                                  >  1:1127264/1‑150 (MQ=255)
                                                                cTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGtct                                                                                  <  2:2118208/150‑1 (MQ=255)
                                                                cTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGtct                                                                                  >  1:3243179/1‑150 (MQ=255)
                                                                cTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACGGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACCTGGGGCCAGtct                                                                                  >  2:3407740/1‑150 (MQ=255)
                                                                  gtgCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCg                                                                                <  2:2602275/150‑1 (MQ=255)
                                                                     cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACAGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGccc                                                                             >  1:2936670/1‑150 (MQ=255)
                                                                      aaGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCct                                                                            >  1:2519890/1‑150 (MQ=255)
                                                                        ggACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCctct                                                                          >  2:2715091/1‑150 (MQ=255)
                                                                         gACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTg                                                                         <  2:1436444/150‑1 (MQ=255)
                                                                           cTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGt                                                                       <  1:131636/150‑1 (MQ=255)
                                                                            ttCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTc                                                                      >  2:3058521/1‑150 (MQ=255)
                                                                            ttCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTc                                                                      <  1:3662433/150‑1 (MQ=255)
                                                                             tCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGt                                                                       >  1:2880617/1‑148 (MQ=255)
                                                                             tCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCg                                                                     <  1:532733/150‑1 (MQ=255)
                                                                             tCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCg                                                                     <  1:940141/150‑1 (MQ=255)
                                                                                gATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCGGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGgcg                                                                  >  1:1875197/1‑150 (MQ=255)
                                                                                gATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGgcg                                                                  <  1:3945146/150‑1 (MQ=255)
                                                                                   cGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCaa                                                               <  2:2240562/150‑1 (MQ=255)
                                                                                   cGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCaa                                                               <  1:2643000/150‑1 (MQ=255)
                                                                                       aaCAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCACGact                                                           >  2:3608278/1‑148 (MQ=255)
                                                                                       aaCAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAAt                                                           >  1:165213/1‑150 (MQ=255)
                                                                                          cgTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCa                                                        <  1:3896303/149‑1 (MQ=255)
                                                                                              gccgccATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTCCATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTg                                                    >  2:3226204/1‑150 (MQ=255)
                                                                                                   cATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgt                                               <  2:1750260/150‑1 (MQ=255)
                                                                                                   cATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgt                                               <  2:1909159/150‑1 (MQ=255)
                                                                                                   cATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgt                                               <  1:617862/150‑1 (MQ=255)
                                                                                                   cATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTAGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgt                                               >  1:541071/1‑150 (MQ=255)
                                                                                                    aTCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgtg                                              <  1:3760338/150‑1 (MQ=255)
                                                                                                    aTCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtgtg                                              >  2:2367756/1‑150 (MQ=255)
                                                                                                     tCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTgg                                             <  2:3791469/150‑1 (MQ=255)
                                                                                                       gggCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGa                                           <  1:3414894/150‑1 (MQ=255)
                                                                                                       gggCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGa                                           <  1:4289394/150‑1 (MQ=255)
                                                                                                       gggCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGa                                           <  2:1589537/150‑1 (MQ=255)
                                                                                                       gggCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGAATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGa                                           <  2:1515749/150‑1 (MQ=255)
                                                                                                         gcgcGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGc                                         >  2:4010819/1‑150 (MQ=255)
                                                                                                            cgATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTc                                      <  2:1108418/150‑1 (MQ=255)
                                                                                                              aTCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTccc                                    >  2:2611625/1‑150 (MQ=255)
                                                                                                                 ccTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATg                                 >  1:4454397/1‑150 (MQ=255)
                                                                                                                   tGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCGAGCTGTTTGTGGGAGCCTCCCATGtt                               >  2:2501196/1‑150 (MQ=255)
                                                                                                                   tGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGtt                               <  1:1679449/150‑1 (MQ=255)
                                                                                                                      aGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGaa                            <  2:4377854/150‑1 (MQ=255)
                                                                                                                        gAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAaca                          >  2:2786582/1‑150 (MQ=255)
                                                                                                                             ggTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCtt                     >  2:3276086/1‑150 (MQ=255)
                                                                                                                                ttACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTa                  <  1:3562416/150‑1 (MQ=255)
                                                                                                                                   cccGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTAccc               >  1:2424070/1‑150 (MQ=21)
                                                                                                                                   cccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGTATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTAccc               >  1:971060/1‑150 (MQ=255)
                                                                                                                                    ccGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCa              >  1:3068232/1‑150 (MQ=21)
                                                                                                                                    ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCa              >  1:3527544/1‑150 (MQ=255)
                                                                                                                                    ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCa              >  2:2068054/1‑150 (MQ=255)
                                                                                                                                    ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCa              >  1:1062026/1‑150 (MQ=255)
                                                                                                                                    ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTGGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCa              >  2:2188629/1‑150 (MQ=255)
                                                                                                                                              gTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGtt    <  1:880328/150‑1 (MQ=255)
                                                                                                                                                cAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCg  >  1:3152425/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTACCGCATCGGCGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCG  >  NC_002947/372502‑372795

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: