Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	852,211	(T)7→6	coding (387/1584 nt)	opgE ←	phosphoethanolamine transferase OpgE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	852,211	0	T	.	100.0%	35.6 / NA	9	coding (387/1584 nt)	opgE	phosphoethanolamine transferase OpgE
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (4/5);  total (4/5)

CCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACGAAATACCGGATAGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCG  >  NC_000913/851977‑852458                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        ccAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACGAAATACCGGATAGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTT‑CGCCTGATTGCTATTAc                                                                                                                                                                                                                                        <  1:53826/251‑1 (MQ=255)                                                                                           aaaTACCGGATAGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTT‑CGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTg                                                                                                                                              >  1:233678/1‑251 (MQ=255)                                                                                              tACCGGATAGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTT‑CGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCa                                                                                                                                            >  2:16506/1‑250 (MQ=255)                                                                                                      aGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCGCTTTGTG‑CGCCTGATTGCGATTACGCCTAGGTGAAAAATAAATAAGCCGACCAGGTAGCGGTACCACGGAAAAAGCGTCAGGAGGTCGGTCGACGCTTCCATATTGGTGGAGTGCAATGCCAGCa                                                                                                                                   <  1:352703/251‑1 (MQ=255)                                                                                                                    aGTTTACATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTT‑CGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAAtt                                                                                                                      >  2:135430/1‑250 (MQ=255)                                                                                                                         aaaTACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTT‑CGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGgcgc                                                                                                                <  2:1013/251‑1 (MQ=255)                                                                                                                                                                               acTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTT‑CGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACAcc                                                          <  2:137666/251‑1 (MQ=255)                                                                                                                                                                                     tACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTT‑CGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGAcc                                                      >  1:24244/1‑249 (MQ=255)                                                                                                                                                                                               aGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTT‑CGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGt                                          <  1:213441/251‑1 (MQ=255)                                                                                                                                                                                                                                  tctcttttttCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATagcaaaagcaaa         >  2:18887/5‑250 (MQ=255)                                                                                                                                                                                                                                         tttCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATAGTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAAGAGCCg  <  1:305810/250‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        CCAGGTATCTTTCATTCCGGAAAGTTTTGCCATCCGTTCTTGCTCTTCAATCACCTCGTTATTATTGACGATGACATCCTTAGCAAAACGAAATACCGGATAGCCAGTATCTTTCAGTTTAAATACGCCACCCCAGGCCAGGTTTTGCACGGGAGCAACAAAAAATGTCGCCACACTGAATACCAGACACAGGCTGTCGAAGGTATTCCAGCGCGCTTTTTCATTCTCTTTTTTTCGCCTGATTGCTATTACGCCGAGCGCAAAAATAAATAAGCCGACCAGGTAGCTGTACCACGGAAAAATCGTCAGGATTTCGGTCGACTCTTCCATATTGGTGGAGTGCAATGCCAGCAATGTATTGAAATTAGGCGCGCCGTAGGCCTGAGCAAAAGGGAAATAACAAGCCGCCACCAGCGAACTGACACCGACCAGAACTTTTTGTACGCGAGGCAATGTTCGCCATAGCAAAAGCAAAATAGCCG  >  NC_000913/851977‑852458

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