Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	173,198	C→T	T579T (ACC→ACT)	fhuB →	iron(III) hydroxamate ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	173,198	0	C	T	100.0%	137.5 / NA	39	T579T (ACC→ACT)	fhuB	iron(III) hydroxamate ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/18);  total (21/18)

ACCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACG  >  NC_000913/173121‑173275                                                                              |                                                                              aCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGTTGTTAGCGGCTTGCCTGACGGCGACtgc                                                                             >  1:37753/1‑80 (MQ=255)  ccGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGcga                                                                           >  1:37741/1‑81 (MQ=255)   cGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGcgac                                                                          >  2:37721/1‑81 (MQ=255)   cGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGcgac                                                                          <  1:37685/81‑1 (MQ=255)    gCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGcgacg                                                                         <  1:37686/81‑1 (MQ=255)      ccGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGAt                                                                       >  1:37733/1‑81 (MQ=255)       cGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATg                                                                      <  2:37687/81‑1 (MQ=255)        gAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTACTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGa                                                                     >  2:223422/1‑81 (MQ=255)           ccGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTa                                                                  <  2:37688/81‑1 (MQ=255)           ccGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTa                                                                  <  2:37689/81‑1 (MQ=255)                  aaTGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGAcc                                                           >  1:37734/1‑81 (MQ=255)                          tGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGtt                                                   >  1:37724/1‑81 (MQ=255)                          tGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGtt                                                   >  2:37737/1‑81 (MQ=255)                            aCGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGtttt                                                 >  2:37740/1‑81 (MQ=255)                             cGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTttg                                                >  1:37727/1‑81 (MQ=255)                                  aCGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGt                                           <  1:37690/81‑1 (MQ=255)                                    gcgAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGtt                                          >  2:37707/1‑80 (MQ=255)                                    gcgAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGtt                                          >  2:37755/1‑80 (MQ=255)                                       aaTTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAAt                                      <  2:37691/81‑1 (MQ=255)                                        aTTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATg                                     <  1:37692/81‑1 (MQ=255)                                         ttGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATgg                                    <  2:37693/81‑1 (MQ=255)                                           gcgcTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGc                                   <  1:37694/80‑1 (MQ=255)                                           gcgcTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCa                                  >  2:37746/1‑81 (MQ=255)                                            cgcTGCTGCTGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCAc                                 >  2:37739/1‑81 (MQ=255)                                             gctgctgctGTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCAcc                                <  1:37695/81‑1 (MQ=255)                                                 ctgctgTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCa                             >  1:37742/1‑80 (MQ=255)                                                 ctgctgTTAGCGGCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCa                             <  1:37696/80‑1 (MQ=255)                                                            ggCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGgat                 >  2:37758/1‑81 (MQ=255)                                                             gCTTGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGgatg                <  2:37697/81‑1 (MQ=255)                                                               ttGCCTGACGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGgatgat              >  1:37749/1‑81 (MQ=255)                                                                      aCGGCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCtt        <  2:37698/80‑1 (MQ=255)                                                                        ggCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTc      >  1:37757/1‑80 (MQ=255)                                                                        ggCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTc      >  2:37730/1‑80 (MQ=255)                                                                        ggCGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCg     <  2:37699/81‑1 (MQ=255)                                                                          cGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATTGCTCCTCCTATTTCGCTGTTTATGGTCTTTCgac   >  1:37752/1‑81 (MQ=255)                                                                          cGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCTCCGCATATTTCGCGGATGATGGGCTTTCgac   >  2:37772/1‑81 (MQ=255)                                                                          cGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCgac   >  2:37759/1‑81 (MQ=255)                                                                          cGACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCgac   <  1:37700/81‑1 (MQ=255)                                                                           gACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCgacg  <  2:37701/81‑1 (MQ=255)                                                                           gACTGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCgacg  <  1:37702/81‑1 (MQ=255)                                                                              |                                                                              ACCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACG  >  NC_000913/173121‑173275

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