Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,515,480	C→T	100%	G400S (GGC→AGC)	fecA ←	ferric citrate outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,515,480	0	C	T	100.0%	42.4 / NA	13	G400S (GGC→AGC)	fecA	ferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/6);  total (7/6)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTG  >  NC_000913/4515227‑4515685                                                                                                                                                                                                                                                              |                                                                                                                                                                                                              aGATAACAGGGGTGCGCCTCGGTGCCGGAACGTGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTATAGTGACGCATTTCATTCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATTCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGgt                                                                                                                                                                                             <  1:214053/272‑1 (MQ=255)                                    aTCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCCTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTAAATACCGCGCACCCAGTAGTTACGCGGAGAGAGGGGGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCGGAATGTTGAATTTATgctggctat                                                                                                                                                          >  1:87922/1‑270 (MQ=255)                                                                   aaCTTCCCGCTTCTTTTGGCGGTGTAGTAACGCATTTCATGCTTTTATTCATTCAAAGAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCAAAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAg                                                                                                                           <  1:240869/272‑1 (MQ=255)                                                                         ccGCTGCTTGTGGCGTTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGc                                                                                                                      <  1:142767/271‑1 (MQ=255)                                                                                                                                                       gcgGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAg                                                                                                                                                                                                      >  1:45229/1‑113 (MQ=255)                                                                                                                                                        cgGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAgg                                                                                                                                                                                                     <  2:45229/113‑1 (MQ=255)                                                                                                                                                                                                          aCCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAg                                                                             >  1:227618/1‑183 (MQ=255)                                                                                                                                                                                                          aCCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAg                                                                             >  1:228384/1‑183 (MQ=255)                                                                                                                                                                                                           cccAGTAGTTACGCGGCGGGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGc                                                                            <  2:227618/183‑1 (MQ=255)                                                                                                                                                                                                           cccAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGc                                                                            <  2:228384/183‑1 (MQ=255)                                                                                                                                                                                                                         ggcgAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGACCCCCTGAATGTTGAATTTATGCTGGCTATCTGGCTGGACCTGGTAGCCCAGGCTCGCCAGCTTGCGCCGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGACCGGCGTCGTCATCCACGCGAGCCAGGCCACCGGGCATTGCGGCCTTCCCGTCGGAATAtttaagaagg     >  1:65465/1‑232 (MQ=255)                                                                                                                                                                                                                         ggcgAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTTGCTGTCTGGCTGGAACTGGTAACCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGACGGCTTCACCGTCGTACTATTTCAGCAGGCTg  >  1:224080/1‑243 (MQ=255)                                                                                                                                                                                                                                                           aGCTGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAg                                                >  1:60933/1‑163 (MQ=255)                                                                                                                                                                                                                                                              |                                                                                                                                                                                                              AGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTGCTCCAGGTAGCCGCTGCGCAGGGTTTGGGTGTAGAACCCCTGAATGTTGAATTTATGCTGGCTGTCTGGCTGGAACTGGTAGCCCAGGCTCGCCAGCTTGCGACGACCCCAGAAGCGGTCATACGGGCGGGTGGATTGCCAGCGATCGGCGTCGTAATCCGCGCGAGACAGGCCACCGGGCATGTCGGCTTCACCGTCGTAATATTGCAGCAGGCTG  >  NC_000913/4515227‑4515685

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