Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,082,943	A→G	Q234R (CAG→CGG)	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,943	0	A	G	100.0%	142.4 / NA	40	Q234R (CAG→CGG)	efeO	inactive ferrous ion transporter EfeUOB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/23);  total (17/23)

CTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGACTGATTGAGGAAGTGGCAGCCAGCAAAATTAGCGGTG  >  NC_000913/1082812‑1083077                                                                                                                                    |                                                                                                                                       ctcCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCt                                                                                                                                 >  2:564956/1‑139 (MQ=255)    cGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGag                                                                                                                              <  1:1013166/139‑1 (MQ=255)    cGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGag                                                                                                                              <  1:771119/139‑1 (MQ=255)    cGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGag                                                                                                                              <  1:564956/139‑1 (MQ=255)    cGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGag                                                                                                                              <  2:742066/139‑1 (MQ=255)    cGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGag                                                                                                                              <  2:391211/139‑1 (MQ=255)      aTCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGagca                                                                                                                            >  2:5169/1‑139 (MQ=255)         tGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCt                                                                                                                         >  2:246212/1‑139 (MQ=255)           gATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCt                                                                                                                         >  1:268452/1‑137 (MQ=255)           gATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCt                                                                                                                         <  2:268452/137‑1 (MQ=255)                                                 aaaaGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGaa                                                                                 >  2:1322046/1‑139 (MQ=255)                                                  aaaGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAAc                                                                                <  2:1210468/139‑1 (MQ=255)                                                       cgccgACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCt                                                                           <  1:316427/139‑1 (MQ=255)                                                                           gTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCaaaaa                                                                                             <  1:137683/101‑1 (MQ=255)                                                                           gTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCaaaaa                                                                                             >  2:137683/1‑101 (MQ=255)                                                                                 aCCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTcggcgg                                                 <  1:246212/139‑1 (MQ=255)                                                                                     tCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTa                                                                                                                                     <  2:480669/51‑1 (MQ=255)                                                                                     tCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTa                                                                                                                                     >  1:480669/1‑51 (MQ=255)                                                                                     tCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCaa                                                                                                >  1:709202/1‑88 (MQ=255)                                                                                     tCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCaa                                                                                                <  2:709202/88‑1 (MQ=255)                                                                                       tGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGc                                           >  1:661381/1‑139 (MQ=255)                                                                                              aGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGACTg                                    <  1:470502/139‑1 (MQ=255)                                                                                               gCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAAc                                                                                >  2:126883/1‑94 (MQ=255)                                                                                               gCATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAAc                                                                                <  1:126883/94‑1 (MQ=255)                                                                                                cATTGTTTGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGACTGAt                                  <  2:555419/139‑1 (MQ=255)                                                                                                     tttGGCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGACTGATTGAgg                             <  2:661381/139‑1 (MQ=255)                                                                                                         gCGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGACTGATTGAGGAAGt                         >  2:723870/1‑139 (MQ=255)                                                                                                          cGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGa                                       <  2:1225598/124‑1 (MQ=255)                                                                                                          cGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGa                                       >  1:1225598/1‑124 (MQ=255)                                                                                                          cGACAACACCACCAAAGGGATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGACTGATTGAGGAAGTg                        >  2:981182/1‑139 (MQ=255)                                                                                                                          gggATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCaaaa                                                                                              <  2:185079/53‑1 (MQ=255)                                                                                                                          gggATGGATCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCaaaa                                                                                              >  1:185079/1‑53 (MQ=255)                                                                                                                           ggatggatCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGACTGATTGAGGAAGTGGCAGCCAGCAAAATTAg       <  1:5169/139‑1 (MQ=255)                                                                                                                             atggatCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGACTGATTGAGGAAGTGGCAGCCAGCAAAATTAGCg     <  2:637898/139‑1 (MQ=255)                                                                                                                               ggatCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACg                                                                                           <  1:946064/51‑1 (MQ=255)                                                                                                                               ggatCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACg                                                                                           >  2:946064/1‑51 (MQ=255)                                                                                                                               ggatCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGACTGATTGAGGAAGTGGCAGCCAGCAAAATTAGCGGt   >  1:552618/1‑139 (MQ=255)                                                                                                                                gatCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGa                                       >  2:1015854/1‑102 (MQ=255)                                                                                                                                gatCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGa                                       <  1:1015854/102‑1 (MQ=255)                                                                                                                                gatCGGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGACTGATTGAGGAAGTGGCAGCCAGCAAAATTAGCGGTg  <  1:139488/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                       CTCCGATCTGGATGGCAGCATTGACGCCCGTGAAGATGATTACGAGCAAAAAGCCGCCGACCCAAAATTCACTGGTTTCCACCGTCTGGAAAAAGCATTGTTTGGCGACAACACCACCAAAGGGATGGATCAGTACGCTGAGCAGCTTTATACCGATGTGGTCGATTTGCAAAAACGCATCAGTGAACTGGCTTTCCCACCTTCAAAAGTGGTCGGCGGCGCAGCCGGACTGATTGAGGAAGTGGCAGCCAGCAAAATTAGCGGTG  >  NC_000913/1082812‑1083077

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