Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,082,896	T→A	R218R (CGT→CGA)	efeO →	inactive ferrous ion transporter EfeUOB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,082,896	0	T	A	100.0%	125.6 / NA	40	R218R (CGT→CGA)	efeO	inactive ferrous ion transporter EfeUOB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (18/22);  total (18/22)

CGCGCCAGCACTATGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAG  >  NC_000913/1082769‑1083020                                                                                                                                |                                                                                                                             cgcgCCAGCACTATGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGc                                                                                                                   >  1:236013/1‑139 (MQ=255) cgcgCCAGCACTATGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGc                                                                                                                   >  1:756383/1‑139 (MQ=255)     ccAGCACTATGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTg                                                                                                               >  2:370779/1‑139 (MQ=255)      cAGCACTATGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGt                                                                                                              <  2:236013/139‑1 (MQ=255)                aGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGacaa                                                                                                    <  2:236061/139‑1 (MQ=255)                      ttGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAAcaccac                                                                                              >  2:1049570/1‑139 (MQ=255)                        gAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGC‑AAAAGCCACCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCaa                                                                                           <  1:267742/139‑1 (MQ=255)                         aaCCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAGAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCaa                                                                                           <  2:350479/139‑1 (MQ=255)                         aaCCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCaa                                                                                           <  1:132981/139‑1 (MQ=255)                         aaCCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCaa                                                                                           <  1:475774/139‑1 (MQ=255)                           ccGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAg                                                                                         <  2:382805/139‑1 (MQ=255)                                          ttCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAAcacc                                                                                                <  2:740440/117‑1 (MQ=255)                                          ttCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAAcacc                                                                                                >  1:740440/1‑117 (MQ=255)                                          ttCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACg                                                                          >  1:400993/1‑139 (MQ=255)                                          ttCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACg                                                                          >  2:918602/1‑139 (MQ=255)                                               cGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCa                                                                               >  1:200406/1‑129 (MQ=255)                                               cGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCa                                                                               <  2:200406/129‑1 (MQ=255)                                               cGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGag                                                                     <  1:918602/139‑1 (MQ=255)                                               cGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGag                                                                     <  1:1049570/139‑1 (MQ=255)                                                    tGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCAt                                                                                                                 <  1:774523/90‑1 (MQ=255)                                                    tGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCAt                                                                                                                 >  2:774523/1‑90 (MQ=255)                                                    tGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGac                                                                                                      <  1:254539/101‑1 (MQ=255)                                                    tGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGac                                                                                                      >  2:254539/1‑101 (MQ=255)                                                    tGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACAccacca                                                                                            >  1:1024375/1‑111 (MQ=255)                                                    tGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACAccacca                                                                                            <  2:1024375/111‑1 (MQ=255)                                                    tGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCt                                                                <  1:419149/139‑1 (MQ=255)                                                    tGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCt                                                                >  2:280148/1‑139 (MQ=255)                                                                          tGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTg                                          >  1:52960/1‑139 (MQ=255)                                                                                   ttACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCa                                                                                                                  <  2:966259/58‑1 (MQ=255)                                                                                   ttACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCa                                                                                                                  >  1:966259/1‑58 (MQ=255)                                                                                    tACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCa                                >  2:418510/1‑139 (MQ=255)                                                                                      cGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACAccacca                                                                                            <  1:827451/77‑1 (MQ=255)                                                                                      cGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACAccacca                                                                                            >  2:827451/1‑77 (MQ=255)                                                                                       gAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCa                             <  1:370779/139‑1 (MQ=255)                                                                                            aaaaGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACg                                  >  1:845700/1‑129 (MQ=255)                                                                                            aaaaGCCGCCGACCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACg                                  <  2:845700/129‑1 (MQ=255)                                                                                                       aCCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTccc             <  2:400993/139‑1 (MQ=255)                                                                                                       aCCCAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTccc             <  2:756383/139‑1 (MQ=255)                                                                                                        cccAAAATTCACTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCa            <  1:855147/139‑1 (MQ=255)                                                                                                                  aCTGGTTTCCACCGACTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAg  >  1:267224/1‑139 (MQ=255)                                                                                                                                |                                                                                                                             CGCGCCAGCACTATGAGCGTATTGAACCGATTGCTGAACTGTTCTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAG  >  NC_000913/1082769‑1083020

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