Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,505,777	T→C	I137V (ATC→GTC)	mntH ←	Nramp family divalent metal transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,505,777	0	T	C	100.0%	89.1 / NA	27	I137V (ATC→GTC)	mntH	Nramp family divalent metal transporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/14);  total (13/14)

CGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCG  >  NZ_CP009273/2505690‑2505854                                                                                        |                                                                              gggCAACAAACAACAGTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  1:83141/89‑1 (MQ=255) cgGCAACAACCAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  1:1182727/90‑1 (MQ=255) cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  2:444298/90‑1 (MQ=255) cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             >  2:216146/1‑90 (MQ=255) cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  1:855002/90‑1 (MQ=255) cgGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  1:487853/90‑1 (MQ=255) cgGCAAAAAACAACAGTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  1:1085649/90‑1 (MQ=255) cgGAAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  2:675079/90‑1 (MQ=255)      acaaacaaCAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTc                                                                        <  1:322071/90‑1 (MQ=255)           caacaGTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAc                                                                   <  1:925200/90‑1 (MQ=255)           caacaGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCAc                                                                   <  2:150315/90‑1 (MQ=255)              caGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACcgc                                                                <  2:136747/90‑1 (MQ=255)                gTACCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCgcgc                                                              <  1:123783/90‑1 (MQ=255)                gTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCgcgc                                                              <  1:257849/90‑1 (MQ=255)                gTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCgcgc                                                              >  1:1172311/1‑90 (MQ=255)                 tAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGcc                                                             >  1:329110/1‑90 (MQ=255)                  aaCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGccc                                                            <  2:1038660/90‑1 (MQ=255)                  aaCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGccc                                                            <  1:216146/90‑1 (MQ=255)                   aCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             >  2:1011834/1‑72 (MQ=255)                   aCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  1:1011834/72‑1 (MQ=255)                   aCCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCt                                                           >  2:609279/1‑90 (MQ=255)                     ccGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGc                                                         >  2:978931/1‑90 (MQ=255)                     ccGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGc                                                         >  1:560228/1‑90 (MQ=255)                     ccGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGc                                                         >  2:147215/1‑90 (MQ=255)                       gCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCgcgc                                                              >  1:95367/1‑83 (MQ=255)                       gCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCgcgc                                                              <  2:95367/83‑1 (MQ=255)                       gCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGcaa                                                       >  1:401272/1‑90 (MQ=255)                              aCTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             <  1:1215057/61‑1 (MQ=255)                              aCTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGAccc                                                                             >  2:1215057/1‑61 (MQ=255)                              aCTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGaa                                                >  1:724657/1‑90 (MQ=255)                                          ggTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATg                                    <  1:204357/90‑1 (MQ=255)                                          ggTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATg                                    <  2:789401/90‑1 (MQ=255)                                          ggTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATg                                    <  1:837209/90‑1 (MQ=255)                                          ggTTTTCGCCCGCGACGTTGCAGCCTTAAAATCAGGAAAGTAGAGACCCCCGTCAGCCCCGCGCCCTGCCACAACGAAACACCAAGAATg                                    <  1:701100/90‑1 (MQ=255)                                                     gcgACGTTGCAGCATTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAAcc                         >  2:999059/1‑90 (MQ=255)                                                                  aTTAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCa            >  1:881130/1‑90 (MQ=255)                                                                    tAAAATCAGGAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAAt          <  2:554207/90‑1 (MQ=255)                                                                            ggAAAGTCGCGACCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCg  >  1:834427/1‑90 (MQ=255)                                                                            ggAAAGTCGCGACCCCCGCCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCg  >  2:698632/1‑90 (MQ=255)                                                                                        |                                                                              CGGCAACAAACAACAGTAACCCGCCAATCACTTTCTCCAGCGGTTTTTGCCCGCGACGTTGCAGCATTAAAATCAGGAAAGTCGCGATCCCCGTCAGCACCGCGCCCTGCAACAACGAAACACCAAGAATGAGTTTAAAACCGATCGCCGCACCAATAAATTCCG  >  NZ_CP009273/2505690‑2505854

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