Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 697,228 C→G 6.7% 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,2280CG 6.7% 238.9 / 6.4 75intergenic (+155/‑596)mmsA‑I/PP_16SDmethylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (36/34);  new base G (2/3);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

GATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGAC  >  NC_002947/697082‑697375
                                                                                                                                                     |                                                                                                                                                   
gtgcacgCGGGGAGGCGGTCTGGCGGCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCcgcgc                                                                                                                                                   <  5:40498/143‑1 (MQ=18)
   aaaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCgcggc                                                                                                                                                >  3:177308/2‑149 (MQ=255)
          ggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGt                                                                                                                                         <  4:39978/149‑1 (MQ=255)
           gAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTTCGCCACGACAGTTTCAGCGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                                        <  7:66246/149‑1 (MQ=18)
            aGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGtct                                                                                                                                       <  4:132362/149‑1 (MQ=255)
            aGGCCGTCAGGCCTCCCCGCTTTGTGTTTGGGGGTTGGTGGTTGCTGGTGGGGAGGAACTGGTTTGAGAGGGGGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGtct                                                                                                                                       <  5:125773/149‑1 (MQ=21)
              gCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTca                                                                                                                                     <  2:243393/149‑1 (MQ=255)
               ccGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTcac                                                                                                                                    <  4:263227/149‑1 (MQ=255)
                          ccccGCTTTGTTTTTGGGGGTTGGTTGGTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGaaa                                                                                                                         <  5:210379/149‑1 (MQ=255)
                           cccGGTTTGTTTTTGGGGGGTGGGTGGTGCTGGTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAAc                                                                                                                        <  1:358869/149‑1 (MQ=25)
                           cccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGCCTCATCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCt                                                                                                                             <  1:307643/147‑1 (MQ=25)
                           cccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGaa                                                                                                                          <  4:177308/147‑1 (MQ=255)
                           cccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAAc                                                                                                                        >  1:129462/1‑149 (MQ=255)
                            ccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACt                                                                                                                       >  2:157277/1‑149 (MQ=255)
                            ccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCACTGcaagg                                                                                                                       <  2:266233/149‑3 (MQ=21)
                               cTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTc                                                                                                                    <  8:177635/149‑1 (MQ=255)
                                  tGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTg                                                                                                                 >  5:304598/1‑149 (MQ=255)
                                      tttGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGt                                                                                                                  >  5:42313/1‑144 (MQ=255)
                                      tttGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGt                                                                                                                  <  6:42313/144‑1 (MQ=255)
                                            ggttggttGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACtt                                                                                                       <  2:98011/149‑1 (MQ=255)
                                             gttggttgTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTg                                                                                                      <  3:29676/149‑1 (MQ=255)
                                              ttggtttTGGCTGTTGGGGAGTAAGTGTTTTGTGTGGGTGTTGGCTGTGATATG‑‑TGTCCGCCAAAACCCTTGCTTTGCCTTTTTTATCGTGGGCCGCCCGCGGGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                                     <  2:129462/149‑1 (MQ=21)
                                              ttggttgTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCAGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                                     <  8:319491/149‑1 (MQ=255)
                                                 gttgttGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGc                                                                                                               >  4:311068/1‑136 (MQ=255)
                                                 gttgttGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGc                                                                                                               <  3:311068/136‑1 (MQ=255)
                                                 gttgttGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAg                                                                                                  >  4:305155/1‑149 (MQ=255)
                                                 gttgttGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAg                                                                                                  >  6:327111/1‑149 (MQ=255)
                                                  ttgttgCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGc                                                                                                 >  8:316634/1‑149 (MQ=255)
                                                  ttgttgCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGc                                                                                                 >  3:302974/1‑149 (MQ=255)
                                                         tGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATa                                                                                          <  6:304598/149‑1 (MQ=255)
                                                                gaggaACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCtt                                                                                   >  3:207679/1‑149 (MQ=255)
                                                                gaggaACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCtt                                                                                   >  4:6108/1‑149 (MQ=255)
                                                                   gaACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTc                                                                                >  8:152363/1‑149 (MQ=255)
                                                                         ttttGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTta                                                                          >  7:261798/1‑149 (MQ=255)
                                                                         ttttGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTta                                                                          >  3:296918/1‑149 (MQ=255)
                                                                             gagaGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATAtt                                                                      >  1:227331/1‑149 (MQ=255)
                                                                             gagaGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATAtt                                                                      >  2:359966/1‑149 (MQ=255)
                                                                               gagGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGgc                                                                                                                                >  4:107787/1‑89 (MQ=255)
                                                                               gagGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGgc                                                                                                                                <  3:107787/89‑1 (MQ=255)
                                                                                agGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCa                                                                   >  7:141593/1‑149 (MQ=255)
                                                                                agGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCa                                                                   >  7:136507/1‑149 (MQ=255)
                                                                                     gtTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTc                                                              >  1:80760/1‑149 (MQ=255)
                                                                                        ggCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCtt                                                            >  5:310253/1‑148 (MQ=255)
                                                                                        ggCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCtt                                                            <  2:227331/148‑1 (MQ=255)
                                                                                            gtgAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAAt                                                       >  2:227668/1‑149 (MQ=255)
                                                                                              gagaTG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTTTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAAtt                                                      >  5:177160/1‑148 (MQ=255)
                                                                                               agaTG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTaa                                                    <  1:359966/149‑1 (MQ=255)
                                                                                                gaTG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTaaa                                                   <  6:310253/149‑1 (MQ=255)
                                                                                                 aTG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAg                                                  >  8:168291/1‑149 (MQ=255)
                                                                                                      tgtCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTCTCCCCTTATATTTCACTTTCTTGAAATTAAAGGtt                                               <  7:316634/149‑1 (MQ=255)
                                                                                                                    ccTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                                     <  6:78821/82‑1 (MQ=255)
                                                                                                                    ccTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                                     >  5:78821/1‑82 (MQ=255)
                                                                                                                    ccTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTg                                   >  1:217749/1‑148 (MQ=255)
                                                                                                                    ccTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTg                                   >  1:98346/1‑148 (MQ=255)
                                                                                                                        gCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATc                               <  7:152363/148‑1 (MQ=255)
                                                                                                                         cAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATccc                             <  4:296918/149‑1 (MQ=255)
                                                                                                                          aGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGaaa                                                        >  8:320024/1‑121 (MQ=255)
                                                                                                                          aGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGaaa                                                        <  7:320024/121‑1 (MQ=255)
                                                                                                                          aGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATcccc                            <  4:302974/149‑1 (MQ=255)
                                                                                                                          aGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATcccc                            <  2:98346/149‑1 (MQ=255)
                                                                                                                          aGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATcccc                            <  4:207679/149‑1 (MQ=255)
                                                                                                                            tGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCtt                          <  8:141593/149‑1 (MQ=255)
                                                                                                                            tGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCtt                          <  8:136507/149‑1 (MQ=255)
                                                                                                                             gCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACt                                                                 <  6:315900/109‑1 (MQ=255)
                                                                                                                             gCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACt                                                                 >  5:315900/1‑109 (MQ=255)
                                                                                                                               cTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTtat                        >  2:223355/1‑148 (MQ=255)
                                                                                                                                     gaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATgcgc                 <  1:223355/149‑1 (MQ=255)
                                                                                                                                     gaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATgcgc                 <  3:305155/149‑1 (MQ=255)
                                                                                                                                      accGAGCGCCGCGCGGGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGACCACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGcc                >  5:348216/3‑149 (MQ=255)
                                                                                                                                      aTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGcc                >  3:173003/1‑149 (MQ=255)
                                                                                                                                       tCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCATTATAATgcgcac               >  4:235208/1‑147 (MQ=255)
                                                                                                                                        cGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCctct                                     <  5:342803/126‑1 (MQ=12)
                                                                                                                                        cGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCctct                                     >  6:342803/1‑126 (MQ=12)
                                                                                                                                        cGAGCGCCGCCCGCGCGGCGATCGGTATCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCGGCTTATGTTTCACTTTGTTGAAATTAAAGGTTGACAGCCTCTTGAATCCCAGTATAATgcgcgc               >  4:232806/1‑146 (MQ=14)
                                                                                                                                         gAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGCATCCCCTTATAATGCGCCCCa             >  6:291251/1‑149 (MQ=255)
                                                                                                                                            cgccgcCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACtt          >  3:24000/1‑149 (MQ=255)
                                                                                                                                                 cccgcGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGc     <  5:327111/149‑1 (MQ=255)
                                                                                                                                                   cgcgcgGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCATCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGa   >  6:371817/1‑149 (MQ=255)
                                                                                                                                                     cgcgGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGAc  >  5:232946/1‑148 (MQ=40)
                                                                                                                                                     |                                                                                                                                                   
GATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATG‑‑TGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGCCCCACTTCCAGCGAC  >  NC_002947/697082‑697375

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