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

CTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCC  >  NC_000913/4510556‑4510818                                                                                                                                       |                                                                                                                                 cTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTAcc                                                                                                                              <  2:645675/139‑1 (MQ=255)          gcGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGc                                                                                                                     >  1:80335/1‑139 (MQ=255)          gcGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGc                                                                                                                     >  1:623085/1‑139 (MQ=255)          gcGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGc                                                                                                                     >  2:448855/1‑139 (MQ=255)              aTGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACAt                                                                                                                 <  1:917979/139‑1 (MQ=255)                  aCGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGc                                                                                                             <  2:425270/139‑1 (MQ=255)                                  aTCTTTCTGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCt                                                                                             <  1:1151462/139‑1 (MQ=255)                                  aTCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCa                                                                                                                           <  1:415820/109‑1 (MQ=255)                                  aTCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCa                                                                                                                           >  2:415820/1‑109 (MQ=255)                                  aTCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCt                                                                                             >  1:726805/1‑139 (MQ=255)                                    cTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTa                                                                                                                                >  2:565193/1‑102 (MQ=255)                                    cTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTa                                                                                                                                <  1:565193/102‑1 (MQ=255)                                    cTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                           >  1:458949/1‑139 (MQ=255)                                               cGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTc                                                                                >  1:338544/1‑139 (MQ=255)                                                      gTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTcc                                                                         <  1:448855/139‑1 (MQ=255)                                                      gTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTcc                                                                         <  2:1121408/139‑1 (MQ=255)                                                      gTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTcc                                                                         <  2:623085/139‑1 (MQ=255)                                                        tCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCAc                                                                       <  2:450069/139‑1 (MQ=255)                                                                aGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATa                                                               <  2:999084/139‑1 (MQ=255)                                                                      aaCAGGCTTGCTGGCACCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTATGCACATCGGcc                                                                                                            >  1:874377/1‑88 (MQ=255)                                                                      aaCAGGCTTGCTGGCACCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGcc                                                                                                            <  2:874377/88‑1 (MQ=255)                                                                          ggCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAg                                                     <  2:514852/139‑1 (MQ=255)                                                                           gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTAcc                                                                                                                              <  2:1211902/65‑1 (MQ=255)                                                                           gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTAcc                                                                                                                              >  1:1211902/1‑65 (MQ=255)                                                                           gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTAc                                                                                                                               <  1:865708/64‑1 (MQ=255)                                                                           gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTAc                                                                                                                               >  2:865708/1‑64 (MQ=255)                                                                           gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGc                                                                                                         >  1:765156/1‑86 (MQ=255)                                                                           gCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGc                                                                                                         <  2:765156/86‑1 (MQ=255)                                                                                cTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCt                                                                    >  1:673419/1‑118 (MQ=255)                                                                                cTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCt                                                                    <  2:673419/118‑1 (MQ=255)                                                                                  ggCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGTGTGTATTTCCGCTTCCACGCTGAATa                                                               <  1:1193800/121‑1 (MQ=255)                                                                                  ggCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATa                                                               >  2:1193800/1‑121 (MQ=255)                                                                                  ggCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTg                                             <  2:80335/139‑1 (MQ=255)                                                                                         gtccagtTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGtcat                                      <  2:458949/139‑1 (MQ=255)                                                                                            cagtTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCAc                                   <  2:338544/139‑1 (MQ=255)                                                                                             agtTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCAcc                                  <  2:726805/139‑1 (MQ=255)                                                                                               tttCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACctc                                >  1:668984/1‑139 (MQ=255)                                                                                                    ttcaGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                           <  1:886816/75‑1 (MQ=255)                                                                                                    ttcaGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCg                                                                                           >  2:886816/1‑75 (MQ=255)                                                                                                    ttcaGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCtctt                              <  1:903114/136‑1 (MQ=255)                                                                                                    ttcaGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCtctt                              >  2:903114/1‑136 (MQ=255)                                                                                                      caGTCGTGGTTTGGTTCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGgt                         <  2:266546/139‑1 (MQ=255)                                                                                                                 tggttCTTACGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTgtgc              <  2:1208792/139‑1 (MQ=255)                                                                                                                         aCGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGTCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACAt      >  1:1179985/1‑139 (MQ=255)                                                                                                                         aCGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACAt      >  2:476486/1‑139 (MQ=255)                                                                                                                         aCGGCCTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACAt      >  2:712279/1‑139 (MQ=255)                                                                                                                             ccTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCt                                                                                             >  1:1097397/1‑48 (MQ=255)                                                                                                                             ccTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCt                                                                                             <  2:1097397/48‑1 (MQ=255)                                                                                                                             ccTGTGCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTcc  >  2:587613/1‑139 (MQ=255)                                                                                                                                gtgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGtt                                                          <  2:1011693/80‑1 (MQ=255)                                                                                                                                gtgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGtt                                                          >  1:1011693/1‑80 (MQ=255)                                                                                                                                 tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTcc                                                                         >  1:1114532/1‑64 (MQ=255)                                                                                                                                 tgCAATTTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTcc                                                                         <  2:1114532/64‑1 (MQ=255)                                                                                                                                       |                                                                                                                                 CTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGTCCGGAAAAGTTCAGCCCCAGTCTGAACAGGCTTGCTGGCGCCAGTCCAGTTTCATTCAGTCGTGGTTTGGTTCTTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCC  >  NC_000913/4510556‑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: ‑