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

CGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTG  >  NC_000913/3387816‑3388272                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                         cGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTc                                                                                                                                                                                                                 >  2:80539/1‑251 (MQ=255)                   aCGAGTTTTCTATTCGTACCAGCGAGGTTACCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAAGGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGc                                                                                                                                                                                               <  2:267538/251‑1 (MQ=255)                            cTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGCTGCGCAACAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCCACAATCCCCCCTGTCAGTCCCGCCCAGCGGGGCGTTTCCAGc                                                                                                                                                                                      >  2:104002/1‑251 (MQ=255)                                    aCCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTg                                                                                                                                                                              <  1:327374/251‑1 (MQ=255)                                                      gTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAgcgc                                                                                                                                                            >  2:86144/1‑251 (MQ=255)                                                                    ccAGATACAGCACACCAGAATCATCAAAAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATcgccg                                                                                                                                                                                                                      <  2:325829/179‑1 (MQ=255)                                                                                         catcaATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAg                                                                                                                         <  1:327191/250‑1 (MQ=255)                                                                                         catcaATAATGGTGCGCGGATCATCGCAAGGATGATCACCGGTCCGGAAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGc                                                                                                                        <  2:107038/251‑1 (MQ=255)                                                                                                    gggcgcGGATCAGCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGAGGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGCGAATATGTTgg                                                                                                             <  2:74994/249‑1 (MQ=255)                                                                                                                tcGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCTCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGga                                                                                                  <  1:254128/250‑1 (MQ=255)                                                                                                                    aaggctgatCACCAGTCCGGCAATACAGCCAATAAATGGGCCGATGATTCGCAAAAAGCCCCTATAGAGAATAGCGCCAGAATACTGTTCACCTCCCGCAGCAAAGGTCGGACCGGCGGCAACTAGCGCCGCTGTCAGTACCGCCCAGCGTTGCGTTTCCAGCGTGAAGTGAAAGCCAACAAACAGCGCCAGTACAAGGGCGGTCTCC‑AGTTTTACCGCAAAGCGAATATTTTGGTTAGCAATGGAGaaaa                                                                                             <  2:167337/246‑1 (MQ=255)                                                                                                                     atgatgATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGCTGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTAACGCCAAACGCATATGtttggttggaatgggggaaa                                                                                             >  1:90764/1‑232 (MQ=255)                                                                                                                            tCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCC‑AGTTTTACCGCAACACGCATATGTTGGTTTGCAAAGGAGAAAATAACCacc                                                                                     >  1:100350/1‑249 (MQ=255)                                                                                                                             cACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATc                                                                                     <  1:86144/250‑1 (MQ=255)                                                                                                                                      ggCAATACAGCCAATAAATGTGCCGCTGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTCCCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTTCCCCCAAGCGAATATGTTTGTTTGCAATTGCGAAAATAACCAACGTGATTTAcc                                                                           >  1:344031/1‑251 (MQ=255)                                                                                                                                              aGCCAATAAATGTGCCGATTAGGCGCCAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGAGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCTCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCAc                                                                   <  2:262540/251‑1 (MQ=255)                                                                                                                                                   aTAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGGCGCCCGGTTTTAC‑GCACAACGGATATGTTGGTTTGCAATGTGGAAAAAAACCACCGTGGTTAACCCAACTCACCCAgg                                                              >  2:55655/1‑251 (MQ=255)                                                                                                                                                              gatgatGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCAt                                                   <  1:265102/251‑1 (MQ=255)                                                                                                                                                                         aaaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGgaaga                                        <  1:55655/251‑1 (MQ=255)                                                                                                                                                                                 cGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGCGAATATGGTGGTTTGCCATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGGGATGGGGCATTTTACGGGAGGACGaaacc                                >  1:181913/1‑248 (MQ=255)                                                                                                                                                                                            aGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAAACAAACTCACGCAGGCGATGGGGCATTTTACGGGAGAAAGAATCCCGGCTTtcgtc                     >  1:323245/1‑251 (MQ=255)                                                                                                                                                                                                             cACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAACGAATATGTTGGTTAGCAATGGGGAAAAAACCCATCGGGATTAACCAAACTCACGCCGGCGAAGGGCCATTTTTAGGGAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGcc      >  1:107424/1‑249 (MQ=255)                                                                                                                                                                                                             cACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGGTGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATTCCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAATAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCa     <  2:198391/250‑1 (MQ=255)                                                                                                                                                                                                               ccTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGCGAATATGTTGGTTa                                                                                                          <  2:197099/147‑1 (MQ=255)                                                                                                                                                                                                               ccTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGCGAATATGTTGGTTa                                                                                                          >  1:197099/1‑147 (MQ=255)                                                                                                                                                                                                                cTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGGTTAACCAAACTCACCCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTTTTTGCCAGTg  >  1:7369/1‑250 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                         CGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCC‑AGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTG  >  NC_000913/3387816‑3388272

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