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%	32.3 / NA	16	T50P (ACG→CCG)	aaeB	aromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (8/8);  total (8/8)

GATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTT  >  NC_000913/3387799‑3388265                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 gATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACAGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGAAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACcggcgg                                                                                                                                                                                                                          >  1:200581/1‑251 (MQ=255)                    gCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAACAACCTGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTc                                                                                                                                                                                                         >  1:254846/1‑249 (MQ=255)                    gCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCTCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAAGGGGCCAATGATACGCAAAAAGCCACGATAGCGAAGAGCGCCAGAATACGGTGCACCTCCCGCAGCAAAGGCCGGAcc                                                                                                                                                                                                                               <  1:287585/227‑1 (MQ=255)                    acaacgccccACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAACCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCTGGACCGGCGGCAACAATCGCCGCTGTc                                                                                                                                                                                                         <  1:94873/243‑1 (MQ=255)                        gCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCAAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTAcc                                                                                                                                                                                                   >  1:240137/1‑251 (MQ=255)                                                                                    cccAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCg                                                                                                                                       >  1:271729/1‑251 (MQ=255)                                                                                            cAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCCATACAGCCCATAAATGTGCCGATGATGCGCAAACAGCAACGCTCGCGAATAGCGCCAGAATACGATTCACCTCCCGCCGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGt                                                                                                                                                                                                          >  1:51556/1‑176 (MQ=255)                                                                                                           atcaATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCAAATAAATGTGCCGCTGATGCGCAACAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATcgccgc                                                                                                                                                                                                             >  1:205648/1‑158 (MQ=255)                                                                                                                            gATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAgc                                                                                                                                                      <  1:11662/196‑1 (MQ=255)                                                                                                                                  ngnaannaagtncACCATTCCGCCAATACAGCCAAANAATGTGCCCATGGTGCGCAAAAAGCCACGAANGNGAATAGCGCCAAAATACGTTTCACCTCCCACAGCAAAGGCCGGACCGGNGGCAACAATCNCCGNNGACANTNCNGCCCATCGTGGCGTTTCCAGCTGGAAGTGNAAGCCAACAAACAGCGCCAGTACNNTGNNNGTNNCCAgc                                                                                                                              <  2:254846/202‑2 (MQ=255)                                                                                                                                                                          gtgCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATgn                                                 <  1:72612/251‑2 (MQ=255)                                                                                                                                                                             ccGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCg                                                                                                                                                                                      <  1:59834/115‑1 (MQ=255)                                                                                                                                                                                 tgaggcgcAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGGGGCGTTTCCAGCTGGAAGTGAAAgn                                                                                                                                                                 <  1:85092/128‑2 (MQ=255)                                                                                                                                                                                                                      aTACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTg     <  1:115369/251‑1 (MQ=255)                                                                                                                                                                                                                         cGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGAGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGGTGGGCCATTTTACGGACGAACGAATCCTGGCTTTCGTCGCGATCTTGttt  >  1:183388/1‑251 (MQ=255)                                                                                                                                                                                                                                          caAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGGTGGAAGGGAAAGCCAAAAAACAGCGCCAGGACAATGGCGGTCGCCAGGTTTACCGCAAAGCGCATATGTTGGTTAGCa                                                                                                >  1:324311/1‑140 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 GATCAGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTT  >  NC_000913/3387799‑3388265

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑