Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,511,130	G→A	R110C (CGT→TGT)	fecE ←	ferric citrate ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,511,130	0	G	A	100.0%	175.2 / NA	47	R110C (CGT→TGT)	fecE	ferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/19);  total (28/19)

CGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCG  >  NC_000913/4510983‑4511275                                                                                                                                                    |                                                                                                                                                  cGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACAccc                                                                                                                                                >  2:710068/1‑151 (MQ=255)     cGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACAccccc                                                                                                                                              >  2:693638/1‑149 (MQ=255)                 cAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGCCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGcc                                                                                                                                  >  2:451319/1‑149 (MQ=255)                    gACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGa                                                                                                                             <  2:242920/151‑1 (MQ=255)                       cATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGAtt                                                                                                                           >  1:554272/1‑150 (MQ=255)                                aaaTGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATAc                                                                                                                 >  2:583096/1‑151 (MQ=255)                                aaaTGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATAc                                                                                                                 >  1:577602/1‑151 (MQ=255)                                    gcgcgcTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGaa                                                                                                              <  1:23524/150‑1 (MQ=255)                                     cgcgcTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAAc                                                                                                            >  1:135370/1‑151 (MQ=255)                                            gcggtgACCCCCGGAAAGCTCCGTTTACCCACCGACCGCCAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCt                                                                                                       <  2:596474/149‑1 (MQ=255)                                                tgACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTgg                                                                                                  >  1:533680/1‑150 (MQ=255)                                                    cgccgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAGCTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACtgt                                                                                             >  2:632651/1‑151 (MQ=255)                                                       cgGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGAt                                                                                          >  1:822880/1‑151 (MQ=255)                                                            aGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCct                                                                                     <  2:197231/151‑1 (MQ=255)                                                                cGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCctct                                                                                   <  2:559255/149‑1 (MQ=255)                                                                            gcaCGGCAAGATGATTGATCCGGGGCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTg                                                                     >  2:466224/3‑151 (MQ=255)                                                                            gaACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTg                                                                     <  2:75390/151‑1 (MQ=255)                                                                                  cAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAgg                                                               <  1:643136/151‑1 (MQ=255)                                                                                   aaGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgc                                                              >  1:784867/1‑151 (MQ=255)                                                                                    aGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgca                                                             <  1:146196/151‑1 (MQ=255)                                                                                      atgatTGATCCGGTTCCATTTCAGGGCGACATTAACTCGTGCATTGTCTTCAACGGACAGACACCCCAAGGGTGACAGCCAGGTATTAGGACCATACGAAACCAGCTCCTGAACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGgcag                                                            <  2:162645/150‑1 (MQ=255)                                                                                           tgatCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCg                                                       >  2:265667/1‑150 (MQ=255)                                                                                           tgatCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCg                                                       >  1:339976/1‑150 (MQ=255)                                                                                             atCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGaa                                                     >  1:166619/1‑150 (MQ=255)                                                                                                  ggTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTg                                               <  1:394512/151‑1 (MQ=255)                                                                                                     cTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCgg                                            >  1:90876/1‑151 (MQ=255)                                                                                                      tGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCggg                                           >  1:7455/1‑151 (MQ=255)                                                                                                       ggTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCgg                                            <  2:64901/149‑1 (MQ=255)                                                                                                       ggTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCg                                             <  1:560626/148‑1 (MQ=255)                                                                                                        gTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGcc                                         >  1:425430/1‑151 (MQ=255)                                                                                                         ttCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGcc                                         >  1:192244/1‑150 (MQ=255)                                                                                                          tCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCa                                        <  2:45053/150‑1 (MQ=255)                                                                                                              ggCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTgg                                   >  2:181004/1‑151 (MQ=255)                                                                                                              ggCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTg                                    <  2:710291/150‑1 (MQ=255)                                                                                                                       aaCTCGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATg                            >  2:127615/1‑149 (MQ=255)                                                                                                                           cGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGagag                        <  2:135805/149‑1 (MQ=255)                                                                                                                           cGTGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGagag                        <  1:280648/149‑1 (MQ=255)                                                                                                                             tGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCat                     >  2:681840/1‑150 (MQ=255)                                                                                                                             tGCATTGTCTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCat                     >  1:748450/1‑150 (MQ=255)                                                                                                                                     cTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATggg            >  2:535785/1‑151 (MQ=255)                                                                                                                                     cTTCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATggg            >  2:760659/1‑151 (MQ=255)                                                                                                                                       tCAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGa           <  2:7433/150‑1 (MQ=255)                                                                                                                                        cAGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGAtt         >  1:896469/1‑151 (MQ=255)                                                                                                                                         aGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGAtt         >  1:67633/1‑150 (MQ=255)                                                                                                                                         aGCGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTa        <  1:812595/151‑1 (MQ=255)                                                                                                                                           cGGAGAGACACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATc      <  2:242985/151‑1 (MQ=255)                                                                                                                                                agaCACCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCg  >  2:468146/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CGGGCGTATTCTGGGCCAGGACCATTGCCAGAAATGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACTGTGATCCCCTCTGGCGTTAAATGGTGCTGAGGCAGCAGCGAAAGCCTGCGGGCCAACTGGCGCGATGAGAGCATATTTATGGGATTATCGCCG  >  NC_000913/4510983‑4511275

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