Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 170,839 T→G 12.5% intergenic (+20/‑550) PP_0159 → / → PP_16SA family III CoA‑transferase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947170,8390TG12.5% 123.8 / 9.7 55intergenic (+20/‑550)PP_0159/PP_16SAfamily III CoA‑transferase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (26/22);  new base G (5/2);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-01

GAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAG  >  NC_002947/170694‑170977
                                                                                                                                                  |                                                                                                                                          
gagaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGAGCTCTAGATGAAAAGCGGGGAGGCCGTCAg                                                                                                                                         <  1:31590/149‑1 (MQ=255)
  gaCGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGc                                                                                                                                       <  8:230049/149‑1 (MQ=255)
     gCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGCGGCCGTCAGGCCTc                                                                                                                                    >  1:172278/1‑149 (MQ=255)
     gCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGGGTTGGGGCTGGATGCTGGCGAGGTACAGCGGCTGCGTGAGGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGGCAGGCCTc                                                                                                                                    >  5:289559/1‑149 (MQ=255)
       cAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTccc                                                                                                                                  <  4:215131/149‑1 (MQ=255)
        aGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGGGGTGCGGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCTTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTcccc                                                                                                                                 <  8:288625/149‑1 (MQ=255)
        aGTGCAGTATCCCCGGCCGCCACCACTGCGGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGGGTTGGGGCGGGATGCGGGCGCGGTACAGCGTCTGCGTGATGCTGGGGGGCTCTAGATGAAACGCGGGGAGGCCTGACGGCCTcccc                                                                                                                                 <  4:65354/149‑1 (MQ=255)
           gtaGTTTCGCGGGGCGCACCCAGTGGTGGGGGAGCATACCGATGGGGTGCTGGGGTATGTGTGGGGGCTGGATGCTGTCGAGGTCCAGCGTCGGGGTGTGGCTGGGGTGCGCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCt                                                                                                                              <  8:197504/147‑1 (MQ=255)
           ggAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCt                                                                                                                              >  2:126095/1‑149 (MQ=255)
                            caccacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTa                                                                                                              >  3:103255/1‑148 (MQ=255)
                             accccggctgGGCGAGCATACCGCCGTGGGGCTGGGGGAAGTGGTGGGGCTGGAGGCGGGCGAGGTACAGCGTCTGCGTGATGGTGGGGGGCGCTAGCTGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATc                                                                                                            <  6:289559/143‑1 (MQ=37)
                                                aCCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCGATCAAAGGTAATTGAAAttt                                                                                           >  3:240489/1‑146 (MQ=255)
                                                  ctaaGTGGTGGTGGGGTATGTGTGGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTTATGCTGGGTTGCTCTATATCAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAgg                                                                                       <  5:203582/147‑1 (MQ=38)
                                                  cGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGttttt                                                                                                                      >  3:229882/1‑118 (MQ=255)
                                                          tgCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGcc                                                                               >  1:51937/1‑149 (MQ=255)
                                                            cTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGTCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGccc                                                                              >  5:113374/1‑148 (MQ=255)
                                                               ggggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTAAGCGGTATGAAAGGTAATTGAAATTTAAGGTTGACGCCctcgc                                                                          >  4:152885/1‑147 (MQ=255)
                                                               ggggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTAAGGCCAACCCGCTTTTTTTTTAAGCGCTATAAAAGGTAATTGAAATTAAAGGTTTACGCCCtgtc                                                                          >  8:56011/1‑149 (MQ=255)
                                                                 ggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGGGGGCGGCAGGCCTCTCCGCGTTGTTTTTGAGCGCGATAAAAGGTTATTGGAATTAAAGGTTTAAGCCCtctct                                                                         >  3:10383/1‑147 (MQ=255)
                                                                 ggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCtctc                                                                          <  8:107704/147‑1 (MQ=255)
                                                                 ggAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCg                                                                         >  7:107704/1‑148 (MQ=255)
                                                                  gAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCaaa                                                                                                         >  1:259286/1‑115 (MQ=255)
                                                                  gAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCaaa                                                                                                         <  2:259286/115‑1 (MQ=255)
                                                                    agggtTGGGGCTGGAGGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGTGGCTCTATATTAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCctct                                                                           <  4:327261/140‑1 (MQ=39)
                                                                    aGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCctct                                                                           >  3:327261/1‑143 (MQ=255)
                                                                      tgtTGGGGCTGTAGGCTGCCGAGGTACAGGGTCTGGGTGAGGCGGGGGTGCTCTAGATCAAACGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcc                                                                    <  4:240489/148‑1 (MQ=38)
                                                                        ttGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCt                                                                 >  3:344152/1‑149 (MQ=255)
                                                                          ggggCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCATTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTTAAGGTTGAAGCGCTCTCGAATCCCCTta                                                               >  4:239520/1‑149 (MQ=255)
                                                                             gCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTtata                                                             >  2:65402/1‑148 (MQ=255)
                                                                             gCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTtata                                                             <  1:65402/148‑1 (MQ=255)
                                                                                  aTGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATgcgc                                                       <  4:103255/149‑1 (MQ=255)
                                                                                   tGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGcc                                                      >  7:113316/1‑149 (MQ=255)
                                                                                     cTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGCGGCCGTCAGGCCTCCCCGCGTTGTTTTTGAGCGCTTTCAAAGGTAATTGAAATTAAAGGTTGACGCCCTATCGAATCCCCTTATAATGCGccc                                                     >  3:61032/1‑148 (MQ=38)
                                                                                      tGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCACCCCGCTTTGTTTTTCAGCGCTATAAAAGGTAATTGAAATTTAACGTTGAAGCCCGCTCGAATCCCCTTATAATGCGccccc                                                   >  8:325115/1‑148 (MQ=39)
                                                                                        gCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACt                                                 <  4:344152/149‑1 (MQ=255)
                                                                                               aCAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATAAAAGGTAATTTTAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGAGCCCCACTTCCAGc                                           >  3:62192/1‑148 (MQ=255)
                                                                                                   cGTCTTCGTGATGGGGGGGTGCCTATTAAGAAAAGCGGTGAGGGCGTGAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACAt                                      <  7:79794/149‑1 (MQ=21)
                                                                                                          gttaTGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGc                                           <  1:341015/134‑1 (MQ=255)
                                                                                                          gTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCGTTGTTTTTAAGGGCTATAAAAGGTAATTTAAATTAAAGGTTGACGCCCGCTCGAATGCCCTTATAATGAGCCCCACATCCAga                                           >  2:341015/1‑136 (MQ=38)
                                                                                                           tGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTATCGAATCCCCGTATAATGCGCCCCACTTCCAGCGACATCGGAACg                               >  5:7547/1‑148 (MQ=255)
                                                                                                             aTGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGGCGTCAGGCCTACCCGCTTTGTTTTTCAGCGCTTTCAAAGGTAATTGTAATTAAAGGTTGACGGCCTGGGGCATGCCCGTATAATGCGGCCCACTTCGCGCGCCATCGGAACgta                             >  6:124868/1‑146 (MQ=255)
                                                                                                                      tGCTCTAGATTAAAAGCGGGCAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTg                   <  3:152885/149‑1 (MQ=255)
                                                                                                                        ctctAGATGAAAAGCGGGGAGGCCGTCATGCCGACCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGCATCCCCGTAAAATGAGCCCCCCTTACAGCGACATCGGAACGAAAAACTCCtt                    >  3:324003/1‑146 (MQ=25)
                                                                                                                           tAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGAt               >  1:250525/1‑148 (MQ=255)
                                                                                                                              gATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCa            <  2:30319/149‑1 (MQ=255)
                                                                                                                              gATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTGTGTTTTTCACCGCAATCAAAGGTAATTGAAAATAAAGGTTTACGCCCACTCGAATCCCCATATAATGCGCCCCACTTCCAGCGACATAGGAACGAAAAACAACTTGAGAttg             >  8:183597/1‑147 (MQ=21)
                                                                                                                               aaacaaaGCGGGGAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGCCGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa           >  6:265872/5‑149 (MQ=18)
                                                                                                                               aTGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa           <  2:250525/149‑1 (MQ=255)
                                                                                                                                ttaaaaGGGGGCAGGCCGTCAGGCCGCCCCGCTTTTTTCATCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa           <  1:126095/146‑1 (MQ=18)
                                                                                                                                tGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAAt          <  7:11399/149‑1 (MQ=255)
                                                                                                                                   aaaGCGGGGAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGt      >  3:66772/1‑149 (MQ=18)
                                                                                                                                   aaaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGa        <  1:114361/148‑1 (MQ=255)
                                                                                                                                    aaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGt      <  3:239520/149‑1 (MQ=255)
                                                                                                                                    aaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGAAGCCCTCTCGAATCCCCTTATAATGCACCCAACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGt      >  2:246827/1‑149 (MQ=255)
                                                                                                                                        ggggAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag  >  4:76201/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
GAGACGCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTGTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAG  >  NC_002947/170694‑170977

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