Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,513,044	T→C	S122S (TCA→TCG)	fecC ←	ferric citrate ABC transporter membrane subunit FecC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,513,044	0	T	C	100.0%	89.7 / NA	32	S122S (TCA→TCG)	fecC	ferric citrate ABC transporter membrane subunit FecC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/15);  total (17/15)

CCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGG  >  NC_000913/4512803‑4513285                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  ccGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAgaga                                                                                                                                                                                                                                           >  2:85249/1‑250 (MQ=255)    gCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAg                                                                                                                                                                                                                                       <  1:110341/251‑1 (MQ=255)                                      ggACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGAGGAACGCCAGAGCATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTggga                                                                                                                                                                                                      >  2:130238/1‑248 (MQ=255)                                               gCCTGCAAGCCCGAAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGAGATGCGGGTAAGGCACATACAAAAGGCCGAGAGCGCGATACCCGCGCGGCACAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCACCGCATACCGCGATGAACGACAGAGAATAACCTGCAATCGGCGGCGGACTCAGCGGGCTGGAAAgagacat                                                                                                                                                                                              >  2:221463/1‑243 (MQ=255)                                                 cTGCCAGCCAGTAACAGATGCGGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATCCCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGCAATCCGCCTCCTGCGGTCCTGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCACCGCGCTGGCAAGCGCCATCGcc                                                                                                                                                                                          >  2:125651/1‑251 (MQ=255)                                                                   ggCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTa                                                                                                                                                                        <  2:104308/250‑1 (MQ=255)                                                                     ccGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAAt                                                                                                                                                                      >  1:225466/1‑251 (MQ=255)                                                                     ccGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCGCACTCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAAt                                                                                                                                                                      <  1:39410/251‑1 (MQ=255)                                                                      cGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCGGTTAATg                                                                                                                                                                     <  1:58516/251‑1 (MQ=255)                                                                        tAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCACATACAAAAGGCCGACAGCGCGATACCCGCGAGGCTCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCCTGACCAGCAGCCAGCTCACGCCGCCCCCACATGCCGCGAGGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCCGCGCCGCGCCGCTGTTACTGcc                                                                                                                                                                   >  2:168793/1‑251 (MQ=255)                                                                                      ggCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGgag                                                                                                                                                     <  1:215412/251‑1 (MQ=255)                                                                                      ggCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGgag                                                                                                                                                     <  2:139916/251‑1 (MQ=255)                                                                                      ggCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGAAAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGCTGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCACCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAgg                                                                                                                                                       >  1:38456/1‑249 (MQ=255)                                                                                          agcagccgGGCGATGCGGGTAACGCCCATACAAAAGGCCGACCGCGCGATACCCGCGAGGATCCGTTTGTTTCTGTCATGGGTATGCCGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGCTGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGc                                                                                                                                                 >  2:31258/1‑251 (MQ=255)                                                                                                     gATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGGTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATTACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCTATGAACGACAGAGAATAGCCTGCAATCTGCGTCGGACTCATCGCGCTGGTAAGCGCCATAGCCAgcaccgcgcc                                                                                                                                                                               <  1:242210/210‑1 (MQ=255)                                                                                                      aTGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTtgtg                                                                                                                                     <  1:234534/251‑1 (MQ=255)                                                                                                                           aagaccGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGt                                                                                                                <  1:90797/247‑1 (MQ=255)                                                                                                                                     gcgcGATACCCGCGAGGATCAGTTTGTTTCTGTCAGGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGGCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTAAAGGAGAGGCGATTGGGTTGTGGGTCAGGGTTTGCAGCTGCGTGCCCGCGAgc                                                                                                      <  2:38456/251‑1 (MQ=255)                                                                                                                                      ggagATACCCGCGAGGATCAATTTGTTACTGTCAGGGGAAAGAAGAAAGCCGCCACCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCACCCGCATGCCGCGAGGAACGACAGAGAAGAGCCGGCAATCGGCGGCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGGCGCGCGGCGGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGTGTCAAGGTGTGCAGGAGCGTGCCGGCGAgcg                                                                                                     <  2:31632/248‑1 (MQ=255)                                                                                                                                        cgctaaccGCGCGGACCAGATTGTTTCTATCATGGGAATGCAGAAATCCGCCTCCTACGGCCATGACCAGCAGCCAGCTCACGCCGCCCCCGCCTGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGGCGCACTCAGCGCGCTGGTAAGCGCAATAGCCAGCGCCGCGCCGCTGTTAATGCCGCGCAGTGAAGGCGAGGCCTTTGGGGTGGGGGTCAGGGTTTGCAGCAGCGAGCCCGCGCGCGcc                                                                                                   >  2:159666/7‑251 (MQ=255)                                                                                                                                           tACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCCTGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAgg                                                                                                >  1:34608/1‑251 (MQ=255)                                                                                                                                              ccGCGAGGATCAGTTTGTTTATGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGACAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCtt                                                                                             >  2:136243/1‑251 (MQ=255)                                                                                                                                                  gAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTgcgc                                                                                         >  2:81045/1‑251 (MQ=255)                                                                                                                                                              gtttCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGACGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCCATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCGGAACgg                                                                             >  2:120806/1‑251 (MQ=255)                                                                                                                                                              gtttCTGTCATGGGTATGACGAAACCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCTGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGCTCAGAACgg                                                                             <  2:184372/251‑1 (MQ=255)                                                                                                                                                                      cATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAgg                                                                     >  1:94291/1‑251 (MQ=255)                                                                                                                                                                          ggTATGACGAGATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGACGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGGTCAGAACGGCGACCAGGCTTc                                                                 >  2:139850/1‑251 (MQ=255)                                                                                                                                                                                        cgccgcCTGCGGTCATGACCAGCATCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCATGGTTGGCAGCAGCGTGCCCGCGAGCTCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAgg                                                   <  1:172859/251‑1 (MQ=255)                                                                                                                                                                                             cctGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGtttt                                               <  1:81045/250‑1 (MQ=255)                                                                                                                                                                                                                       cACGCCGCCCCCGCATGCCGCGATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCg                    <  2:171925/251‑1 (MQ=255)                                                                                                                                                                                                                                   gCATGCCGCGATGAACGACAGAGAATAGCATGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGgc         >  1:198626/1‑250 (MQ=255)                                                                                                                                                                                                                                         cgcgATGAACGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAggg  >  1:86071/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  CCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTCGGCCAGCAGCAGGGTGATGCGGGTAAGGCCCATACAAAAGGCCGACAGCGCGATACCCGCGAGGATCAGTTTGTTTCTGTCATGGGTATGACGAAATCCGCCTCCTGCGGTCATGACCAGCAGCCAGCTCACGCCGCCCCCGCATGCCGCGATGAATGACAGAGAATAGCCTGCAATCGGCGTCGGACTCAGCGCGCTGGTAAGCGCCATAGCCAGCGCCGCGCCGCTGTTAATGCCGAGCAGTGAAGGAGAGGCCATTGGGTTGTGGGTCAGGGTTTGCAGCAGCGTGCCCGCGAGCGCCAGGCTTGCGCCGATCAGAACGGCGACCAGGCTTCGTGGCAAACGAAGGTTTTGCACCAGCGCTTCTGGTAGCGTTGGCGTGTGTCCAGGCAGCAGGG  >  NC_000913/4512803‑4513285

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