Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,216,313	T→C	Q60R (CAA→CGA)	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,313	0	T	C	100.0%	69.5 / NA	21	Q60R (CAA→CGA)	nfeF	NADPH‑dependent ferric chelate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/8);  total (13/8)

CATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGG  >  NC_000913/3216200‑3216441                                                                                                                  |                                                                                                                                 cATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCt                                                                                                               <  2:430268/133‑1 (MQ=255)         aGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGt                                                                                                       >  1:385487/1‑133 (MQ=255)           ttCGTCATACAGCGGCGTATAGTCACGCGCCGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGACAGGTCGTGGAAAGAAGAGTTTGCTGTGATCGTCa                                                                                                     >  2:480393/1‑133 (MQ=255)                cATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGGGTTTGCTGTGATCTTCAAAGcc                                                                                                >  1:175795/1‑133 (MQ=255)                   aCAGCGGCATATAGTTACCCGACGGTGGGCCTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCAcg                                                                                             <  1:144500/133‑1 (MQ=255)                        ggcgTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGATGTGATCGTCAAAGCCACGCGATg                                                                                        <  2:299652/133‑1 (MQ=255)                         gcgTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGt                                                                                       >  1:490607/1‑133 (MQ=255)                          cgTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTa                                                                                      >  2:352085/1‑133 (MQ=255)                                 tCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGt                                                                               >  2:94920/1‑133 (MQ=255)                                      cgACGGTGGGCGTGGTCCTTCCGGGCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAgc                                                                          >  2:347051/1‑133 (MQ=255)                                            tGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTcg                                                                    <  2:490607/133‑1 (MQ=255)                                                    gTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAg                                                            >  1:269438/1‑133 (MQ=255)                                                                       gCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTGGCCTCCGAGGACGATGCCCTGAAAACCg                                         >  2:182092/1‑133 (MQ=255)                                                                       gCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCg                                         >  2:443671/1‑133 (MQ=255)                                                                            ctTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCt                                    >  2:1020/1‑133 (MQ=255)                                                                                cGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATg                                >  2:402080/1‑133 (MQ=255)                                                                                    aCCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCt                            <  1:352085/133‑1 (MQ=255)                                                                                    aCCGTTGGCGGCACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCt                            >  1:390140/1‑133 (MQ=255)                                                                                      cGTTGGCGGCGCAAATTGAGCGTCGGGTTGAGGAAAGAAGACTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCa                          <  2:219646/133‑1 (MQ=255)                                                                                          ggcggcACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAAcgc                      <  2:385487/133‑1 (MQ=255)                                                                                          ggcggcACAAAGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAAcgc                      >  2:90131/1‑133 (MQ=255)                                                                                                    aGTGAGCGTCAGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACAcagt            <  2:290404/133‑1 (MQ=255)                                                                                                              aGGTCGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACgg  <  2:390140/133‑1 (MQ=255)                                                                                                                  |                                                                                                                                 CATGGCGTAGTTCGTCATACAGCGGCGTATAGTCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATCGTCAAAGCCACGCGATGTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGG  >  NC_000913/3216200‑3216441

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