Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,388,726	A→G	V134A (GTT→GCT)	aaeA ←	aromatic carboxylic acid efflux pump membrane fusion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,388,726	0	A	G	100.0%	121.8 / NA	35	V134A (GTT→GCT)	aaeA	aromatic carboxylic acid efflux pump membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/21);  total (14/21)

TAATAAACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAACCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCAGTACCTGATAATAAGCAA  >  NC_000913/3388598‑3388846                                                                                                                                 |                                                                                                                         taataaACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGt                                                                                                                      <  1:601583/133‑1 (MQ=255) taataaACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGt                                                                                                                      <  1:478575/133‑1 (MQ=255)   ataaACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGttt                                                                                                                    >  1:443038/1‑133 (MQ=255)        cTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTgtagt                                                                                                               <  1:608311/133‑1 (MQ=255)         tCACCGGTCTAGACGATGTGGTCGGTCACCTAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTgtagta                                                                                                              <  1:185862/133‑1 (MQ=255)         tCACCGGTATCGACGCTGAGGTTGCTCAACCAGCCATCTGCCGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTgtagta                                                                                                              <  2:178864/133‑1 (MQ=255)         tCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTgtagta                                                                                                              <  2:487452/133‑1 (MQ=255)         tCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTgtagta                                                                                                              <  2:483073/133‑1 (MQ=255)         tCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTgtagta                                                                                                              <  1:370192/133‑1 (MQ=255)         tCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTgtagta                                                                                                              <  1:443418/133‑1 (MQ=255)           aCCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACg                                                                                                            >  2:145530/1‑133 (MQ=255)           aCCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTATCGGTAATTGATGCAGAGCCGTTTGTAGTACg                                                                                                            >  2:185862/1‑133 (MQ=255)              ggTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACgttg                                                                                                         <  2:319458/133‑1 (MQ=255)                   aGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGc                                                                                                    <  1:579034/133‑1 (MQ=255)                         tGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTc                                                                                              >  1:608623/1‑133 (MQ=255)                         tGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTc                                                                                              >  1:121289/1‑133 (MQ=255)                          gAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTGACTGATGCTGAGCCTTTTGTAGTACGTTGTTCGCCTGGttg                                                                                             >  2:299556/1‑131 (MQ=255)                          gAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCg                                                                                             >  1:36092/1‑133 (MQ=255)                            ggTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGTCGTTTGTGGTAAGTGGTTGGCGTGGTCGAt                                                                                           >  1:174654/1‑133 (MQ=255)                            ggTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGAt                                                                                           <  1:145530/133‑1 (MQ=255)                                           cATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACAt                                                                            <  2:370956/133‑1 (MQ=255)                                               tgctgGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGcc                                                                        <  1:409063/133‑1 (MQ=255)                                                gctgGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTATTACGTTGTAGGCCTGGAAGATCTCTACGCGATACATCGcgt                                                                       >  2:296307/1‑131 (MQ=255)                                                       cgcgGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACAcc                                                                >  1:6732/1‑133 (MQ=255)                                                                      gTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTacgacg                                                 >  1:23405/1‑133 (MQ=255)                                                                          aaGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCg                                             <  2:608623/133‑1 (MQ=255)                                                                                cAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGcctcc                                       >  1:502410/1‑133 (MQ=255)                                                                                        ccAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGttt                               >  1:625065/1‑133 (MQ=255)                                                                                            aTCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTcc                           <  2:23405/133‑1 (MQ=255)                                                                                            aTCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTcc                           <  2:443038/133‑1 (MQ=255)                                                                                               gcgGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCtgt                        <  1:60622/133‑1 (MQ=255)                                                                                                 gggTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGc                      <  1:148425/133‑1 (MQ=255)                                                                                                   gTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCa                    <  2:121289/133‑1 (MQ=255)                                                                                                                   cTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCAGTACCTGATAATAAGc    <  2:625065/133‑1 (MQ=255)                                                                                                                     aaCTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCAGTACCTGATAATAAGCaa  >  2:350853/1‑133 (MQ=255)                                                                                                                                 |                                                                                                                         TAATAAACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAACCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCAGTACCTGATAATAAGCAA  >  NC_000913/3388598‑3388846

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