Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,285,630	T→C	N581S (AAC→AGC)	acs ←	acetyl‑CoA synthetase (AMP‑forming)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,285,630	0	T	C	100.0%	73.8 / NA	22	N581S (AAC→AGC)	acs	acetyl‑CoA synthetase (AMP‑forming)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (7/15);  total (7/15)

AGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTGACCTTTAATATTGTGCGGA  >  NC_000913/4285495‑4285722                                                                                                                                        |                                                                                             aGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCg                                                                                          <  1:631908/140‑1 (MQ=255)   gTTGCTGGTTTCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGa                                                                                        <  2:87309/140‑1 (MQ=255)                  cgccgcAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGtt                                                                         >  2:149433/1‑140 (MQ=255)                               cgcAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTGTTTACGCACCCAGCTGGGGACTTCTTCGTATAGTTCTGGTGACGGTtc                                                            >  2:177705/1‑140 (MQ=255)                                   gAATACGGCGCATAACTTTGCCTGACCGGCTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTtcctcc                                                        <  1:177705/140‑1 (MQ=255)                                              aTAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTATGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAg                                             >  2:415821/1‑140 (MQ=255)                                              aTAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAg                                             >  1:70075/1‑140 (MQ=255)                                                     tGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACg                                      <  1:149433/140‑1 (MQ=255)                                                                    tAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGcc                       <  2:351861/140‑1 (MQ=255)                                                                      ggCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTATGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTg                     <  1:415821/140‑1 (MQ=255)                                                                      ggCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTg                     <  2:562348/140‑1 (MQ=255)                                                                        cAGGGAGTCGGTCCAGTGCAGCCCGTCTGACGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTGAc                   <  1:303772/140‑1 (MQ=255)                                                                        cAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTGAc                   <  2:535341/140‑1 (MQ=255)                                                                        cAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTGAc                   <  2:70075/140‑1 (MQ=255)                                                                        cAGGGAGTCGGTCCAGTCCAGCTCGTCTGGCTTCGCCAGTGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTGAc                   <  2:173383/140‑1 (MQ=255)                                                                           ggAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTGAc                   >  1:527618/1‑137 (MQ=255)                                                                           ggAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTGAc                   <  2:527618/137‑1 (MQ=255)                                                                           ggAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTGACCtt                >  1:425078/1‑140 (MQ=255)                                                                                cGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTGACCTTTAata           <  2:482612/140‑1 (MQ=255)                                                                                    ccAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTGACCTTTAATATtgt       <  1:569930/140‑1 (MQ=255)                                                                                       gTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTGACCTTTAATATTGTGCg    >  2:648901/1‑140 (MQ=255)                                                                                         gcagcaCGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGCTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTGACCTTTAATATTGTGCGGa  <  1:569700/140‑1 (MQ=255)                                                                                                                                        |                                                                                             AGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTAAGCGTGACGTAGGCGTAGATCGCCTGACCTTTAATATTGTGCGGA  >  NC_000913/4285495‑4285722

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