Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,106,950	A→G	N161S (AAC→AGC)	fieF →	Zn(2(+))/Fe(2(+))/Cd(2(+)) exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,106,950	0	A	G	100.0%	87.3 / NA	25	N161S (AAC→AGC)	fieF	Zn(2(+))/Fe(2(+))/Cd(2(+)) exporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/15);  total (10/15)

CGGGGTTATCGTGACAATTGTGGCGCTAATTTGTACGATTATCCTTGTCTCGTTTCAGCGTTGGGTGGTGCGGCGGACGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAACGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCTATATTTTATATAGCGCGTTACGCATGGGATATGAGGCGGTACAGT  >  NC_000913/4106816‑4107078                                                                                                                                       |                                                                                                                                 cGGGGTTATCGTGACAATTGTGGCGCTAATTTGTACGATTATCCTTGTCTCGTTTCAGCGTTGGGTGGTGCGGCGGACGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAgcggcg                                                                                                                             >  2:244864/1‑140 (MQ=255)               cAATTGTGGCGCTAATTTGTACGATTATCCTTGTCTCGTTTCAGCGTTGGGTGGTGCGGCGGACGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGgcg                                                                                                               <  2:74625/140‑1 (MQ=255)                aaTTGTGGCGCTAATTTGTACGATTATCCTTGTCTCGTTTCAGCGTTGGGTGGTGCGGCGGACGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGgcgc                                                                                                              <  1:244864/140‑1 (MQ=255)                          cTAATTTGTACGATTATCCTTGTCTCGTTTCAGCGTTGGGTGGTGCGGCGGACGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCAGCGCAATTCTGCTGGCGCTGGGGTTGTc                                                                                                    >  2:222579/1‑140 (MQ=255)                                                       tCAGCGTTGGGTGGTGCGGCGGACGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGc                                                                         <  2:204768/138‑1 (MQ=255)                                                       tCAGCGTTGGGTGGTGCGGCGGACGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGc                                                                         >  1:204768/1‑138 (MQ=255)                                                          gCGTTGGGTGGTGCGGCGGACGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGt                                                                    <  1:321492/140‑1 (MQ=255)                                                                  tggtgCGGCGGACGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGAATTg                                                            <  2:282928/140‑1 (MQ=255)                                                                     tgCGGCGGACGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGa                                                         >  1:190058/1‑140 (MQ=255)                                                                       cggcggACGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAAt                                                       >  1:62961/1‑140 (MQ=255)                                                                          cggACGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCgg                                                    >  1:42664/1‑140 (MQ=255)                                                                             aCGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCAt                                                 <  1:87976/140‑1 (MQ=255)                                                                             aCGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCAt                                                 <  2:27446/140‑1 (MQ=255)                                                                               gCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCt                                               >  1:384998/1‑140 (MQ=255)                                                                               gCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCt                                               <  1:158853/140‑1 (MQ=255)                                                                               gCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCt                                               <  2:190058/140‑1 (MQ=255)                                                                                                      gagacagaCATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCTATATTTTATATAGCGCGTTACGc                        <  1:276389/133‑1 (MQ=255)                                                                                                       tatGCTACATTACCAGTCTGATGTTATGATGAGCGGTGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCTATATTTTATATAGCGCGTTACGCa                       <  1:44914/140‑1 (MQ=255)                                                                                                        atGCTACATTACCAGTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCTATATTTTATATAGCGCGTTACGCAt                      <  2:62961/140‑1 (MQ=255)                                                                                                                      gTCTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCTATATTTTATATAGCGCGTTACGCATGGGATATGAGGCgg        >  2:236106/1‑140 (MQ=255)                                                                                                                        cTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCTATATTTTATATAGCGCGTTACGCATGGGATATGAGGCGGTa      >  2:220084/1‑140 (MQ=255)                                                                                                                        cTGATGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCTATATTTTATAGAGCGTCGTACGCATGGGATATGactgtctc      >  2:276389/1‑133 (MQ=255)                                                                                                                            tGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCTATATTTTATATAGCGCGTTACGCATGGGATATGAGGCGGTACAGt  <  2:336436/140‑1 (MQ=255)                                                                                                                            tGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCTATATTTTATATAGCGCGTTACGCATGGGATATGAGGCGGTACAGt  <  2:390495/140‑1 (MQ=255)                                                                                                                            tGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCTATATTTTATATAGCGCGTTACGCATGGGATATGAGGCGGTACAGt  <  2:42664/140‑1 (MQ=255)                                                                                                                            tGTTATGATGAGCGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCTATATTTTATATAGCGCGTTACGCATGGGATATGAGGCGGTACAGt  >  1:390495/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                 CGGGGTTATCGTGACAATTGTGGCGCTAATTTGTACGATTATCCTTGTCTCGTTTCAGCGTTGGGTGGTGCGGCGGACGCAAAGCCAGGCGGTGCGGGCTGATATGCTACATTACCAGTCTGATGTTATGATGAACGGCGCAATTCTGCTGGCGCTGGGGTTGTCCTGGTACGGCTGGCATCGCGCCGATGCTCTGTTTGCATTGGGAATCGGCATCTATATTTTATATAGCGCGTTACGCATGGGATATGAGGCGGTACAGT  >  NC_000913/4106816‑4107078

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: ‑