Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,576,911	C→A	G237C (GGC→TGC)	yddB ←	putative TonB‑dependent receptor YddB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,576,911	0	C	A	100.0%	36.4 / NA	13	G237C (GGC→TGC)	yddB	putative TonB‑dependent receptor YddB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (6/7);  total (6/7)

AGACTGCGGGAAGGTGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGC  >  NC_000913/1576693‑1577155                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     aGACTGCGGGAAGGTGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAAACGGCGGTAACGCAGAAGTTATCAGCGAGTTCCTGATTAAACGAc                                                                                                                                                                                                                       >  1:5401/1‑250 (MQ=255)                                     ggAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCAGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGtaa                                                                                                                                                                                   <  1:178218/250‑1 (MQ=255)                                                        gTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCAGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGcc                                                                                                                                                               >  1:131125/1‑251 (MQ=255)                                                                                            gTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATCATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCAGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTCAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATg                                                                                                                           >  1:131778/1‑251 (MQ=255)                                                                                                       gaaTGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGGGACAATGCCGTCATTCGAAACATAAGCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCAGAAGTTATCAGCGATTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTGCTGCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAAc                                                                                                                 <  1:101082/249‑1 (MQ=255)                                                                                                             aGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCAGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCg                                                                                                          >  1:90958/1‑251 (MQ=255)                                                                                                                                                     aTGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCAGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTatcatc                                                                  <  1:91239/251‑1 (MQ=255)                                                                                                                                                                      agaaTCCGCGCGGGTGATATCAGACTGGCGTCGCGATAAACCGGCGGTAACGCAGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCTTCAGCGTTGAAGCGTTTGa                                                 <  2:157591/249‑1 (MQ=255)                                                                                                                                                                                   gTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCAGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAAt                                    <  1:73600/251‑1 (MQ=255)                                                                                                                                                                                            gACTGGCGGCGCGAGAAACCTGCGGTAACGCAGAATTTATCAGCGAGTTCCTGATTAAACGAAAAGGTATAAAATTTCGTTTTAAAATCTGGGGAGTAATAAGTACTTCCTTAAGAGCCTTGATTAAATTCGCTCTTTTTATCCTCATCGATATGCGATGTTATCCAGTACTAACGCTTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTAAAGCGTTTGATCTTTGTATCAGTTACCCCGcc                           <  2:182809/251‑1 (MQ=255)                                                                                                                                                                                                          tAAACCGGCGGTACCGCAGAAGATATCAGCGAGGTCCTGATTAAACGACAAGGTATAACAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGCAAAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACCCGTAGTGCGATACCCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGCTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCa              >  2:100078/1‑250 (MQ=255)                                                                                                                                                                                                          tAAACCGACGGTAACGCAGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCa              >  2:57098/1‑250 (MQ=255)                                                                                                                                                                                                                      aaCGCAGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCCCCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGc  <  1:57098/250‑1 (MQ=255)                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     AGACTGCGGGAAGGTGCTGGTATTATAATCACGGCTGGAGCCGGTATATTTTAAGGTTAAATCGTGGGTAAAGCGGTCGCTGGCAAACCAGGTAAATTTGCTTAATGCAGTATCGATAACGTTTTTATACTGTGCCCGACCGGCGACAATGCCGTCATTCGAAACATAATCCGCGCGGGTGATATCAGACTGGCGGCGCGATAAACCGGCGGTAACGCCGAAGTTATCAGCGAGTTCCTGATTAAACGACAAGGTATAAAAGTTCTTTTTAAAATCTGGGGAGTAATAAGTACTTCCTGAAGAGCCTTGATTAAATGCGCTCTTGTTATTCTCATCGATATGCGATGTTAACCAGTCCGAACGCGTAGTGCGATAACCCAATTTCACCTTGCTATCATCAGCGTTGAAGCGTTTGATCTTTGCATCAATTACCCCGCCATTGAAGCGACCAAATTCAACCGGC  >  NC_000913/1576693‑1577155

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