Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 697,228 C→G 12.0% intergenic (+152/‑599) 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,2280CG12.0% 400.9 / 37.6 133intergenic (+152/‑599)mmsA‑I/PP_16SDmethylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (68/49);  new base G (9/7);  total (77/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 145 total aligned reads displayed.

GATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGA  >  NC_002947/697085‑697377
                                                                                                                                               |                                                                                                                                                     
gATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGcgggcg                                                                                                                                                 <  2:1788554/150‑1 (MQ=255)
 aTGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCgcgg                                                                                                                                                <  1:4956671/150‑1 (MQ=255)
       cGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                                                          >  1:3458872/1‑150 (MQ=255)
       cGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCg                                                                                                                                          <  2:3907671/150‑1 (MQ=255)
        ggggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCgg                                                                                                                                         <  2:2919915/150‑1 (MQ=255)
        ggggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGGTGCTGTTGGGGAGGCACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCga                                                                                                                                         <  2:3041709/150‑2 (MQ=21)
         gggCGGCCGACAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGt                                                                                                                                        <  1:3426402/150‑1 (MQ=255)
         gggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGTGATCGAGCGCCGCCCGCGCGGCGCTCGGt                                                                                                                                        >  2:4852856/1‑150 (MQ=255)
         gggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGt                                                                                                                                        >  1:2391090/1‑150 (MQ=255)
         gggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGt                                                                                                                                        >  1:5119747/1‑150 (MQ=255)
          ggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGtc                                                                                                                                       <  2:937996/150‑1 (MQ=255)
           gAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGtct                                                                                                                                      <  2:5436782/150‑1 (MQ=255)
           gAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                      <  2:778283/150‑1 (MQ=255)
              gCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTcac                                                                                                                                   <  2:437853/150‑1 (MQ=255)
               ccGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTcaca                                                                                                                                  <  1:2727297/150‑1 (MQ=255)
                 gTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAgg                                                                                                                                >  1:5173768/1‑150 (MQ=255)
                  gacggCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGgc                                                                                                                               >  1:2235765/4‑150 (MQ=255)
                    aGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGgcgc                                                                                                                             <  2:2983137/150‑1 (MQ=255)
                        cTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGaa                                                                                                                         >  2:4677135/1‑150 (MQ=255)
                         tCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGaaa                                                                                                                        >  1:5082501/1‑150 (MQ=255)
                         tCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGaaa                                                                                                                        >  1:4003664/1‑150 (MQ=255)
                         gccccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGaaa                                                                                                                        <  1:5031806/149‑1 (MQ=255)
                           cccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACt                                                                                                                      >  1:3086135/1‑150 (MQ=255)
                           cccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACt                                                                                                                      >  2:1361583/1‑150 (MQ=255)
                           cccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACt                                                                                                                      <  1:4685679/150‑1 (MQ=255)
                            ccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGcgggcg                                                                                                                                                 <  2:2530984/122‑1 (MQ=255)
                            ccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGcgggcg                                                                                                                                                 >  1:2530984/1‑122 (MQ=255)
                             cGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGt                                                                                                                    >  2:3717081/1‑150 (MQ=255)
                                tttGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                      >  1:3255260/1‑129 (MQ=255)
                                tttGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                                      <  2:3255260/129‑1 (MQ=255)
                                 ttGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTg                                                                                                                <  2:1936550/150‑1 (MQ=255)
                                       ttGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAac                                                                                                          >  1:859534/1‑150 (MQ=255)
                                          ggggTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACt                                                                                                       >  1:3854623/1‑150 (MQ=255)
                                           gggTTGGTTGTTGCTGTTGGGGAGGATCTGTTGTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACtt                                                                                                      <  2:197552/150‑1 (MQ=255)
                                           gggTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACtt                                                                                                      <  2:5066367/150‑1 (MQ=255)
                                             gttggttgTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                                    >  1:5224467/1‑150 (MQ=255)
                                              ttggttgTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGc                                                                                                   >  1:2270293/1‑150 (MQ=255)
                                              ttggttgTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGc                                                                                                   >  2:4680013/1‑150 (MQ=255)
                                                ggttgTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAg                                                                                                 <  1:4409605/150‑1 (MQ=255)
                                                  ttgttgCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGcc                                                                                               >  1:1568402/1‑150 (MQ=255)
                                                                    aaCTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTccc                                                                             >  1:2819356/1‑150 (MQ=255)
                                                                    aaCTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTccc                                                                             <  1:3285149/150‑1 (MQ=255)
                                                                     aCTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTcccc                                                                            <  2:329042/150‑1 (MQ=255)
                                                                     aCTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCgat                                                                                                                                        >  2:4262219/1‑88 (MQ=17)
                                                                     aCTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCgat                                                                                                                                        <  1:4262219/90‑3 (MQ=17)
                                                                         ttttGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTtat                                                                        <  1:896404/150‑1 (MQ=255)
                                                                              agagGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTc                                                                   >  2:2557534/1‑150 (MQ=255)
                                                                               gagGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGccgcgcgg                                                                                                                                                    >  1:4991487/1‑67 (MQ=255)
                                                                                agGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCAc                                                                 >  1:226736/1‑150 (MQ=255)
                                                                                    tgtTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTc                                                             >  1:2929902/1‑150 (MQ=255)
                                                                                      ttGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCtt                                                           >  2:4847868/1‑150 (MQ=255)
                                                                                       tGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTATTTTCCCCTTATATTTCACTTTCTTg                                                          <  2:771647/150‑1 (MQ=255)
                                                                                         gCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCtt                                                                          <  2:4182343/132‑1 (MQ=255)
                                                                                         gCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCtt                                                                          >  1:4182343/1‑132 (MQ=255)
                                                                                         gCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGaa                                                        <  1:1582139/150‑1 (MQ=255)
                                                                                            gtgAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAAtt                                                     >  1:2239684/1‑150 (MQ=255)
                                                                                            gtgAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAAtt                                                     >  2:3306303/1‑150 (MQ=255)
                                                                                               agctgtgTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTaaa                                                  <  2:577302/147‑1 (MQ=255)
                                                                                               agaTGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTaaa                                                  <  2:1434888/150‑1 (MQ=255)
                                                                                                gaTGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAg                                                 <  2:2393542/150‑1 (MQ=255)
                                                                                                gaTGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAg                                                 <  2:5082501/150‑1 (MQ=255)
                                                                                                     gtCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGa                                            <  1:4529581/150‑1 (MQ=255)
                                                                                                      tCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTaa                                                   >  2:4990275/1‑142 (MQ=255)
                                                                                                      tCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGac                                           >  1:4569236/1‑150 (MQ=255)
                                                                                                         gCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACAcc                                        >  2:2254578/1‑150 (MQ=255)
                                                                                                         gCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAAGGCAAGGTTGACAcc                                        >  2:1775328/1‑150 (MQ=255)
                                                                                                         gCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGACTCACAGGCGCAGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCATATATTTCACTTTCTTGAAATTAATGGTTGACAcc                                        >  2:4010915/1‑150 (MQ=255)
                                                                                                                   ttGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCGGCGATACGGTCTTTTCGCCTTATATTTCAGTTTGTTGAAATTATAGGTTGTCTCGGTCTTGTATc                              >  2:3787054/1‑150 (MQ=17)
                                                                                                                   ttGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATc                              >  2:4514930/1‑150 (MQ=255)
                                                                                                                   ttGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATc                              <  2:5004174/150‑1 (MQ=255)
                                                                                                                    tGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATcc                             <  2:2819356/150‑1 (MQ=255)
                                                                                                                    tGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATcc                             <  2:226736/150‑1 (MQ=255)
                                                                                                                      cAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATcccc                           <  1:4429964/150‑1 (MQ=255)
                                                                                                                       aGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCt                          >  1:4119919/1‑150 (MQ=255)
                                                                                                                         tGCCTGTGTGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTta                        <  1:4852856/150‑1 (MQ=255)
                                                                                                                         tGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTta                        <  2:1592593/150‑1 (MQ=255)
                                                                                                                          gCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTGATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCGtat                       >  1:518945/1‑150 (MQ=11)
                                                                                                                           ccTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTtata                      <  1:2557534/150‑1 (MQ=255)
                                                                                                                           ccTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTtata                      >  2:913087/1‑150 (MQ=255)
                                                                                                                            cTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATaa                     >  1:4978024/1‑150 (MQ=255)
                                                                                                                             tgtgAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAAt                    >  1:223986/1‑150 (MQ=255)
                                                                                                                             tgtgAGACCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAAt                    >  1:3312371/1‑150 (MQ=255)
                                                                                                                               tgAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATgc                  >  1:2071974/1‑150 (MQ=255)
                                                                                                                                gagaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATgcg                 <  2:2391090/150‑1 (MQ=255)
                                                                                                                                gagaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATgcg                 <  2:5119747/150‑1 (MQ=255)
                                                                                                                                gagaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTATTGAAATTAAAGGTTGACACCCTCTTGAATCACATTATAATgcg                 >  1:1700691/1‑150 (MQ=255)
                                                                                                                                 agaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATgcgc                >  1:2868163/1‑150 (MQ=255)
                                                                                                                                 agaCCGAGCGCCGCGCGGGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATgcgc                <  2:518945/150‑1 (MQ=16)
                                                                                                                                  gaccGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGcc               >  2:5031183/4‑150 (MQ=255)
                                                                                                                                   accGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGccc              >  1:2718291/3‑150 (MQ=255)
                                                                                                                                   aTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGccc              <  1:3237249/150‑1 (MQ=255)
                                                                                                                                    tCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGcccc             >  2:5206686/1‑150 (MQ=255)
                                                                                                                                       aGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACt          >  2:3050195/1‑150 (MQ=255)
                                                                                                                                        gcgcCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACtt         >  1:140217/1‑150 (MQ=255)
                                                                                                                                          gccgccCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTcc       >  2:2654041/1‑150 (MQ=255)
                                                                                                                                          gccgccCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTcc       >  1:2953506/1‑150 (MQ=255)
                                                                                                                                           ccgcgcggGCGGCGCTCGATCTCACAGGCACTGCAAGGGTTGTGGCGGACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCa      >  2:2836826/1‑150 (MQ=11)
                                                                                                                                           ccgcccgcGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCa      >  1:1872940/1‑150 (MQ=255)
                                                                                                                                           ccgcccgcGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCa      >  2:4952825/1‑150 (MQ=255)
                                                                                                                                               ccgcGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGa  >  2:4732244/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGA  >  NC_002947/697085‑697377

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: