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

CCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTT  >  NC_000913/4510555‑4510821                                                                                                                                        |                                                                                                                                    ccTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTAc                                                                                                                                  <  2:746052/139‑1 (MQ=255)           gcGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGc                                                                                                                        >  1:653117/1‑139 (MQ=255)                   aCGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGc                                                                                                                <  2:33366/139‑1 (MQ=255)                         aGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCa                                                                                                          >  1:191168/1‑139 (MQ=255)                           tATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGa                                                                                                        >  2:438218/1‑139 (MQ=255)                             ttGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATa                                                                                                      >  2:577252/1‑139 (MQ=255)                                     cTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                              >  1:722353/1‑139 (MQ=255)                                     cTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                              >  2:662113/1‑139 (MQ=255)                                        tCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTa                                                                                                                                   >  1:893796/1‑99 (MQ=255)                                        tCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCgggg                                                                                           >  1:526061/1‑139 (MQ=255)                                          agagTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGgtg                                                                                         >  2:830954/1‑139 (MQ=255)                                           gagTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGgtgt                                                                                        <  2:1059694/139‑1 (MQ=255)                                           gagTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGgtgt                                                                                        <  1:859171/139‑1 (MQ=255)                                               ccGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTAttt                                                                                    >  1:268488/1‑139 (MQ=255)                                                cGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTc                                                                                   >  1:463255/1‑139 (MQ=255)                                                  gAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCAGCCTGCCAGATACCGGCTCGGGGTGTATTTCCg                                                                                 <  1:424534/139‑1 (MQ=255)                                                       gTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTcc                                                                            <  1:699241/139‑1 (MQ=255)                                                        ttCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTc                                                                                               >  2:7085/1‑119 (MQ=255)                                                        ttCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTc                                                                                               <  1:7085/119‑1 (MQ=255)                                                         tCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCAc                                                                          <  1:362737/139‑1 (MQ=255)                                                                 aGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATa                                                                  <  2:268488/139‑1 (MQ=255)                                                                  gTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATAc                                                                 >  1:936204/1‑139 (MQ=255)                                                                   tCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACt                                                                <  1:130244/139‑1 (MQ=255)                                                                      ggaCAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCACTTCCACGCTGAATACTGtt                                                             <  2:653117/137‑1 (MQ=255)                                                                      gAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTCTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGtt                                                             <  2:722353/139‑1 (MQ=255)                                                                           ggCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAg                                                        <  2:936204/139‑1 (MQ=255)                                                                           ggCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAg                                                        >  1:367282/1‑139 (MQ=255)                                                                            gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTa                                                                                                                                   >  1:265739/1‑63 (MQ=255)                                                                                           tccagtTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATa                                                                                                      <  2:593535/77‑1 (MQ=255)                                                                                           tccagtTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATa                                                                                                      >  1:593535/1‑77 (MQ=255)                                                                                           tccagtTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtcatc                                        <  1:577252/139‑1 (MQ=255)                                                                                           tccagtTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtcatc                                        <  1:438218/139‑1 (MQ=255)                                                                                             cagtTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCAc                                      <  2:367282/139‑1 (MQ=255)                                                                                                            gTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGcc                       >  1:678946/1‑139 (MQ=255)                                                                                                            gTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGcc                       <  2:586959/139‑1 (MQ=255)                                                                                                              ggtttggttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCtt                     <  2:526061/139‑1 (MQ=255)                                                                                                                  tggttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTgcgc                 <  2:267636/139‑1 (MQ=255)                                                                                                                          aCGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACAt         <  1:662113/139‑1 (MQ=255)                                                                                                                              ccTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTA‑TTCCGCTTCCACg                                                                         <  2:172007/70‑1 (MQ=255)                                                                                                                              ccTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTA‑TTCCGCTTCCACg                                                                         >  1:172007/1‑70 (MQ=255)                                                                                                                                 gtgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGtt  >  2:520525/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                    CCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTT  >  NC_000913/4510555‑4510821

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