Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 170,839 T→G 16.1% 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,8390TG16.1% 67.2 / 8.1 37intergenic (+20/‑550)PP_0159/PP_16SAfamily III CoA‑transferase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (17/14);  new base G (4/2);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

GCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGT‑‑GTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGG  >  NC_002947/170699‑170986
                                                                                                                                              |                                                                                                                                                   
gCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGTAAAGTGGGGATGGATGATGGCGAGGTAAAGCGTATGGGTGATGATGGGGTGATATAGATGAAAAGCGGGGAGGCCGTCAGGcg                                                                                                                                               >  1:307343/1‑148 (MQ=255)
     tGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGT‑‑GTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGc                                                                                                                                        >  1:67039/1‑149 (MQ=255)
          tATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGT‑‑GTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTccc                                                                                                                                           >  7:151923/1‑141 (MQ=255)
          tATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGT‑‑GTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCTGACGGCCTccc                                                                                                                                           <  8:151923/141‑1 (MQ=255)
          tATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGT‑‑GTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGGCACGCCTCCCCCCGTTGt                                                                                                                                   >  2:162438/1‑149 (MQ=255)
                       caccacTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGT‑‑GTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAATAGCGGGGAGGACGTCAGGCCTGCCCGCTTTTTTTTTAAGGGCTa                                                                                                                       >  7:355816/1‑148 (MQ=255)
                                      agaaTACCGATGTGGTGGTGGGGGATGT‑‑GTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGaaa                                                                                                       <  2:220030/146‑1 (MQ=255)
                                           aCCGAAGTGGTGCTGGGGGAAGT‑‑GTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTGAGGGCTCCCCGCGTTGTTTTTGAGCGCGTTCAAATGTTATTGATATTa                                                                                                    >  8:332352/1‑147 (MQ=255)
                                           aCCGAAGTGGTGCTGGGGGAAGT‑‑GTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCg                                                                                                                                         >  1:209707/1‑110 (MQ=255)
                                           aCCGAAGGGGTGCGGGGGGAAGT‑‑GTGGGGGCTGGATGCTGGCGAGGTACAGCGTCGGCGTGATGCGGGGGGGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCg                                                                                                                                         <  2:209707/110‑1 (MQ=255)
                                             cGAAGTGGTGCTGGGGGAAGT‑‑GTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAgg                                                                                                >  6:85273/1‑149 (MQ=255)
                                                  gggtgCGGGGGCATGG‑‑GTTGGGGCTGGATGCTGGCGCGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGAc                                                                                           <  7:290063/148‑1 (MQ=39)
                                                            ggAAGT‑‑GTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGcc                                                                                        <  1:123735/142‑1 (MQ=255)
                                                            ggAAGT‑‑GTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGcc                                                                                        >  2:123735/1‑142 (MQ=255)
                                                                 t‑‑gtgGGGGCTGTAGGGTGGCGCGGTACAGCGTCTGCGGGATGCTGGGGTGGTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATccc                                                                            <  8:355816/149‑1 (MQ=38)
                                                                 t‑‑gtTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATccc                                                                            >  6:134142/1‑149 (MQ=255)
                                                                     ttGGGGGTGCATGCTGGGGAGTAACAGCGTCTGGGTTATGCGGGGGGGCTCTAGATGAAACGGGGGCAGCCCGTCTGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcccc                                                                           <  1:150653/148‑1 (MQ=255)
                                                                     tgggggCTTGATGCTGGCGAGGTACAGCTGCTGCGGGATGCTGGGGTGCTCTATATGAAAAGCGGGTAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATcccc                                                                           <  2:307343/146‑1 (MQ=255)
                                                                         ggCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTAAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTtat                                                                       >  4:116938/1‑148 (MQ=255)
                                                                           cTGGATGCTGGCGAGGTCCTGCGTCTGCGTTATGCTGGGGGGCTATAGATTAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATccc                                                                            <  5:37720/141‑1 (MQ=39)
                                                                           cTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGGGGCCGTCAGGCCTCCCCGCGTTTTTTTTAAGCGCGATAAAAGGTAATTTTAATTAAAAGTTTACGCCCACTCGAATaccccaaacg                                                                     >  6:37720/1‑142 (MQ=255)
                                                                                   tGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCACGCCTCCCCGCTTTGTTTTTAAGCGCTATCAAAGGTAATTTAAATTAAAAGTTGACGCCCTCTCGAATCCCCTTATAATGCGaccc                                                             >  4:143739/1‑148 (MQ=39)
                                                                                                      cGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATACCCTTATAATGCGCCCCAATTACAGCGACATCGGaa                                          >  5:94487/1‑148 (MQ=255)
                                                                                                          aTGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGaaa                                     <  3:116938/149‑1 (MQ=255)
                                                                                                            gCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGaaaa                                    >  7:180924/1‑148 (MQ=255)
                                                                                                             cTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAAc                                  >  4:83469/1‑149 (MQ=255)
                                                                                                                  gggcactAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCtt                             <  5:293450/144‑1 (MQ=255)
                                                                                                                   tGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGACAAACTCCTTg                            <  5:134142/149‑1 (MQ=255)
                                                                                                                        tAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCGCTCGACTCCCCGTATAATGCGCCCCACTTGCAGCGACATCGGAACGAAAAACTGCTTGAGAtt                       >  5:207019/1‑149 (MQ=25)
                                                                                                                            ttaaaaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCaa                    <  7:332352/146‑1 (MQ=255)
                                                                                                                                aaGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGt               >  4:91691/1‑149 (MQ=255)
                                                                                                                                   cGGGGAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGtt              >  7:295580/1‑147 (MQ=255)
                                                                                                                                   cGGGGAGGCCTGACGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGtt              <  8:295580/147‑1 (MQ=255)
                                                                                                                                    ggggAGGCCTGACGGCCTCCCCGCTTTGTTTTTAAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTACCGCGACATAGGAACGAAAAACTCCTTGAGATTAAATGAGTTAag           >  2:185976/1‑149 (MQ=18)
                                                                                                                                    ggggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAag           <  3:83469/149‑1 (MQ=255)
                                                                                                                                           ccGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGtt      >  4:348459/1‑147 (MQ=255)
                                                                                                                                             gTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAgg  <  5:85273/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GCCAGTGGAGTATCGCCGGGCGCCACCACTGCTGGGCGAGCATACCGAAGTGGTGCTGGGGGAAGT‑‑GTTGGGGCTGGATGCTGGCGAGGTACAGCGTCTGCGTGATGCTGGGGTGCTCTAGATGAAAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTCAGCGCTATCAAAGGTAATTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCAGCGACATCGGAACGAAAAACTCCTTGAGATTCAATGAGTTAAGTAGTTAAGG  >  NC_002947/170699‑170986

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