Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,511,516	A→G	V300A (GTG→GCG)	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,516	0	A	G	100.0%	80.2 / NA	25	V300A (GTG→GCG)	fecD	ferric citrate ABC transporter membrane subunit FecD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (9/16);  total (9/16)

TGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTAC  >  NC_000913/4511392‑4511642                                                                                                                             |                                                                                                                               tGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAAcc                                                                                                                        >  1:409084/1‑133 (MQ=255)   aaaCGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCgg                                                                                                                      <  2:458056/133‑1 (MQ=255)    aaCGTCGTTAAGTACATTGTCTGTGCCGTAACAGACCGTCAAATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCggg                                                                                                                     <  2:203738/133‑1 (MQ=255)         cgtTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGACAATGGCGGACAGCGCGCCAACCCGGAGccc                                                                                                                >  1:233479/1‑131 (MQ=255)               gTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTgg                                                                                                          <  1:402108/133‑1 (MQ=255)                      gtctgtcCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTg                                                                                                           <  2:485807/125‑1 (MQ=255)                      gtctgtcCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTg                                                                                                           >  1:485807/1‑125 (MQ=255)                       tctgtcCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAAGGGGGGGTGaa                                                                                                  >  1:13170/1‑133 (MQ=255)                         tgtcCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGaata                                                                                                >  2:321894/1‑133 (MQ=255)                          gtcCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGaataa                                                                                               >  1:456222/1‑133 (MQ=255)                                tAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTTAACGCGCCCAACCGGATCTCCCGTGGGGGGTGGATAATTCTCg                                                                                         >  2:228264/1‑133 (MQ=255)                                tAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAAGCGGAAGCACCAGTGGGAGATGGATAATTCTCg                                                                                         >  2:147789/1‑133 (MQ=255)                                     gACCGTCACATTTTCTGTCCGTCAAGTCATCTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCcagc                                                                                    <  1:147789/133‑1 (MQ=255)                                      aCCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCcaaca                                                                                   >  1:61462/1‑133 (MQ=255)                                                 ttCAGTTCCTAAAGTAATTTATTGCATTCTGACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAacc                                                                        <  1:228264/133‑1 (MQ=255)                                                        cGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAGCCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACa                                                                 >  1:73349/1‑133 (MQ=255)                                                           aaaGTCATATATCGCCTTATCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCa                                                              <  1:150500/133‑1 (MQ=255)                                                                         cATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAgg                                                >  1:22396/1‑133 (MQ=255)                                                                            tctcACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCCACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGcc                                             >  1:125893/1‑133 (MQ=255)                                                                               cacaAGCAACCAGACGAACCACGGGGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCACCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACAGACAGCAGCTGTCCTTTCTGGGCCGaa                                          <  2:233479/133‑1 (MQ=255)                                                                               cacaAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGaa                                          <  2:409084/133‑1 (MQ=255)                                                                                   aGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCCGCAGATCGGCAACCACCGACAGCAAAGCACCTGTCAGGGCCGAAActg                                      >  1:177089/1‑131 (MQ=255)                                                                                         cAGACAAACCTCGGCTCACCGATAATGGCGGTGAGCGCGCCTACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGgcagca                                <  2:137292/133‑1 (MQ=255)                                                                                              aaaCCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTg                           <  2:30541/133‑1 (MQ=255)                                                                                                aCCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCg                         <  1:321894/133‑1 (MQ=255)                                                                                                aCCACGGCGCACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCg                         <  2:456222/133‑1 (MQ=255)                                                                                                      gcgcACCGATCATGGAGGTCAGCGCGCCAACCAGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACg                   <  2:177089/133‑1 (MQ=255)                                                                                                      gcgcACCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACg                   >  1:325404/1‑133 (MQ=255)                                                                                                          aCCGATAATGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCAGAAACAGTCAGTAGCCTGCGGAGACGTCCa               >  1:216561/1‑133 (MQ=255)                                                                                                                 aTGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGa        <  1:282238/133‑1 (MQ=255)                                                                                                                  tGGCGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGAt       <  1:116307/133‑1 (MQ=255)                                                                                                                     cGGTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCt    <  2:325404/133‑1 (MQ=255)                                                                                                                      ggTCAGCGCGCCAACCGGGAGCTCCAGTGGGGGGTGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCTAAACAGGCGGCACTCTGCGGTGACGGCCACCCGTGATGCTa   >  2:136577/1‑133 (MQ=255)                                                                                                                       gTCAGCGCGCGAACCTGGAGCTCCAGTGAGGGTTGAATAATGCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTAc  <  2:216561/133‑1 (MQ=255)                                                                                                                             |                                                                                                                               TGAAACGTCGTTAAGTACCTTGTCTGTCCCGTAACTGACCGTCAGATTTTCAGTTCGTAAAGTCATTTATCGCATTCTCACAAGCAACCAGACAAACCACGGCGCACCGATAATGGCGGTCAGCACGCCAACCGGGAGCTCCAGTGGGGGATGAATAATTCTCGCCAGCAGATCGGCAACCACCAACAGCAACGCACCTGTCAGGGCCGAAACAGGCAGCAGTCTGCGGTGACGTCCACCGGTGATGCTAC  >  NC_000913/4511392‑4511642

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