Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	172,269	(G)5→6	coding (808/1983 nt)	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	172,264	1	.	G	100.0%	116.2 / NA	34	L268R (CTG→CGG)	fhuB	iron(III) hydroxamate ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (16/18);  total (16/18)

AATCTCGGGCTGGCCTTGTCGCTTGCGCGTCTGGCAGCGCTGTCGCTGGCGATTGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCT‑GGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCA  >  NC_000913/172131‑172392                                                                                                                                       |                                                                                                                                 aaTCTCGGGCTGGCCTTGTCGCTTGCGCGTCTGGCAGCGCTGTCGCTGGCGATTGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCT‑gggggg                                                                                                                            >  2:81730/1‑139 (MQ=255)       gggCTGGCCTTGTCGCTTGCGCGTCTGGCAGCGCTGTCGCTGGCGATTGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCgcggc                                                                                                                       <  2:323072/140‑1 (MQ=255)             gCCTTGTCGCTTGCGCGTCTGGCAGCGCTGTCGCTGGCGATTGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTctgc                                                                                                                 >  1:481340/1‑140 (MQ=255)                  gTCGCTTGCGCGTCTGGCAGCGCTGTCGCTGGCGATTGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGCCTTTATCGGGTTGTTCGCTCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCa                                                                                                            <  2:649838/140‑1 (MQ=255)                      cTTGCGCGTCTGGCAGCGCTGTCGCTGGCGATTGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGAc                                                                                                        <  1:369443/140‑1 (MQ=255)                        tGCGCGTCTGGCAGCGCTGTCGCTGGCGATTGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTg                                                                                                      >  1:467156/1‑140 (MQ=255)                            cgTCTGGCAGCGCTGTCGCTGGCGATTGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGc                                                                                                  >  1:446287/1‑140 (MQ=255)                               cTGGCAGCGCTGTCGCTGGCGATTGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTgg                                                                                               <  2:122061/140‑1 (MQ=255)                                 ggCAGCGCTGTCGCTGGCGATTGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGcg                                                                                             <  1:714198/140‑1 (MQ=255)                                         tGTCGCTGGCGATTGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGcg                                                                                             <  2:369925/132‑1 (MQ=255)                                         tGTCGCTGGCGATTGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGcg                                                                                             >  1:369925/1‑132 (MQ=255)                                                  cGATTGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCt                                                                            <  1:234660/140‑1 (MQ=255)                                                      tGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATc                                                                        <  2:446287/140‑1 (MQ=255)                                                        tcatcaGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCct                                                                      >  2:173747/1‑140 (MQ=255)                                                               tGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCtt                                                               >  2:548528/1‑140 (MQ=255)                                                                 cgcTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGC‑AAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTcc                                                            >  2:202828/1‑140 (MQ=255)                                                                   ctgctgGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGC‑AAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGa                                                          <  1:202828/140‑1 (MQ=255)                                                                    tgctgGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGa                                                          >  1:84891/1‑140 (MQ=255)                                                                                cTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAatcatc                                                >  1:128090/1‑138 (MQ=255)                                                                                cTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAatcatc                                                <  2:128090/138‑1 (MQ=255)                                                                                cTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCct                                              >  2:712436/1‑140 (MQ=255)                                                                                cTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTCCCACGACTGATGCCGGCGTCGTTGATTGGTGCGCTGCTCCTCTGGCTTTCCGATCAAATCATCct                                              <  2:481340/140‑1 (MQ=255)                                                                                 tgcggggATTATCGGCTTTATCGGGTTGTTCGCCCCGCTGCTGTCAAAAATTCTGGGGGGCGCCGCGTCTGCGGCGACGACTGATGGTGGCGTCGTTGATTGGGGCGCTGATCCTCTGGCTTTCCGATCAACTCATCctc                                             >  2:455649/4‑140 (MQ=255)                                                                                   tgGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTg                                           >  2:439625/1‑140 (MQ=255)                                                                                         ttATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGAc                                     >  2:396534/1‑140 (MQ=255)                                                                                           aTCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTc                                   >  1:136938/1‑140 (MQ=255)                                                                                                  ttATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGtgtg                            >  2:622843/1‑140 (MQ=255)                                                                                                         gttgttCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGaa                     <  1:712436/140‑1 (MQ=255)                                                                                                             ttCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGaa                     <  2:348229/136‑1 (MQ=255)                                                                                                             ttCGCGCCGCTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGaa                     >  1:348229/1‑136 (MQ=255)                                                                                                                   ccgcTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCg           >  2:144354/1‑140 (MQ=255)                                                                                                                    cgcTGCTGGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCgg          <  1:622843/140‑1 (MQ=255)                                                                                                                       tgctgGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTc       <  2:467156/140‑1 (MQ=255)                                                                                                                          tgGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGt    <  2:136938/140‑1 (MQ=255)                                                                                                                          tgGCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGt    <  1:171292/140‑1 (MQ=255)                                                                                                                           ggCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTc   <  2:20370/140‑1 (MQ=255)                                                                                                                            gCAAAAATGCTGGGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCa  <  1:439625/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                 AATCTCGGGCTGGCCTTGTCGCTTGCGCGTCTGGCAGCGCTGTCGCTGGCGATTGTCATCAGTGCGCTGCTGGTGAACGCTGTGGGGATTATCGGCTTTATCGGGTTGTTCGCGCCGCTGCTGGCAAAAATGCT‑GGGGGCGCGGCGTCTGCTGCCACGACTGATGCTGGCGTCGTTGATTGGTGCGCTGATCCTCTGGCTTTCCGATCAAATCATCCTCTGGCTGACTCGCGTGTGGATGGAAGTGTCCACCGGTTCGGTCA  >  NC_000913/172131‑172392

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