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

AGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAA‑AGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAA‑CC‑AAACTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCA  >  NC_000913/3387803‑3388279                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   aGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGcaac                                                                                                                                                                                                                                        >  1:158928/1‑251 (MQ=255)         ataACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATcgcc                                                                                                                                                                                                                                 <  1:66344/250‑1 (MQ=255)                              aTACGAGTTTTCTATTCGTACAAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCggg                                                                                                                                                                                                           >  2:263845/1‑249 (MQ=255)                                aCGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGc                                                                                                                                                                                                         <  2:32722/251‑1 (MQ=255)                                            ttCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCCGCTgg                                                                                                                                                                                             >  2:64512/1‑251 (MQ=255)                                            ttCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTgg                                                                                                                                                                                             >  2:309408/1‑251 (MQ=255)                                                         ggATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCca                                                                                                                                                                                 >  1:235888/1‑250 (MQ=255)                                                                           aCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACaa                                                                                                                                                              <  1:118058/251‑1 (MQ=255)                                                                                        cAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCa                                                                                                                                                                                                  >  2:77379/1‑202 (MQ=255)                                                                                        cAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCa                                                                                                                                                                                                  <  1:77379/202‑1 (MQ=255)                                                                                                   aatcagaaataatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCaa                                                                                                                                      <  2:76594/244‑1 (MQ=255)                                                                                                    atcatcAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCaa‑a                                                                                                                                    <  2:133709/251‑1 (MQ=255)                                                                                                      catcaATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAA‑Ag                                                                                                                                   <  2:128064/250‑1 (MQ=255)                                                                                                                      cgGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAA‑AGCGAATATGTTGGTTAGc                                                                                                                  <  1:227833/251‑1 (MQ=255)                                                                                                                                               gTCCGGCAATAGAGCCAATAAAGGTGACGATGCTGCGCAAAAAGCCACGATAGCGAATAGAGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGAGTTTCCAGCTGGAAGTGAACGCCAACAAACAGCGCCCGTACAATGGCGGAACCCAGTTTTACCGCAACACCGA‑TATGTTGGTTAACAATGGGGAACATaacccaacgtatt                                                                                         >  2:197161/1‑239 (MQ=255)                                                                                                                                                 ccGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAA‑AGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTa                                                                                        <  1:64512/250‑1 (MQ=255)                                                                                                                                                   ggCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAA‑AGCGAATATGTTTGTTAGCAATGGAGAAAATACCCCTCGTGATTAA‑cc                                                                                    >  2:69231/1‑251 (MQ=255)                                                                                                                                                            gCCAATCAATGTGCAGATGATGCGCAAAAAGCCCCGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGCCGCCAGGTTTTCCCCCA‑ACAGAATATGTTGGGTAGCCATGGGGAAAATACCCCTCGTGATTAA‑CC‑AAAATCAcc                                                                          >  2:129984/1‑250 (MQ=255)                                                                                                                                                               aaTAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGGTTTACCGCAA‑AACGAAAATGGTGGTTAGCAATGGGGGAAAAAACCATCGTTATTAAACC‑AAA‑TCACGc                                                                         >  2:238465/1‑249 (MQ=255)                                                                                                                                                                 aaaaTGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAA‑AGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAA‑CC‑AAACTCACGCAGGc                                                                     <  2:77705/250‑1 (MQ=255)                                                                                                                                                                   aaTGTGCCGGTGATGCGCAACAAGACACGATAGCGAATAGCGCCAGAATACGGTTCCCCTACCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTgca                                                                                                                                                                                      >  2:317776/1‑137 (MQ=255)                                                                                                                                                                   aaTGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCATAATACGGGTCACCTCCCGCAGCAAAGGCCGGACGGGCGGCAACAATCTCCGCTGTCAGTGGCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGaa                                                                                                                                                                                      <  1:317776/139‑1 (MQ=255)                                                                                                                                                                   aaTGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAA‑AGCGAATATGTTGGGTAGCAATGGAGAAAATACCCATCGTTGTTAA‑CC‑AAACTCACGCCGGCg                                                                    >  1:83955/1‑250 (MQ=255)                                                                                                                                                                        gCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAA‑AGCTAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAA‑CC‑AAACTCACGCAGGCGATGGG‑c                                                             <  1:263845/251‑1 (MQ=255)                                                                                                                                                                                     aaaaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAA‑AGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAA‑CC‑AAACTCACGCAGGCGATGGG‑CCATTTTACGgaag                                                <  1:285091/251‑1 (MQ=255)                                                                                                                                                                                         aGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAA‑AGCGGATATGGTTGTTAGCAATGGGGAAAATACCCATCGTGATTAA‑AC‑AAACTCACGCAGGCGATGGG‑CCATTTTACGgaagaaac                                            >  1:88608/1‑249 (MQ=255)                                                                                                                                                                                           ccACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAA‑AGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAA‑CC‑AAACTCACGCAGGCGGTGGG‑CCATTTTACGGAAGAACGa                                           >  1:209683/1‑250 (MQ=255)                                                                                                                                                                                               gATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCCCCCGGTTTACCGCCA‑AGCGAAAATGTTGGGTAGCCATGGAGAAAATACCCATCGTGATTAA‑CC‑AAACTCACGCAGGCGCTGGG‑GCCTTTTACGgaaaaaagactcct                                      >  2:140079/1‑240 (MQ=255)                                                                                                                                                                                                   gCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAA‑AGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAA‑CC‑AAACTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCt                                  <  1:69231/251‑1 (MQ=255)                                                                                                                                                                                                     gAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAACGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCCGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCCCAC‑AGCGAATATGTTGGTTAGCAATGGAGGAAATACCCAACGTGATTAA‑CCCAAAC‑CAACCAGGCGAGTGG‑CCCTTTTACgggcgaactacatccgtcgtt                                >  2:58472/1‑232 (MQ=255)                                                                                                                                                                                                                  aTACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGGTTTTCCGCAA‑AACGAATATGTTTGTTAGACATGGAGAAAATACCCCTCGTGATTAAACC‑AAAC‑CACGCAGGCGAGGGG‑GCCTATTACGGAAGAAACAATCCAGGCTTtcttcgagtgcttg                   >  2:317275/1‑243 (MQ=255)                                                                                                                                                                                                                     cGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAA‑AGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAA‑CC‑AAACTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGttt                <  1:308112/251‑1 (MQ=255)                                                                                                                                                                                                                           aCCTCCCGCAGCAAAGGCCGGACCGGCGGCAACCCTCGACGCTGTCAGTACCGCCCCGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACCATGGCGGTCGCCAGCTTTACCGCAA‑AGCGAATATGTTGGTTAGCCATGGAGAAAATACCCCTCGCGATTAA‑CC‑AAACTCACGCAGGCGATGGG‑CCCTTTTACGGCACCACGAATCCTGGCCTTCCTCGCGATCTTGTTTGCCAGt          >  1:141712/1‑251 (MQ=255)                                                                                                                                                                                                                                ccGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGGTTTACCCCAA‑AGCGAATATGTTGGGTAGCAATGGGGAAAAAACCCCTCGTGATTAA‑CC‑CAACTCACCCAGGCGATGGGGCCTTTTTA‑GGAAGAACGAATCCCGGCTTTCGTCGCGATCTTGTTTGCCAGTTaccc     >  2:128250/1‑250 (MQ=255)                                                                                                                                                                                                                                   cagcaAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGGTTTTCCGCAA‑AGCGAATATGTTGGTTTGCCATGGAGGAAATACCCACCGTGATTAA‑CC‑AAACTCACCCAGGCGATGGG‑CCCTTTTACGGAACAACGAAACCTGGGTTTCGTCGCGATCtttttttcccgtggcccccc  >  2:228340/1‑237 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   AGCGCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAA‑AGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAA‑CC‑AAACTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCA  >  NC_000913/3387803‑3388279

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