Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,776,100	G→A	E172K (GAA→AAA)	ydiF →	putative acetate‑CoA transferase YdiF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,776,100	0	G	A	100.0%	57.3 / NA	18	E172K (GAA→AAA)	ydiF	putative acetate‑CoA transferase YdiF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/10);  total (8/10)

CACTGGGGACAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATA  >  NC_000913/1775869‑1776331                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        cACTGGGGACAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAATACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAg                                                                                                                                                                                                                       <  2:96564/250‑1 (MQ=255)                              aCTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTcgcg                                                                                                                                                                                         <  2:234920/251‑1 (MQ=255)                                               aaatTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGACTGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGt                                                                                                                                                                        >  2:103126/1‑251 (MQ=255)                                                                 aCCCACAAGGTGTACTTACACAAACCTTACGCGCAGCCGCAGACCACAAGCCTGGTATTATTAGTGATATTGGCATCGGGCCATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAACGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACtttt                                                                                                                                                      >  1:191560/1‑251 (MQ=255)                                                                    cacaAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGaa                                                                                                                                                   <  1:173466/251‑1 (MQ=255)                                                                                 tACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTGGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTTACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCTCTACCACCTGCGACAGTTAAGGCTACGCCACTTTTGAAGATGAGGTGATg                                                                                                                                       <  1:290028/250‑1 (MQ=255)                                                                                    acaAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTAtc                                                                                                                                   <  2:287448/251‑1 (MQ=255)                                                                                                              accaGCCCGGAACTATTAGTGATATTTGCATCGGAACATATGTCGATCCACGCCAGCAAGGCGCCAAACTGAATGAAGTCACTAAAGAAGACCTGCTTAAACTGGTCGCGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCACGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATATCGACGCATTGGTTATTGCCCAGGCg                                                                                                         >  2:184310/1‑251 (MQ=255)                                                                                                               ccaGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTacc                                                                                                                                                                                    >  1:61344/1‑175 (MQ=255)                                                                                                               ccaGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTacc                                                                                                                                                                                    <  2:61344/175‑1 (MQ=255)                                                                                                                           tattaGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGACAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGATTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGTCGGTGCACAATAAcg                                                                                            >  2:52445/1‑251 (MQ=255)                                                                                                                                        ggCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAAATGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGtgat                                                                               <  1:302135/251‑1 (MQ=255)                                                                                                                                                      tGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCGCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGAGGTATCTCGACGCATTGGGTATTTCCCAGGCGGTGCACAAGAACGGCGGTATTGTGATGATGCAGGTGCAg                                                                  <  2:191560/250‑1 (MQ=255)                                                                                                                                                                        aGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACtt                                                                                                                                                        >  2:222438/1‑146 (MQ=255)                                                                                                                                                                        aGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACtt                                                                                                                                                        <  1:222438/146‑1 (MQ=255)                                                                                                                                                                                      atgaaGTCACTAAAGAAGCCCTGATTAAACTGGTCGAGTTTGATAACAAAAAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGGTTGTGCAGGGGCAGAAAATTGATAAGaaaacaccgcagcactccaa                                  >  2:159343/1‑233 (MQ=255)                                                                                                                                                                                                                 aaCTGGTCGAGTTTGATAACAAAAAATATCTCTATTAAAACGCGATTGAACCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGATACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCAAAATAACGGCGGTATTGTGATGATGGAGGGGCAGGAAAAGGGTAAGAAAACAAAGCTGCCTCCTAAAACTGTCAGAAATCCGGGGTATCtggtg      >  1:294266/1‑251 (MQ=255)                                                                                                                                                                                                                     ggTCGAGTTTGATAACAAAAAATAGCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGata  <  2:1763/251‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        CACTGGGGACAATCGCCGCGTATTTCTGAACTCGCAGAACAAAATAAAATTATTGCTTATAACTACCCACAAGGTGTACTTACACAAACCTTACGCGCCGCCGCAGCCCACCAGCCTGGTATTATTAGTGATATTGGCATCGGGACATTTGTCGATCCACGCCAGCAAGGCGGCAAACTGAATGAAGTCACTAAAGAAGACCTGATTAAACTGGTCGAGTTTGATAACAAAGAATATCTCTATTACAAAGCGATTGCGCCAGATATTGCCTTCATTCGCGCTACCACCTGCGACAGTGAAGGCTACGCCACTTTTGAAGATGAGGTGATGTATCTCGACGCATTGGTTATTGCCCAGGCGGTGCACAATAACGGCGGTATTGTGATGATGCAGGTGCAGAAAATGGTTAAGAAAGCCACGCTGCATCCTAAATCTGTCCGTATTCCGGGTTATCTGGTGGATA  >  NC_000913/1775869‑1776331

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