Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,149 A→G 14.2% intergenic (+103/‑210) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1490AG14.2% 121.4 / 9.0 53intergenic (+103/‑210)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (22/20);  new base G (3/4);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT  >  NC_002947/955003‑955294
                                                                                                                                                     |                                                                                                                                                  
aCAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                    <  8:284700/147‑1 (MQ=255)
aCAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                    >  7:284700/1‑147 (MQ=255)
aCAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAgg                                                                                                                                                  >  4:266175/1‑149 (MQ=255)
    ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                              >  7:219986/1‑149 (MQ=255)
      aGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                             <  4:442/149‑1 (MQ=255)
               cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                     <  1:46782/147‑1 (MQ=255)
               cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                   <  6:342508/149‑1 (MQ=255)
               cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                   <  5:103400/149‑1 (MQ=255)
                ggATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                  <  2:25950/149‑1 (MQ=37)
                     tgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAc                                                                                                                            <  5:283484/149‑1 (MQ=255)
                          atgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAcc                                                                                                                           <  5:135980/145‑1 (MQ=255)
                          atgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAcc                                                                                                                           >  6:135980/1‑145 (MQ=255)
                          atgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCcg                                                                                                                         <  3:266175/147‑1 (MQ=255)
                               cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCAACAAAAAGGTTCCTGTAGGAGCGGGTTTAACCGCGAagag                                                                                                                  >  6:127291/1‑149 (MQ=39)
                                 acgcacgcTCACCCTGCTGCGTGTTGCACGAAAACCCCGGCCCCTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTGCGCGGGTAAACCCGCTCCTACACGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                <  8:219986/149‑1 (MQ=37)
                                     acgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAg                                                                                                             >  5:252928/1‑148 (MQ=255)
                                       gcTGAAACTGCTGCCTGTTTCACGTAAACCCCGGCCACTGGGCCGGGGTTTTTTTTTGGGGGTTTTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTGCTACAAGGGTCCTTTAGGCGGGGGTTTACCCCGCGAAGAGGCCAGca                                                                                                           <  4:285913/149‑1 (MQ=37)
                                         tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGGGAAGAGGCCAGca                                                                                                           >  8:237616/1‑146 (MQ=255)
                                                    ccTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGGTAAACCCCGGCGTCAACGGTGCTGTTGGGGGGGGGTTTACCGCGGAGGGGGCCGCCCAGgaaaataaaaa                                                                                             >  3:119694/1‑149 (MQ=255)
                                                    ccTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTTAACCCGCACCGTCAAGGTTACGGTAAGAGCGGGGTTTACCGCGAAGCCGCCAGGACAGaaaaaaaaaa                                                                                              >  6:85887/1‑148 (MQ=255)
                                                    ccTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTTCAAGGTTACTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaataaaaa                                                                                             >  1:65502/1‑149 (MQ=255)
                                                      tGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGt                                                                                                                                 <  4:335190/112‑1 (MQ=21)
                                                      tGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGt                                                                                                                                 >  3:335190/1‑112 (MQ=21)
                                                              aaaaaCCCCGGCACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGCCGACAAGGTTACGGTGGGGGCGGGGTTGACCGCGCAGAGGCCCGCACAGCAAAATAAAAATAACTGGttt                                                                                   >  4:239102/2‑149 (MQ=25)
                                                                 aaCCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCAACTAAAAGGTTACTGTATGAGCGGGTTTAACCGCGAAGAGGCCAGCACAGAAAAATAAAAAAAACAGAtttta                                                                                 >  6:125696/1‑147 (MQ=255)
                                                                 aaCCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTAGCGGGTTAAACCCCAACAAAAAGGGTACGGGATGAGCCGGTTTTACCGCCAAGAGGGCACCACAGAAAAATAAAAATAAATGAattta                                                                                 >  6:186162/1‑147 (MQ=255)
                                                                 aaCCCCGGCACAGTGGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTATTCGCGGGTAAACCCACAACAAAAAAGTTACAGTAGGAACCGGTTTAACCCCCAAGCAGCCAAAACACAAAAAAAAAAAAAACTGAattt                                                                                  >  6:275271/1‑148 (MQ=255)
                                                                   ccccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                <  5:325392/115‑1 (MQ=18)
                                                                      cGGCCACTGTGCGGGGGTTTTTTTTTGGGGTTATTTTTCGTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACATGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                                                 <  4:256942/111‑1 (MQ=25)
                                                                             tgtgCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTACTACAAGGTTACTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATAACTGATTTTCAAAAAATAGAtat                                                                    >  1:12795/1‑147 (MQ=255)
                                                                                         tttttttttGGGTTTATTTTTCCTGTGCCGGCCTATTCGCGGTAAAACCCGCTCATAAAATTTTACTTTAGTAGCGTTTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                                                           <  3:185613/148‑1 (MQ=255)
                                                                                                 tGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAAAAGGTTCCTTTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATaaaa                                                   <  7:237616/148‑1 (MQ=255)
                                                                                                  gggggTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAAc                                                  <  6:252928/148‑1 (MQ=255)
                                                                                                       tATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa                                                                           <  5:271721/118‑1 (MQ=255)
                                                                                                                  tgCCGGCCTCTTCGCGGTTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTa                                  <  2:12795/148‑1 (MQ=255)
                                                                                                                  tgCCGGCCTATTCGCTGTTAAACACTATCATAAAAGTTTCCTTTAGGAGCTTGTTTAACCGCGAAGAGGACAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTa                                  <  5:275271/148‑1 (MQ=255)
                                                                                                                    ccGGCCTCTTCGCGGTAAAACCCGATCAAAAAAGTTACCTTTAGTAGGGGTTTAACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACa                                <  5:186162/148‑1 (MQ=255)
                                                                                                                             tCGGGGGTAAACCCGCTCCAACATGTTTCCTTTATGTGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTg                       <  2:307597/148‑1 (MQ=255)
                                                                                                                                    gTAAACCCGCTCCTACAAGGTTCCTTTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGc                <  5:127291/149‑1 (MQ=255)
                                                                                                                                    gTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCgg                 >  7:14301/1‑148 (MQ=255)
                                                                                                                                       aaCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACt             <  8:14301/149‑1 (MQ=255)
                                                                                                                                        aCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGc                             >  4:91699/1‑132 (MQ=255)
                                                                                                                                        aCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGc                             <  3:91699/132‑1 (MQ=255)
                                                                                                                                         cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAt           >  5:323248/1‑149 (MQ=37)
                                                                                                                                         cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAt           >  1:197576/1‑149 (MQ=255)
                                                                                                                                          ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa                                                                           <  7:28453/84‑1 (MQ=17)
                                                                                                                                          ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCaaaaaa                                                                           >  8:28453/1‑84 (MQ=17)
                                                                                                                                          ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc          >  2:237664/1‑149 (MQ=255)
                                                                                                                                            gCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtctt        >  5:309989/1‑149 (MQ=255)
                                                                                                                                             cTCCTACATGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtcttc       <  5:85887/149‑1 (MQ=255)
                                                                                                                                              tCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTAc                                 <  6:337757/122‑1 (MQ=255)
                                                                                                                                              tCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTAc                                 >  5:337757/1‑122 (MQ=255)
                                                                                                                                              tCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATc                          <  6:237122/129‑1 (MQ=255)
                                                                                                                                              tCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATc                          >  5:237122/1‑129 (MQ=255)
                                                                                                                                              tCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCAGCACTAtctact      >  8:294176/1‑149 (MQ=255)
                                                                                                                                                    aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt  >  1:236839/1‑147 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT  >  NC_002947/955003‑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: