Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,149 A→G 13.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,1490AG13.2% 122.3 / 9.5 55intergenic (+103/‑210)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (20/26);  new base G (2/5);  total (22/33)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

AGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGA  >  NC_002947/955005‑955292
                                                                                                                                                     |                                                                                                                                                
aGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGTTCCTACAAGGt                                                                                                                                               >  7:109271/1‑148 (MQ=255)
 gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                             <  2:359155/149‑1 (MQ=255)
 gCCAGACCGCTGCGGATGTGCGGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                             <  8:247704/149‑1 (MQ=255)
  ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTcc                                                                                                                                            >  7:178431/1‑149 (MQ=255)
  ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTTCAAGGTTcc                                                                                                                                            >  4:208409/1‑149 (MQ=255)
   ctatttgCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                           <  1:309103/143‑1 (MQ=255)
   cAGACCGCTGCGGCTGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGTAAACCCCGGCCACTGGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                           <  6:361979/149‑1 (MQ=255)
       ccGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAg                                                                                                                                       <  8:51803/149‑1 (MQ=255)
        cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                           >  2:309103/1‑144 (MQ=255)
             cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGTTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                 <  8:109271/149‑1 (MQ=255)
             cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                 <  5:156814/149‑1 (MQ=38)
             cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                 <  1:337936/149‑1 (MQ=38)
             cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCg                                                                                                                                  <  1:9285/148‑1 (MQ=38)
             cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                                                                                                   <  7:161293/147‑1 (MQ=255)
              gtaTGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACTAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                <  1:233788/147‑1 (MQ=37)
              ggATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCggg                                                                                                                                <  4:61479/149‑1 (MQ=37)
               gATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGt                                                                                                                               <  7:139948/149‑1 (MQ=255)
                     cTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAccc                                                                                                                         <  4:14946/149‑1 (MQ=255)
                      tgatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCAACAGGAACCATGTAGGATCGGGTTTACCcg                                                                                                                        >  2:41298/1‑149 (MQ=37)
                       gatgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGACAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAGTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCcgc                                                                                                                       <  5:263995/149‑1 (MQ=255)
                        atgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCcgc                                                                                                                       <  1:178304/148‑1 (MQ=255)
                            cttACCCACGCTGAAAGTGCTGCCTTTTTCACTAAAACCCGGGCCAGTTTGGGGGGTTTTTTTTTTGGGGTTATTTTTCGTGTGCCGGCCTCTTGGGGGTAAAACCCGCTCCTACAAGGTTCCTTTTGTTGGGGTTTTACCCGCGAaga                                                                                                                  <  6:43139/147‑1 (MQ=255)
                                     gcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAAACGCTACTAAAATGTTAATTTAGGAGCGGGTTTACCCGCGAAGAGGCCAGc                                                                                                           >  5:72804/1‑148 (MQ=37)
                                       tGAAACTGCTGCCTGTTTCACGAAAACCCCGCCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCAATTCGCGGGTTAACCCCCTCCTACAAGGTTACAGTAGGAGCGCGTTTAACCGCGAAACGGCCAGca                                                                                                          >  1:38584/1‑146 (MQ=255)
                                           aCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAGGTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGAACCTACAAGGTTACTTTATGAGCTGCTTTACAAGCCAAAACGCCAGCAAAGaa                                                                                                     >  1:103185/1‑148 (MQ=25)
                                            ctgctgCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCGGCTTCAAGGTTTCTTTTGGAGCGGGGTTTCCCGGGAAGAGGGCCGCACAGaaaaa                                                                                                  >  3:243710/1‑149 (MQ=38)
                                               ctgcCTGTTTCACCAAAACCCCGGCCACTGGGCGGGGTTTTTTTTTTTGGGGTTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGTAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaa                                                                                                  <  6:72804/147‑1 (MQ=25)
                                                   cTGTTTCACGAAACCCCGGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCGGGCGTCTTGGCGTGTAAACCCGCTCCAACATGGTTCGTTTTGGTGCGGGTTTACCCGCGAATAGGCCAGCACAGaaaaataaaaat                                                                                           <  2:38584/149‑1 (MQ=255)
                                                       aaaaaaaaaaCCCCGGCACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTTAACCCGCTGCTACAAGGTTGCGGTTGGGGCGGGGTTAACCGGGAAGAGGCCCGCCCAGAAAAATAAAAATACCt                                                                                       >  1:81758/7‑149 (MQ=25)
                                                        tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTTAACCCCCGGCGACCAGGGTGGGGTTGGAGGGGGGTTTCCCGGGAAGCGGGCCGGCCAGAAAAATAAAAATTCCt                                                                                       >  7:62319/1‑149 (MQ=255)
                                                        tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCt                                                                                       <  8:178431/148‑1 (MQ=255)
                                                        tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                 <  4:235071/106‑1 (MQ=255)
                                                        tCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                 >  3:235071/1‑106 (MQ=255)
                                                                  cccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCACCGACAAGGTTGCTGTATGAGCGGGTTTACCCGCGAAGAGGCCCGCACAGAAAAATTAAAAAAACTGATTTTGaaa                                                                             >  8:64559/1‑148 (MQ=37)
                                                                     ggCCCCTGGGCGGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTGCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAAt                                                                         <  6:300256/149‑1 (MQ=255)
                                                                      gCCACTGTGCCGGGGTTTTTTTTTGGGGTAATTTTTCCTGTGCCGGCCTCTTCGCGGTAAAACCCGCTCCTACAAGTTTCCTGTAGTAGCGTGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATa                                                                        <  7:303171/149‑1 (MQ=38)
                                                                          cTGTGCCGGGTTTTTTTTTTTGGGGTATTTTTCCTTTGCCGGCCTCTTCTCGGTTAAACCCTCTCCTACAATGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATc                                                                    <  1:34599/149‑1 (MQ=37)
                                                                             cctctctggTTTTTTTTTGGGGGTATTTTTCCTTTGCCGGCCTCTTCTCGGGTAAACCCGCTCATACAAGTATCCTGTAGTAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGAt                                                                     <  3:30829/139‑1 (MQ=25)
                                                                              ggCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGCCCTCTTCGCGGTAAAACCCGCTCCAACAAGTTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCt                                                                   <  3:208409/146‑1 (MQ=38)
                                                                                  ggggTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCACTTCGCGGGTAAACACGAGAACTACAAGGTTGGTGTATGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaata                                                                                                >  1:99251/1‑115 (MQ=25)
                                                                                         ttttttttGGGGGTATTTTTCCTGTCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCt                                                        <  8:62319/148‑1 (MQ=255)
                                                                                                             tcttgggCCGGCCTCTTCGCTTTTAAACCCTCTCATAAAATTTTCTTTTATTAGCGTTTTTACCCGCTAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATc                                     <  5:226772/142‑1 (MQ=255)
                                                                                                                 gtgCCGGCCTCTTCGCTGTTAAACCCGCTCCTACATTTTTCCTGTAGTTGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATc                                     <  1:117341/144‑1 (MQ=255)
                                                                                                                         ctctTCGCGGGAAAACCCGCTCCTACAAGGTTCCTTTTGGAGGGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATc                        <  3:57329/149‑1 (MQ=39)
                                                                                                                            ttCGCGGTTAAAACCGCTCCAAAAAGGTTCCTTTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAACAGTTCTAATAAAACAGGAAATTATcc                                    <  1:189013/135‑1 (MQ=38)
                                                                                                                            ttCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTTCAGGCATGCTg                     >  7:105174/1‑149 (MQ=255)
                                                                                                                            ttCGCGGGTAAAACCCCGCCGACAAGGTTGCTGTTGGAGCGGGGTTTCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATAACGGATTTTAAAAAAAAAGATCTGAAAAGTTGTAATAAAAAAGGAAATTAtgc                                    >  2:189013/1‑132 (MQ=255)
                                                                                                                                     aaaCCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCAc            >  1:186656/1‑149 (MQ=255)
                                                                                                                                        cccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTa          >  7:199091/1‑148 (MQ=255)
                                                                                                                                         ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                               >  7:96856/1‑110 (MQ=255)
                                                                                                                                         ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                               <  8:96856/110‑1 (MQ=255)
                                                                                                                                           gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct       >  2:100438/1‑148 (MQ=38)
                                                                                                                                                cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGa  >  3:300419/1‑149 (MQ=255)
                                                                                                                                                   cAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCtt                                  <  2:40785/114‑1 (MQ=255)
                                                                                                                                                   cAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCtt                                  >  1:40785/1‑114 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGA  >  NC_002947/955005‑955292

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: