Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,510,690	C→T	*256* (TAG→TAA)	fecE ←	ferric citrate ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,510,690	0	C	T	100.0%	128.2 / NA	41	*256* (TAG→TAA)	fecE	ferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/20);  total (21/20)

AAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCC  >  NC_000913/4510558‑4510818                                                                                                                                     |                                                                                                                                 aaaGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGCGGTTTGGTTCTTACGACCTGTGCAATTTACCTc                                                                                                                            >  2:310912/1‑139 (MQ=255)             tGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATc                                                                                                                >  1:411275/1‑139 (MQ=255)                      aGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCa                                                                                                       >  1:151944/1‑139 (MQ=255)                      aGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCa                                                                                                       >  1:583011/1‑139 (MQ=255)                          ttGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATa                                                                                                   <  2:583011/139‑1 (MQ=255)                                aTCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCt                                                                                             >  1:170548/1‑139 (MQ=255)                                aTCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTGGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCt                                                                                             >  2:953380/1‑139 (MQ=255)                                  cTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                           >  2:370794/1‑139 (MQ=255)                                  cTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                           <  1:953380/139‑1 (MQ=255)                                  cTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                           <  2:564836/139‑1 (MQ=255)                                  cTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                           <  1:370794/139‑1 (MQ=255)                                     tCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTa                                                                                                                                >  1:197426/1‑99 (MQ=255)                                     tCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTa                                                                                                                                <  2:197426/99‑1 (MQ=255)                                           tCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGc                                                                                              <  1:1176070/127‑1 (MQ=255)                                           tCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGc                                                                                              >  2:1176070/1‑127 (MQ=255)                                            ccGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTAttt                                                                                 >  1:162933/1‑139 (MQ=255)                                             cGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTAcc                                                                                                                              <  1:214549/93‑1 (MQ=255)                                             cGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTAcc                                                                                                                              >  2:214549/1‑93 (MQ=255)                                             cGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTc                                                                                <  2:411275/139‑1 (MQ=255)                                             cGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTc                                                                                >  2:98209/1‑139 (MQ=255)                                                      tCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCAc                                                                       <  2:511181/139‑1 (MQ=255)                                                       cAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACg                                                                      <  1:190428/139‑1 (MQ=255)                                                                        ggCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCGCGCTGAATACTGTTCTCAg                                                     <  2:966840/139‑1 (MQ=255)                                                                         gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGc                                                    >  1:1164506/1‑139 (MQ=255)                                                                              cTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                           >  1:347504/1‑95 (MQ=255)                                                                              cTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                           <  2:347504/95‑1 (MQ=255)                                                                               tGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCt                                              <  1:359831/139‑1 (MQ=255)                                                                                 gcgcCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTgg                                            >  2:413706/1‑139 (MQ=255)                                                                                   gcCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAg                                                     <  1:1147908/128‑1 (MQ=255)                                                                                   gcCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAg                                                     >  2:1147908/1‑128 (MQ=255)                                                                                        tccagtTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtcatc                                     <  1:239376/139‑1 (MQ=255)                                                                                               tcattcaGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTAc                                                                                                                               <  2:1076269/42‑1 (MQ=255)                                                                                               tcattcaGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTAc                                                                                                                               >  1:1076269/1‑42 (MQ=255)                                                                                                          tggtttggttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCt                   <  2:162933/139‑1 (MQ=255)                                                                                                          tggtttggttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCt                   <  2:753795/139‑1 (MQ=255)                                                                                                                 gttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGACTTGCGCCa            >  2:796226/1‑139 (MQ=255)                                                                                                                      tACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATa                                                               <  1:709013/83‑1 (MQ=255)                                                                                                                      tACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATa                                                               >  2:709013/1‑83 (MQ=255)                                                                                                                           ccTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTcc  >  1:919125/1‑139 (MQ=255)                                                                                                                               tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtc                                        <  2:466825/97‑1 (MQ=255)                                                                                                                               tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtc                                        >  1:466825/1‑97 (MQ=255)                                                                                                                                     |                                                                                                                                 AAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCC  >  NC_000913/4510558‑4510818

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