Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 115,909 G→A 12.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,9090GA12.5% 76.8 / 2.8 32R71Q (CGG→CAG) PP_0110CyoE‑like protoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base G (13/15);  new base A (2/2);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AAGGTGGTGGTGCTGATGCTGATCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCA  >  NC_002947/115761‑116040
                                                                                                                                                    |                                                                                                                                   
aaGGTGGTGGTGCTGATGCTGATCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCgg                                                                                                                                    >  2:555630/1‑150 (MQ=255)
           gctgatgctgatCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATgc                                                                                                                         <  1:210830/150‑1 (MQ=255)
             tgatgctgatCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATgcgc                                                                                                                       <  1:555630/150‑1 (MQ=255)
                          cTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGANCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGAGCGGCGCATTGATGCGCTGATTgcacgcac                                                                                                          >  1:238586/1‑150 (MQ=255)
                             gCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCca                                                                                                       <  2:116397/150‑1 (MQ=255)
                                  cGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAAc                                                                                                  <  1:443695/150‑1 (MQ=255)
                                   gggCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCAGCGCATCAATgcgccga                                                                                                                    >  2:578688/1‑128 (MQ=255)
                                   gggCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCAGCGCATCAATGCGCCGa                                                                                                                    <  1:578688/131‑1 (MQ=255)
                                     gCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGc                                                                                               >  1:56351/1‑150 (MQ=255)
                                     gCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGc                                                                                               >  2:11579/1‑150 (MQ=255)
                                               ggCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCa                                                                                     >  1:206109/1‑150 (MQ=255)
                                                 ccACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAg                                                                                   <  2:206109/150‑1 (MQ=255)
                                                    cccGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTGGGCCAAAGGcc                                                                                >  2:309330/1‑150 (MQ=255)
                                                     ccGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCg                                                                               >  2:367696/1‑150 (MQ=255)
                                                      cgggcgggCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGa                                                                              <  2:525468/150‑1 (MQ=255)
                                                          cgggcgTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTGGGCCAAAGGCCgtttgt                                                                          <  1:309330/150‑6 (MQ=255)
                                                             gcgTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGc                                                                       <  1:216907/150‑1 (MQ=255)
                                                                             gCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGGATTGATGCGCGGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCtgtt                                                       <  1:23576/150‑1 (MQ=255)
                                                                             gCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTg                                                                  <  1:276287/139‑1 (MQ=255)
                                                                             gCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTg                                                                  >  2:276287/1‑139 (MQ=255)
                                                                                       ggCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACACACCCACAAACGGCCTTTGGCCTAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTg                                             >  1:349713/1‑150 (MQ=255)
                                                                                          aaCCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCtgttg                                                      <  1:380204/138‑1 (MQ=255)
                                                                                          aaCCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCtgttg                                                      >  2:380204/1‑138 (MQ=255)
                                                                                             cTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTg                                       >  1:385072/1‑150 (MQ=255)
                                                                                             cTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTg                                       >  1:410191/1‑150 (MQ=255)
                                                                                              tGGGCATCGGCTTGTGTGCGGGTGGTGCGGCAGTGGTTAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTgg                                      <  2:144435/150‑1 (MQ=255)
                                                                                              tGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTgg                                      <  2:385072/150‑1 (MQ=255)
                                                                                                       gCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTgg                             >  1:535657/1‑150 (MQ=255)
                                                                                                                 gggTGGGGCGGCAGTGGTCACTCATGTGGTGGCTCGGCGCCTTGATGCGCTGATGGCACGCTCCAACACACGGCCTTTGCCCCAAGGCTGAGTCGAGCCCCTGCCGGCGCTCTTGTTCGCGTTGGCCCTGGCGCTGCTGGGCCTGGTCct                   <  2:238586/150‑1 (MQ=255)
                                                                                                                           gCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCACGGCCCAGTCGAGAACCCGTCGTCGGTTCGGGTCTTCTGGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGgtg         <  1:11579/150‑1 (MQ=255)
                                                                                                                                  tCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTAGCGTGCCGAAACCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCa  >  2:16132/1‑150 (MQ=255)
                                                                                                                                                 aTCAGCGCATCAATGCGCCGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgctg            <  2:364287/125‑1 (MQ=255)
                                                                                                                                                 aTCAGCGCATCAATGCGCCGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgctg            >  1:364287/1‑125 (MQ=255)
                                                                                                                                                    |                                                                                                                                   
AAGGTGGTGGTGCTGATGCTGATCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCA  >  NC_002947/115761‑116040

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