Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,568,337	A→G	G61G (GGT→GGC)	mdtM ←	multidrug efflux pump/bile salt:H(+) antiporter/Na(+):H(+) antiporter/K(+):H(+) antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,568,337	0	A	G	100.0%	118.5 / NA	34	G61G (GGT→GGC)	mdtM	multidrug efflux pump/bile salt:H(+) antiporter/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/23);  total (11/23)

GCGATAAGAAACTGTGTCATAGACGTTGTGAACATTGTCGCGGCGCAGGCAAGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATACCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCA  >  NC_000913/4568200‑4568451                                                                                                                                          |                                                                                                                   gcgATAAGAAACTGTGTCATAGACGTTGTGAACATTGTCGCGGCGCAGGCAAGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGcc                                                                                                                  >  2:89590/1‑140 (MQ=255)   gATAAGAAACTGTGTCGTAGACATTGTGAACTTAGCCGCGGCGCAGGCAAGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCGCTGTAACGCCATGCCGc                                                                                                                <  2:299762/140‑1 (MQ=255)      aaGAAACTGTGTCATAGACGTTGTGAACATTGTCGCGGCGCAGGCAAGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAg                                                                                                             <  1:135849/140‑1 (MQ=255)               tgtCATAGACGTTGTGAACATTGTCGCGGCGCAGGCAAGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATaga                                                                                                    >  1:31931/1‑140 (MQ=255)               tgtCATAGACGTTGTGAACATTGTCGCGGCGCAGGCAAGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATaga                                                                                                    >  2:45302/1‑140 (MQ=255)               tgtCATAGACGTTGTGAACATTGTCGCGGCGCAGGCAAGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATaga                                                                                                    >  1:487901/1‑140 (MQ=255)                      gACGTTGTGAACATTGTCGCGGCGCAGGCAAGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGAc                                                                                             >  2:482977/1‑140 (MQ=255)                           tgtgAACATTGTCGCGGCGCAGGCAAGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGgcag                                                                                        >  1:429151/1‑140 (MQ=255)                             tgAACATTGTCGCGGCGCAGGCAAGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGgcagca                                                                                      <  1:581799/140‑1 (MQ=255)                                  aTTGTCGCGGCGCAGGCAAGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGc                                                                                 <  2:55956/140‑1 (MQ=255)                                       cgcgGCGCAGGCAAGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCGCGATAGAGCCTGAGCGCAGCAGAGGCCAGAc                                                                            >  2:321022/1‑140 (MQ=255)                                         cggcgCAGGCAAGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTg                                                                          >  2:353714/1‑140 (MQ=255)                                                   aaGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCAt                                                                <  2:31931/140‑1 (MQ=255)                                                           aaaTTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCa                                                        <  1:482977/140‑1 (MQ=255)                                                           aaaTTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCa                                                        >  2:166235/1‑140 (MQ=255)                                                             aTTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACg                                                      <  1:353714/140‑1 (MQ=255)                                                               tAGCGCCCCGGTAATCAGCACCTGCCTGCGGCCAATTCTGTCGGAAAGTGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTa                                                    <  2:282959/140‑1 (MQ=255)                                                               tAGCGCCCCGGTAATCAGCACCGGCCTGCGGTCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGCCCAGACTGACATCGGCATTAAAATCACGTa                                                    <  2:413673/140‑1 (MQ=255)                                                               tAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTa                                                    <  1:87004/140‑1 (MQ=255)                                                                aGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTAc                                                   <  1:1087/140‑1 (MQ=255)                                                                aGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTAc                                                   <  2:347731/140‑1 (MQ=255)                                                                aGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTAc                                                   <  1:166235/140‑1 (MQ=255)                                                                              cAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATcc                                     <  2:13983/140‑1 (MQ=255)                                                                              cAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATcc                                     <  2:562726/140‑1 (MQ=255)                                                                              cAGCACCGGCCTGCGGCAGATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATcc                                     <  2:662871/140‑1 (MQ=255)                                                                                  aCCGGCCTGCGGCCAATTCTGTCGGAAAGCGTCCCAAGCAGCCACTGTTACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAgg                                 <  2:256909/140‑1 (MQ=255)                                                                                   ccGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGTCATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGc                                >  2:617066/1‑140 (MQ=255)                                                                                   ccGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGGTCCCAGGc                                <  1:45302/140‑1 (MQ=255)                                                                                   ccGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGc                                <  2:487901/140‑1 (MQ=255)                                                                                      gCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCCgg                             <  2:429151/140‑1 (MQ=255)                                                                                        cTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGAt                           <  1:614306/140‑1 (MQ=255)                                                                                           cGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCCGGATCAg                        >  1:62490/1‑140 (MQ=255)                                                                                           cGGCCAATTCTGTCGGAAAGCGG‑CCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCAAATAAATGATCCCAGGCTGGATCAGa                       >  1:325196/1‑140 (MQ=255)                                                                                              ccAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATc                     <  1:617066/140‑1 (MQ=255)                                                                                                                 gCCCCAGCAGCCACTGTAACGCCATGCCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACgcagca  <  1:379868/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                   GCGATAAGAAACTGTGTCATAGACGTTGTGAACATTGTCGCGGCGCAGGCAAGGGTAAAAATTAGCGCCCCGGTAATCAGCACCGGCCTGCGGCCAATTCTGTCGGAAAGCGGCCCCAGCAGCCACTGTAACGCCATACCGCCAGCAAGATAGAGACTGACGGCAGCAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCA  >  NC_000913/4568200‑4568451

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