Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 524,352 G→C 9.8% intergenic (+216/‑600) tyrS → / → PP_16SC tyrosine‑‑tRNA ligase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947524,3520GC 9.8% 152.4 / 7.6 51intergenic (+216/‑600)tyrS/PP_16SCtyrosine‑‑tRNA ligase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (22/24);  new base C (2/3);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

GGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCA  >  NC_002947/524209‑524491
                                                                                                                                               |                                                                                                                                           
ggtTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCgcgggcg                                                                                                                                       <  2:434325/150‑1 (MQ=255)
    gCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGgcgc                                                                                                                                   <  1:667274/150‑1 (MQ=255)
      tgttgtAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                 >  2:523628/1‑150 (MQ=255)
         tgtAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGt                                                                                                                              <  2:813803/150‑1 (MQ=255)
                ggagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAGACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATAg                                                                                                                       <  2:1187081/150‑1 (MQ=25)
                ggagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAGACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATAg                                                                                                                       <  2:583693/150‑1 (MQ=25)
                ggagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAg                                                                                                                       >  2:632897/1‑150 (MQ=255)
                    ggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGgcgc                                                                                                                   >  1:737038/1‑150 (MQ=255)
                      ggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCAAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTg                                                                                                                 >  2:744913/1‑150 (MQ=255)
                       gCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGa                                                                                                                >  1:1110312/1‑150 (MQ=255)
                       gCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGa                                                                                                                >  1:849114/1‑150 (MQ=255)
                         tGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGaaa                                                                                                              <  1:1033975/150‑1 (MQ=255)
                              tGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTc                                                                                                         <  1:632897/150‑1 (MQ=255)
                                ggagGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGTTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCATGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGt                                                                                                       <  1:743490/149‑1 (MQ=255)
                                agagGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGt                                                                                                       <  2:817875/150‑1 (MQ=255)
                                        ttGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGtctc                                                                                                                           >  2:1155184/1‑122 (MQ=255)
                                        ttGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGtctc                                                                                                                           <  1:1155184/122‑1 (MQ=255)
                                        ttGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAaca                                                                                               >  1:749938/1‑150 (MQ=255)
                                        ttGGCGGTGAGATGTGTCCGCCACAACCCTCGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGCCGGCTCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAaca                                                                                               <  1:802488/150‑1 (MQ=255)
                                         tGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTATCGAGGctgccaa                                                                                              <  1:523628/150‑7 (MQ=255)
                                            cGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTg                                                                                           <  2:737038/150‑1 (MQ=255)
                                             ggTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCAAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                          <  1:744913/150‑1 (MQ=255)
                                                   aTGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCt                                                                                    >  2:675995/1‑150 (MQ=255)
                                                      tgtCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGa                                                                                 <  1:213612/150‑1 (MQ=255)
                                                        tCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATa                                                                               >  1:307117/1‑150 (MQ=255)
                                                           gCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACgg                                                                            >  1:110953/1‑150 (MQ=255)
                                                            ccACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAacac                                                                                              <  1:1162818/131‑1 (MQ=255)
                                                            ccACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAacac                                                                                              >  2:1162818/1‑131 (MQ=255)
                                                            ccACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTc                                                                                   >  2:929185/1‑142 (MQ=255)
                                                            ccACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTc                                                                                   <  1:929185/142‑1 (MQ=255)
                                                                    cTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTccc                                                                   >  2:173282/1‑150 (MQ=255)
                                                                    cTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTGccc                                                                   <  2:352966/150‑1 (MQ=255)
                                                                                  gagaTCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACtt                                                     >  1:1239973/1‑150 (MQ=255)
                                                                                    gaTCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTc                                                   <  2:110953/150‑1 (MQ=255)
                                                                                     aTCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCt                                                  >  2:1255631/1‑150 (MQ=255)
                                                                                     aTCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCt                                                  <  1:1255631/150‑1 (MQ=255)
                                                                                           cgccgcCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTCCCCCTTATATTTCACTTTCTTGAAAt                                            <  1:675995/150‑1 (MQ=255)
                                                                                                cccgcGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTaa                                         <  2:1140195/148‑1 (MQ=255)
                                                                                                cccgcGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTaa                                         >  1:1140195/1‑148 (MQ=255)
                                                                                                           cTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCct                            >  1:353201/1‑150 (MQ=255)
                                                                                                                gaCATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCctc                           <  1:709334/146‑1 (MQ=255)
                                                                                                                gaCATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCctc                           >  2:709334/1‑146 (MQ=255)
                                                                                                                 aCATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTcccc                                                                  <  2:721995/106‑1 (MQ=255)
                                                                                                                 aCATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTcccc                                                                  >  1:721995/1‑106 (MQ=255)
                                                                                                                        ccGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCtt               <  2:849114/150‑1 (MQ=255)
                                                                                                                           ccTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTatat                                                            <  1:712285/102‑1 (MQ=255)
                                                                                                                           ccTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTatat                                                            >  2:712285/1‑102 (MQ=255)
                                                                                                                                 aaaCCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATgcgc      >  2:68871/1‑150 (MQ=255)
                                                                                                                                     cGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCa  <  2:749938/150‑1 (MQ=255)
                                                                                                                                     cGAGCGCCGCCCGCGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCa  >  1:713515/1‑150 (MQ=16)
                                                                                                                                      gAGCGCCGCCCGCGCGGCGCTCGGTTTCATAGGCGGTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCctct                          <  2:610024/125‑1 (MQ=11)
                                                                                                                                      gAGCGCCGCCCGCGCGGCGCTCGGTTTCATAGGCGGTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCctct                          >  1:610024/1‑125 (MQ=11)
                                                                                                                                               |                                                                                                                                           
GGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCA  >  NC_002947/524209‑524491

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