Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	173,121	A→G	T554A (ACC→GCC)	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,121	0	A	G	100.0%	155.1 / NA	46	T554A (ACC→GCC)	fhuB	iron(III) hydroxamate ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/23);  total (23/23)

GACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCG  >  NC_000913/172988‑173255                                                                                                                                      |                                                                                                                                       gACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGccc                                                                                                                                  <  1:858850/140‑1 (MQ=255)  aCCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCg                                                                                                                                 <  1:796833/140‑1 (MQ=255)  aCCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCg                                                                                                                                 <  1:560991/140‑1 (MQ=255)            ccGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGa                                                                                                                       >  1:697896/1‑140 (MQ=255)                        aCAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACg                                                                                                           <  1:472036/140‑1 (MQ=255)                              cgACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGAcg                                                                                                     >  2:472676/1‑140 (MQ=255)                               gACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGAcgc                                                                                                    >  1:550995/1‑140 (MQ=255)                                aCCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGAcgcg                                                                                                   >  1:687566/1‑140 (MQ=255)                                   gATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTCCCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAAt                                                                                                <  1:980392/140‑1 (MQ=255)                                    aTGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAAtt                                                                                               >  2:266465/1‑140 (MQ=255)                                            gTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTGCCGCTGGGTGGTGATGCCCCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCgctgct                                                                                       <  1:942580/140‑1 (MQ=255)                                             tCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCgctgctg                                                                                      <  1:630230/140‑1 (MQ=255)                                                      ccGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCgct                                                                                          <  1:250853/127‑1 (MQ=255)                                                      ccGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCgct                                                                                          >  2:250853/1‑127 (MQ=255)                                                      ccGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGtt                                                                                 <  2:1055981/136‑1 (MQ=255)                                                      ccGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGtt                                                                                 >  1:1055981/1‑136 (MQ=255)                                                                 aTGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGAc                                                                  >  1:1035776/1‑140 (MQ=255)                                                                 aTGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGAc                                                                  >  2:1001611/1‑140 (MQ=255)                                                                    gTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGc                                                               >  1:1010966/1‑140 (MQ=255)                                                                        ttttGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGAcc                                                           <  1:266465/140‑1 (MQ=255)                                                                              tGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGcgacg                                                     <  2:1010966/140‑1 (MQ=255)                                                                                gCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGAt                                                   <  1:477373/140‑1 (MQ=255)                                                                                gCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGAt                                                   <  1:471489/140‑1 (MQ=255)                                                                                 cGATTACACCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATg                                                  >  2:691219/1‑140 (MQ=255)                                                                                          cGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGGGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTATCGGCTTTCCTGACGGCGACCGCGACGCTGGCTAtttgg                                         >  2:527068/1‑137 (MQ=255)                                                                                            cTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGAcc                                       >  2:352308/1‑140 (MQ=255)                                                                                              gtgCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGt                                     >  1:633970/1‑140 (MQ=255)                                                                                                 ccgccgCTGGCTGACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGa                                  >  1:561587/1‑140 (MQ=255)                                                                                                      gctggctgACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGtttt                             >  2:566679/1‑140 (MQ=255)                                                                                                      gctggctgACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGtttt                             >  2:688159/1‑140 (MQ=255)                                                                                                       ctggctgACCATTTTCCCGCTGGGTGGTTAGGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTttg                            <  1:1194151/140‑1 (MQ=255)                                                                                                        tggctgACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTttgt                           <  1:352308/140‑1 (MQ=255)                                                                                                         ggctgACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGGCCGCGACGATGACTTTTGGCCCGTTCAGTTttgtt                          >  2:1010189/1‑140 (MQ=255)                                                                                                         ggctgACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTttgtt                          <  2:633970/140‑1 (MQ=255)                                                                                                         ggctgACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTttgtt                          >  1:415457/1‑140 (MQ=255)                                                                                                           ctgACCATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTgg                        <  2:561587/140‑1 (MQ=255)                                                                                                              aCCATTTTACCGGTGGGTGGTTATGCCGCCCGAACCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCGGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGGTCTGGACTATTGGACCGTTGAGTGTTTTTGGttt                     >  2:1194151/1‑140 (MQ=255)                                                                                                                cATTTTACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTaa                   <  2:415457/140‑1 (MQ=255)                                                                                                                    ttACCGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGc               <  1:1001611/140‑1 (MQ=255)                                                                                                                    tgaCCGCGGGGGGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGc               <  1:1010189/138‑1 (MQ=255)                                                                                                                      cccGCTGGGTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCAc             <  2:1082829/139‑1 (MQ=255)                                                                                                                            gggTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCata       <  2:391194/140‑1 (MQ=255)                                                                                                                            gggTGGTGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCata       <  1:566679/140‑1 (MQ=255)                                                                                                                             ggtggtGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATg                 <  2:145481/129‑1 (MQ=255)                                                                                                                             ggtggtGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATg                 >  1:145481/1‑129 (MQ=255)                                                                                                                             ggtggtGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCatat      >  1:423362/1‑140 (MQ=255)                                                                                                                             ggtggtGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCAGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTaa                   <  1:580645/127‑1 (MQ=255)                                                                                                                             ggtggtGATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCAGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTaa                   >  2:580645/1‑127 (MQ=255)                                                                                                                                 gtgATGCCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTgcg  >  2:860139/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                       GACCTGGATTTCCGGTTCGACCTACAACGCGACCGATGCGCAGGTCTGGCGCACCGGAATTGTGATGGTGATTTTGCTGGCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCG  >  NC_000913/172988‑173255

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