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%	145.6 / NA	35	coding (2219/2325 nt)	fecA	ferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (16/19);  total (16/19)

AGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCAT  >  NC_000913/4514381‑4514535                                                                               |                                                                             aGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGttt                                                                            >  1:935230/1‑81 (MQ=255)        gcgGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACAcc                                                                    <  1:935188/81‑1 (MQ=255)        gcgGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACAcc                                                                    <  2:935189/81‑1 (MQ=255)          ggCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGATGTAGTCCTGGTCGAAGATG‑TTTTCACACCg                                                                   >  2:935237/1‑80 (MQ=255)              ggCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGc                                                              >  2:935238/1‑81 (MQ=255)               gCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGcc                                                             <  1:935191/81‑1 (MQ=255)               gCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGcc                                                             <  2:935192/81‑1 (MQ=255)               gCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGcc                                                             <  2:935190/81‑1 (MQ=255)               gCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGc                                                              >  1:935241/1‑80 (MQ=255)                ccTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCa                                                            <  2:935193/81‑1 (MQ=255)                    cATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGt                                                         <  1:935194/80‑1 (MQ=255)                     aTAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTc                                                       <  2:935195/81‑1 (MQ=255)                                 gttgttgtCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATc                                           >  2:935231/1‑81 (MQ=255)                                  ttgttgtCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCt                                          <  2:935196/81‑1 (MQ=255)                                      tgtCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCgg                                      >  1:935236/1‑81 (MQ=255)                                       gtcgtcATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCgg                                      <  1:935197/80‑1 (MQ=255)                                        tcgtcATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCggg                                     <  2:935198/80‑1 (MQ=255)                                          gtcATAAGGGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCa                                  <  1:935199/81‑1 (MQ=255)                                          gtcATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCa                                  <  2:935200/81‑1 (MQ=255)                                            cATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCa                                  >  2:935253/1‑79 (MQ=255)                                             aTAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCATGTTCAGATCTGCCATCTGCGGGCCAAAg                               >  1:935247/1‑81 (MQ=255)                                             aTAAGAGCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAg                               <  1:935201/81‑1 (MQ=255)                                                  agCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAACGTCATa                          >  2:935242/1‑81 (MQ=255)                                                   gCGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATAc                         <  2:935202/81‑1 (MQ=255)                                                    cGGATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACg                        <  1:935203/81‑1 (MQ=255)                                                     ggATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGc                       >  1:935266/1‑81 (MQ=255)                                                     ggATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGc                       >  2:935239/1‑81 (MQ=255)                                                      gATGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGc                       >  1:935244/1‑80 (MQ=255)                                                        tGAAGTAGTCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCAc                    >  1:935240/1‑81 (MQ=255)                                                           agtagtCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCAcgcg                 >  2:935255/1‑81 (MQ=255)                                                            gtagtCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGt                <  2:935204/81‑1 (MQ=255)                                                              agtCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTgc              >  2:935252/1‑81 (MQ=255)                                                               gtCCTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTgc              <  2:935205/80‑1 (MQ=255)                                                                  cTGGTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGccc          >  1:935250/1‑81 (MQ=255)                                                                     gTCGAAGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGccc          >  1:935249/1‑78 (MQ=255)                                                                          aGATG‑TTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCAt  <  1:935206/81‑1 (MQ=255)                                                                               |                                                                             AGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCAT  >  NC_000913/4514381‑4514535

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