Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,227,656	G→A	intergenic (‑184/‑25)	pliG ← / → ycgL	inhibitor of g‑type lysozyme/PF05166 family protein YcgL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,227,656	0	G	A	100.0%	42.2 / NA	14	intergenic (‑184/‑25)	pliG/ycgL	inhibitor of g‑type lysozyme/PF05166 family protein YcgL
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (4/10);  total (4/10)

AGCAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATGAACTGTTATAATATAACAATCCCTAACGGGAGAATTCCAGTGGCTTTCCTGCGTTCGGTAAAATCACGTTGCGGCATTTTTGCCGGGGGGATTACTTATAAAGGGGGCATG‑GACTCAACAAACCGTAGAAATCTACAGCTGGCTGGCTTGTGGCGCGGGTTTCATTGATGTTAGCTTATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTAA  >  NC_000913/1227438‑1227895                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                agcaATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATGAACTGTTATAATATAACAATCCCTAACGGGAGAATTCCAGTGGCTTTCCTGCGTTCGGTAAAATCACGTTGCGGCATTTTTGCCGGGGGGATTACTTATAAAGGGGGCATG‑GACTCAACAAACCGTAGAAATCTACAGCTGGCTGGCTTGTGACGCGGGTTTCATTGATGTTAGCTTATGCCaa                                                                                                                                                                                                                  >  1:3680/1‑250 (MQ=255)           ccTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATGAACTGTTATAATATAACAATCCCTAACGGGAGAATTCCAGTGGCTTTCCTGCGTTCGGTAAAATCACGTTGCGGCATTTTTGCCGGGGGGATTACTTATAAAGGGGGCATG‑GACTCAACAAACCGTAGAAATCTACAGCTGGCTGGCTTGTGACGCGGGTTTCATTGATGTTAGCTTATGCCAAAACCTGGCATc                                                                                                                                                                                                       <  1:77176/251‑1 (MQ=255)                                    tGGTAGTCTCCTGTTCAGAATGAGTACGATGAACTGTTATAATATAACAATCCCTAACGGGAGAATTCCAGTGGCTTTCCTGCGTTCGGTAAAATCACGTTGCGGCATTTTTGCCGGGGGGATTACTTATAAAGGGGGCATG‑GACTCAACAAACCGTAGAAATCTACAGCTGGCTGGCTTGTGACGCGGGTTTCATTGATGTTAGCTTATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTtgtg                                                                                                                                                                              >  1:239763/1‑251 (MQ=255)                                           ctcCTGTTCAGAATGAGTACGATGAACTGTTATAATATAACAATCCCTAACGGGAGAATTCCAGTGGCTTTCCTGCGTTCGGTAAAATCACGTTGCGGCATTTTTGCCGGGGGGATTACTTATAAAGGGGGCATG‑GACTCAACAAACCGTAGAAATCTACAGCTGGCTGGCTTGTGACGCGGGTTTc                                                                                                                                                                                                                                       <  1:139447/187‑1 (MQ=255)                                           ctcCTGTTCAGAATGAGTACGATGAACTGTTATAATATAACAATCCCTAACGGGAGAATTCCAGTGGCTTTCCTGCGTTCGGTAAAATCACGTTGCGGCATTTTTGCCGGGGGGATTACTTATAAAGGGGGCATG‑GACTCAACAAACCGTAGAAATCTACAGCTGGCTGGCTTGTGACGCGGGTTTc                                                                                                                                                                                                                                       >  2:139447/1‑187 (MQ=255)                                                                  gAACTGTTATAATATAACAATCCCTAACGGGAGAATTCCAGTGGCTTTCCTGCGTTCGGTAAAATCACGTTGCGGCATTTTTGCCGGGGGGATTACTTATAAAGGGGGCATG‑GACTCAACAAACCGTAGAAATCTACAGCTGGCTGGCTTGTGACGCGGGTTTCATTGATGTTAGCTTATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCCAGCGTGACCAGa                                                                                                                                                <  2:13043/251‑1 (MQ=255)                                                                        ttATAATATAACAATCCCTAACGGGAGAATTCCAGTGGCTTTCCTGCGTTCGGTAAAATCACGTTGCGGCATTTTTGCCGGGGGGATTACTTATAAAGGGGGCATG‑GACTCAACAAACCGTAGAAATCTACAGCTGGCTGGCTTGTGACGCGGGTTTCATTGATGTTAGCTTATGTCAAAGCCTGGCATCGTAAAGAGGAAAAGGATGTTTTGTGTGATGTATCGAagcaggaatcttggccagtcctat                                                                                                                                           <  2:44845/250‑20 (MQ=255)                                                                               ataaCAATCCCTAACGGGAGAATTCCAGTGGCTTTCCTGCGTTCGGTAAAATCACGTTGCGGCATTTTTGCCGGGGGGATTACTTATAAAGGGGGCATG‑GACTCAACAAACCGTAGAAATCTACAGCTGGCTGGCTTGTGACGCGGGTTTCATTGATGTTAGCTTATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATg                                                                                                                                    <  1:184461/250‑1 (MQ=255)                                                                                           aaCGGGAGAATTCAAGTGGCTTTCCTGCGTTCGGTAAAATCACGTTGCGGCATTTTTGCCGGGGGGATTACTTATAAAGGGGGCATG‑GACTCAACAAACCGTAGAAATCTACAGCTGGCTGGCTTGTGACGCGGGTTTCATTGATGTTAGCTTATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGaaaaaaa                                                                                                                          >  1:82611/1‑248 (MQ=255)                                                                                                   aagtCCAGTGGCTTTCCTTCGTTCGGTAAAATCACGTTGCGGCATTTTTGCCGGGGGGATTACTTATAAAGGGGGCATG‑GACTCAACAAACCGTAGAAATCTACAGCTGGCTGGCTTGTGACGCGGGTTTCATTGATGTTAGCTTATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAAGACGAttttt                                                                                                               <  2:38273/248‑1 (MQ=255)                                                                                                                                                 ttttccgggggattgtcTTATAAAGGGGGCATGCGA‑TAAACTAACCGTTGAAATCTACAGCTGGCTGGCTTGTGACGCGGTTTTCATTGATGTTAGCTTATGCCAAAACCTGGACTCTTAAAGAGTAAAAGGATGTTTTGTGAGATTTATCGAAGCAGCAAGCGTGCCCAGACCTATTAATATGTCGAAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGACCTGATGAAAGGATTTTGTCGTCCTCAgc                                                                  <  2:7328/234‑2 (MQ=255)                                                                                                                                                                       aaGGGGGCATG‑GACTCAACAAACCGTAGAAATCTACAGCTGGCTTGCTTGTGACGCGGGTTTCATTGATGTTAGCTTATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAAGACGATTTTTCGCGTGTTCCTGAGTAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGac                                           <  2:53284/251‑2 (MQ=255)                                                                                                                                                                                                gTAGAAATCTACAGCTGGCTGGCTTGTGACGCGGGTTTCATTGATGTTAGCTTATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGcc                   <  2:3680/251‑1 (MQ=255)                                                                                                                                                                                                                 gctggctTGTGACTCGGGTTTCATTGATGTTAGCTTATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATTGTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATTGCGAAAAAAAAGACGAGGTTTCGCGTGTTCCTGAGGGACTGATGAAATTTTTGGGTCAGCGTTAGTTATCGATGATTCTTCCGCTGGGTGGGCGGGAGAAACTGGTCAATGTCGATCTTGTAATGGTTaa  <  2:258355/251‑1 (MQ=255)                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                AGCAATACAGCCTTCCTGATGCTCTTGATTTTCATTGGTAGTCTCCTGTTCAGAATGAGTACGATGAACTGTTATAATATAACAATCCCTAACGGGAGAATTCCAGTGGCTTTCCTGCGTTCGGTAAAATCACGTTGCGGCATTTTTGCCGGGGGGATTACTTATAAAGGGGGCATG‑GACTCAACAAACCGTAGAAATCTACAGCTGGCTGGCTTGTGGCGCGGGTTTCATTGATGTTAGCTTATGCCAAAACCTGGCATCTTAAAGAGTAAAAGTATGTTTTGTGTGATTTATCGAAGCAGCAAGCGTGACCAGACCTATTTATATGTCGAAAAAAAAGACGATTTTTCGCGTGTTCCTGAGGAACTGATGAAAGGTTTTGGTCAGCCTCAGTTAGCGATGATTCTGCCGCTGGATGGGCGGAAGAAACTGGTCAATGCCGATATTGAAAAGGTTAA  >  NC_000913/1227438‑1227895

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: ‑