Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,461,176	A→G	D95G (GAC→GGC)	paaJ →	3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,461,176	0	A	G	100.0%	60.9 / NA	18	D95G (GAC→GGC)	paaJ	3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/8);  total (10/8)

GTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAGTGCATTTTCTCGT  >  NC_000913/1461046‑1461300                                                                                                                                   |                                                                                                                             gtgATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTgg                                                                                                                      >  2:330028/1‑139 (MQ=255)         cGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGcc                                                                                                              >  2:256672/1‑139 (MQ=255)               tgCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACgg                                                                                                        >  2:246629/1‑139 (MQ=255)                            ggAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCgg                                                                                           >  2:98692/1‑139 (MQ=255)                                             cGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACCACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTg                                                                          >  1:115734/1‑139 (MQ=255)                                                aGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGCGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGGTTTGCTGATc                                                                       >  1:12040/1‑139 (MQ=255)                                                                    tGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAAt                                                   <  1:98692/139‑1 (MQ=255)                                                                         gggCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCAc                                              >  2:34549/1‑139 (MQ=255)                                                                          ggCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACg                                             <  1:256672/139‑1 (MQ=255)                                                                           gCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACgg                                            >  1:115518/1‑139 (MQ=255)                                                                            cTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACggg                                           <  2:115734/139‑1 (MQ=255)                                                                                  cAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGt                                     <  1:330028/139‑1 (MQ=255)                                                                                                      cATTAGCCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGcc                 <  1:472352/139‑1 (MQ=255)                                                                                                               cTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGt                                     >  1:460783/1‑110 (MQ=255)                                                                                                               cTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGt                                     <  2:460783/110‑1 (MQ=255)                                                                                                               cTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAGTGCAttt        >  2:48039/1‑139 (MQ=255)                                                                                                                    gtgGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAGTGCATTTTCTCg   <  1:246629/139‑1 (MQ=255)                                                                                                                     tgGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAGTGCATTTTCTCGt  <  1:34549/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                             GTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAGTGCATTTTCTCGT  >  NC_000913/1461046‑1461300

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