Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,514,423	T→C	D752G (GAC→GGC)	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,514,423	0	T	C	100.0%	122.9 / NA	35	D752G (GAC→GGC)	fecA	ferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/17);  total (18/17)

AGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTAC  >  NC_000913/4514286‑4514558                                                                                                                                          |                                                                                                                                        aGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCa                                                                                                                                       >  1:172111/1‑140 (MQ=255) aGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCCCGGCTGGACTGCATGGATGCCATTGTTGTTGTCGCCa                                                                                                                                       >  1:366595/1‑140 (MQ=255)      aCGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAaga                                                                                                                                  <  2:637997/140‑1 (MQ=255)             aaTGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATg                                                                                                                           >  2:623529/1‑140 (MQ=255)                 ccAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAagt                                                                                                                       >  2:718682/1‑140 (MQ=255)                             caTTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGa                                                                                                           <  1:623529/140‑1 (MQ=255)                                 ccAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGAt                                                                                                       <  1:795750/140‑1 (MQ=255)                                        aaaGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTc                                                                                                >  1:151288/1‑140 (MQ=255)                                               ggCAAGCCGGGCGTTAGTACATTAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGa                                                                                         <  2:366595/140‑1 (MQ=255)                                               ggCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGa                                                                                         <  2:172111/140‑1 (MQ=255)                                                  aaGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATg                                                                                      <  1:763005/140‑1 (MQ=255)                                                    gCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCTTAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGcc                                                                                    <  2:791619/140‑1 (MQ=255)                                                       gggCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAgg                                                                                 <  2:627925/140‑1 (MQ=255)                                                           gTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCa                                                                             <  1:175293/140‑1 (MQ=255)                                                              aaCACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGa                                                                                         <  1:184505/125‑1 (MQ=255)                                                              aaCACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGa                                                                                         >  2:184505/1‑125 (MQ=255)                                                              aaCACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGAt                                                                          >  1:512873/1‑140 (MQ=255)                                                                             cAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGc                                                           >  1:450426/1‑140 (MQ=255)                                                                               aCGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCa                                                         <  1:718682/140‑1 (MQ=255)                                                                               aCGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCa                                                         <  2:512873/140‑1 (MQ=255)                                                                                            ataCAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCa                                            >  1:414348/1‑140 (MQ=255)                                                                                                 gCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGt                                       >  1:447686/1‑140 (MQ=255)                                                                                                  cGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTg                                      >  2:703058/1‑140 (MQ=255)                                                                                                   gTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTgc                                     >  2:204289/1‑140 (MQ=255)                                                                                                   gTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTgc                                     >  1:821330/1‑140 (MQ=255)                                                                                                    tGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGt                                                     <  2:434395/123‑1 (MQ=255)                                                                                                    tGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGt                                                     >  1:434395/1‑123 (MQ=255)                                                                                                    tGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTgcg                                    <  2:39498/140‑1 (MQ=255)                                                                                                      cgcgGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGcc                                  >  1:106031/1‑140 (MQ=255)                                                                                                      cgcgGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGcc                                  >  1:16573/1‑140 (MQ=255)                                                                                                       gcgGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGccc                                 <  2:414348/140‑1 (MQ=255)                                                                                                          gCTGGCCTGCATAGACGCCTTTTTTCTTGCCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCag                              <  1:211630/140‑1 (MQ=255)                                                                                                                   cATAGATGCCTTTGTTGTTGTCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGACCAAAGTGATACGCCAAGCGTGCGCCCCAGAGCATGAAg                     >  2:211630/1‑140 (MQ=255)                                                                                                                               tgttgttgtCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCgg         >  1:31077/1‑140 (MQ=255)                                                                                                                                     tgtCGCCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTa   <  1:703058/140‑1 (MQ=255)                                                                                                                                      gtcgccATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTAc  >  2:504617/1‑140 (MQ=255)                                                                                                                                      gtcgccATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGAGCGCCCCCGAGAATGAAGCCGGGAATGCGGCCGGTAc  >  1:362889/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        AGAAAACGGATAAATGCCAACATAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGTGCGCCCCAGAGCATGAAGCCGGGAATGCGGCCGGTAC  >  NC_000913/4514286‑4514558

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