Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	955,085:1	+G	14.9%	intergenic (+36/‑277)	cyoE → / → alaC	protoheme IX farnesyltransferase/aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	955,085	1	.	G	14.9%	160.1 / 25.5	47	intergenic (+36/‑277)	cyoE/alaC	protoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base . (18/22);  new base G (3/4);  total (21/26)
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.94e-01

CCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAG  >  NC_002947/954944‑955227                                                                                                                                               |                                                                                                                                               ccGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGt                                                                                                                                        >  1:1093190/1‑150 (MQ=255) ccGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGCATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGt                                                                                                                                        >  2:115527/1‑150 (MQ=255)       ggTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGGttttcgt                                                                                                                                  <  2:910128/150‑4 (MQ=255)                 ttCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGATAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCATGAAAATTCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTAt                                                                                                                        <  2:85752/150‑1 (MQ=255)                     ccATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTAttttt                                                                                                                    >  2:1130164/1‑150 (MQ=255)                                        ccTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGc                                                                                                         >  1:620829/1‑142 (MQ=255)                                        ccTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGc                                                                                                         <  2:620829/142‑1 (MQ=255)                                         cTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGGCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGGttt                                                                                                                                      >  1:288311/1‑112 (MQ=255)                                         cTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGGCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGGttt                                                                                                                                      <  2:288311/112‑1 (MQ=255)                                          tGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCgg                                                                                               <  1:1287037/150‑1 (MQ=255)                                             gCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTa                                                                                            >  2:712099/1‑150 (MQ=255)                                                tgatgaTGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAAc                                                                                         <  1:712099/150‑1 (MQ=255)                                                    gatgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCtgt                                                                                                               >  1:127660/1‑124 (MQ=255)                                                    gatgGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCtgt                                                                                                               <  2:127660/124‑1 (MQ=255)                                                         gtgtGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTa                                                                                <  2:1174669/150‑1 (MQ=255)                                                             ggACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAg                                                                            <  2:1313260/150‑1 (MQ=255)                                                             ggACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAg                                                                            <  1:390092/150‑1 (MQ=255)                                                                   ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                      >  1:1242685/1‑150 (MQ=255)                                                                   ccAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                      <  2:1093190/150‑1 (MQ=255)                                                                    cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTg                                                                     <  2:1242685/150‑1 (MQ=255)                                                                         cGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTGCTGTAggg                                                                >  1:1389913/1‑149 (MQ=255)                                                                             gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTTTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                            <  1:886472/150‑1 (MQ=255)                                                                                  tgtgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTa                                                       <  2:1323575/150‑1 (MQ=255)                                                                                          tgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGaa                                               <  1:1130164/150‑1 (MQ=255)                                                                                            aCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGGttttcgt                                                                                                                                  <  2:1129306/65‑4 (MQ=18)                                                                                            aCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGG‑CACAGTGGCCGGGGttttcgt                                                                                                                                  >  1:1129306/1‑62 (MQ=18)                                                                                              cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                              <  2:107290/131‑1 (MQ=255)                                                                                              cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGc                                                              >  1:107290/1‑131 (MQ=255)                                                                                                        tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCgg                                                                                               <  2:412245/88‑1 (MQ=255)                                                                                                        tGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCgg                                                                                               >  1:412245/1‑88 (MQ=255)                                                                                                          aaaCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                        >  2:1317833/1‑109 (MQ=255)                                                                                                          aaaCTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                        <  1:1317833/109‑1 (MQ=255)                                                                                                               cttgcCTTTTTCACGAAACCCCCGGCCCCTGT‑GCGGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaataa                          <  2:1389913/148‑1 (MQ=255)                                                                                                                     tGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTgg                                                                                                                              <  1:1066443/44‑1 (MQ=255)                                                                                                                     tGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTgg                                                                                                                              >  2:1066443/1‑44 (MQ=255)                                                                                                                          cACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGATCGGGTTTACCTGCTAAGCGGCCAGCACAGAAAAATAAAAATACCTGAt               >  2:241742/1‑150 (MQ=255)                                                                                                                           aCGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGc                                                                                     <  2:1041035/79‑1 (MQ=255)                                                                                                                           aCGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGc                                                                                     >  1:1041035/1‑79 (MQ=255)                                                                                                                               aaaCCCCGG‑CACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCgg                                                                                               >  2:247594/1‑65 (MQ=25)                                                                                                                               aaaCCCCGG‑CACAGTGGCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCgg                                                                                               <  1:247594/65‑1 (MQ=25)                                                                                                                                  ccccGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCgg                                                                                               <  2:107708/62‑1 (MQ=255)                                                                                                                                  ccccGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCgg                                                                                               >  1:107708/1‑62 (MQ=255)                                                                                                                                   cccGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAc                                                      >  2:1143651/1‑102 (MQ=255)                                                                                                                                   cccGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAc                                                      <  1:1143651/102‑1 (MQ=255)                                                                                                                                   cccGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAgg                                           <  1:949616/113‑1 (MQ=255)                                                                                                                                   cccGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAgg                                           >  2:949616/1‑113 (MQ=255)                                                                                                                                       gCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTTTCCTGTAGGAGCGGGTTTTCCCGCGAAAAGGCCAGCACCGAAAAATAAAAATACCTGATTTTCAAAAAATaa  >  2:1358630/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                               CCGCAAGGTGTTCGGCTTCTCCATCCTCACCGTCACTGCCCTGAGCGTGATGATGGGTGTGGACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGT‑GCCGGGGTTTTTTTTTGGGGGTATTTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAG  >  NC_002947/954944‑955227

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