Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 115,909 G→A 7.5% R71Q (CGG→CAG)  PP_0110 → CyoE‑like protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947115,9090GA 7.5% 149.4 / 3.1 53R71Q (CGG→CAG) PP_0110CyoE‑like protoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base G (20/29);  new base A (2/2);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

CAAGGTGGTGGTGCTGATGCTGATCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGC  >  NC_002947/115760‑116058
                                                                                                                                                     |                                                                                                                                                     
cAAGGTGGTGGTGCTGATGCTGATCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCg                                                                                                                                                       >  1:89032/1‑150 (MQ=255)
  aGGTGGTGTTGCTGATGCTGATCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGgc                                                                                                                                                     <  1:221519/150‑1 (MQ=255)
  aGGTGGTGGTGCTGATGCTGATCACCTCGCTGGCGGGCATATTCCTGGCCAACTGTGCGGGCGTCAGCTGGAGCGTACTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGgc                                                                                                                                                     <  1:731232/150‑1 (MQ=255)
         ggtgCTGATGCTGATCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGa                                                                                                                                              <  1:726608/150‑1 (MQ=255)
         ggtgCTGATGCTGATCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGa                                                                                                                                              >  1:588380/1‑150 (MQ=255)
          gtgCTGATGCTGATCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGAt                                                                                                                                             <  1:713149/150‑1 (MQ=255)
          gtgCTGATGCTGACCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGAt                                                                                                                                             >  2:216444/1‑150 (MQ=255)
           tgctgatgctgaTCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATg                                                                                                                                            >  2:409525/1‑150 (MQ=255)
                         aCCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGgcacgc                                                                                                                              >  2:654166/1‑150 (MQ=255)
                           cTCGCTGGCGGGCATGTTCCTGGCCACCC‑GGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCAcc                                                                                                                           <  1:37116/150‑1 (MQ=255)
                                     ggCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCAGCGCATCAATgcgccga                                                                                                                                      >  1:749640/1‑127 (MQ=255)
                                     ggCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCAGCGCATCAATGCGCCGa                                                                                                                                      <  2:749640/130‑1 (MQ=255)
                                      gCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGc                                                                                                                 <  2:89032/150‑1 (MQ=255)
                                            tCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTgg                                                                                                           <  1:471627/150‑1 (MQ=255)
                                            tCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTgg                                                                                                           <  1:290219/150‑1 (MQ=255)
                                             ccTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCCCCCACAAACGGCCTTTGGc                                                                                                          <  1:35580/150‑1 (MQ=255)
                                              cTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGcc                                                                                                         <  1:654166/150‑1 (MQ=255)
                                               gggCCACCGGGGCGGGCGTCAGCGGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGGGGGTGGTGCGGCAGTGGTAAATCATGCGGTGGATCGGCGCATTGAGGCGCTGATGGCACGCACCCACAAACGGCTTTTGGccc                                                                                                        <  2:485732/149‑1 (MQ=255)
                                                 gCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCaa                                                                                                      >  2:737078/1‑150 (MQ=255)
                                                  ccACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAg                                                                                                     <  1:120354/150‑1 (MQ=255)
                                                     cccGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGcc                                                                                                  <  2:423825/150‑1 (MQ=255)
                                                     cccGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGcc                                                                                                  >  2:359028/1‑150 (MQ=255)
                                                            gggcgTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCAGCGCATCAATgcgccga                                                                                                                                      >  2:752591/1‑104 (MQ=255)
                                                            gggcgTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCAGCGCATCAATGCGCCGa                                                                                                                                      <  1:752591/107‑1 (MQ=255)
                                                                   aGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTg                                                                                    >  2:241285/1‑150 (MQ=255)
                                                                                 gttgttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGgcac                                                                                                                                >  1:496570/1‑92 (MQ=255)
                                                                                 gttgttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGgcac                                                                                                                                <  2:496570/92‑1 (MQ=255)
                                                                                    gttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCtgtt                                                                         <  1:438635/144‑1 (MQ=255)
                                                                                    gttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCtgtt                                                                         >  2:438635/1‑144 (MQ=255)
                                                                                     ttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGcc                                                                  >  1:513107/1‑150 (MQ=255)
                                                                                       cGGCAGCCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCtt                                                                <  2:604444/150‑1 (MQ=255)
                                                                                       cGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCtt                                                                <  2:320562/150‑1 (MQ=255)
                                                                                       cGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCCGTTGTTCGCCtt                                                                >  1:149748/1‑150 (MQ=255)
                                                                                              cTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTg                                                         >  1:316069/1‑150 (MQ=255)
                                                                                               tGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTgg                                                        <  1:67044/150‑1 (MQ=255)
                                                                                                gggCATCGGCATGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGgc                                                       >  2:184233/1‑150 (MQ=255)
                                                                                                 ggCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGCTGCGCTGATGGCACGCACCCACAAACGGCCTGTTGCCCGAGTCCTAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGgcg                                                      >  2:273472/1‑150 (MQ=255)
                                                                                                                gCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGTCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTc                                       <  1:359028/150‑1 (MQ=255)
                                                                                                                gCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTc                                       <  2:513107/150‑1 (MQ=255)
                                                                                                                  gggTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCct                                     >  2:530469/1‑150 (MQ=255)
                                                                                                                   ggtggtGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATgcgc                                                                                                                                         >  2:529556/1‑49 (MQ=255)
                                                                                                                   ggtggtGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATgcgc                                                                                                                                         <  1:529556/49‑1 (MQ=255)
                                                                                                                     tggtgCGGCAGTGGTCAATCATGTGGTGGATCGGCGGATTGATGCGCTGATGGCACGCACCCACAATCGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCttgct                                  <  2:693712/150‑1 (MQ=255)
                                                                                                                                tGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCa                       <  1:530469/150‑1 (MQ=255)
                                                                                                                                 ggTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCAc                      <  1:723282/150‑1 (MQ=255)
                                                                                                                                  gTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCgctgc                                                  >  2:408067/1‑121 (MQ=255)
                                                                                                                                  gTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCgctgc                                                  <  1:408067/121‑1 (MQ=255)
                                                                                                                                   tCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgctg                              >  1:340445/1‑140 (MQ=255)
                                                                                                                                   tCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgctg                              <  2:340445/140‑1 (MQ=255)
                                                                                                                                   tCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCa                    >  1:601349/1‑150 (MQ=255)
                                                                                                                                             ggtggATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCAc          <  1:234744/150‑1 (MQ=255)
                                                                                                                                                ggATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgc                                   <  2:407100/122‑1 (MQ=255)
                                                                                                                                                ggATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgc                                   >  1:407100/1‑122 (MQ=255)
                                                                                                                                                   tCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTg    <  2:316069/150‑1 (MQ=255)
                                                                                                                                                     ggCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGc  <  1:737078/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAGGTGGTGGTGCTGATGCTGATCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCTGGC  >  NC_002947/115760‑116058

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