Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,511,496	C→A	A307S (GCG→TCG)	fecD ←	ferric citrate ABC transporter membrane subunit FecD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,511,496	0	C	A	100.0%	108.7 / NA	31	A307S (GCG→TCG)	fecD	ferric citrate ABC transporter membrane subunit FecD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/15);  total (16/15)

CCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCA  >  NC_000913/4511419‑4511573                                                                              |                                                                              cccGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAAcc                                                                            >  2:831321/1‑81 (MQ=255)   cGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGa                                                                          >  2:831330/1‑81 (MQ=255)   cGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGa                                                                          <  2:831291/81‑1 (MQ=255)    gTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGa                                                                          >  2:831334/1‑80 (MQ=255)      aaCTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATaa                                                                       >  2:831328/1‑81 (MQ=255)      aaCTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATaa                                                                       >  1:831325/1‑81 (MQ=255)       aCTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAAt                                                                      <  1:831292/81‑1 (MQ=255)              gTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTc                                                               <  2:831293/81‑1 (MQ=255)                 aGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAg                                                             >  2:831342/1‑80 (MQ=255)                 aGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGc                                                            <  2:831294/81‑1 (MQ=255)                  gATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCa                                                           <  2:831295/81‑1 (MQ=255)                       tCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGcc                                                       >  2:831331/1‑80 (MQ=255)                         aGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCaa                                                     <  2:831296/80‑1 (MQ=255)                          gTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAAcc                                                   <  1:831297/81‑1 (MQ=255)                          gTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAAcc                                                   <  2:831298/81‑1 (MQ=255)                           ttCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCg                                                  <  1:831299/81‑1 (MQ=255)                              gTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGa                                               <  1:831300/81‑1 (MQ=255)                                 aaGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGAGc                                             <  1:831301/80‑1 (MQ=255)                                 aaGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGAGCt                                            >  2:831344/1‑81 (MQ=255)                                  aGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTc                                           <  1:831302/81‑1 (MQ=255)                                      aTTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGt                                       >  1:831348/1‑81 (MQ=255)                                         tATCGCATTCTCACAAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTgg                                     <  2:831303/80‑1 (MQ=255)                                                     acaAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCAAGTGGGGGATGAATAAtt                        >  2:831308/1‑81 (MQ=255)                                                      caAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATtc                       >  1:831335/1‑81 (MQ=255)                                                      caAGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATtc                       >  2:831333/1‑81 (MQ=255)                                                        aGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATtctc                     <  1:831304/81‑1 (MQ=255)                                                        aGCAACCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATtctc                     >  1:831350/1‑81 (MQ=255)                                                           aaCCAGACAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGcc                  <  1:831305/81‑1 (MQ=255)                                                                  cAAACCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCcag                >  1:831343/1‑76 (MQ=255)                                                                   aaaCCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCcag                >  2:831312/1‑75 (MQ=255)                                                                    aaCCACGGCGAACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCcag                >  1:831347/1‑74 (MQ=255)                                                                           gcgaACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAacca  >  1:831337/1‑81 (MQ=255)                                                                              |                                                                              CCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCA  >  NC_000913/4511419‑4511573

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