Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,320,398	G→A	A89V (GCG→GTG)	phnI ←	carbon‑phosphorus lyase core complex subunit PhnI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,320,398	0	G	A	100.0%	51.3 / NA	19	A89V (GCG→GTG)	phnI	carbon‑phosphorus lyase core complex subunit PhnI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (6/13);  total (6/13)

CCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCG  >  NC_000913/4320156‑4320628                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       ccTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCATGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGGCTTGATAAACTCCGGAGAGACGGCGTTCGAGACGCAGCCCGGTGGTGTCGAGCGGCTCGCTTACCACCAGCTTc                                                                                                                                                                                                                                <  1:137431/251‑1 (MQ=255)      tCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGATGAACACGTGCGGCGACGGGCGCCGTTCGCCAGCGAAGGAGGAGCGGGTAGGCGATTCGCCGATTGACAGCAGGGTCAAATCGAGCCGGCGATGGGGGTAGTCGAAGGGTGGGCAAAGCAGCTGGCCGCCGGGAAAGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGAAACCCGGAGGTGGCGAGCGGCGCGCTTACAACCAGCGTCGCCa                                                                                                                                                                                                                            >  1:164893/1‑250 (MQ=255)      tCAAACGTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCGGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCTATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTGTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCACCAGCTTCGc                                                                                                                                                                                                                              <  1:207141/248‑1 (MQ=255)                    ccccTGCCGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGTGATACGGCGTTCGAGACGCATCGCGGTGGTGTCGAGCGGCTCGCTTACCACCAGCTTCGCCAACGTGGTGCGGTAgn                                                                                                                                                                                                             <  1:57316/251‑2 (MQ=255)                                     ggCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGCGAAGGCGCTGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATCCGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCACCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATc                                                                                                                                                                                            >  1:50447/1‑251 (MQ=255)                                                                                           gcgcTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCACCAGCTTCGCCAACgtgn                                                                                                                                                                                                                      <  1:34739/171‑2 (MQ=255)                                                                                                   ccGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGCGCCAAGCAGCTGGCCGCAGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGCGACGCATCCCGGTGGAGTCGAGCGGCTCGCTTACCACCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAAt                                                                                                                              >  1:55222/1‑251 (MQ=255)                                                                                                                                       ttgtAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAAGGTCTTTATAAACGGCGGAGAGACTGCGGTCGAGACTCATCCCGGTGGTGTCTGGCGGCTCGCTTACCACCAGCTTCTCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTcag                                                                                                                                                            <  1:180870/184‑1 (MQ=255)                                                                                                                                         gtAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGAAGATACGGCGTTCGATACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCACCAGCTTCGCCAACGTGGTGCGGTAGGCGCGcagcag                                                                                                                                                                                                   <  1:327230/144‑1 (MQ=255)                                                                                                                                              cGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCACCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAgctgn                                                                                   <  1:195516/251‑2 (MQ=255)                                                                                                                                                           aaGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCACCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTggg                                                                      <  1:249154/251‑1 (MQ=255)                                                                                                                                                                gcTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCACCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACg                                                                 >  1:290698/1‑251 (MQ=255)                                                                                                                                                                        ccgGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGAGTCGAGCGGCTCACTTACCACCACCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCACCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTc                                                         >  1:285928/1‑251 (MQ=255)                                                                                                                                                                                                     aCGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCACCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTTCCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCgg                             <  1:274597/250‑1 (MQ=255)                                                                                                                                                                                                                    cagcccGGTGGTGTCGAGCGGCTCGCTTACCACCAGCTTCGCCAACGTGGGGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCgg                                              <  1:279291/215‑1 (MQ=255)                                                                                                                                                                                                                        ccGGTGGTGTCGAGCGGCTCGCTTACCACCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTcagcgc                                                                                                                                                         <  1:284262/107‑1 (MQ=255)                                                                                                                                                                                                                          ggtggtGTCGAGCGGCTCGCTTACCACCAGTTTCGCCAACGTGGTGCTGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGTCAGTTTGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCATTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGggcg        <  1:309711/250‑1 (MQ=255)                                                                                                                                                                                                                                gtCGAGCGGCTCGCTTACCACCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCg  >  1:175829/1‑250 (MQ=255)                                                                                                                                                                                                                                  cGAGCGGCTCGCTTACCACCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTgg                                                                       <  1:171229/179‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       CCTCTTCAAACTTCGCCAGCCCCTGACGCGCCAGCAGGCTGAAAACGTGCGGCGACGGCTGCTGTTCGCTGTCGGCGGTGGTCAGCGTCGGCGCTTCGCCGTTTGCCAGCAGGGTAAAATCGAGCAGGCGATGGGTGTAGTCGTAGGTTGGGCCAAGCAGCTGGCCGCCGGGAATGTCTTTATAAACGGCGGAGATACGGCGTTCGAGACGCATCCCGGTGGTGTCGAGCGGCTCGCTTACCGCCAGCTTCGCCAACGTGGTGCGGTAGGCGCGCAGCAGGAAAATCGCTTCAACGTTATCGCCGCTGGCCTGTTTCAGCGCCAGCGCCGCCAGTTCGCGGTCGGCAATGCCGCCTTCGGTCATCACGCGATCTACCGCGAGGTTAAGCTGCTGTTCAATCTGGGCGACGCTCAGTTCGGGCAAATCGGTATCGCCTCGGCGTCGGCTCTCTTGCAGGGCGTGGGCGGCGTCG  >  NC_000913/4320156‑4320628

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