Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,911,840	C→G	12.8%	R129P (CGT→CCT)	ccmH ←	cytochrome c heme lyase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,911,840	0	C	G	12.8%	105.2 / 8.9	39	R129P (CGT→CCT)	ccmH	cytochrome c heme lyase subunit
Reads supporting (aligned to +/- strand):  ref base C (19/15);  new base G (2/3);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

CGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCGATGATCTGC  >  NC_002947/4911691‑4911989                                                                                                                                                       |                                                                                                                                                      cGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCAc                                                                                                                                                       <  1:430973/150‑1 (MQ=255)      gcagcaAAAAGCTGAGGGCAGCGAGCTGCAGGAGGCCTTCGCTAAGCCAGAATTCAATCATGGGTTCATAATCGCTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGT‑GCAGGgcgc                                                                                                                                                  <  1:1040638/150‑1 (MQ=255)      gcagcaAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGgcgc                                                                                                                                                  >  1:905067/1‑150 (MQ=255)           aaaaaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGT‑GCAGGGCGCCgatg                                                                                                                                             >  2:1109405/1‑147 (MQ=255)            aaaaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGT‑GCAGGGCGCCgatga                                                                                                                                            <  1:1109405/150‑5 (MQ=255)              aaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACATATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTAc                                                                                                                                          >  1:289659/1‑150 (MQ=255)                 cTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGAt                                                                                                                                       <  1:771936/150‑1 (MQ=255)                                    gAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCAc                                                                                                                                                       >  2:420008/1‑115 (MQ=255)                                    gAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCAc                                                                                                                                                       <  1:420008/115‑1 (MQ=255)                                                     cAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGAt                                                                                                   >  2:890048/1‑150 (MQ=255)                                                          ttCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCAAAGCTCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCg                                                                                              <  1:850972/150‑1 (MQ=255)                                                             aaTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGc                                                                                           >  2:890616/1‑150 (MQ=255)                                                             aaTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGc                                                                                           <  1:696581/150‑1 (MQ=255)                                                                      tCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACcaca                                                                                  >  2:900458/1‑150 (MQ=255)                                                                          ttCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagcag                                                                              >  1:934882/1‑150 (MQ=255)                                                                              ttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAg                                                                          <  1:890616/150‑1 (MQ=255)                                                                              ttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAg                                                                          <  1:1059301/150‑1 (MQ=255)                                                                                    agcagTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGc                                                                    >  2:231989/1‑150 (MQ=255)                                                                                              cGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAgcg                                                          >  1:546204/1‑150 (MQ=255)                                                                                               gagaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGc                                                                                                                             >  1:804340/1‑82 (MQ=255)                                                                                               gagaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGc                                                                                                                             <  2:804340/82‑1 (MQ=255)                                                                                                agaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGTgcggt                                                        <  2:546204/150‑3 (MQ=255)                                                                                                 gaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCg                                                       >  1:962995/1‑150 (MQ=255)                                                                                                        gcgTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTa                                                >  1:663461/1‑150 (MQ=255)                                                                                                                 ggCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGaa                                       >  2:320303/1‑150 (MQ=255)                                                                                                                  gCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAg                                      >  2:178927/1‑150 (MQ=255)                                                                                                                   ctgACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGt                                     <  1:308053/148‑1 (MQ=255)                                                                                                                      aCAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTCGCCCACTaatgcag                                  >  2:308053/1‑145 (MQ=255)                                                                                                                        aaCTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTa                                                >  1:458744/1‑134 (MQ=255)                                                                                                                        aaCTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTa                                                <  2:458744/134‑1 (MQ=255)                                                                                                                          cTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCg                                                                                              >  2:753877/1‑86 (MQ=255)                                                                                                                          cTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCg                                                                                              <  1:753877/86‑1 (MQ=255)                                                                                                                           tCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAg                             >  1:338553/1‑150 (MQ=255)                                                                                                                               cGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGt                         <  1:681888/150‑1 (MQ=255)                                                                                                                                ggCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTc                        >  1:605393/1‑150 (MQ=255)                                                                                                                                 gcgcCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCg                       <  2:962995/150‑1 (MQ=255)                                                                                                                                    ccgtgctaCGGCCGT‑GCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCg                       >  2:597737/8‑147 (MQ=255)                                                                                                                                    ccgtgctaCGGCCGT‑GCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCg                       <  1:597737/140‑1 (MQ=255)                                                                                                                                     ccTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGAACACCGCCGGCACCACGAAGCCACCGGCCAGCAAGATCCGCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGAccc                   >  1:933388/1‑149 (MQ=255)                                                                                                                                                cGTAGCCCGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCgat         <  1:490197/150‑1 (MQ=255)                                                                                                                                                cGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGTCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCgat         <  2:605393/150‑1 (MQ=255)                                                                                                                                                    gCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCGATGATc     <  1:231989/150‑1 (MQ=255)                                                                                                                                                       cGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCGATGATctgc  <  1:890048/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CGGGAGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGC‑ACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGACCATGAAGTCGATGATCTGC  >  NC_002947/4911691‑4911989

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: ‑