Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,348,587	C→T	A25V (GCG→GTG)	yfaE →	ferredoxin involved with ribonucleotide reductase diferric‑tyrosyl radical (Y*) cofactor maintenance

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,348,587	0	C	T	100.0%	171.0 / NA	49	A25V (GCG→GTG)	yfaE	ferredoxin involved with ribonucleotide reductase diferric‑tyrosyl radical (Y*) cofactor maintenance
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/31);  total (18/31)

GGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTT  >  NC_000913/2348439‑2348733                                                                                                                                                     |                                                                                                                                                     ggaagtCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGgtgc                                                                                                                                                    >  1:498860/1‑151 (MQ=255)   aagtCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTg                                                                                                                                                  >  2:831294/1‑151 (MQ=255)    agtCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTgg                                                                                                                                                 <  2:158721/151‑1 (MQ=255)     gtCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGa                                                                                                                                                <  2:851842/151‑1 (MQ=255)     gtCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGa                                                                                                                                                <  1:501393/151‑1 (MQ=255)              tATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCca                                                                                                                                          <  2:670932/148‑1 (MQ=255)              tATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCca                                                                                                                                          <  2:461306/148‑1 (MQ=255)                    gTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTgg                                                                                                                                  >  1:33281/1‑150 (MQ=255)                    gTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGc                                                                                                                                 <  2:250876/151‑1 (MQ=255)                     tCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCg                                                                                                                                >  1:293569/1‑151 (MQ=255)                       gggCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGt                                                                                                                              <  1:402719/151‑1 (MQ=255)                                    cGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAgt                                                                                                                 <  1:735965/149‑1 (MQ=255)                                     ggAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTc                                                                                                              <  1:160984/151‑1 (MQ=255)                                       aaGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTcgc                                                                                                            <  2:568263/151‑1 (MQ=255)                                               acCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGGTGG‑GGA‑CCAATGGCGCGGAgg                                                                                                       >  1:413690/1‑148 (MQ=255)                                                 cgacgaTTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACt                                                                                                  <  1:73010/151‑1 (MQ=255)                                                  gacgaTTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGGTGC‑GTA‑CCCGTGTTGGGAAgggttct                                                                                                  >  1:108887/1‑145 (MQ=255)                                                    cgaTTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCg                                                                                               <  1:359405/151‑1 (MQ=255)                                                              aaCTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTTGCGAAGGTTACTGCGGCTCCTGTc                                                                                      <  2:546172/150‑1 (MQ=255)                                                              aaCTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTc                                                                                      <  2:83744/150‑1 (MQ=255)                                                                    cAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTAT‑GTA‑ACCACGTTGCGGAGGGTACTGCGGTTCCTGTCGcac                                                                                  >  1:700580/1‑148 (MQ=255)                                                                     aGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGG‑GTA‑CCCATGTCGCGAAAGGTACTGGGGGGCCCGTCGcaca                                                                                 >  2:816013/1‑148 (MQ=255)                                                                           gATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTg                                                                           <  1:110250/148‑1 (MQ=255)                                                                           gATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGt                                                                         <  2:413722/150‑1 (MQ=255)                                                                            aTGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTg                                                                       <  2:817323/151‑1 (MQ=255)                                                                                  cgcgTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTc                                                                 <  2:242829/151‑1 (MQ=255)                                                                                        aCCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGGTGG‑GGA‑ACAGTGTCGCGAGGGATACTGCGGGCACTGTCGCCCCCCCCTggcttcccggcaaggt                                                            >  1:710832/1‑136 (MQ=255)                                                                                         cccTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGa                                                          <  2:529711/151‑1 (MQ=255)                                                                                           cTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAAGGTCGCGAAGGGTTCTTCGGCCCCCGTCGCACCCGCCTGGGTGTAGGGCAAggt                                                            >  1:339407/1‑145 (MQ=255)                                                                                              cgcATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑ACAGGGTCGCCAAGGGTAATTCGGGCCCCTTCGCCCACCCCTGGGTGCAggggcaggtgacagg                                                      >  1:806273/1‑137 (MQ=255)                                                                                               gcATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTTA‑GTA‑CCCGAGTCGCGGAGGGTTATGCGGCCCCCGTCGCCCACGCCCGGGTGCAGGGCAAGTTGAATgg                                                      >  2:53908/1‑149 (MQ=255)                                                                                               gcATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGGTGA‑GGA‑CCAGTGTCGCGGAGGGTACTGCGGCTCCTGTCGCCCACGCCCCGGTGCCGGGCCAGGTGACTgg                                                      >  2:817026/1‑149 (MQ=255)                                                                                                cATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTTA‑GGA‑CCAGGGTCGCGAAGGGTAATGCGGCCCCTGTCGCACACCCCCGGTTGCGGGGCAAATTTACTggggt                                                   >  1:618528/1‑148 (MQ=255)                                                                                                cATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGAtt                                                   <  2:629150/151‑1 (MQ=255)                                                                                                   cACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTTAGGTA‑CC‑GTGTCGGGGAGGGGAATGGGGCCCCTTTCGCACACCCCCCGTTGGAGGGCCAGTTGagcgggttg                                                  >  2:686758/1‑141 (MQ=255)                                                                                                     cTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGGTGG‑GGACCCAGGTTCCGGAGGGTTACTGC‑GCGCCCGTCGCACCCGCCTGGTTTCAGGGCAAGTTGGCTGGGTTgtc                                                >  1:732081/1‑147 (MQ=255)                                                                                                         cacacaACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCg                                          <  1:537356/151‑1 (MQ=255)                                                                                                              aaCTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGc                                     <  2:337289/151‑1 (MQ=255)                                                                                                                  gctgTGCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTTGAATCCCACAATGTGCCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTGACTGCGTCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCttt                                 <  1:388465/151‑1 (MQ=255)                                                                                                                     gtgCCAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTAtt                              <  1:636463/151‑1 (MQ=255)                                                                                                                        ccAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAg                           <  2:318021/151‑1 (MQ=255)                                                                                                                         cAGGATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGGGTTGCGAAGGGTACTGCGGCTCCTGTCGCACACCCCTGGTTGGAGGGCAAGTTGACTGGGTTGCCGAACCGTTAGCCTTTTTTTAGc                          >  2:347219/1‑151 (MQ=255)                                                                                                                          aagattaataacctgcccttcTGGCGGTGCTGGAATCCCACAATGTCGCGGTTGA‑GTA‑CCAGGGTCCCGAAGGTTACTGGGGTGCTTGTCGCACACTCCTGGTTGTAGGTCAAGTTTACTGGATTGCCTAACCGTTAGCCTTTATTCAGc                          <  1:795540/139‑1 (MQ=255)                                                                                                                           ggATGAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTTA‑GTA‑CCACTGTCGCGGAGGGTAACGCGGCCCCCGTCGCCCACGCCCGGGTGGCGGGCCAGTTTGCTGGGTTTCCGCACCCTTAACCtttttttcggc                         >  1:616407/1‑145 (MQ=255)                                                                                                                               gAACACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCgggg                     <  2:223849/150‑1 (MQ=255)                                                                                                                                aaCACCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGa                   <  1:816970/151‑1 (MQ=255)                                                                                                                                 acacCCTTCCCTTCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGaa                  <  2:558587/151‑1 (MQ=255)                                                                                                                                    cccttcccttCTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAAtt               <  2:304577/151‑1 (MQ=255)                                                                                                                                      cttcccttcTGGCGGTGCTGGAATCCCACAATGTGGCGGTTTT‑GTA‑CCAGGGTCGCGAAGGGTACTGCGGGTCCTTTCGCACACGCCTGGGTTCCGGGCAAGTTGAATGGGTTGCCGAACCGGTAGCCTTTATTCTGCCGGGGGGAAtttt             >  2:673404/1‑151 (MQ=255)                                                                                                                                             tcTGGCGGTGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCt        <  2:496330/149‑1 (MQ=255)                                                                                                                                                ggcggtGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCtgt      <  1:833903/148‑1 (MQ=255)                                                                                                                                                ggcggtGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCtgt      <  2:732124/148‑1 (MQ=255)                                                                                                                                                 gcggtgCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTgttgtt  <  1:303285/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GGAAGTCAGTTCTTATCTGGTCGGGCAGATTGACTCGGAAGTGGACACCGACGATTTGAGTAACTTCCAGCTCTGATGGCCCGCGTTACCCTGCGCATCACTGGCACACAACTGCTGTGCCAGGATGAACACCCTTCCCTTCTGGCGGCGCTGGAATCCCACAATGTGGCGGTTGA‑GTA‑CCAGTGTCGCGAAGGTTACTGCGGCTCCTGTCGCACACGCCTGGTTGCAGGTCAAGTTGACTGGATTGCCGAACCGTTAGCCTTTATTCAGCCGGGGGAAATTTTGCCCTGTTGTT  >  NC_000913/2348439‑2348733

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