Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,325,644	C→T	S33N (AGC→AAC)	yjdN ←	PF06983 family protein YjdN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,325,644	0	C	T	100.0%	34.0 / NA	13	S33N (AGC→AAC)	yjdN	PF06983 family protein YjdN
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/6);  total (7/6)

GATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCC  >  NC_000913/4325426‑4325871                                                                                                                                                                                                                           |                                                                                                                                                                                                                                    naTTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGCTTTTGGCATTTCGCCGAAGTTGATGTTATAGAGCAgatacg                                                                                                                                                                                                                >  1:208871/2‑235 (MQ=255)                                   ttGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGt                                                                                                                                                                               <  1:294236/239‑1 (MQ=255)                                                 ccAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATAATGCTGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGCTGGCGGTATCGGGAACTTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCaaa                                                                                                                                                            >  1:289182/1‑244 (MQ=255)                                                         aaGCCGGAGCCGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGAACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTg                                                                                                                                                                                                >  1:76806/1‑200 (MQ=255)                                                                                           nnnnnntcnnnnNGATNTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCcg                                                                                                                                                                                                                <  2:208871/144‑1 (MQ=255)                                                                                                                   cGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCACTCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAAttt                                                                                    >  1:344153/1‑250 (MQ=255)                                                                                                                                        gggCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGGTGGTGACCTCAGTTAAGAGAAGc                                                                                                           >  1:98988/1‑206 (MQ=255)                                                                                                                                           cGATGGCGGTATCGGGAAATTGCATTTCGGAAGGGCAGTTCTCGGCGCTGGCCTGAGCTGATTTTGGCATTTCGCCGACGTTGATTTTATCGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGACGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGAATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGCGTAGTTCGAAATTTATACAGATGAGAGACGGCCTGATa                                                            >  1:130704/1‑250 (MQ=255)                                                                                                                                                gCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGccc                                                                                                                                                                                                           <  1:94316/102‑1 (MQ=255)                                                                                                                                                  ggTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACGGATGAGAGGCGGCCTGATAAGAcgcgc                                                    <  1:100688/251‑1 (MQ=255)                                                                                                                                                   gTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGAcgcgcn                                                   <  1:104280/251‑2 (MQ=255)                                                                                                                                                                                  tctcGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGc                     <  1:108931/250‑1 (MQ=255)                                                                                                                                                                                                    cTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGCGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATAGCGTCGGCACAGTTACCGGCAAAAGCGAGGTAGGGACTTAACGGCCTGATGGTGACCTCCGTTAAGAGAAGCCAGTTAAGAGTAGTTCGCAATTTATACAGATGAGAGGCGGCCTGATAAGCCGCGCGAGCGTCGCATCAGGCCGTCGGCACTGTTGCCGGATGCGGCGTAAACGcc  >  1:355752/1‑251 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                                    GATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCC  >  NC_000913/4325426‑4325871

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