Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,152 T→A 13.3% intergenic (+106/‑207) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1520TA13.3% 118.8 / 11.6 53intergenic (+106/‑207)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (20/23);  new base A (3/4);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

CCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT  >  NC_002947/955007‑955294
                                                                                                                                                    |                                                                                                                                               
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                                                                                                                           >  6:135980/1‑145 (MQ=255)
                      atgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑cc                                                                                                                           <  5:135980/145‑1 (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                                                                                                                                 <  4:335190/112‑1 (MQ=21)
                                                  tGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGt                                                                                                                                 >  3:335190/1‑112 (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                             >  4:91699/1‑132 (MQ=255)
                                                                                                                                    aCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGc                             <  3:91699/132‑1 (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)
                                                                                                                                                    |                                                                                                                                               
CCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGG‑TTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGC‑GGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTA‑CCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT  >  NC_002947/955007‑955294

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