Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 697,231 C→G 5.6% intergenic (+155/‑596) mmsA‑I → / → PP_16SD methylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947697,2310CG 5.6% 481.4 / 14.7 142intergenic (+155/‑596)mmsA‑I/PP_16SDmethylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (77/57);  new base G (4/4);  total (81/61)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 158 total aligned reads displayed.

TAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGAC  >  NC_002947/697083‑697378
                                                                                                                                                    |                                                                                                                                                   
tAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGGTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCAGTGAGATCGAGCGCCGCCcgcg                                                                                                                                                    >  2:1714913/1‑150 (MQ=255)
 ataTGACGGGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCcgcgc                                                                                                                                                   <  1:1972853/148‑1 (MQ=255)
    tGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCgcggc                                                                                                                                                <  2:3803066/150‑1 (MQ=255)
       cgCGGGGAGGCCGTCCGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCt                                                                                                                                             <  1:2420879/149‑1 (MQ=255)
        gCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                            >  1:2344630/1‑150 (MQ=255)
        gCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                            >  1:1838831/1‑150 (MQ=255)
          ggggAGGCCTGACGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGGTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCgg                                                                                                                                          >  2:3009654/1‑150 (MQ=18)
          ggggAGGCCTGACGGCCACCCCGCTTTGTTTTTGGGGGGTGGTTGTTGCTGTTTGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCgg                                                                                                                                          >  1:860062/1‑150 (MQ=18)
          ggggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGGTGGGGAGGACCTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCcgcgcggcgcgcgg                                                                                                                                          >  2:1999079/1‑150 (MQ=255)
          ggggAGGCCGTCAGGCCTCCCCGCTTTGGTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCgg                                                                                                                                          >  2:4405971/1‑150 (MQ=255)
            ggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                                        <  1:2830128/150‑1 (MQ=18)
             gAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGtct                                                                                                                                       <  2:2336901/150‑1 (MQ=255)
               ggCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTca                                                                                                                                     <  1:3746627/150‑1 (MQ=255)
               ggCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTca                                                                                                                                     <  2:2498545/150‑1 (MQ=255)
                gCCGTCTGGCCTCCCCGCTTTGTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTcaca                                                                                                                                   <  2:737298/150‑1 (MQ=255)
                gCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTcac                                                                                                                                    >  2:859797/1‑150 (MQ=255)
                  ctgacggCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAg                                                                                                                                  >  2:3379175/6‑150 (MQ=255)
                   gTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAgg                                                                                                                                 >  2:2221814/1‑150 (MQ=255)
                   gTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAgg                                                                                                                                 >  2:2221791/1‑150 (MQ=255)
                   gTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAgg                                                                                                                                 >  1:1002889/1‑150 (MQ=255)
                        gCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGATAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTTCAGTGCATTTGAGATAGATCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTg                                                                                                                            >  1:37428/1‑150 (MQ=25)
                        gCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTg                                                                                                                            >  1:4506945/1‑150 (MQ=255)
                        gCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTg                                                                                                                            >  1:794257/1‑150 (MQ=255)
                          cTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAAGTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCACTgca                                                                                                                          <  1:3131530/150‑3 (MQ=21)
                          cTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGaa                                                                                                                          <  1:4439036/150‑1 (MQ=255)
                          cTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGaa                                                                                                                          <  1:4141527/150‑1 (MQ=255)
                           tCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGaaa                                                                                                                         >  2:3598356/1‑150 (MQ=255)
                             cccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACt                                                                                                                       >  2:3525407/1‑150 (MQ=255)
                             cccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACt                                                                                                                       >  1:4375729/1‑150 (MQ=255)
                             cccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACt                                                                                                                       >  1:2048171/1‑150 (MQ=255)
                             cccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACt                                                                                                                       >  1:1699180/1‑150 (MQ=255)
                              ccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTg                                                                                                                      >  2:2617426/1‑150 (MQ=255)
                                gCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTc                                                                                                                    >  1:3169814/1‑150 (MQ=255)
                                  tttGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGt                                                                                                                  <  2:1847941/150‑1 (MQ=255)
                                    tGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTgg                                                                                                                >  1:458719/1‑150 (MQ=255)
                                        tttGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGaa                                                                                                            <  2:2740324/150‑1 (MQ=255)
                                             gggTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACtt                                                                                                       >  2:1411757/1‑150 (MQ=255)
                                               gttggttgTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                                     <  1:1561509/150‑1 (MQ=255)
                                                ttggttgTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCAGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGc                                                                                                    >  2:2335108/1‑150 (MQ=255)
                                                   gttgttGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGc                                                                                                 >  1:660683/1‑150 (MQ=255)
                                                   gttgttGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGc                                                                                                 >  1:665156/1‑150 (MQ=255)
                                                   gttgttGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGc                                                                                                 >  2:1009189/1‑150 (MQ=255)
                                                              tGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGt                                                                                      >  2:1957327/1‑150 (MQ=255)
                                                              tGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGt                                                                                      <  1:126919/150‑1 (MQ=255)
                                                                gggAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATATGGTCt                                                                                    <  2:3150728/150‑1 (MQ=255)
                                                                   aggaACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCtttt                                                                                 <  2:1495929/150‑1 (MQ=255)
                                                                      tagtGTTTTGAGAGGGTGTTGGGTGTGAGTTGTGTCCGCCACAACCCTTGCTGTGGCTGTGAGATCGTGCGGGGCCCGCGCGGCGCTCGGTGTGTGAGGGGGTGAATGTGTCGTGGCGTACACTTGGCAGCCCCGATACGGTCTTTTccg                                                                              >  2:58896/4‑149 (MQ=21)
                                                                      aaCTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTccc                                                                              >  2:53725/1‑150 (MQ=255)
                                                                         tGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCtt                                                                           >  2:2767427/1‑150 (MQ=255)
                                                                          gTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTta                                                                          >  1:4213828/1‑150 (MQ=255)
                                                                           ttttGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTtat                                                                         <  1:1013309/150‑1 (MQ=255)
                                                                           ttttGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTtat                                                                         >  1:3581210/1‑150 (MQ=255)
                                                                           ttttGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTtat                                                                         >  1:3683118/1‑150 (MQ=255)
                                                                            tttGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTtata                                                                        <  1:3334280/150‑1 (MQ=255)
                                                                                agagGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTc                                                                    <  2:3555961/150‑1 (MQ=255)
                                                                                 gagGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCa                                                                   >  1:2132143/1‑150 (MQ=255)
                                                                                 gagGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCa                                                                   >  1:4005561/1‑150 (MQ=255)
                                                                                      tgtTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTc                                                              >  2:210721/1‑150 (MQ=255)
                                                                                       gtTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCt                                                             <  2:619523/150‑1 (MQ=255)
                                                                                          ggCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGa                                                          <  1:2229655/150‑1 (MQ=255)
                                                                                            cTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACGGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGaaa                                                        <  2:1800068/150‑1 (MQ=255)
                                                                                            cTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATCTCACTTTCTTGaaa                                                        <  2:4047513/150‑1 (MQ=255)
                                                                                            cTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCACGCGCGGCGCTCGGTCTCACGGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGaaa                                                        <  2:3081567/150‑1 (MQ=255)
                                                                                             tgtgAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAAt                                                       >  2:3978496/1‑150 (MQ=255)
                                                                                              gtgAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAAtt                                                      >  1:3649590/1‑150 (MQ=255)
                                                                                               tgAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTa                                                     >  2:1260895/1‑150 (MQ=255)
                                                                                               tgAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTa                                                     >  2:1901798/1‑150 (MQ=255)
                                                                                                   aTGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAgg                                                 >  2:982419/1‑150 (MQ=255)
                                                                                                     gtgtCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGtt                                               <  2:2129664/150‑1 (MQ=255)
                                                                                                     gtgtCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGtt                                               <  1:2226562/150‑1 (MQ=255)
                                                                                                     gtgtCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTAACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTACCCTTATATTTCACTTTCTTGAAATTAAAGGtt                                               >  1:1546169/1‑150 (MQ=255)
                                                                                                          cGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGacac                                          >  1:4322974/1‑150 (MQ=255)
                                                                                                           gCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACAcc                                         <  1:934736/150‑1 (MQ=255)
                                                                                                           gCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACAcc                                         >  2:1777309/1‑150 (MQ=255)
                                                                                                           gCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACAcc                                         >  2:1510126/1‑150 (MQ=255)
                                                                                                                     ttGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATc                               >  1:2907963/1‑150 (MQ=255)
                                                                                                                       gCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTATTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATccc                             <  1:3866835/150‑1 (MQ=255)
                                                                                                                        cAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATcccc                            <  2:458719/150‑1 (MQ=255)
                                                                                                                        cAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATcccc                            <  2:4180975/150‑1 (MQ=255)
                                                                                                                        cAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATcccc                            <  1:4369299/150‑1 (MQ=255)
                                                                                                                             ccTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTtata                       >  2:500402/1‑150 (MQ=255)
                                                                                                                               tgtgAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAAt                     >  2:255864/1‑150 (MQ=255)
                                                                                                                               tgtgAGATCGAGCGCCGCACGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGAGAACACTTGGCAGCACCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAAt                     >  2:257079/1‑150 (MQ=255)
                                                                                                                               tgtgAGACCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCACTGCAAGGGTTGTGGCGGACACATGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAAt                     >  1:2727909/1‑150 (MQ=11)
                                                                                                                                   agaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATgcgc                 <  1:4424538/150‑1 (MQ=255)
                                                                                                                                    gaccGAGCGCCGCGCGGGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGACTTTTCCCCTTATATTTCACTTTCTTGAAATGAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGcc                >  2:1806463/4‑150 (MQ=11)
                                                                                                                                     aTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGccc               >  2:4099930/1‑150 (MQ=255)
                                                                                                                                     aTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGccc               >  1:383239/1‑150 (MQ=255)
                                                                                                                                      tCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGcccc              >  1:647609/1‑150 (MQ=255)
                                                                                                                                       cGAGCGCCGCGCGGGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCa             >  2:4395729/1‑150 (MQ=21)
                                                                                                                                          gcgcCGCGCGGGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACtt          >  2:4140140/1‑150 (MQ=21)
                                                                                                                                              cgcccgcGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAACTTAAAGGTTGACACCCTCTTGAATCCCCTGATAATGCGCCCCACTTCCAg      >  2:2420195/1‑150 (MQ=255)
                                                                                                                                              cgcccgcGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAg      >  1:4208429/1‑150 (MQ=255)
                                                                                                                                              cgcccgcGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAg      >  2:1446394/1‑150 (MQ=255)
                                                                                                                                              cgcccgcGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGGTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAg      >  1:3735426/1‑150 (MQ=255)
                                                                                                                                               gcccgcGCGTCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGc     <  2:37428/150‑1 (MQ=255)
                                                                                                                                               gcccgcGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTAGATTTCACTTTCTTGAAATTAAAGGTTGCCACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGc     <  1:3525407/150‑1 (MQ=255)
                                                                                                                                                cccgcGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCg    >  2:3370764/1‑150 (MQ=255)
                                                                                                                                                cccgcGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCg    >  1:4081252/1‑150 (MQ=255)
                                                                                                                                                 ccgcGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGa   <  1:1330265/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGAC  >  NC_002947/697083‑697378

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