Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 372,657 C→G 6.1% 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,6570CG 6.1% 210.8 / 3.5 66intergenic (+77/‑102)dgcA/dgcBdimethylglycine dehydrogenase subunit/dimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base G (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

CGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGG  >  NC_002947/372511‑372802
                                                                                                                                                  |                                                                                                                                                 
cGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTg                                                                                                                                                 <  6:3406/149‑1 (MQ=255)
      cGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGc                                                                                                                                             <  4:141228/147‑1 (MQ=255)
      cGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCt                                                                                                                                           <  8:56065/149‑1 (MQ=255)
         ggCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAc                                                                                                                                        <  3:146663/149‑1 (MQ=255)
          gCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAcc                                                                                                                                       <  8:148769/149‑1 (MQ=255)
          gCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGAc                                                                                                                                        <  7:169619/148‑1 (MQ=255)
           cccAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCg                                                                                                                                      >  8:83674/1‑149 (MQ=255)
             cAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGa                                                                                                                                    >  1:269292/1‑149 (MQ=255)
             cAGCGCAACGTGCATGCGGCGATCAACGACGCCTTGCGCCTGTGAAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGa                                                                                                                                    >  7:228080/1‑149 (MQ=255)
              aGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGAc                                                                                                                                   <  1:142176/149‑1 (MQ=255)
                      gTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTcgc                                                                                                                           <  3:120149/149‑1 (MQ=21)
                         cATGCGTCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCTTGTAGGAGCGGGTTTACCCGAGAATGCGTCAGCACTGAGAACACTGTTGCCTGACCGGACGCATTCGCggg                                                                                                                        <  8:280954/149‑1 (MQ=255)
                         cATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCgg                                                                                                                         <  2:269292/148‑1 (MQ=255)
                              ggcgAGCTACGACGCCTTGCGCCTGTGAAAGGACTTCTCATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAAc                                                                                                                   <  8:81072/149‑1 (MQ=255)
                                 gATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCg                                                                                                                >  5:269412/1‑149 (MQ=255)
                                      acgacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGctct                                                                                                            >  3:341053/1‑148 (MQ=255)
                                        gacgCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTAc                                                                                                         <  5:147696/149‑1 (MQ=255)
                                             cTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATggg                                                                                                    >  3:33601/1‑149 (MQ=255)
                                             cTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATggg                                                                                                    >  8:292720/1‑149 (MQ=255)
                                             cTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATggg                                                                                                    <  1:105560/149‑1 (MQ=255)
                                                gcgcCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGc                                                                                                  <  4:79465/148‑1 (MQ=255)
                                                  gcCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAg                                                                                               <  7:118551/149‑1 (MQ=255)
                                                   ccTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATCCGTCAGCACTGACAACACTGATGCCTGACCGGACGCATTCGCGCGTCAACACGCTCTTACATGGGGACAGt                                                                                              >  3:240322/1‑149 (MQ=255)
                                                      gtgCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGtctc                                                                                           >  5:292131/1‑149 (MQ=255)
                                                       tgCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGtctc                                                                                           <  8:131217/148‑1 (MQ=255)
                                                       tgCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGtctc                                                                                           <  2:238985/148‑1 (MQ=255)
                                                         cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGcc                                                                                        <  2:1261/149‑1 (MQ=255)
                                                         cAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGcc                                                                                        >  2:182235/1‑149 (MQ=255)
                                                            ggACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCctc                                                                                     <  4:33601/149‑1 (MQ=255)
                                                                ttCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGt                                                                                 >  5:195117/1‑149 (MQ=255)
                                                                     aTCCTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCAGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACTCATTCGCGGGTTAACCCGCTCTTACATGTGGCCAGTCTCGCCCTCTGGTCg                                                                               <  8:228080/146‑1 (MQ=255)
                                                                        gTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGgcgc                                                                           >  4:191147/1‑147 (MQ=255)
                                                                        gTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGgcgc                                                                           >  6:80727/1‑147 (MQ=255)
                                                                         ttAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGc                                                                                                                                 >  4:87113/1‑92 (MQ=255)
                                                                         ttAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGc                                                                                                                                 <  3:87113/92‑1 (MQ=255)
                                                                         ttAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGgcgc                                                                           >  1:94126/1‑146 (MQ=255)
                                                                                tCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAg                                                                 <  7:292720/149‑1 (MQ=255)
                                                                                  gccgccATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCt                                                               <  6:292131/149‑1 (MQ=255)
                                                                                  gccgccATCGGGCGCGATCCCTGTAGGAGCGAATTTACCCGCGAATACATCAGCACTGACAACCCAGTTGCCTGACCGAAAGCATTCGCGGGTGAACCCGCGCTTACAAGGGGCCAGACACGCCCGCTGGGCGGCGCAAGAAACCacgg                                                               >  8:59205/1‑146 (MQ=255)
                                                                                       cATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTtg                                                          >  4:50385/1‑149 (MQ=255)
                                                                                           gggCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCtt                                                                                                           >  4:318169/1‑96 (MQ=255)
                                                                                           gggCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCtt                                                                                                           <  3:318169/96‑1 (MQ=255)
                                                                                             gcgcGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAg                                                    <  7:59205/149‑1 (MQ=255)
                                                                                              cgcgATCCCTGTAGGAGCGGGTTTACCGGCTAATGCGTCCACCTTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGc                                                   <  1:130854/149‑1 (MQ=255)
                                                                                                         tAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGtt                                         <  2:94126/148‑1 (MQ=255)
                                                                                                           ggAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGaa                                      <  3:318436/149‑1 (MQ=255)
                                                                                                             aGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAaca                                    <  1:182235/149‑1 (MQ=255)
                                                                                                               cGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACAcc                                  >  3:189286/1‑149 (MQ=255)
                                                                                                                   tttACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCtgc                              >  4:214280/1‑147 (MQ=255)
                                                                                                                       cccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTAcc                          >  3:105176/1‑149 (MQ=255)
                                                                                                                        ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGCACACCCTTCTAcc                          >  5:207245/1‑148 (MQ=255)
                                                                                                                        ccGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTAcc                          >  6:237075/1‑148 (MQ=255)
                                                                                                                         cgcgAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCGGGTAAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCa                        >  2:64936/1‑149 (MQ=21)
                                                                                                                          gcgAATGCGTCCGGTCAGGCAACAGTGTTGTCAGTGCTGACGCATTCGCGGGTAAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCAt                       >  6:364603/1‑149 (MQ=21)
                                                                                                                            gAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTccc                                              >  6:336090/1‑124 (MQ=255)
                                                                                                                            gAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTccc                                              <  5:336090/124‑1 (MQ=255)
                                                                                                                              aTGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCctg                   <  4:341053/149‑1 (MQ=255)
                                                                                                                               tGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCctgc                  >  5:253123/1‑149 (MQ=255)
                                                                                                                                   tCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGtt              >  5:387108/1‑149 (MQ=255)
                                                                                                                                     aGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCg            <  4:189286/149‑1 (MQ=255)
                                                                                                                                       cACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCt          >  6:292509/1‑149 (MQ=255)
                                                                                                                                         cTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGGTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGc        >  8:296333/1‑149 (MQ=255)
                                                                                                                                         cTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGc        <  5:237075/149‑1 (MQ=255)
                                                                                                                                          tGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGcc       >  2:227502/1‑149 (MQ=255)
                                                                                                                                           ggcaacaGTGTTGTCAGTGCTGACGCATTCGCGGGTAAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGccc      <  1:64936/147‑1 (MQ=37)
                                                                                                                                             caacaCTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCtt    <  6:195117/149‑1 (MQ=255)
                                                                                                                                               acacTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTgg  >  6:52372/1‑149 (MQ=255)
                                                                                                                                                cacTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTgg  >  8:251509/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGACTGCGTGGCCCAGCGCAACGTGCATGCGGCGATCTACGACGCCTTGCGCCTGTGCAAGGACTTCTGATCGTTAACAGTCGCCGCCATCGGGCGCGATCCCTGTAGGAGCGGGTTTACCCGCGAATGCGTCAGCACTGACAACACTGTTGCCTGACCGGACGCATTCGCGGGTGAACCCGCTCTTACATGGGGCCAGTCTCGCCCTCTGGTCGGCGCAAGAATCCAGCTGTTTGTGGGAGCCTCCCATGTTGAACACCCTTCTACCCATCCTGCTGTTCGCTGCCCTTGG  >  NC_002947/372511‑372802

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: