Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,514,875	G→A	S601S (TCC→TCT)	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,514,875	0	G	A	100.0%	67.8 / NA	21	S601S (TCC→TCT)	fecA	ferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/4);  total (17/4)

CTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAG  >  NC_000913/4514753‑4515002                                                                                                                           |                                                                                                                                cTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTagag                                                                                                                              >  2:465131/1‑126 (MQ=255) cTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTagag                                                                                                                              <  1:465131/126‑1 (MQ=255) cTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGtt                                                                                                                >  1:506912/1‑140 (MQ=255)    gCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTaaa                                                                                                             >  2:776002/1‑140 (MQ=255)       tAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGtt                                                                                                          >  1:115571/1‑140 (MQ=255)           tGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATc                                                                                                      <  2:19091/140‑1 (MQ=255)                   ttATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGTCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAg                                                                                              >  2:606470/1‑140 (MQ=255)                    tATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGc                                                                                             <  2:760056/140‑1 (MQ=255)                                gCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGCATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTccgcc                                                                                 >  1:773917/1‑140 (MQ=255)                                                aTCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGTGCGCcgtcgt                                                                 >  2:114672/1‑140 (MQ=255)                                                aTCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCcgtcgt                                                                 >  1:627205/1‑140 (MQ=255)                                                 tCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCcgtcgtc                                                                >  1:168514/1‑140 (MQ=255)                                                        gTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCgg                                                         >  1:395806/1‑140 (MQ=255)                                                           ccTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTa                                                      >  2:410035/1‑140 (MQ=255)                                                            cTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTAc                                                     >  2:126214/1‑140 (MQ=255)                                                                            gTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTcg                                     >  1:651577/1‑140 (MQ=255)                                                                                  cgcgTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCttt                               >  1:669721/1‑140 (MQ=255)                                                                                  cgcgTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCttt                               >  2:694234/1‑140 (MQ=255)                                                                                                       aCGGTGTCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGTCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCt          <  2:438602/140‑1 (MQ=255)                                                                                                          gtgtCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTg       >  1:205695/1‑140 (MQ=255)                                                                                                             tCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGcac    >  2:2463/1‑140 (MQ=255)                                                                                                             tCGTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGcac    <  1:2463/140‑1 (MQ=255)                                                                                                               gTTGGTCTGGTTAGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAg  >  1:595946/1‑140 (MQ=255)                                                                                                                           |                                                                                                                                CTGGCGTAGATGGAAACGTTATCAAGCGTTGGCGTTAGCGTACCCAGATCGTAACGTGCCTGCGTTTCCAGCCCGGTATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCAACATTGCCGCTTTGCACAG  >  NC_000913/4514753‑4515002

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑