Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,285,533	A→G	R613R (CGT→CGC)	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,533	0	A	G	100.0%	64.6 / NA	19	R613R (CGT→CGC)	acs	acetyl‑CoA synthetase (AMP‑forming)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/9);  total (10/9)

GGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTC  >  NC_000913/4285411‑4285664                                                                                                                           |                                                                                                                                    ggTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGcc                                                                                                                    >  2:109079/1‑140 (MQ=255)           gCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGtt                                                                                                          <  1:303070/140‑1 (MQ=255)           gCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGtt                                                                                                          >  1:83857/1‑140 (MQ=255)            cATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGttt                                                                                                         >  2:21633/1‑140 (MQ=255)                 cgATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGc                                                                                                    <  1:21633/140‑1 (MQ=255)                   aTAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAg                                                                                                  <  1:7736/140‑1 (MQ=255)                        cTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGt                                                                                             <  1:109079/140‑1 (MQ=255)                                ctTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGACGAGGCATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAg                                                                                     >  2:84082/1‑140 (MQ=255)                                          cTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGt                                                                           <  2:83857/140‑1 (MQ=255)                                                    aCGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGc                                                                 >  1:122152/1‑140 (MQ=255)                                                      gCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCa                                                               >  2:72732/1‑140 (MQ=255)                                                               ggCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCaa                                                      >  2:191446/1‑140 (MQ=255)                                                                                cGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAg                                     <  1:45785/140‑1 (MQ=255)                                                                                cGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAg                                     <  2:122152/140‑1 (MQ=255)                                                                                       tttGCTGTTTTCGCCCGCCCCAACTTTGCGCAGACTGAGACGGATAATTTTTCCGGAGCCGGTTTTAGGCCGGGAGTCGGTCCAGTGGTGCCCGTTTGCTGTATCCTGCGGGCCAATCTCTTTACGCACCTAGTTGCGGa                              <  2:352692/139‑1 (MQ=255)                                                                                       gTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGa                              >  1:352692/1‑140 (MQ=255)                                                                                                        ccgcAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCt             >  1:3163/1‑140 (MQ=255)                                                                                                         cgcAATTTTGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTg            >  2:286945/1‑140 (MQ=255)                                                                                                              ttttGCGCAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGAc       <  2:333662/140‑1 (MQ=255)                                                                                                                   cgcAGAATGCGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTtc  <  1:72732/140‑1 (MQ=255)                                                                                                                           |                                                                                                                                    GGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTC  >  NC_000913/4285411‑4285664

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