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%	170.5 / NA	37	coding (1172/2244 nt)	fhuA	ferrichrome outer membrane transporter/phage receptor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (19/18);  total (19/18)

CGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCG  >  NC_000913/168586‑168727                                                                      |                                                                         cGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCa                                                              <  1:33663/81‑1 (MQ=255)    aTAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCagag                                                           >  2:33710/1‑81 (MQ=255)    aTAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCagag                                                           >  2:33700/1‑81 (MQ=255)     tAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCagag                                                           >  2:33701/1‑80 (MQ=255)     tAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGc                                                          <  1:33664/81‑1 (MQ=255)           ccATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGttt                                                    >  1:33711/1‑81 (MQ=255)           ccATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGttt                                                    <  2:33666/81‑1 (MQ=255)           ccATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGttt                                                    <  2:33665/81‑1 (MQ=255)            cATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTg                                                   <  1:33667/81‑1 (MQ=255)                aTCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCa                                                <  2:33669/80‑1 (MQ=255)                aTCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCa                                                <  2:33668/80‑1 (MQ=255)                aTCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCAc                                               >  1:33715/1‑81 (MQ=255)                     gCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTg                                             <  2:33670/78‑1 (MQ=255)                     gCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGc                                           <  2:195187/80‑1 (MQ=255)                      cACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGa                                         <  1:33671/81‑1 (MQ=255)                       aCGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGa                                         >  1:33722/1‑80 (MQ=255)                        cGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGata                                       >  1:33712/1‑81 (MQ=255)                         gTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGatat                                      >  1:33708/1‑81 (MQ=255)                              tACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGAcc                                 <  1:33672/81‑1 (MQ=255)                              tACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGAc                                  >  1:33714/1‑80 (MQ=255)                                 gtcgtTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACcaca                              <  1:33673/81‑1 (MQ=255)                                   cgtTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACAcc                            >  2:33706/1‑81 (MQ=255)                                   cgtTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACAcc                            >  1:33705/1‑81 (MQ=255)                                    gtTGATGATGAGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACAccc                           <  2:33674/81‑1 (MQ=255)                                      tgatgatgaGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCctg                         >  1:33721/1‑81 (MQ=255)                                         tgatgaGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCctgctg                      <  2:33675/81‑1 (MQ=255)                                           atgaGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGAc                    >  1:33717/1‑81 (MQ=255)                                           atgaGAAGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGAc                    >  2:33718/1‑81 (MQ=255)                                                 aGCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGgtgt              <  2:33676/81‑1 (MQ=255)                                                  gCTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTc             >  1:33724/1‑81 (MQ=255)                                                   cTGCAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCg            <  1:33677/81‑1 (MQ=255)                                                      cAAAACTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACt         <  2:33678/81‑1 (MQ=255)                                                       aaaaCTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACtt        >  1:33716/1‑81 (MQ=255)                                                       aaaaCTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACtt        >  2:33727/1‑81 (MQ=255)                                                           cTTCTCCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTtatg    >  1:33720/1‑81 (MQ=255)                                                              ctcCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTtatgcg  >  2:33723/1‑80 (MQ=255)                                                              ctcCGTTG‑TACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTtatgc   <  1:33679/79‑1 (MQ=255)                                                                      |                                                                         CGGATAAAGGCCATTATCTGGCACGTAAATACGTCGTTGATGATGAGAAGCTGCAAAACTTCTCCGTTGATACCCAGTTGCAGAGCAAGTTTGCCACTGGCGATATCGACCACACCCTGCTGACCGGTGTCGACTTTATGCG  >  NC_000913/168586‑168727

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