Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	168,921	Δ1 bp	coding (1438/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,920	0	G	.	100.0%	152.0 / NA	36	coding (1437/2244 nt)	fhuA	ferrichrome outer membrane transporter/phage receptor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (23/13);  total (23/13)

CCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGT  >  NC_000913/168841‑168995                                                                                |                                                                            ccGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTgt                                                                            >  2:36714/1‑80 (MQ=255)  cGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTgt                                                                            >  1:36712/1‑79 (MQ=255)      ccTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCAccc                                                                      >  1:36698/1‑81 (MQ=255)       cTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCt                                                                     <  2:36675/81‑1 (MQ=255)         tACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAg                                                                   >  2:36700/1‑81 (MQ=255)          aCCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAgg                                                                  >  2:36709/1‑81 (MQ=255)          aCCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAg                                                                   <  2:36676/80‑1 (MQ=255)           ccGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGc                                                                 >  1:36708/1‑81 (MQ=255)           ccGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGc                                                                 <  2:36677/81‑1 (MQ=255)             gCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCgg                                                               <  2:36678/81‑1 (MQ=255)               aTTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTc                                                             >  1:36706/1‑81 (MQ=255)               aTTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTc                                                             >  1:36718/1‑81 (MQ=255)                     aaTAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGtt                                                          <  2:36679/78‑1 (MQ=255)                     aaTAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTAt                                                        >  2:36713/1‑80 (MQ=255)                     aaTAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTAt                                                        >  1:36719/1‑80 (MQ=255)                        aaaCAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACt                                                    >  2:36707/1‑81 (MQ=255)                          aCAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTgg                                                  >  2:36722/1‑81 (MQ=255)                                 caaacGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATc                                           <  2:36680/81‑1 (MQ=255)                                   aacGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATc                                           >  2:36727/1‑79 (MQ=255)                                   aacGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCaa                                         <  2:36681/81‑1 (MQ=255)                                       ggCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAAt                                     <  1:36682/81‑1 (MQ=255)                                              aTGTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTaa                              >  1:36728/1‑81 (MQ=255)                                                gTTCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGAACAAGAATCTCTTAAcc                            >  2:36716/1‑81 (MQ=255)                                                  tCAGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACcgc                          >  1:36731/1‑81 (MQ=255)                                                    aGGATCAGGCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGGATCTCTTAACCGCGt                        <  2:36683/81‑1 (MQ=255)                                                           ggCGCAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCgg                  >  1:36717/1‑80 (MQ=255)                                                              gcAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGgacg              <  1:36684/81‑1 (MQ=255)                                                              gcAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGgac               >  1:36710/1‑80 (MQ=255)                                                               cAGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGACCAAGAATCTCTTAACCGCGTTGCCGGgacga             >  2:36720/1‑81 (MQ=255)                                                                aGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGgacgac            >  1:36738/1‑81 (MQ=255)                                                                aGTGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGgacgac            >  1:36715/1‑81 (MQ=255)                                                                  tGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCg          <  2:36686/81‑1 (MQ=255)                                                                  tGGGATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCg          <  2:36685/81‑1 (MQ=255)                                                                   gggATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGa         >  1:36740/1‑81 (MQ=255)                                                                   gggATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGa         >  1:36729/1‑81 (MQ=255)                                                                     gATAAAGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATa       <  2:36687/81‑1 (MQ=255)                                                                        aaaGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAAc    >  2:36737/1‑81 (MQ=255)                                                                          aGTGCT‑GTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGt  >  2:36735/1‑81 (MQ=255)                                                                                |                                                                            CCGGCCCTTACCGCATTCTGAATAAACAGAAACAAACGGGCGTTTATGTTCAGGATCAGGCGCAGTGGGATAAAGTGCTGGTCACCCTAGGCGGTCGTTATGACTGGGCAGATCAAGAATCTCTTAACCGCGTTGCCGGGACGACCGATAAACGT  >  NC_000913/168841‑168995

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