Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,511,280	A→G	F60L (TTT→CTT)	fecE ←	ferric citrate ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,511,280	0	A	G	100.0%	121.4 / NA	34	F60L (TTT→CTT)	fecE	ferric citrate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/20);  total (14/20)

TGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGA  >  NC_000913/4511204‑4511357                                                                             |                                                                              tGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACg                                                                           >  1:911836/1‑81 (MQ=255)          ggCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGAc                                                                  <  1:911778/81‑1 (MQ=255)             agcagcGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGc                                                               >  1:911842/1‑81 (MQ=255)                   gAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCAtt                                                         <  2:911779/81‑1 (MQ=255)                    aaaGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTa                                                        <  1:911780/81‑1 (MQ=255)                      aGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTaa                                                       >  2:911838/1‑80 (MQ=255)                       gCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTaaaa                                                     <  1:911781/81‑1 (MQ=255)                              ggCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGa                                              <  1:911782/81‑1 (MQ=255)                               gCCAACTGGCGCGATGAGGGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAGGCCGCGaa                                             <  2:911783/81‑1 (MQ=255)                               gCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGaa                                             <  2:911784/81‑1 (MQ=255)                               gCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGa                                              >  1:911833/1‑80 (MQ=255)                                ccAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGaaa                                            <  1:911785/81‑1 (MQ=255)                                       gcgcGATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAg                                       <  2:911786/79‑1 (MQ=255)                                        cgcgATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGttt                                    <  2:911787/81‑1 (MQ=255)                                        cgcgATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGtt                                     >  2:911840/1‑80 (MQ=255)                                          cgATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTaa                                  >  1:911829/1‑81 (MQ=255)                                           gATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTaa                                  >  1:911791/1‑80 (MQ=255)                                           gATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAAc                                 <  1:911790/81‑1 (MQ=255)                                           gATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAAc                                 <  2:911788/81‑1 (MQ=255)                                           gATGAGAGCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAAc                                 <  2:911789/81‑1 (MQ=255)                                               agagCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCg                             >  1:911835/1‑81 (MQ=255)                                               agagCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCg                             <  1:911792/81‑1 (MQ=255)                                                  gCATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCg                          <  1:911793/81‑1 (MQ=255)                                                   cATATTTATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGa                         <  1:911794/81‑1 (MQ=255)                                                         tATGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTcc                    <  1:911795/80‑1 (MQ=255)                                                          aTGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCg                  >  2:911844/1‑81 (MQ=255)                                                           tGGGATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCg                  <  2:911796/80‑1 (MQ=255)                                                            gggATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGc                 >  1:911834/1‑80 (MQ=255)                                                              gATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAAc              <  1:911798/81‑1 (MQ=255)                                                              gATTATCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAAc              <  1:911797/81‑1 (MQ=255)                                                                  aTCGCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGt           >  1:911849/1‑80 (MQ=255)                                                                     gCCGAGAAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAg        >  2:911852/1‑80 (MQ=255)                                                                          gaAGTATGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGCTTTCCCGCAACCGTTAGGACCGa  >  1:911839/1‑81 (MQ=255)                                                                          gaAGTACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGa  >  1:911841/1‑81 (MQ=255)                                                                             |                                                                              TGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCGAGAAATACGGTGCCAGACTGCGGCATTAAAAGCCGCGAAAAACAGTTTAACAGCGTCGATTTCCCGCAACCGTTAGGACCGA  >  NC_000913/4511204‑4511357

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