Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,513,844	T→C	A155A (GCA→GCG)	fecB ←	ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,513,844	0	T	C	100.0%	175.5 / NA	51	A155A (GCA→GCG)	fecB	ferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/29);  total (22/29)

CCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCT  >  NC_000913/4513699‑4513990                                                                                                                                                  |                                                                                                                                                   ccAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTg                                                                                                                                               <  2:646139/151‑1 (MQ=255) ccAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTg                                                                                                                                               <  2:488430/151‑1 (MQ=255)  cAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGc                                                                                                                                              <  2:706928/151‑1 (MQ=255)                  cTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTAc                                                                                                                              >  1:424918/1‑151 (MQ=255)                   tGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTAcc                                                                                                                             >  1:259011/1‑151 (MQ=255)                    gCCGGTCCAGGTCTCCTGAGTATCCAGGTTGAATTGCTGTTGGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAGGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTAc                                                                                                                              <  1:49723/149‑1 (MQ=255)                       ggTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACccac                                                                                                                          <  2:609502/150‑1 (MQ=255)                         tCCAGTTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACccacca                                                                                                                        <  1:296896/150‑1 (MQ=255)                          ccAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTCATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCAtt                                                                                                                      >  1:530723/1‑151 (MQ=255)                             ggTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCAtctc                                                                                                                    <  1:273374/150‑1 (MQ=255)                                 tcCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCg                                                                                                                >  2:819091/1‑150 (MQ=255)                                  ccTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgat                                                                                                              <  1:847673/151‑1 (MQ=255)                                     gAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatgat                                                                                                           <  2:347494/151‑1 (MQ=255)                                      aGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatgat                                                                                                           >  2:407979/1‑150 (MQ=255)                                      aGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATa                                                                                                          <  1:417446/151‑1 (MQ=255)                                       gTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAg                                                                                                         <  1:485815/151‑1 (MQ=255)                                        tATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGc                                                                                                        >  2:354839/1‑151 (MQ=255)                                              ggTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGc                                                                                                     >  2:339147/1‑148 (MQ=255)                                                  gAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTg                                                                                               <  2:712907/150‑1 (MQ=255)                                                     ttGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCaa                                                                                            >  2:320109/1‑150 (MQ=255)                                                     ttGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCaa                                                                                            >  2:209440/1‑150 (MQ=255)                                                         tGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAAtttt                                                                                       >  2:219538/1‑151 (MQ=255)                                                          gTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAAtttt                                                                                       <  2:116602/150‑1 (MQ=255)                                                             cgcgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGc                                                                                   >  1:697/1‑151 (MQ=255)                                                                  gATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAg                                                                               <  2:834222/150‑1 (MQ=255)                                                                     gtgCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGttt                                                                           >  2:209707/1‑151 (MQ=255)                                                                         cAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGtt                                                                       <  2:459714/151‑1 (MQ=255)                                                                          aaaGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGtt                                                                       <  1:374958/150‑1 (MQ=255)                                                                          aaaGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTg                                                                      <  2:539457/151‑1 (MQ=255)                                                                                 aCGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACt                                                               <  1:321112/151‑1 (MQ=255)                                                                                       gtCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAagc                                                          <  2:603965/150‑1 (MQ=255)                                                                                                   aGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGgcg                                             >  2:440938/1‑151 (MQ=255)                                                                                                          tggcCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTtg                                      >  2:316083/1‑151 (MQ=255)                                                                                                           ggcCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTtg                                      <  1:316059/150‑1 (MQ=255)                                                                                                            gcCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTtgct                                    <  1:690454/151‑1 (MQ=255)                                                                                                             cccACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTtgctg                                   <  2:127085/151‑1 (MQ=255)                                                                                                                   gcgcCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAgg                              <  1:780613/150‑1 (MQ=255)                                                                                                                   gcgcCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAgg                              <  2:142249/150‑1 (MQ=255)                                                                                                                    cgcCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCg                            <  2:123873/151‑1 (MQ=255)                                                                                                                     gcCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCg                            >  2:454996/1‑150 (MQ=255)                                                                                                                     gcCATCCTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCg                            >  2:701118/1‑150 (MQ=255)                                                                                                                          ccTCTCTTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTa                       <  1:613773/150‑1 (MQ=255)                                                                                                                             ctctTTATGTTGTTCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAAcc                   >  2:796961/1‑151 (MQ=255)                                                                                                                                tttATGTTGTTCCAGACGCGCCTGCGTCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAAccccc                >  1:414735/1‑151 (MQ=255)                                                                                                                                     gttgttCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATg           >  1:262237/1‑151 (MQ=255)                                                                                                                                        gttCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCg        <  2:419772/151‑1 (MQ=255)                                                                                                                                          tCCAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGAc      >  1:479648/1‑151 (MQ=255)                                                                                                                                           ccAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGAct     >  1:304759/1‑151 (MQ=255)                                                                                                                                           ccAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGAct     <  2:821104/151‑1 (MQ=255)                                                                                                                                           ccAGACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGAct     >  2:825937/1‑151 (MQ=255)                                                                                                                                              gACGCGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGActgct  <  1:454960/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCT  >  NC_000913/4513699‑4513990

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