Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,516,058	T→C	Y207C (TAC→TGC)	fecA ←	ferric citrate outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,516,058	0	T	C	100.0%	161.9 / NA	45	Y207C (TAC→TGC)	fecA	ferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/21);  total (24/21)

TGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCT  >  NC_000913/4515924‑4516188                                                                                                                                       |                                                                                                                                   tGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAggg                                                                                                                               >  2:192169/1‑140 (MQ=255) tGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAggg                                                                                                                               >  2:712955/1‑140 (MQ=255)    cACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCg                                                                                                                            <  2:589695/140‑1 (MQ=255)       gggTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAg                                                                                                                         <  2:116246/140‑1 (MQ=255)         gTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAggg                                                                                                                       <  1:889636/140‑1 (MQ=255)           aaCAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGa                                                                                                                     <  1:149715/140‑1 (MQ=255)              aaaGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACgg                                                                                                                  <  2:584751/140‑1 (MQ=255)                    caccacGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGc                                                                                                            <  2:818727/140‑1 (MQ=255)                          gccgccCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTcc                                                                                                       <  1:72944/139‑1 (MQ=255)                          gccgccCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTcc                                                                                                       >  2:72944/1‑139 (MQ=255)                                cACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGa                                                                                                <  1:370148/140‑1 (MQ=255)                                cACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGa                                                                                                <  1:192169/140‑1 (MQ=255)                                      ctGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCg                                                                                          >  1:194318/1‑140 (MQ=255)                                         cGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGc                                                                                       >  1:258968/1‑140 (MQ=255)                                              cGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTg                                                                                  >  2:859804/1‑140 (MQ=255)                                              cGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTg                                                                                  >  1:816625/1‑140 (MQ=255)                                                   cgcaccgcacCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAg                                                                             >  1:178605/1‑140 (MQ=255)                                                      accgcaccACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGgcg                                                                          <  2:105318/140‑1 (MQ=255)                                                      accgcaccACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGgcg                                                                          >  1:907583/1‑140 (MQ=255)                                                       ccgcaccACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGgcgc                                                                         <  2:584848/140‑1 (MQ=255)                                                        cgcaccACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGcgcg                                                                        >  2:69012/1‑140 (MQ=255)                                                        cgcaccACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGcgcg                                                                        >  2:450631/1‑140 (MQ=255)                                                         gcaccACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCgg                                                                       <  2:360963/140‑1 (MQ=255)                                                            ccaccaccGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGtt                                                                    <  2:194318/140‑1 (MQ=255)                                                               ccaccGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAg                                                                 >  1:458668/1‑140 (MQ=255)                                                                  ccGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGcc                                                              >  2:757476/1‑140 (MQ=255)                                                                    gcgTACCACGTCAATGGCATCCATGTTACCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGcccc                                                            <  1:885078/140‑1 (MQ=255)                                                                     cgTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCg                                                           >  2:893887/1‑140 (MQ=255)                                                                     cgTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGAGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCg                                                           >  2:34004/1‑140 (MQ=255)                                                                       tACCACTTCAATGGCATCCATGTCGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGa                                                         >  1:672730/1‑140 (MQ=255)                                                                       tACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGa                                                         <  2:286005/140‑1 (MQ=255)                                                                                aaTGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGt                                                <  2:515601/140‑1 (MQ=255)                                                                                         cATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCa                                       <  2:258968/140‑1 (MQ=255)                                                                                               gCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGt                                 >  2:43429/1‑140 (MQ=255)                                                                                                   aGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCt                             >  2:417538/1‑140 (MQ=255)                                                                                                   aGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCt                             >  2:293899/1‑140 (MQ=255)                                                                                                     cGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGc                           <  2:448406/140‑1 (MQ=255)                                                                                                          cGGGAGCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTg                      <  2:590587/140‑1 (MQ=255)                                                                                                               gCCAGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCgtt                 >  1:690658/1‑140 (MQ=255)                                                                                                                  aGTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCgttgtt              >  2:914641/1‑140 (MQ=255)                                                                                                                   gTGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGttt             >  2:818736/1‑140 (MQ=255)                                                                                                                    tGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGtttt            >  1:757900/1‑140 (MQ=255)                                                                                                                    tGAAAGCTGCGGCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGtttt            >  1:192372/1‑140 (MQ=255)                                                                                                                              ggCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCt  <  2:458668/140‑1 (MQ=255)                                                                                                                              ggCTGACCGCAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCt  <  2:178605/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                   TGGCACGGGTAACAAAGTTCACCACGCCGCCCACGCTCTGCGGTCCGTAACGCACCGCACCACCACCGCGTACCACGTCAATGGCATCCATGTTGCCGAGCGAAACGGGAGCCAGTGAAAGCTGCGGCTGACCGTAAGGGGCGAAGGGGACGGGGATGCCGTCCATCAGGACGGTCGAGCGGCTGGCGAGGCGCGGGTTCAGGCCCCGGATGCCAAAGTTCATCGCCAGGTCGTGGCTGCCGGTGCCGTTGTTTTCCGGCGCGCT  >  NC_000913/4515924‑4516188

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: ‑