Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,007,245	G→A	M328I (ATG→ATA)	ygeW →	putative carbamoyltransferase YgeW

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,007,245	0	G	A	100.0%	33.8 / NA	14	M328I (ATG→ATA)	ygeW	putative carbamoyltransferase YgeW
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/5);  total (9/5)

GGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAGTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACC  >  NC_000913/3007020‑3007454                                                                                                                                                                                                                                  |                                                                                                                                                                                                                  ggTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAGTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATACACTGCCTGCCAGCTGATATCAgcg                                                                                                                                                                                          >  1:133014/1‑251 (MQ=255)  gTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAGTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGCTCACGACGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATACACTGCCTGCCAGCTGatcggcccg                                                                                                                                                                                          >  1:150334/1‑243 (MQ=255)      ccAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAGTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATACACTGCCTGCCAGCTGATATCAGCGGCGt                                                                                                                                                                                      >  1:13120/1‑250 (MQ=255)                                gCAGAAATCGTTTATCCGAAGTCATGGGCACCTTAAAAAGTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGACGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACAGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATACACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTg                                                                                                                                                            >  1:130463/1‑250 (MQ=255)                                                                                               ctgcGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATACACTGCCTGCCAGCTGTTATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATatcat                                                                                            >  1:166011/1‑251 (MQ=255)                                                                                                      ngAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATACACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCcgcg                                                                                     >  1:15093/2‑251 (MQ=255)                                                                                                                             gCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGTAACTGACCCGTGATGGCGAAGCCCTGTACATACACTGCCTTCCAGCGTATATCAGCGGCGTATCCTGTAAAGAAGGGGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGc                                                                                                      <  1:83934/211‑1 (MQ=255)                                                                                                                               aCTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGCCTGGCATTGCACTGAAGCGATGATGGCACTGACCCGTGATGGCGAAGCCCTGTACATACACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGcc                                                             >  1:55096/1‑250 (MQ=255)                                                                                                                                           gtgtCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGCTGGCACTGACCCGTGATGGCGAAGCCCTGTACATACACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAg                                                      >  1:168259/1‑245 (MQ=255)                                                                                                                                           gtgtCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATACACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCan                                                <  1:29670/251‑2 (MQ=255)                                                                                                                                                 ggCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATACActgcctgc                                                                                                                                                                                                         <  1:78956/92‑1 (MQ=255)                                                                                                                                                          cgcgCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATACACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCTATGGTCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgc                                 <  1:110129/251‑1 (MQ=255)                                                                                                                                                                                         gagaTGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATACACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCATATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAAcc  <  1:170672/251‑1 (MQ=255)                                                                                                                                                                                                       gACCCGTGATGGCGAAGCCCTGTACCTACACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAACAATACCGTATCGCTAACTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTCAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGCAcg              >  1:63104/1‑225 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                  GGTCACCAGCATGGAAGAAGCCTTCAAAGACGCAGACATCGTTTATCCGAAGTCATGGGCACCTTACAAAGTGATGGAAGAGCGTACTGAATTGCTGCGTGCGAACGATCACGAAGGCTTAAAAGCACTGGAAAAACAGTGTCTGGCACAGAACGCGCAACACAAAGACTGGCATTGTACTGAAGAGATGATGGAACTGACCCGTGATGGCGAAGCCCTGTACATGCACTGCCTGCCAGCTGATATCAGCGGCGTATCCTGTAAAGAAGGTGAAGTGACTGAAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACC  >  NC_000913/3007020‑3007454

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