Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 373,393 C→G 5.3% T212S (ACC→AGC)  dgcB → dimethylglycine dehydrogenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947373,3930CG 5.3% 245.9 / 3.8 76T212S (ACC→AGC) dgcBdimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base C (38/34);  new base G (2/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

ATTCTGGGCCTGGGCGTGCTGTGGGGTGTGTCGGAGCTGTTCTTTGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCTTTGCCGCTGGC  >  NC_002947/373251‑373541
                                                                                                                                              |                                                                                                                                                    
aTTCTGGGCCTGGGCGTGCTGTGGGGTGTGTCGGAGCTGTTCTTTGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCt                                                                                                                                                >  3:326328/1‑149 (MQ=255)
    tgggcctgggcGTGCTGTGGGGTGTGTCGGAGCTGTTCTTTGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAg                                                                                                                                            >  3:49940/1‑149 (MQ=255)
         ctgggcGTGCTGTGGGGTGTGTCGGAGCTGTTCTTTGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGGTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCt                                                                                                                                       <  4:206892/149‑1 (MQ=255)
            ggcGTGCTGTGGGGTGTGTCGGAGCTGTTCTTTGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGa                                                                                                                                    >  8:11752/1‑149 (MQ=255)
             gcGTGCTGTGGGGTGTGTCGGAGCTGTTCTTTGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGAc                                                                                                                                   <  5:44058/149‑1 (MQ=255)
                tGCTGTGGGGTGTGTCGGAGCTGTTCTTTGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCAGCGGCTTGAGGCCGGTGGagcgg                                                                                                                                <  4:369332/149‑5 (MQ=255)
                 gCTGTGGGGTGTGTCGGAGCTGTTCTTTGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTg                                                                                                                                >  6:195082/1‑148 (MQ=255)
                    gtgGGGTGTGTCGGAGCTGTTCTTTGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTgg                                                                                                                               >  5:126877/1‑146 (MQ=255)
                        ggTGTGTCGGAGCTGTTCTTTGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGAccc                                                                                                                        <  3:427652/149‑1 (MQ=255)
                         gtgtgtCGGAGCTGTTCTTTGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCg                                                                                                                       <  6:126877/149‑1 (MQ=255)
                            tgtCGGAGCTGTTCTTTGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCAGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAAc                                                                                                                    >  5:310982/1‑149 (MQ=255)
                                ggAGCTGTTCTTTGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACgcgc                                                                                                                >  7:210209/1‑149 (MQ=255)
                                            tGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTg                                                                                                     <  7:199151/148‑1 (MQ=255)
                                                aTGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGaaaa                                                                                                <  2:81278/149‑1 (MQ=255)
                                                aTGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGaaa                                                                                                 >  2:271993/1‑148 (MQ=255)
                                                aTGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGaaa                                                                                                 <  1:271993/148‑1 (MQ=255)
                                                     cTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCgg                                                                                           >  8:166511/1‑149 (MQ=255)
                                                       ggggCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCgg                                                                                           >  2:238391/1‑147 (MQ=255)
                                                       ggggCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCgg                                                                                           <  1:238391/147‑1 (MQ=255)
                                                         ggcggcCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGa                                                                                       >  4:223956/1‑149 (MQ=255)
                                                          gcggcCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGAc                                                                                      >  3:312816/1‑149 (MQ=255)
                                                             gcCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACtt                                                                                    >  2:282891/1‑148 (MQ=255)
                                                              cccGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCa                                                                                  >  2:372270/1‑149 (MQ=255)
                                                               ccGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCAc                                                                                 >  2:58778/1‑149 (MQ=255)
                                                                  aTGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTg                                                                              <  8:210209/149‑1 (MQ=255)
                                                                   tGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTgg                                                                             <  1:282891/149‑1 (MQ=255)
                                                                        cACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAAc                                                                          >  6:260439/1‑147 (MQ=255)
                                                                         aCGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCCCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAg                                                                       >  3:281087/1‑149 (MQ=255)
                                                                              ttCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAgctgct                                                                  >  5:121762/1‑149 (MQ=255)
                                                                              ttCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAgctgct                                                                  >  1:411972/1‑149 (MQ=255)
                                                                              ttCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAgctgct                                                                  >  4:170212/1‑149 (MQ=255)
                                                                                 gCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCAGCGGCTTGAGGCCGGTGGa                                                                                                                                    >  3:209834/1‑80 (MQ=38)
                                                                                 gCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCAGCGGCTTGAGGCCGGTGGa                                                                                                                                    <  4:209834/80‑1 (MQ=38)
                                                                                 gCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTggg                                                               >  1:244073/1‑149 (MQ=255)
                                                                                 gCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTggg                                                               >  2:388261/1‑149 (MQ=255)
                                                                                 gCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTgg                                                                <  3:51815/148‑1 (MQ=255)
                                                                                 gCCAGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCTCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTgg                                                                <  6:310982/148‑1 (MQ=255)
                                                                                     gTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTc                                                           >  5:234864/1‑149 (MQ=255)
                                                                                     gTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTc                                                           <  5:195082/149‑1 (MQ=255)
                                                                                      tGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCCCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCg                                                          <  4:281087/149‑1 (MQ=255)
                                                                                        ccTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGAt                                                        <  4:49940/149‑1 (MQ=255)
                                                                                         cTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATg                                                       >  6:33078/1‑149 (MQ=255)
                                                                                         cTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATg                                                       <  7:166511/149‑1 (MQ=255)
                                                                                              aCCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGc                                                  <  2:244073/149‑1 (MQ=255)
                                                                                                ctggcctggcACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAg                                                                       >  6:258707/1‑126 (MQ=255)
                                                                                                ctggcctggcACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAg                                                                       <  5:258707/126‑1 (MQ=255)
                                                                                                        gcACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCg                                        >  6:248501/1‑149 (MQ=255)
                                                                                                         cACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGCACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCgg                                       >  2:225676/1‑149 (MQ=255)
                                                                                                             gccgTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTaaa                                   >  8:14373/1‑149 (MQ=255)
                                                                                                             gccgTGCCGAACGCCTCTGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAACCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTaaa                                   <  2:118315/149‑1 (MQ=255)
                                                                                                               cgTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAAtg                                 >  8:123481/1‑149 (MQ=255)
                                                                                                                gTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAAtgt                                >  3:330882/1‑149 (MQ=255)
                                                                                                                 tGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAAtgtg                               <  7:117607/149‑1 (MQ=255)
                                                                                                                  gCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAAtgtg                               <  4:56439/148‑1 (MQ=255)
                                                                                                                  gCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAAtgtg                               >  3:56439/1‑148 (MQ=255)
                                                                                                                    cGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAg                            >  4:122969/1‑149 (MQ=255)
                                                                                                                        cGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCAt                        <  5:33078/149‑1 (MQ=255)
                                                                                                                        cGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCAt                        <  3:286704/149‑1 (MQ=255)
                                                                                                                        cGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCAt                        <  6:234864/149‑1 (MQ=255)
                                                                                                                         gCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCAtg                       >  3:31083/1‑149 (MQ=255)
                                                                                                                           ttCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCAtgtg                     >  4:185215/1‑149 (MQ=255)
                                                                                                                             cggcggcggCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCAtgtg                     <  7:123481/147‑1 (MQ=255)
                                                                                                                               gcggcggcCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGACCGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACAAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCg                 >  3:3570/1‑149 (MQ=255)
                                                                                                                                    ggcCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCtt            >  3:258416/1‑149 (MQ=255)
                                                                                                                                     gcCGCTGCACCGGCCTCAAGGCGGGGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCttt           <  1:296977/149‑1 (MQ=255)
                                                                                                                                     gcCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCttt           <  4:166265/149‑1 (MQ=255)
                                                                                                                                     gcCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCttt           <  1:58778/149‑1 (MQ=255)
                                                                                                                                      ccGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCTTTg          >  6:65970/1‑149 (MQ=255)
                                                                                                                                      ccGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCTTTg          <  1:388261/149‑1 (MQ=255)
                                                                                                                                       cGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCTTTGc         <  5:248501/149‑1 (MQ=255)
                                                                                                                                         cTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCTTTGCCg       >  7:192993/1‑149 (MQ=255)
                                                                                                                                           ccAGCGGCTTGAGGCCGGTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCTTTGCCGCt     >  8:228627/1‑149 (MQ=255)
                                                                                                                                            cACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCTTTGCCGCt     <  2:411972/148‑1 (MQ=255)
                                                                                                                                            cACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCTTTGCCGCTg    <  6:121762/149‑1 (MQ=255)
                                                                                                                                            cACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCTTTGCCGCTg    <  4:31083/149‑1 (MQ=255)
                                                                                                                                             aCCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCTTTGCCGCTgg   <  5:260439/149‑1 (MQ=255)
                                                                                                                                              ccGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCTTTGCCGCTGGc  <  7:14373/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
ATTCTGGGCCTGGGCGTGCTGTGGGGTGTGTCGGAGCTGTTCTTTGGCATGACCTGGGGCGGCCCGATGAAGCACGCCTTCGCCGGTGCCTTGCACCTGGCCTGGCACCGCCGTGCCGAACGCTTCGGCGGCGGCCGCTCCACCGGCCTCAAGCCGCTGGACCTGGAAGACCCGAACGCGCCCCTTGGCGTGGAAAAGCCGGTGGACTTCACCTGGAACCAGCTGCTGGGCTTCGATGCCTGCGTGCAGTGCGGTAAATGTGAAGCCATGTGCCCGGCCTTTGCCGCTGGC  >  NC_002947/373251‑373541

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