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

CTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTT  >  NC_000913/4510556‑4510822                                                                                                                                       |                                                                                                                                     cTAAAGCAAGCGGATGGACGATGAGTATTGGTAACCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTAcc                                                                                                                                  <  1:824889/139‑1 (MQ=255)          gcGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGc                                                                                                                         >  2:719519/1‑139 (MQ=255)                                   tCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTc                                                                                                >  2:407594/1‑139 (MQ=255)                                    cTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                               >  2:1049368/1‑139 (MQ=255)                                    cTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                               >  1:492213/1‑139 (MQ=255)                                       tCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCgggg                                                                                            >  2:139998/1‑139 (MQ=255)                                         agagTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGgtg                                                                                          <  1:139998/139‑1 (MQ=255)                                                      gTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTcc                                                                             <  2:449295/139‑1 (MQ=255)                                                       ttCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATa                                                                                                       <  2:492996/112‑1 (MQ=255)                                                       ttCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATa                                                                                                       >  1:492996/1‑112 (MQ=255)                                                        tCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCAc                                                                           <  1:719519/139‑1 (MQ=255)                                                                      aaCAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCAc                                                                           <  2:799358/125‑1 (MQ=255)                                                                      aaCAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCAc                                                                           >  1:799358/1‑125 (MQ=255)                                                                          ggCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCTGCTTCCACGCTGAATACTGTTCTCAg                                                         <  1:407594/139‑1 (MQ=255)                                                                           gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGc                                                        >  1:330331/1‑139 (MQ=255)                                                                           gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGc                                                        >  2:973816/1‑139 (MQ=255)                                                                                caggCGCCAGTCCAGTTTCTTTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACAt                                                                                                                     <  2:994845/71‑1 (MQ=255)                                                                                cTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACAt                                                                                                                     >  1:994845/1‑73 (MQ=255)                                                                                     gcCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGc                                                                                                                 <  1:499371/72‑1 (MQ=255)                                                                                     gcCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGc                                                                                                                 >  2:499371/1‑72 (MQ=255)                                                                                            cagtTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCAc                                       <  2:533621/139‑1 (MQ=255)                                                                                            cagtTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCAc                                       <  1:400384/139‑1 (MQ=255)                                                                                            cagtTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCAc                                       <  2:330331/139‑1 (MQ=255)                                                                                                   attcaGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCtcttct                                >  2:578158/1‑139 (MQ=255)                                                                                                                   gttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGc                                                                                                                 <  1:453848/42‑1 (MQ=255)                                                                                                                   gttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGc                                                                                                                 >  2:453848/1‑42 (MQ=255)                                                                                                                   gttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCt                                                                                                 >  2:107423/1‑58 (MQ=255)                                                                                                                   gttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCt                                                                                                 <  1:107423/58‑1 (MQ=255)                                                                                                                         aCGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACAt          <  1:578158/139‑1 (MQ=255)                                                                                                                             ccTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAg                                                         <  2:475239/88‑1 (MQ=255)                                                                                                                             ccTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAg                                                         >  1:475239/1‑88 (MQ=255)                                                                                                                                 tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTcc                                                                             >  1:856017/1‑64 (MQ=255)                                                                                                                                 tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTcc                                                                             <  2:856017/64‑1 (MQ=255)                                                                                                                                 tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACg                                                                          <  1:722307/67‑1 (MQ=255)                                                                                                                                 tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACg                                                                          >  2:722307/1‑67 (MQ=255)                                                                                                                                 tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAg                                                         <  2:133846/84‑1 (MQ=255)                                                                                                                                 tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAg                                                         >  1:133846/1‑84 (MQ=255)                                                                                                                                 tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtc                                            <  2:152956/97‑1 (MQ=255)                                                                                                                                 tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtc                                            >  1:152956/1‑97 (MQ=255)                                                                                                                                 tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGttt  >  2:13475/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                     CTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTT  >  NC_000913/4510556‑4510822

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