Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,517,263	T→C	Y152C (TAC→TGC)	fecR ←	anti‑sigma transmembrane signal transducer for ferric citrate transport; periplasmic FecA‑bound ferric citrate sensor and cytoplasmic FecI ECF sigma factor activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,517,263	0	T	C	100.0%	159.9 / NA	46	Y152C (TAC→TGC)	fecR	anti‑sigma transmembrane signal transducer for ferric citrate transport; periplasmic FecA‑bound ferric citrate sensor and cytoplasmic FecI ECF sigma factor activator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/24);  total (22/24)

TGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGC  >  NC_000913/4517119‑4517407                                                                                                                                                 |                                                                                                                                                 tgctgctgcACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCat                                                                                                                                             <  1:520734/150‑2 (MQ=255)   ctgctgcACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGc                                                                                                                                          >  2:732122/1‑151 (MQ=255)   ctgctgcACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGc                                                                                                                                          >  1:561775/1‑151 (MQ=255)       tgcACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCgg                                                                                                                                       >  1:356932/1‑150 (MQ=255)          aCGTCAAGCTGCGTGATATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACgg                                                                                                                                   <  2:727273/151‑1 (MQ=255)           cGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGt                                                                                                                                  <  1:496902/151‑1 (MQ=255)            gTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTg                                                                                                                                 >  1:654803/1‑151 (MQ=255)               aaGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTgcgc                                                                                                                              >  2:172842/1‑151 (MQ=255)                          aaTTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATc                                                                                                                   >  2:555701/1‑151 (MQ=255)                           aTTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCa                                                                                                                  <  1:457022/151‑1 (MQ=255)                              aTCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCaaa                                                                                                                >  2:252433/1‑150 (MQ=255)                                ccTGCTGGCGGACGGTAAATTCTGACCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACg                                                                                                             <  1:560013/151‑1 (MQ=255)                                   gctgGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTCGCCCTGACGGGTCAGGACACGGCAGGGGCGTTGCAGGGCATCTTTCGCGGAAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGcac                                                                                                          >  1:478800/1‑151 (MQ=255)                                      ggCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGcaca                                                                                                         <  1:43181/149‑1 (MQ=255)                                        cGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCGCCAGAGGCGGACGGTGCGCTGATGCGCATCAAAACGCACATccg                                                                                                     <  1:12186/151‑1 (MQ=255)                                         ggACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATccgc                                                                                                    >  1:301094/1‑151 (MQ=255)                                          gACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATccgcc                                                                                                   <  2:295777/151‑1 (MQ=255)                                             ggTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCcgcg                                                                                                >  2:213513/1‑151 (MQ=255)                                                aaaTTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCt                                                                                              <  2:251210/150‑1 (MQ=255)                                                      tGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGt                                                                                        >  1:421822/1‑150 (MQ=255)                                                       gTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTAt                                                                                      >  1:260066/1‑151 (MQ=255)                                                             aaaGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCg                                                                                <  1:30349/151‑1 (MQ=255)                                                                   gTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCa                                                                          <  2:413633/151‑1 (MQ=255)                                                                        cTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAggg                                                                       >  1:91664/1‑149 (MQ=255)                                                                           gCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCa                                                                  <  2:657718/151‑1 (MQ=255)                                                                              cTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCt                                                               >  2:677855/1‑151 (MQ=255)                                                                                gACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCt                                                               <  1:7284/149‑1 (MQ=255)                                                                                gACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTc                                                             >  1:119279/1‑151 (MQ=255)                                                                                  cGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTcc                                                            <  2:181010/150‑1 (MQ=255)                                                                                  cGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCa                                                           <  2:30994/151‑1 (MQ=255)                                                                                   gggTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCa                                                           <  1:177832/150‑1 (MQ=255)                                                                                         ggACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTg                                                    <  1:732048/151‑1 (MQ=255)                                                                                                        gTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGcc                                     <  2:723355/151‑1 (MQ=255)                                                                                                             aGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTgg                                >  1:417229/1‑151 (MQ=255)                                                                                                              gggCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGgc                               <  2:109290/151‑1 (MQ=255)                                                                                                                 cATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGgcgg                             >  1:194946/1‑150 (MQ=255)                                                                                                                      ttCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGgcggtg                           <  1:253604/147‑1 (MQ=255)                                                                                                                          cgGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATc                   >  2:483055/1‑151 (MQ=255)                                                                                                                          cgGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATc                   >  2:508430/1‑151 (MQ=255)                                                                                                                            gtagtaATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTg                 >  1:410499/1‑151 (MQ=255)                                                                                                                             tagtaATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGc                >  2:226206/1‑151 (MQ=255)                                                                                                                                 aaTGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCg             <  1:686425/150‑1 (MQ=255)                                                                                                                                 aaTGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCg             <  1:555649/150‑1 (MQ=255)                                                                                                                                   tGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGc            >  2:655878/1‑149 (MQ=255)                                                                                                                                           tCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCg   <  2:71453/150‑1 (MQ=255)                                                                                                                                            cACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGc  <  1:154048/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                 TGCTGCTGCACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCGC  >  NC_000913/4517119‑4517407

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