Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,388,041	T→G	T50P (ACG→CCG)	aaeB ←	aromatic carboxylic acid efflux pump subunit AaeB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,388,041	0	T	G	100.0%	78.7 / NA	25	T50P (ACG→CCG)	aaeB	aromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (18/7);  total (18/7)

GAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAA  >  NC_000913/3387857‑3388239                                                                                                                                                                                         |                                                                                                                                                                                                       gAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACTCCAGAAGCATCAATAATGGTGCGCGGATCATCGCAATGATTATCACCAGTCCGGCAATACAGCCAATAAATGTGCCTATGATGCGCAAAAATCCACGATAGCGAATAGCGCCAGAAGACGGTTCACCGCCCGCAGCAAAGGCCGGACCTGCGGCAACAATCTCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGcc                                                                                                                                      <  1:383599/251‑1 (MQ=255)            ggTACAAAAACCGGCCCAGATGCAGCACACCAGAATCATCCATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCACTACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAgcg                                                                                                                           >  1:487935/1‑251 (MQ=255)             gTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATcgc                                                                                                                                                                                      <  1:22990/191‑1 (MQ=255)             gTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATAGCACTGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGCTGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAgcgc                                                                                                                          >  1:80933/1‑251 (MQ=255)                         gCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGCTGATCACCAGTCCGGCAATACAGCCAATAAATGTGACGATGATGCGCAACAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGc                                                                                                              >  1:484240/1‑251 (MQ=255)                                                   caATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGa                                                                                     >  1:20422/1‑250 (MQ=255)                                                                  gATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGCTGCGCAAACAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAAtag                                                                                                                                                                                       >  1:491662/1‑135 (MQ=255)                                                                  gATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGcaaca                                                                                                                                                                                           >  1:230658/1‑133 (MQ=255)                                                                                      ccAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTg                                                 >  1:112444/1‑251 (MQ=255)                                                                                                   aCAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGANNGGCGGcnnc                                                                                                                                                                                            >  2:7112/1‑99 (MQ=255)                                                                                                    cAGCCTATAAATGTGCCGATCATGCGCAAAAATCCACGATAGCGAATAGCGCCAGCATACGGTTCACCTCCCGCAGCAAAGGCCGGACNNGCGGcanna                                                                                                                                                                                           >  2:35537/1‑99 (MQ=255)                                                                                                        cAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCGCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTgn                                                                                                                                                              <  1:61336/124‑2 (MQ=255)                                                                                                                   ccGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTa                                                                         >  1:226440/1‑198 (MQ=255)                                                                                                                     gatgatGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTAAGTACCGCCCAGCGTGGCGTTTCCAGCTGGCAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGCTTAACCAAACTCACGCAGGCGATGGGCCAt                  >  1:290006/1‑251 (MQ=255)                                                                                                                        gatgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCAtttt               <  1:7112/251‑1 (MQ=255)                                                                                                                             gcAAAAAGCCACGATAGCGAATAGCGCCAGACTACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCAGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAAACAAACTCACGCAGGCGATGGGCCATTTTACgg           >  1:209586/1‑250 (MQ=255)                                                                                                                              cgaaaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGGGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGga          >  1:33897/3‑250 (MQ=255)                                                                                                                                      cACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCTAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCCTCGTGATTAACCAAACTCACGCAGGCGATGGGCCCTTTTACGGAAGAACGaa  >  1:488727/1‑250 (MQ=255)                                                                                                                                      cACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGACCGGACCGGAGGCAACCCTCGCCGCTGTCAGTCACGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCTCGTTTTACCGCAAAACGAATCTGTTGGTTAGCAATGGAGAAAATACCCATCGATATTAACCAAACTCACGCAGGCGATGGGCCCTTTTACGGCAGAACGaa  >  1:488706/1‑250 (MQ=255)                                                                                                                                            aNNGAATAGNNCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCt                                                                                                                                                   <  2:9086/99‑1 (MQ=255)                                                                                                                                                       ccAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACg    <  1:21123/231‑1 (MQ=255)                                                                                                                                                         aGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGATGTAAGTACCGCCCAGCGTGGCGTTTCAAGCTGGACGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTAACGCAGGCGATGGGCCCTTTTACGGAAGAACg    >  1:284139/1‑229 (MQ=255)                                                                                                                                                                      cctcctGCAGCAAAGGCCGGACCGGCGGCAACACTCGCCGCTGTCAGTAACGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGc                                                                                                              >  1:402634/1‑110 (MQ=255)                                                                                                                                                                         cccGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCgt                                                                                                        <  1:211420/113‑2 (MQ=255)                                                                                                                                                                           ngcagcaAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCAAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAAAAGCGCCAg                                                                                                                       >  1:99556/2‑96 (MQ=255)                                                                                                                                                                                         |                                                                                                                                                                                                       GAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAA  >  NC_000913/3387857‑3388239

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