Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	167,651	(G)5→4	coding (168/2244 nt)	fhuA →	ferrichrome outer membrane transporter/phage receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	167,647	0	G	.	100.0%	127.9 / NA	30	coding (164/2244 nt)	fhuA	ferrichrome outer membrane transporter/phage receptor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (15/15);  total (15/15)

TTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGC  >  NC_000913/167569‑167721                                                                               |                                                                           tttATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATg                                                                            <  1:38685/79‑1 (MQ=255)   tATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATggggc                                                                        >  1:38723/1‑80 (MQ=255)   tATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATgggg                                                                         <  2:38686/80‑1 (MQ=255)       cacaGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCctgc                                                                   <  2:38687/81‑1 (MQ=255)              aGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACt                                                            <  2:38688/81‑1 (MQ=255)                 ggTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTAtt                                                         <  1:38689/81‑1 (MQ=255)                 ggTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTAtt                                                         >  1:38742/1‑81 (MQ=255)                 ggTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTAtt                                                         >  1:38734/1‑81 (MQ=255)                 ggTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTAtt                                                         >  1:38732/1‑81 (MQ=255)                     gAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTgcgg                                                     >  1:38746/1‑81 (MQ=255)                      aaCCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCGAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTgcggc                                                    <  2:38690/81‑1 (MQ=255)                      aaCCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTgcggc                                                    >  1:38724/1‑81 (MQ=255)                       aCCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTgcggcg                                                   >  1:38739/1‑81 (MQ=255)                          gAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGa                                                >  2:38737/1‑81 (MQ=255)                           aaaGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGAc                                               <  2:38691/81‑1 (MQ=255)                            aagaagACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACa                                              <  1:38692/81‑1 (MQ=255)                             agaagaCACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAg                                             <  1:38693/81‑1 (MQ=255)                              gaagaCACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGt                                            >  2:38726/1‑81 (MQ=255)                                        tCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCg                                   <  1:38694/80‑1 (MQ=255)                                          aCCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGc                                 <  1:38695/80‑1 (MQ=255)                                             gTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACt                              >  2:38729/1‑80 (MQ=255)                                               tACCGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTaa                            >  1:38743/1‑80 (MQ=255)                                                  cGCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAAcc                        >  1:38740/1‑81 (MQ=255)                                                   gCTGCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAAcc                        >  1:38747/1‑80 (MQ=255)                                                      gCACCTGCGCCGCAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATa                    >  2:38727/1‑81 (MQ=255)                                                             cgccgcAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCg               <  2:38696/79‑1 (MQ=255)                                                             cgccgcAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGAt             <  2:38697/81‑1 (MQ=255)                                                              gccgcAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGAtt            <  2:38698/81‑1 (MQ=255)                                                              gccgcAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGAtt            >  2:38763/1‑81 (MQ=255)                                                               ccgcAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGAtt            >  1:38733/1‑80 (MQ=255)                                                                cgcAAGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCa          <  1:38699/81‑1 (MQ=255)                                                                    aGAAAGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCaaaaa      <  1:38700/81‑1 (MQ=255)                                                                        aGCGCAT‑GGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGc  <  1:807443/81‑1 (MQ=255)                                                                               |                                                                           TTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGGGGGCCTGCTGCAACTATTGCGGCGCGACAGTCTGCTACCGGCACTAAAACCGATACGCCGATTCAAAAAGTGC  >  NC_000913/167569‑167721

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