Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	167,915	Δ1 bp	coding (432/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,915	0	G	.	100.0%	149.4 / NA	33	coding (432/2244 nt)	fhuA	ferrichrome outer membrane transporter/phage receptor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (18/15);  total (18/15)

ACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTT  >  NC_000913/167845‑167990                                                                       |                                                                            accaccTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACt                                                                   <  1:35936/80‑1 (MQ=255) accaccTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACt                                                                   <  1:35935/80‑1 (MQ=255)    accTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCta                                                               <  1:35937/81‑1 (MQ=255)        gATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAAc                                                           >  2:35983/1‑81 (MQ=255)        gATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAAc                                                           >  1:35979/1‑81 (MQ=255)         aTCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATaa                                                            <  1:35938/79‑1 (MQ=255)         aTCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATaa                                                            <  2:35939/79‑1 (MQ=255)              tCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCg                                                     >  2:36069/1‑81 (MQ=255)                             aGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGtat                                      <  1:35940/81‑1 (MQ=255)                               aaGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTatat                                    >  1:35989/1‑81 (MQ=255)                               aaGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTatat                                    >  2:35977/1‑81 (MQ=255)                                  gCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCt                                 >  2:35982/1‑81 (MQ=255)                                  gCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCt                                 <  2:35941/81‑1 (MQ=255)                                    cAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTgg                               <  2:35942/81‑1 (MQ=255)                                    cAAAGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTg                                >  2:36005/1‑80 (MQ=255)                                       aGCCAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAAc                            <  1:35943/81‑1 (MQ=255)                                         ccAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAAcgc                          >  1:35998/1‑81 (MQ=255)                                          cAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAAcgcg                         <  2:35944/81‑1 (MQ=255)                                          cAGAATAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAAcgcg                         >  2:35991/1‑81 (MQ=255)                                              aTAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTg                      >  1:35981/1‑80 (MQ=255)                                              aTAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGa                     <  1:35946/81‑1 (MQ=255)                                              aTAACTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGa                     <  2:35945/81‑1 (MQ=255)                                                aaCTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGaaa                   >  1:35993/1‑81 (MQ=255)                                                 aCTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAAt                  <  2:35947/81‑1 (MQ=255)                                                 aCTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAAt                  >  1:35994/1‑81 (MQ=255)                                                 aCTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAAt                  >  2:35978/1‑81 (MQ=255)                                                 aCTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAAt                  >  1:35987/1‑81 (MQ=255)                                                  cTATCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAAtt                 <  2:35948/81‑1 (MQ=255)                                                     tCTGAATGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTAt               <  1:35949/80‑1 (MQ=255)                                                         aaTGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTg          >  2:35990/1‑81 (MQ=255)                                                          aTGGCCTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTgg         >  1:35986/1‑81 (MQ=255)                                                              ccTGAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCg     >  2:35999/1‑81 (MQ=255)                                                                 gAAGTT‑CAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGtt  >  1:35988/1‑81 (MQ=255)                                                                       |                                                                            ACCACCTGATCATTCGCGGCTTTGCGGCAGAAGGCCAAAGCCAGAATAACTATCTGAATGGCCTGAAGTTGCAGGGCAACTTCTATAACGATGCGGTCATTGACCCGTATATGCTGGAACGCGCTGAAATTATGCGTGGCCCGGTT  >  NC_000913/167845‑167990

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