Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,155 C→T 13.2% intergenic (+109/‑204) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1550CT13.2% 125.0 / 9.2 53intergenic (+109/‑204)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (21/25);  new base T (3/4);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

CCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCC  >  NC_002947/955007‑955300
                                                                                                                                                       |                                                                                                                                                  
ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                                    >  7:219986/1‑149 (MQ=255)
  aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                   <  4:442/149‑1 (MQ=255)
           cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                           <  1:46782/147‑1 (MQ=255)
           cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         <  6:342508/149‑1 (MQ=255)
           cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                         <  5:103400/149‑1 (MQ=255)
            ggATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                        <  2:25950/149‑1 (MQ=37)
                 tgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑c                                                                                                                                  <  5:283484/149‑1 (MQ=255)
                      atgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑cc                                                                                                                                 <  5:135980/145‑1 (MQ=255)
                      atgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑cc                                                                                                                                 >  6:135980/1‑145 (MQ=255)
                      atgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCcg                                                                                                                               <  3:266175/147‑1 (MQ=255)
                           cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCAACAAAAAGGTTCCTGTAGGAGCGGGTTTA‑ACCGCGAagag                                                                                                                        >  6:127291/1‑149 (MQ=39)
                             acgcacgcTCACCCTGCTGCGTGTTGCACGAAAACCCCGGCCCCTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTGCGC‑GGGTAAACCCGCTCCTACACGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGc                                                                                                                      <  8:219986/149‑1 (MQ=37)
                                 acgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAg                                                                                                                   >  5:252928/1‑148 (MQ=255)
                                   gcTGAAACTGCTGCCTGTTTCACGTAAACCCCGGCCACTGG‑GCCGGGG‑TTTTTTTTTGGGGGTTTTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTGCTACAAGGGTCCTTTAGGCGGGGGTTTACCCCGCGAAGAGGCCAGca                                                                                                                 <  4:285913/149‑1 (MQ=37)
                                     tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGGGAAGAGGCCAGca                                                                                                                 >  8:237616/1‑146 (MQ=255)
                                                ccTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGTAAACCC‑CGGCGTCAACGGTGCTGTTGGGGGGGGGTTT‑ACCGCGGAGGGGGCCGCCCAGgaaaataaaaa                                                                                                   >  3:119694/1‑149 (MQ=255)
                                                ccTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTTAACCCGCACCGTCAAGGTTACGGTAAGAGCGGGGTTT‑ACCGCGAAGCCGCCAGGACAGaaaaaaaaaa                                                                                                    >  6:85887/1‑148 (MQ=255)
                                                ccTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTTCAAGGTTACTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGaaaaataaaaa                                                                                                   >  1:65502/1‑149 (MQ=255)
                                                  tGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGt                                                                                                                                       >  3:335190/1‑112 (MQ=21)
                                                  tGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGt                                                                                                                                       <  4:335190/112‑1 (MQ=21)
                                                          aaaaaCCCCGG‑CACAGTGGCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCGCCGACAAGGTTACGGTGGGGGCGGGGTTG‑ACCGCGCAGAGGCCCGCACAGCAAAATAAAAATAACTGGttt                                                                                         >  4:239102/2‑149 (MQ=25)
                                                             aaCCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAAACCCCAACTAAAAGGTTACTGTATGAGCGGGTTTA‑ACCGCGAAGAGGCCAGCACAGAAAAATAAAAAAAACAGAtttta                                                                                       >  6:125696/1‑147 (MQ=255)
                                                             aaCCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTAGC‑GGGTTAAACCCCAACAAAAAGGGTACGGGATGAGCCGGTTTT‑ACCGCCAAGAGGGCACCACAGAAAAATAAAAATAAATGAattta                                                                                       >  6:186162/1‑147 (MQ=255)
                                                             aaCCCCGG‑CACAGTGGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTATTCGC‑GGGTAAACCCACAACAAAAAAGTTACAGTAGGAACCGGTTTA‑ACCCCCAAGCAGCCAAAACACAAAAAAAAAAAAAACTGAattt                                                                                        >  6:275271/1‑148 (MQ=255)
                                                               ccccGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTCCAGGAACCTTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGc                                                                                                                      <  5:325392/115‑1 (MQ=18)
                                                                  cGGCCACTGT‑GCGGGGG‑TTTTTTTTTGGGGTTATTTTTCGTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACATGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAgg                                                                                                                       <  4:256942/111‑1 (MQ=25)
                                                                         tgt‑gCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTACTACAAGGTTACTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATAACTGATTTTCAAAAAATAGAtat                                                                          >  1:12795/1‑147 (MQ=255)
                                                                                     tttttttttGGGTTTATTTTTCCTGTGCCGGCCTATTCGC‑GGTAAAACCCGCTCATAAAATTTTACTTTAGTAGCGTTTTAA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                                                                 <  3:185613/148‑1 (MQ=255)
                                                                                             tGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTAAAAGGTTCCTTTAGGAGCGGTTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATaaaa                                                         <  7:237616/148‑1 (MQ=255)
                                                                                              gggggTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAAc                                                        <  6:252928/148‑1 (MQ=255)
                                                                                                   tATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa                                                                                 <  5:271721/118‑1 (MQ=255)
                                                                                                              tgCCGGCCTCTTCGC‑GGTTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTa                                        <  2:12795/148‑1 (MQ=255)
                                                                                                              tgCCGGCCTATTCGC‑TGTTAAACACTATCATAAAAGTTTCCTTTAGGAGCTTGTTTA‑ACCGCGAAGAGGACAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTa                                        <  5:275271/148‑1 (MQ=255)
                                                                                                                ccGGCCTCTTCGC‑GGTAAAACCCGATCAAAAAAGTTACCTTTAGTAGGGGTTTAA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACa                                      <  5:186162/148‑1 (MQ=255)
                                                                                                                         tCGG‑GGGTAAACCCGCTCCAACATGTTTCCTTTATGTGCGGTTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTg                             <  2:307597/148‑1 (MQ=255)
                                                                                                                                gTAAACCCGCTCCTACAAGGTTCCTTTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGc                      <  5:127291/149‑1 (MQ=255)
                                                                                                                                gTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg                       >  7:14301/1‑148 (MQ=255)
                                                                                                                                   aaCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACt                   <  8:14301/149‑1 (MQ=255)
                                                                                                                                    aCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGc                                   <  3:91699/132‑1 (MQ=255)
                                                                                                                                    aCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGc                                   >  4:91699/1‑132 (MQ=255)
                                                                                                                                     cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAt                 >  5:323248/1‑149 (MQ=37)
                                                                                                                                     cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAt                 >  1:197576/1‑149 (MQ=255)
                                                                                                                                      ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa                                                                                 <  7:28453/84‑1 (MQ=17)
                                                                                                                                      ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa                                                                                 >  8:28453/1‑84 (MQ=17)
                                                                                                                                      ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc                >  2:237664/1‑149 (MQ=255)
                                                                                                                                        gCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtctt              >  5:309989/1‑149 (MQ=255)
                                                                                                                                         cTCCTACATGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc             <  5:85887/149‑1 (MQ=255)
                                                                                                                                          tCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTAc                                       <  6:337757/122‑1 (MQ=255)
                                                                                                                                          tCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTAc                                       >  5:337757/1‑122 (MQ=255)
                                                                                                                                          tCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATc                                <  6:237122/129‑1 (MQ=255)
                                                                                                                                          tCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATc                                >  5:237122/1‑129 (MQ=255)
                                                                                                                                          tCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCAGCACTAtctact            >  8:294176/1‑149 (MQ=255)
                                                                                                                                                aaGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt        >  1:236839/1‑147 (MQ=255)
                                                                                                                                                     tCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGgcc  <  3:239102/148‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
CCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATTCAGGCC  >  NC_002947/955007‑955300

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: ‑