Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,306,168	C→T	G144S (GGC→AGC)	yjcS ←	linear primary‑alkylsulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,306,168	0	C	T	100.0%	17.9 / NA	8	G144S (GGC→AGC)	yjcS	linear primary‑alkylsulfatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (3/5);  total (3/5)

CAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATCACGC  >  NC_000913/4305973‑4306370                                                                                                                                                                                                    |                                                                                                                                                                                                           cAGGGGCAATCACCTGAACTTTGCCGGTTTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCTTTTCTCg                                                                                                                                                                                                     <  1:251376/203‑1 (MQ=255)                                                                                gtgGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCTTTTCTCGCTCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTccn                                                                      <  1:60603/251‑2 (MQ=255)                                                                                                                     tGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCTTTTCTCGCTCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCcg                                  <  1:121592/250‑1 (MQ=255)                                                                                                                                                     gtcttgtGCGGCATGGGGCGTCACCAGCTGGTCGATAACAATAATGCTTGTCTCGCTCTCAACGAACGTAATGTTAGAGATATCCTGGc                                                                                                                                                                   <  1:324867/83‑1 (MQ=255)                                                                                                                                                     ggCTTTTGCGGCAGGCGGCGTCCCCAGCGGGTCGATAACAATAATGCTTTTCTCGCTCTCAACGCACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTCATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATCACGc  >  1:118260/1‑250 (MQ=255)                                                                                                                                                                cAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCTTTTCTCGCTCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTgn                                                                                                        <  1:94663/137‑2 (MQ=255)                                                                                                                                                                   ncggcgTCACCAGCGGGTCGATAACAATAATGCTTTTCTCGCTCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATAc                                                                                                                                                     >  1:79298/2‑89 (MQ=255)                                                                                                                                                                                  ggTCGATAACAATAATGCTTTTCTCGCTCTCAACGAACGTAATGTTAGAGATATACTGGCCGCGCACCTGATACCTTTTCTCGGTGACTTTGAACAGGCCAGAAATACCGTTGCTTTGCGACTGACGACACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCAt       >  1:137002/1‑216 (MQ=255)                                                                                                                                                                                                    |                                                                                                                                                                                                           CAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTGAACAGGCCAGAAATACCGTTGATTTGCGACTGACGCCACAGGCTGGGGTTTACGGTTTCCGGCGCTGCGGCATTAATATCAAATTTGTAATCGTCCGCGCGGTAGTAAACTTTACCGTTCGCATCACGC  >  NC_000913/4305973‑4306370

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