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%	43.2 / NA	13	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 (5/8);  total (5/8)

TGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAGTGCATTTTCTCG  >  NC_000913/1461042‑1461299                                                                                                                                       |                                                                                                                            tgatgTGTTCCTCAGCTGTGCTAATCAGACGGGAGAAGACGACCGTCACGTAGCGCGGATGGCGACTTTACTGGGGCGACTGCCCAAGGGTGTCTTCTGCACAACCATTACCGCCGTTGGTGGTTCCGGGCTGGgcgcc                                                                                                                         >  2:300387/1‑138 (MQ=255)   atgTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCa                                                                                                                         <  1:385580/137‑1 (MQ=255)   atgTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCa                                                                                                                         >  2:385580/1‑137 (MQ=255)                                   gaagaTAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGggcg                                                                                       <  2:365543/139‑1 (MQ=255)                                                        cGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCg                                                                  <  1:64021/139‑1 (MQ=255)                                                               cGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACg                                                                                                        <  1:278073/94‑1 (MQ=255)                                                               cGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACg                                                                                                        >  2:278073/1‑94 (MQ=255)                                                                             gggCCGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTTATCGCCGGTTGCGTGGAGGCAATGTCAc                                             <  1:300387/139‑1 (MQ=255)                                                                                                     acaACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGtt                                                                                                                 <  1:63833/47‑1 (MQ=255)                                                                                                     acaACCATTAACCGCTTGTGTGGTTCCGGGCTGGGCGCACTGGGGtt                                                                                                                 >  2:63833/1‑47 (MQ=255)                                                                                                                   cTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGAt                                                                       <  1:417209/75‑1 (MQ=255)                                                                                                                   cTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGAt                                                                       >  2:417209/1‑75 (MQ=255)                                                                                                                   cTTGTGTGGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAGTGCAttt       >  1:126313/1‑139 (MQ=255)                                                                                                                        gtgGTTCCGGGCTGGGCGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAGTGCATTTTCTCg  <  1:387205/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                            TGATGTGATCCTCGGCTGTGCTAATCAGGCGGGAGAAGATAACCGTAACGTAGCCCGGATGGCGACTTTACTGGCGGGGCTGCCGCAGAGTGTTTCCGGCACAACCATTAACCGCTTGTGTGGTTCCGGGCTGGACGCACTGGGGTTTGCCGCACGGGCGATTAAAGCGGGCGATGGCGATTTGCTGATCGCCGGTGGCGTGGAGTCAATGTCACGGGCACCGTTTGTTATGGGCAAGGCAGCCAGTGCATTTTCTCG  >  NC_000913/1461042‑1461299

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