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%	182.6 / NA	51	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 (19/32);  total (19/32)

CACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCG  >  NC_000913/4517127‑4517406                                                                                                                                         |                                                                                                                                                cACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACg                                                                                                                                   <  1:762481/151‑1 (MQ=255)   cGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGt                                                                                                                                 <  1:637685/151‑1 (MQ=255)      cAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTgcg                                                                                                                              <  1:846046/151‑1 (MQ=255)        aGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCt                                                                                                                            <  2:679037/151‑1 (MQ=255)        aGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCt                                                                                                                            <  2:618541/151‑1 (MQ=255)         gCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTg                                                                                                                           >  2:859494/1‑151 (MQ=255)            gCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGa                                                                                                                          <  1:214161/149‑1 (MQ=255)              gTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGTCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATgc                                                                                                                       <  2:353962/150‑1 (MQ=255)                gAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATgcgc                                                                                                                     >  1:467694/1‑150 (MQ=255)                 aaaTTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCAt                                                                                                                   <  2:260251/151‑1 (MQ=255)                   aTTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCa                                                                                                                 <  1:468715/151‑1 (MQ=255)                     tATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCaa                                                                                                                <  1:70599/150‑1 (MQ=255)                      aTCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCaaaa                                                                                                              <  2:905357/151‑1 (MQ=255)                             tgGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCCTCAAAACGcaca                                                                                                        >  2:471675/1‑150 (MQ=255)                               gCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATcc                                                                                                     <  2:171510/151‑1 (MQ=255)                                  gACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATccgcc                                                                                                  >  1:54829/1‑151 (MQ=255)                                  gACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATccgcc                                                                                                  >  1:255078/1‑151 (MQ=255)                                   aCGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATccgccg                                                                                                 >  1:606649/1‑151 (MQ=255)                                   aCGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATccgccg                                                                                                 <  2:686388/151‑1 (MQ=255)                                      gTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCgcgc                                                                                              >  2:106876/1‑151 (MQ=255)                                        aaaTTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCtt                                                                                            <  2:567523/151‑1 (MQ=255)                                        aaaTTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCtt                                                                                            >  2:90174/1‑151 (MQ=255)                                          aTTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTg                                                                                          <  1:846604/151‑1 (MQ=255)                                            tCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTgg                                                                                         <  2:753497/150‑1 (MQ=255)                                             cTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTggg                                                                                        >  1:497203/1‑150 (MQ=255)                                                tCCCTAAAGCCGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGTGTTGCAGGGCATCTTTCGCGGTAGTAAGGGCTATTTCGCCGCACAAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTAt                                                                                     <  2:790214/150‑1 (MQ=255)                                                    tAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGc                                                                                >  1:158114/1‑151 (MQ=255)                                                            tGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGc                                                                          >  1:438595/1‑149 (MQ=255)                                                                 tGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGc                                                                   >  2:49412/1‑151 (MQ=255)                                                                  ggCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGcc                                                                  >  1:273263/1‑151 (MQ=255)                                                                   gCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCa                                                                 <  2:367045/151‑1 (MQ=255)                                                                      cggACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCt                                                              <  2:110040/149‑1 (MQ=255)                                                                        gACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCt                                                              <  2:542972/149‑1 (MQ=255)                                                                         aCGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTcc                                                           >  1:709365/1‑151 (MQ=255)                                                                          cGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCa                                                          <  1:335674/151‑1 (MQ=255)                                                                            ggTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGt                                                        <  1:499860/151‑1 (MQ=255)                                                                                   acacGGAAGGGGGGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGAc                                                 <  2:798782/151‑1 (MQ=255)                                                                                       ggAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCt                                             <  1:198248/151‑1 (MQ=255)                                                                                       ggAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCt                                             <  2:757220/151‑1 (MQ=255)                                                                                        gAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTg                                            <  1:365885/151‑1 (MQ=255)                                                                                                gTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGc                                     <  2:731977/150‑1 (MQ=255)                                                                                                  tGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCt                                   <  2:273237/150‑1 (MQ=255)                                                                                                   gCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCttt                                 <  2:797029/151‑1 (MQ=255)                                                                                                      gggCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGgc                              >  2:406683/1‑151 (MQ=255)                                                                                                       ggCATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGgcg                             <  2:599436/151‑1 (MQ=255)                                                                                                         cATCTTTCGCGGTAGTAATGGCGATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGgcgg                            >  1:220323/1‑150 (MQ=255)                                                                                                                               gATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCt     >  1:234317/1‑151 (MQ=255)                                                                                                                               gATTTCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCt     >  1:103487/1‑151 (MQ=255)                                                                                                                                 tttCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTc   >  2:390382/1‑151 (MQ=255)                                                                                                                                  ttCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTc   <  2:136155/150‑1 (MQ=255)                                                                                                                                   tCACCGCACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCg  <  1:893841/150‑1 (MQ=255)                                                                                                                                         |                                                                                                                                                CACGTCAAGCTGCGTGAAATTATCCTGCTGGCGGACGGTAAATTCTGTCCCTAAAGCAGTGAGCTGGCCCTGACGGGTCAGGACACGGAAGGGGCGTTGCAGGGCATCTTTCGCGGTAGTAATGGCGATTTCACCGTACCAGAGCCGGACGGTGCGCTGATGCGCATCAAAACGCACATCCGCCGCGCTTTGGGTATTCAGCGTGAGCAGGGAGCCATCTTCCAGTTGCTGACGGCTGACTGTGCCTTTGGCGGTGCGGTAATCTGCCCGCAGACCTTCG  >  NC_000913/4517127‑4517406

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