Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,732,931	A→T	19.0%	M245K (ATG→AAG)	sucD ←	succinyl‑CoA synthetase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,732,931	0	A	T	19.0%	42.1 / 5.9	21	M245K (ATG→AAG)	sucD	succinyl‑CoA synthetase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (7/10);  new base T (2/2);  total (9/12)
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.14e-01

AGCTCGGCCAGGGCCTTGCCGATGTCAGCCAGGGAACGCACGGTTTTCACACCCGCGTCCTGCAGGGCGGCGAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGACGTGGGCCTTGATGTAGGCCGCAGCTTCTTCTTCAGCCGAACCGCCGATCTCACCGATCATGACGATCGCT  >  NC_002947/4732805‑4733059                                                                                                                               |                                                                                                                                 aGCTCGGCCAGGGCCTTGCCGATGTCAGCCAGGGAACGCACGGTTTTCACACCCGCGTCCTGCAGGGCGGCGAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGc                                                                                                                         <  2:129016/136‑1 (MQ=255)      ggCCAGGGCCTTGCCGATGTCAGCCAGGGAACGCACGGTTTTCACACCCGCGTCCTGCAGGGCGGCGAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCa                                                                                                                    >  2:9867/1‑136 (MQ=255)                   ccGATGTCAGCCAGGGAACGCACGGTTTTCACACCCGCGTCCTGCAGGGCGGCGAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAacac                                                                                                       <  1:94833/136‑1 (MQ=255)                     gATGTCAGCCAGGGAACGCACGGTTTTCACACCCGCGTCCTGCAGGGCGGCGAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCg                                                                                                     >  1:88055/1‑136 (MQ=255)                     gATGTCAGCCAGGGAACGCACGGTTTTCACACCCGCGTCCTGCAGGGCGGCGAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCg                                                                                                     >  2:399486/1‑136 (MQ=255)                          cAGCCAGGGAACGCACGGTTTTCACACCCGCGTCCTGCAGGGCGGCGAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCTTGCGCTTGCCCGCAGGTGCAGTAACACCGGCGAt                                                                                                >  2:295475/1‑136 (MQ=255)                          cAGCCAGGGAACGCACGGTTTTCACACCCGCGTCCTGCAGGGCGGCGAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCTTGCGCTTGCCCGCAGGTGCAGTAACACCGGCGAt                                                                                                >  2:438362/1‑136 (MQ=255)                          cAGCCAGGGAACGCACGGTTTTCACACCCGCGTCCTGCAGGGCGGCGAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCGGCGAt                                                                                                >  2:288419/1‑136 (MQ=255)                                                 acacCCGCGTCCTGCAGGGCGGCGAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGa                                                                         <  2:48044/136‑1 (MQ=255)                                                  cacCCGCGTCCTGCAGGGCGGCGAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGAc                                                                        <  1:399486/136‑1 (MQ=255)                                                           ccTGCAGGGCGGCGAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGACGTGGGCCtt                                                               >  1:4947/1‑136 (MQ=255)                                                           ccTGCAGGGCGGCGAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGACGTGGGCCtt                                                               >  1:93175/1‑136 (MQ=255)                                                                        gAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGACGTGGGCCTTGATGTAGGCCGCa                                                  <  1:430712/136‑1 (MQ=255)                                                                        gAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGACGTAGGCCTTGATGTAGGCCGCa                                                  <  1:9867/136‑1 (MQ=255)                                                                                     gCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGACGTGGGCCTTGATGTAGGCCGCAGCTTCTTCTTCAg                                     <  2:498757/136‑1 (MQ=255)                                                                                       aGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGACGTGGGCCTTGATGTAgg                                                       <  1:384259/116‑1 (MQ=255)                                                                                       aGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGACGTGGGCCTTGATGTAgg                                                       >  2:384259/1‑116 (MQ=255)                                                                                                    ccgGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGACGTGGGCCTTGATGTAGGCCGCAGCTTCTTCTTCAGCCGAACCGCCGAtct                      <  1:348087/136‑1 (MQ=255)                                                                                                                        gTGGCCCTTGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGACGTGGGCCTTGATGTAGGCCGCAGCTTCTTCTTCAGCCGAACCGCCGATCTCACCGATCATGACGATCGCt  <  1:438362/136‑1 (MQ=255)                                                                                                                        gTGGCCCTTGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGACGTGGGCCTTGATGTAGGCCGCAGCTTCTTCTTCAGCCGAACCGCCGATCTCACCGATCATGACGATCGCt  <  1:295475/136‑1 (MQ=255)                                                                                                                        gTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGACGTGGGCCTTGATGTAGGCCGCAGCTTCTTCTTCAGCCGAACCGCCGATCTCACCGATCATGACGATCGCt  <  1:288419/136‑1 (MQ=255)                                                                                                                               |                                                                                                                                 AGCTCGGCCAGGGCCTTGCCGATGTCAGCCAGGGAACGCACGGTTTTCACACCCGCGTCCTGCAGGGCGGCGAACTTCTCGTCCGCAGTGCCCTTGCCGCCGGAGATGATGGCGCCAGCGTGGCCCATGCGCTTGCCCGCAGGTGCAGTAACACCGGCGATGTAGGAAACGACAGGCTTGGTGACGTGGGCCTTGATGTAGGCCGCAGCTTCTTCTTCAGCCGAACCGCCGATCTCACCGATCATGACGATCGCT  >  NC_002947/4732805‑4733059

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