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%	113.8 / NA	36	*256* (TAG→TAA)	fecE	ferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/20);  total (16/20)

CCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTG  >  NC_000913/4510555‑4510823                                                                                                                                        |                                                                                                                                      ccTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTAc                                                                                                                                    >  2:44722/1‑139 (MQ=255)                   aCGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGc                                                                                                                  <  1:379053/139‑1 (MQ=255)                         aGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCa                                                                                                            >  1:681094/1‑139 (MQ=255)                                   aTCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCt                                                                                                  <  2:690669/139‑1 (MQ=255)                                   aTCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCt                                                                                                  >  2:665042/1‑139 (MQ=255)                                   aTCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCt                                                                                                  >  1:297107/1‑139 (MQ=255)                                     cTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                                <  2:652634/139‑1 (MQ=255)                                     cTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                                >  2:40044/1‑139 (MQ=255)                                        tCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCt                                                                                                  <  1:80262/134‑1 (MQ=255)                                        tCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCt                                                                                                  >  2:80262/1‑134 (MQ=255)                                         cAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGgt                                                                                            >  1:361772/1‑139 (MQ=255)                                               ccGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTAttt                                                                                      >  1:117397/1‑139 (MQ=255)                                                          cAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACg                                                                           <  1:43861/139‑1 (MQ=255)                                                          cAGCCCCAGTCTGAACAGCCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACg                                                                           <  1:517983/139‑1 (MQ=255)                                                            gCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATTGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCt                                                                         <  2:545198/139‑1 (MQ=255)                                                                 aGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATa                                                                    <  1:596665/139‑1 (MQ=255)                                                                 aGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATa                                                                    <  2:221923/139‑1 (MQ=255)                                                                           ggCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAg                                                          <  1:44722/139‑1 (MQ=255)                                                                            gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGc                                                                                                   <  1:143185/97‑1 (MQ=255)                                                                            gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGc                                                                                                   >  2:143185/1‑97 (MQ=255)                                                                            gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                                <  1:553943/100‑1 (MQ=255)                                                                            gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                                >  2:553943/1‑100 (MQ=255)                                                                            gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCCGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGc                                                         <  2:544989/139‑1 (MQ=255)                                                                             cTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGcc                                                                                                             <  2:373585/86‑1 (MQ=255)                                                                             cTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGcc                                                                                                             >  1:373585/1‑86 (MQ=255)                                                                                      gcCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCa                                                                                                            <  2:511851/78‑1 (MQ=255)                                                                                      gcCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCa                                                                                                            >  1:511851/1‑78 (MQ=255)                                                                                                               gtttggttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTg                      <  2:326964/139‑1 (MQ=255)                                                                                                                  tggttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtcatc                                          <  1:767603/116‑1 (MQ=255)                                                                                                                  tggttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtcatc                                          >  2:767603/1‑116 (MQ=255)                                                                                                                    gttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCa                 >  2:135997/1‑139 (MQ=255)                                                                                                                    gttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCa                 <  1:269999/139‑1 (MQ=255)                                                                                                                                tgtgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGt     >  1:271973/1‑139 (MQ=255)                                                                                                                                  tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCggg                                                                                              >  1:509583/1‑48 (MQ=255)                                                                                                                                  tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCggg                                                                                              <  2:509583/48‑1 (MQ=255)                                                                                                                                   gCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTg  <  1:135997/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      CCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTG  >  NC_000913/4510555‑4510823

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