Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,513,455	G→A	A285V (GCT→GTT)	fecB ←	ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,513,455	0	G	A	100.0%	31.8 / NA	12	A285V (GCT→GTT)	fecB	ferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (6/6);  total (6/6)

CACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATCGCGGTCATTTCACAACGGTAAGCGGCTGATGGTGGAAGATTTTTACCGTGTCAGCGGCAATACGCTCTGCAGCAAAAATACCGCGCATCCGCGCCCAGGTGTTACTGTCGACCGAAGCAACCTGCTGCTTCTGCGCGGCGGTTAACATCTGCCAGAGCGGATCTTGTTGCCAGCGTTTAACAATGCTCTCTTCGCGATAGTGGGCAACCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCG  >  NC_000913/4513242‑4513652                                                                                                                                                                                                                      |                                                                                                                                                                                                      cACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCCGCAGGGCGCGGGTTGCATATGCTCCGGAAAAAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATCGCGGTCATTTCACAACGGTAAGCGGCTGATGGTGGAAGATTTTTACCGTGTCAACGGCAATACGCACTGCAGCAAAAATACCGCGCATccc                                                                                                                                                                  >  1:120616/1‑250 (MQ=255)        gcTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCGTCTGCTCCGGAAACAGGAATGGCCGAATAGCAAACCTGTCTCCGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCCCCGGGTGTTGTATCGCGGTCATTTCACAACGGTAAGCGGCTGATGGTGGAAGATTGTTCCCGTGTCAACGGCAATACGCTCTGCAGCAAAAATACCGCGCATCCGCGCCCAg                                                                                                                                                           <  1:109742/251‑1 (MQ=255)           tCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATCGCGGTCATTTCACAACGGTAAGCGGCTGATGGTGGAAGATTTTTACCGTGTCAACGGCAATACGCTCTGCAGCAAAAATACCGCGCATCCGCGCCCAgn                                                                                                                                                          <  1:52319/249‑2 (MQ=255)                 aGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATCGCGGTCATTTCACAACGGTAAGCGGCTGATGGTGGAAGATTTTTACCGTGTCAACGGCAATACGTTCTGCAGCAAAAATACCGCGCATCCGCGCCCAGGTGTTACTg                                                                                                                                                  <  1:31986/251‑1 (MQ=255)                            gtgtCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATCGCGGTCATTTCACAACGGTAAGCGGCTGATGGTGGAAGATTTTTACCGTGTCAACGGCAATa                                                                                                                                                                                               >  1:56552/1‑195 (MQ=255)                                                                 gAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATCGCGGTCATTTCACAACGGTAAGCGGCTGATGGTGGAAGATTTTTATCGTGTCAACGGCAATACGCTCTGCAGCAAAAATACCGCGCATCCGCGCCCAGGTGTTACTGTCGACCGAAGCAACCTGCTGCTTCTGCGCGGCGGTTAACATCTGCCag                                                                                                  <  1:204536/251‑1 (MQ=255)                                                                               ngAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATCGCGGTCATTTCACAACGGTAAGCGGCTGATGGTGGAAGATTTTTACCGTGTCAACGGCAATACGCTCTGCAGCAAAAATACCGCGCATCCGCGCCCAGGTGTTACTGTCGACCGAAGCAACCTGCTGCTTCTGCGCGGCGGTTAACATCTGCCAGAGCGGATCttgttg                                                                                    >  1:54349/2‑251 (MQ=255)                                                                                    aGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATCGCGGTCATTTCACAACGGTAAGCGGCTGATGGTGGAAGATTTTTACCGTGTCAACGGCAATACGCTCTGCAGCAAAAATACCGCGCATCCGCGCCCAGGTGTTACTGTCGACCGAAGCAACCTGCTGCTTCTGCGCGGCGGTTAACATCTGCCAGAGCGGATCTTGTTGCCAgn                                                                               <  1:53221/251‑2 (MQ=255)                                                                                                                                     cccccACAGCAGCACCGGGTGTTTTATCGCGGTCATTTCACAACGGTAAGCGGCTGCTGGTGGAAGATTTTTACCGTGTCAACGGCAATACGATCTGCAGCAAAAATACCGCGCATCCGCGCCCAGGTGTTAc                                                                                                                                                    >  1:144415/1‑133 (MQ=255)                                                                                                                                                                 gcgGTCATTTCACAACGGTAAGCGGCTGATGGTGGAAGCATTTTACCGTGTCAACGGCAATCCGCTCTGCAGCAACAATACAGCGCATCCGCGACCAGGTGTTACTGTCGACCGAAGCAACCTGCTACTTCTGCGCGGCGGTTAACATCTGCCAGAGCGGATCTTGTTGACAGCGTTTAACCATGCTCTCTTCGCGATAGTGGGCAACCAGCAGCCAGGCAGGATTGACAGCCAGCAGTTGCTCCAGGCCg  >  1:159070/1‑251 (MQ=255)                                                                                                                                                                                                                gtgtCAACGGCAATACGCTCTGCAGCAAAAATACcg                                                                                                                                                                          <  1:2492/36‑1 (MQ=255)                                                                                                                                                                                                                gtgtCAACGGCAATACGCTCTGCAGCAAAAATACcg                                                                                                                                                                          >  2:2492/1‑36 (MQ=255)                                                                                                                                                                                                                      |                                                                                                                                                                                                      CACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGGCGCGGGTTGCATCTGCTCCGGAAACAGGAATGGCCGAGTAGCAAAACAGACTCAGCCAGAAAATAATAATAAGTGCTGCAACGGGAAGCCCCCACAGCAGCACCGGGTGTTTTATCGCGGTCATTTCACAACGGTAAGCGGCTGATGGTGGAAGATTTTTACCGTGTCAGCGGCAATACGCTCTGCAGCAAAAATACCGCGCATCCGCGCCCAGGTGTTACTGTCGACCGAAGCAACCTGCTGCTTCTGCGCGGCGGTTAACATCTGCCAGAGCGGATCTTGTTGCCAGCGTTTAACAATGCTCTCTTCGCGATAGTGGGCAACCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCG  >  NC_000913/4513242‑4513652

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