Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,216,481	(G)6→7	coding (11/765 nt)	nfeF ←	NADPH‑dependent ferric chelate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,216,475	1	.	G	100.0%	127.0 / NA	36	R6P (CGC→CCC)	nfeF	NADPH‑dependent ferric chelate reductase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (15/21);  total (15/21)

TCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATAT  >  NC_000913/3216337‑3216593                                                                                                                                            |                                                                                                                       tCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAGCACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑g                                                                                                                       <  2:201712/140‑1 (MQ=255)   gTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑ggg                                                                                                                     >  1:712006/1‑140 (MQ=255)    tCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑gggg                                                                                                                    <  1:219251/140‑1 (MQ=255)    tCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑gggg                                                                                                                    <  2:712006/140‑1 (MQ=255)     cAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑ggggg                                                                                                                   >  2:689996/1‑140 (MQ=255)         gCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑gggggggtg                                                                                                               <  2:817668/140‑4 (MQ=255)             cgcgATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTAt                                                                                                            >  2:132771/1‑140 (MQ=255)             cgcgATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTAt                                                                                                            >  1:774840/1‑140 (MQ=255)             cgcgATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTAt                                                                                                            >  2:619342/1‑140 (MQ=255)               cgATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTc                                                                                                          <  1:286681/140‑1 (MQ=255)                gATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCa                                                                                                         <  1:572809/140‑1 (MQ=255)                  tGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCAtt                                                                                                       >  1:390717/1‑140 (MQ=255)                    tAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGt                                                                                                     >  2:617471/1‑140 (MQ=255)                           gTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGcc                                                                                              <  1:132771/140‑1 (MQ=255)                           gTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGcc                                                                                              <  1:112551/140‑1 (MQ=255)                           gTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGcc                                                                                              <  2:390717/140‑1 (MQ=255)                            tCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCt                                                                                             >  1:748860/1‑140 (MQ=255)                               aGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCg                                                                                          <  2:308609/140‑1 (MQ=255)                                  gcCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGa                                                                                       <  1:617471/140‑1 (MQ=255)                                    cTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATg                                                                                     <  1:432051/140‑1 (MQ=255)                                      cgccgcCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGt                                                                                   >  2:736373/1‑140 (MQ=255)                                      cgccgcCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGt                                                                                   <  1:619342/140‑1 (MQ=255)                                        ccgccgAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTaa                                                                                 <  2:492451/140‑1 (MQ=255)                                             gAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGa                                                                            <  2:774840/140‑1 (MQ=255)                                                 aCAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGatata                                                                        <  1:522118/140‑1 (MQ=255)                                                     tGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCta                                                                    >  1:283684/1‑140 (MQ=255)                                                                                 aaCGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGt                                        <  1:765391/140‑1 (MQ=255)                                                                                   cgcgTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTa                                      >  2:140604/1‑140 (MQ=255)                                                                                      gTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATc                                   >  2:477593/1‑140 (MQ=255)                                                                                             agtcagtTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATaa                            >  1:928320/1‑140 (MQ=255)                                                                                             agtcagtTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATaa                            >  2:722916/1‑140 (MQ=255)                                                                                             agtcagtTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATaa                            >  1:120026/1‑140 (MQ=255)                                                                                               tcagtTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAAt                          >  1:249104/1‑140 (MQ=255)                                                                                                      aCGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCg                   <  1:689996/140‑1 (MQ=255)                                                                                                        ggAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGc                        >  2:472971/1‑133 (MQ=255)                                                                                                        ggAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGc                        <  1:472971/133‑1 (MQ=255)                                                                                                               cAGATCATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGa          <  1:736373/140‑1 (MQ=255)                                                                                                                    cATTGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGa     <  1:140604/140‑1 (MQ=255)                                                                                                                      ttGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGata   <  2:249104/140‑1 (MQ=255)                                                                                                                      ttGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGata   <  2:283684/140‑1 (MQ=255)                                                                                                                       tGCGAACGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGatat  <  1:737482/140‑1 (MQ=255)                                                                                                                           aaCGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTa            >  2:521668/1‑126 (MQ=255)                                                                                                                           aaCGCGCTGCGGGTAGCGGGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTa            <  1:521668/126‑1 (MQ=255)                                                                                                                                            |                                                                                                                       TCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAGC‑GGGGGGTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATAT  >  NC_000913/3216337‑3216593

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