Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,514,463	(T)5→4	coding (2215/2325 nt)	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,459	0	T	.	100.0%	142.9 / NA	34	coding (2219/2325 nt)	fecA	ferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (16/18);  total (16/18)

CAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAG  >  NC_000913/4514380‑4514532                                                                                |                                                                          cAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGt                                                                           <  2:864939/80‑1 (MQ=255)       gcgcGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGttttcaca                                                                    >  2:864988/1‑77 (MQ=255)        cgcgGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCacac                                                                  <  2:864940/81‑1 (MQ=255)           ggCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCg                                                                >  2:864971/1‑80 (MQ=255)                gCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGc                                                           <  2:864941/80‑1 (MQ=255)                  cTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAg                                                        >  1:864978/1‑81 (MQ=255)                   tGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAgg                                                       >  1:864986/1‑81 (MQ=255)                      aTAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGtt                                                     >  2:864975/1‑80 (MQ=255)                      aTAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTc                                                    <  1:864944/81‑1 (MQ=255)                      aTAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTACACACCGAATGCCAGGTTc                                                    <  1:864943/81‑1 (MQ=255)                       tAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCa                                                   <  2:864944/81‑1 (MQ=255)                        aGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGtt                                                     >  1:864990/1‑78 (MQ=255)                             ccTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCt                                              <  2:864945/80‑1 (MQ=255)                             ccTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCt                                              <  1:864947/80‑1 (MQ=255)                                 tgttgttgtCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCa                                          <  1:864948/80‑1 (MQ=255)                                   ttgttgtCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATc                                        <  1:864949/80‑1 (MQ=255)                                   ttgttgtCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCt                                       <  1:864950/81‑1 (MQ=255)                                    tgttgtCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCTCACCGAATGCCAGGTTCAGATCTGCCATCTg                                      >  2:864996/1‑81 (MQ=255)                                    tgttgtCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTg                                      >  2:864982/1‑81 (MQ=255)                                      ttgtCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCg                                    >  1:865005/1‑81 (MQ=255)                                        gtcgtcATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCgg                                   <  2:864950/80‑1 (MQ=255)                                         tcgtcATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCggg                                  <  2:864951/80‑1 (MQ=255)                                           gtcATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGcc                                >  2:864994/1‑80 (MQ=255)                                           gtcATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGcc                                <  2:864953/80‑1 (MQ=255)                                           gtcATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCa                               <  1:864953/81‑1 (MQ=255)                                                       gATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGc                    >  2:865001/1‑80 (MQ=255)                                                        aTGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCa                  >  2:864984/1‑81 (MQ=255)                                                         tGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCAc                 <  2:864954/81‑1 (MQ=255)                                                         tGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCAc                 <  1:864956/81‑1 (MQ=255)                                                           aaGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCAcgc               >  2:865007/1‑81 (MQ=255)                                                               agtCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTgc           >  2:864983/1‑81 (MQ=255)                                                                gtCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTgcg          >  2:865006/1‑81 (MQ=255)                                                                    tGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGccc       >  2:865012/1‑80 (MQ=255)                                                                      gTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCag    <  2:864956/81‑1 (MQ=255)                                                                      gTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCag    <  1:864958/81‑1 (MQ=255)                                                                         gAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCagag  >  1:865000/1‑80 (MQ=255)                                                                                |                                                                          CAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAG  >  NC_000913/4514380‑4514532

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: ‑