Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,514,908	C→T	L590L (CTG→CTA)	fecA ←	ferric citrate outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,514,908	0	C	T	100.0%	150.4 / NA	42	L590L (CTG→CTA)	fecA	ferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/19);  total (23/19)

ATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCA  >  NC_000913/4514830‑4514983                                                                               |                                                                            aTGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAg                                                                           >  2:935342/1‑81 (MQ=255) aTGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAg                                                                           >  1:935371/1‑81 (MQ=255) aTGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAg                                                                           >  2:935377/1‑81 (MQ=255)  tGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGc                                                                          >  2:935348/1‑81 (MQ=255)      gcgTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCa                                                                      >  2:935357/1‑81 (MQ=255)        gTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCAtt                                                                    <  2:935306/81‑1 (MQ=255)        gTTTTGCCACGGGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCAtt                                                                    >  1:935384/1‑81 (MQ=255)         ttttGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCAttt                                                                   <  1:935308/81‑1 (MQ=255)         ttttGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCAttt                                                                   <  2:935307/81‑1 (MQ=255)           ttGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTcc                                                                 >  2:935375/1‑81 (MQ=255)            tGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTcc                                                                 <  2:935309/80‑1 (MQ=255)               cACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTccgccg                                                             <  2:935310/81‑1 (MQ=255)                aCGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGt                                                            >  1:935389/1‑81 (MQ=255)                 cGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTc                                                           <  1:935312/81‑1 (MQ=255)                 cGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTc                                                           <  2:935311/81‑1 (MQ=255)                   tGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAg                                                         >  2:935379/1‑81 (MQ=255)                   tGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAg                                                         <  1:935313/81‑1 (MQ=255)                           cGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCcgt                                                 >  2:935382/1‑81 (MQ=255)                           cGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCcgt                                                 >  1:935383/1‑81 (MQ=255)                            ggTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCcgtc                                                <  1:935314/81‑1 (MQ=255)                               gtCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCcgtcgtc                                             >  2:935387/1‑81 (MQ=255)                                tCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCcgtcgtc                                             <  2:935315/80‑1 (MQ=255)                                  gTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCcgtcgtcg                                            >  2:935400/1‑79 (MQ=255)                                   ttGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAg                                         >  2:935390/1‑81 (MQ=255)                                        cTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCggg                                     >  2:935398/1‑80 (MQ=255)                                           gTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTAcc                                 <  2:935317/81‑1 (MQ=255)                                           gTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTAc                                  <  2:935316/80‑1 (MQ=255)                                              ggAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAg                              >  1:935388/1‑81 (MQ=255)                                              ggAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAg                              >  2:935396/1‑81 (MQ=255)                                               gAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAg                              <  1:935318/80‑1 (MQ=255)                                               gAGTCGTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGt                             >  2:935385/1‑81 (MQ=255)                                                    gTACTGATTGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTcc                          >  1:935408/1‑79 (MQ=255)                                                           ttGTTAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTcgc                 <  2:935319/81‑1 (MQ=255)                                                              ttAAAGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCt              >  1:935421/1‑81 (MQ=255)                                                                aaaGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCttt            >  1:935399/1‑81 (MQ=255)                                                                aaaGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCttt            >  2:935409/1‑81 (MQ=255)                                                                aaaGTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCtt             <  1:935320/80‑1 (MQ=255)                                                                   gTTAATCAGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTc         <  2:935321/81‑1 (MQ=255)                                                                          aGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTc   >  1:935401/1‑80 (MQ=255)                                                                          aGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCa  <  2:935324/81‑1 (MQ=255)                                                                          aGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCa  <  2:935323/81‑1 (MQ=255)                                                                          aGGAATAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCa  <  2:935322/81‑1 (MQ=255)                                                                               |                                                                            ATGGCGCGTTTTGCCACGTGCAGTGACGGTGTCGTTGGTCTGGTTGGAGTCGTACTGATTGTTAAAGTTAATCAGGAACAGCCCCATTTCCGCCGTCAGCGCGCCGTCGTCGTAGCGGGTACCGAGTTCCCAGGTTCGCGCTTTTTCCGGTTCA  >  NC_000913/4514830‑4514983

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