Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,576,634	A→G	R79R (CGA→CGG)	BW25113_RS17825 →	YrhA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,576,634	0	A	G	100.0%	32.4 / NA	11	R79R (CGA→CGG)	BW25113_RS17825	YrhA family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/2);  total (9/2)

TAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAA  >  NZ_CP009273/3576556‑3576719                                                                               |                                                                                      tAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCct                                                                            >  1:393710/1‑90 (MQ=255) tAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCct                                                                            >  2:337856/1‑90 (MQ=255)    tATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCtctc                                                                         <  2:15/90‑1 (MQ=255)                                        tttAAAAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGt                                     >  1:172294/1‑90 (MQ=255)                                                       tggatTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAg                      >  2:411278/1‑90 (MQ=255)                                                         gatTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGaa                    >  2:465285/1‑90 (MQ=255)                                                                   tAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAAtgatga          >  1:243250/1‑90 (MQ=255)                                                                   tAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAAtgatga          >  2:225942/1‑90 (MQ=255)                                                                   tAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAAtgatga          >  2:73588/1‑90 (MQ=255)                                                                          ttACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATa   <  1:465285/90‑1 (MQ=255)                                                                           tACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATaa  >  2:27160/1‑90 (MQ=255)                                                                               |                                                                                      TAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAA  >  NZ_CP009273/3576556‑3576719

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