Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 955,076 Δ1 bp 10.0% intergenic (+30/‑283) cyoE → / → alaC protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,0760C.10.0% 151.6 / 11.2 50intergenic (+30/‑283)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (20/25);  new base . (2/3);  total (22/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 = 6.98e-01

AGTGGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATA  >  NC_002947/954934‑955223
                                                                                                                                               |                                                                                                                                                    
aGTGGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGCTGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACtgt                                                                                                                                                >  3:243308/1‑149 (MQ=255)
 gTGGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACtgt‑g                                                                                                                                              >  4:266869/1‑149 (MQ=255)
       ccGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCgggg                                                                                                                                        <  8:177346/149‑1 (MQ=255)
        cGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGGt                                                                                                                                       <  1:231420/149‑1 (MQ=255)
         gCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATTGGGCGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACATACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGt                                                                                                                                       >  1:101736/1‑149 (MQ=255)
              gtgtTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGGttttcgt                                                                                                                                 <  7:292922/149‑4 (MQ=255)
                    ggCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTgggg                                                                                                                           >  1:177903/1‑149 (MQ=255)
                     gCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGtt                                                                                                                                      <  1:226061/137‑1 (MQ=255)
                     gCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGtt                                                                                                                                      >  2:226061/1‑137 (MQ=255)
                             aTCCTCACCGTCACTGCCCTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACAAAAACCCCGG‑CACAGTGGCCGGGGttttcgt                                                                                                                                 <  8:267983/134‑4 (MQ=39)
                             aTCCTCACCGTCACTGCCCTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACAAAAACCCCGG‑CACAGTGGCCGGGGttttcgt                                                                                                                                 >  7:267983/1‑131 (MQ=39)
                             aTCCTCACCGTCACTGCCCTGAGCGTGATGA‑TGGGCGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTc                                                                                                                  <  2:101736/149‑1 (MQ=255)
                                cTCACCGTCACTGCCCTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCtg                                                                                                               >  2:355330/1‑149 (MQ=255)
                                              ccTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCttag                                                                                                 >  8:1354/1‑147 (MQ=255)
                                               cTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCTGGCCTCTTcgc                                                                                                <  3:202206/149‑1 (MQ=255)
                                               cTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTcgc                                                                                                <  6:364447/149‑1 (MQ=255)
                                                  aGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCgg                                                                                              <  8:190513/148‑1 (MQ=255)
                                                  aGCGAGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCgg                                                                                              <  7:185942/148‑1 (MQ=255)
                                                   gCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGt                                                                                            <  2:177903/149‑1 (MQ=255)
                                                         tga‑tgGGTGGGTACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCAGCCTGTTTCACGAAAACCCCGGCCACTGG‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAccc                                                                                      <  4:371622/149‑1 (MQ=255)
                                                          ga‑tgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCg                                                                                     >  2:267671/1‑149 (MQ=255)
                                                           a‑tgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGc                                                                                    <  1:267671/149‑1 (MQ=255)
                                                           a‑tgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGc                                                                                    <  3:266869/149‑1 (MQ=255)
                                                                gtgtGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTcgc                                                                                                <  5:409299/133‑1 (MQ=255)
                                                                gtgtGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTcgc                                                                                                >  6:409299/1‑133 (MQ=255)
                                                                    ggACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCTGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTcc                                                                                                                 >  7:246858/1‑112 (MQ=255)
                                                                        aGCCAGACCGCTGCGGATGTGGTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGt                                                                                            <  4:157037/129‑1 (MQ=255)
                                                                           cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                     <  2:269519/149‑1 (MQ=255)
                                                                           cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                     >  1:269519/1‑149 (MQ=255)
                                                                               ccGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTc                                                                                                  >  5:51671/1‑116 (MQ=255)
                                                                               ccGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTc                                                                                                  <  6:51671/116‑1 (MQ=255)
                                                                                cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCgg                                                                                              >  5:323658/1‑119 (MQ=255)
                                                                                cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCgg                                                                                              <  6:323658/119‑1 (MQ=255)
                                                                                    gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTT‑CTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                           <  3:65512/149‑1 (MQ=25)
                                                                                     cgtaTGTGCTGATGACCAACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                           <  8:177015/146‑1 (MQ=255)
                                                                                     cGGATGTGCTGATGACCTACGCAAGCTGAAACTGCTTCCTGTTTAACGAAAACCACGGCCACTGT‑GCCGGGGTTTTTTTTTGGTGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                           <  6:193968/149‑1 (MQ=25)
                                                                                       gATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGt                                                         <  1:355330/149‑1 (MQ=255)
                                                                                            gCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAc                                                                                        <  5:378986/113‑1 (MQ=255)
                                                                                            gCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAc                                                                                        >  6:378986/1‑113 (MQ=255)
                                                                                                atgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCggg                                                                                             <  4:6893/104‑1 (MQ=255)
                                                                                                atgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCggg                                                                                             >  3:6893/1‑104 (MQ=255)
                                                                                                           acgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa                                                                  <  1:29133/120‑1 (MQ=255)
                                                                                                           acgcTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa                                                                  >  2:29133/1‑120 (MQ=255)
                                                                                                                             gTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTa                                                                               <  6:45833/89‑1 (MQ=255)
                                                                                                                             gTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTa                                                                               >  5:45833/1‑89 (MQ=255)
                                                                                                                              tttCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTATTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTGCTACAAGGTTGCGGTAGGAGCGGGGTTACCCGCGGAGTGGCCAGCACAGAAAAATAAAAATTcc                  >  1:59425/1‑149 (MQ=37)
                                                                                                                                    gAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                     >  2:246120/1‑92 (MQ=255)
                                                                                                                                    gAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                     <  1:246120/92‑1 (MQ=255)
                                                                                                                                       aaCCCCGG‑CACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCACTTAGCGGGTAAACCCGCTACTACAAGGTTGCGGTAGGAGCGGGTTCACCCGGGAAGAGGGCACCACAGAAAAATAAAAATAACTGAATTTc          >  8:152452/1‑148 (MQ=25)
                                                                                                                                              gCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGGTACTACAAGGTTCCTGTAGGAGCCGGTTTAACAGCGAAAAGGCCAGCACAGAAAAATAAGAATACATGATTTTCAAAAAATa  >  5:371418/1‑149 (MQ=37)
                                                                                                                                               |                                                                                                                                                    
AGTGGGCCCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGA‑TGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATA  >  NC_002947/954934‑955223

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