Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	168,655	Δ1 bp	coding (1172/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	168,655	0	A	.	100.0%	124.4 / NA	27	coding (1172/2244 nt)	fhuA	ferrichrome outer membrane transporter/phage receptor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (11/16);  total (11/16)

GCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGC  >  NC_000913/168581‑168726                                                                           |                                                                        gcCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAg                                                                  >  2:37395/1‑81 (MQ=255)   cAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGtt                                                                <  2:37348/81‑1 (MQ=255)       ggATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCag                                                            <  1:37349/81‑1 (MQ=255)            aaGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGc                                                         >  1:37387/1‑79 (MQ=255)              ggCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGt                                                     <  2:37350/81‑1 (MQ=255)                 cATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTg                                                  <  2:222615/81‑1 (MQ=255)                       cTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTg                                            >  1:37402/1‑81 (MQ=255)                          gCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCg                                         <  2:37351/81‑1 (MQ=255)                           cACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGa                                        >  2:37404/1‑81 (MQ=255)                           cACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGa                                        >  1:37398/1‑81 (MQ=255)                            aCGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGat                                       <  1:37352/81‑1 (MQ=255)                             cGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGata                                      <  1:37353/81‑1 (MQ=255)                              gTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGatat                                     >  2:37399/1‑81 (MQ=255)                                   tACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGAcc                                <  2:37354/81‑1 (MQ=255)                                     cgtcgtTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACcac                              <  2:37355/81‑1 (MQ=255)                                      gtcgtTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACcaca                             <  2:37356/81‑1 (MQ=255)                                        cgtTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACAcc                           >  2:37407/1‑81 (MQ=255)                                         gtTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACAcc                           <  1:37357/80‑1 (MQ=255)                                          ttGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCct                         <  2:37358/81‑1 (MQ=255)                                                atgaGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGAc                   >  2:37416/1‑81 (MQ=255)                                                     aaGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGgtg              >  2:37401/1‑81 (MQ=255)                                                         tGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGa          >  1:37408/1‑81 (MQ=255)                                                           cAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGa          <  2:37359/79‑1 (MQ=255)                                                           cAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACt        >  2:37414/1‑81 (MQ=255)                                                             aaaCTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACttt      <  2:37360/81‑1 (MQ=255)                                                               aCTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTtat    <  1:37361/81‑1 (MQ=255)                                                                   ctcCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTtatgc  <  1:37362/79‑1 (MQ=255)                                                                           |                                                                        GCCAGCGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGC  >  NC_000913/168581‑168726

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