Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,464,489	G→A	*317* (TAG→TAA)  D7N (GAC→AAC)	paaX → paaY →	transcriptional repressor of phenylacetic acid degradation paa operon, phenylacetyl‑CoA inducer thioesterase required for phenylacetic acid degradation, trimeric, phenylacetate regulatory and detoxification protein, hexapeptide repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,464,489	0	G	A	100.0%	154.9 / NA	50	*317* (TAG→TAA)  D7N (GAC→AAC)	paaX paaY	transcriptional repressor of phenylacetic acid degradation paa operon, phenylacetyl‑CoA inducer thioesterase required for phenylacetic acid degradation, trimeric, phenylacetate regulatory and detoxification protein, hexapeptide repeat protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/27);  total (23/27)

CTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAGACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGCCAAATGCCAGCCTGCGTGGCGATTTTGGTCGTATC  >  NC_000913/1464356‑1464620                                                                                                                                      |                                                                                                                                    ctgctgCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGt                                                                                                                                <  2:252185/139‑1 (MQ=255) ctgctgCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGt                                                                                                                                <  2:363815/139‑1 (MQ=255)  tgctgcTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTc                                                                                                                               >  1:225670/1‑139 (MQ=255)              gttcgttAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTtg                                                                                                                   >  1:235102/1‑139 (MQ=255)               ttcgttAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGt                                                                                                                     <  2:414230/136‑1 (MQ=255)               ttcgttAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGt                                                                                                                     >  1:414230/1‑136 (MQ=255)               ttcgttAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTtgt                                                                                                                  >  1:256088/1‑139 (MQ=255)                 cgttAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGc                                                                                                                >  1:266127/1‑139 (MQ=255)                 cgttAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGc                                                                                                                <  2:174566/139‑1 (MQ=255)                 cgttAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGc                                                                                                                >  2:719561/1‑139 (MQ=255)                             aGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGtt                                                                                                    <  2:235102/139‑1 (MQ=255)                                       tcggtcggTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATcc                                                                                          <  1:509966/139‑1 (MQ=255)                                         ggtcggtGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGa                                                                                        >  1:233040/1‑139 (MQ=255)                                         ggtcggtGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGa                                                                                        <  2:256088/139‑1 (MQ=255)                                         ggtcggtGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGa                                                                                        <  1:442856/139‑1 (MQ=255)                                          gtcggtGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGAc                                                                                       <  2:692532/139‑1 (MQ=255)                                          gtcggtGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGAc                                                                                       <  2:266127/139‑1 (MQ=255)                                            cggtGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAg                                                                                     <  2:24319/139‑1 (MQ=255)                                            cggtGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAg                                                                                     <  2:734340/139‑1 (MQ=255)                                             ggtGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGc                                                                                    <  2:145941/139‑1 (MQ=255)                                                     gcctgcGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGa                                                                            >  1:374252/1‑139 (MQ=255)                                                     gcctgcGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGa                                                                            >  1:784068/1‑139 (MQ=255)                                                     gcctgcGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGa                                                                            >  2:105389/1‑139 (MQ=255)                                                         gcgcCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCgg                                                                        <  1:719561/139‑1 (MQ=255)                                                         gcgcCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCgg                                                                        <  2:374252/139‑1 (MQ=255)                                                              gggAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATg                                                                   <  1:105389/139‑1 (MQ=255)                                                                           tttCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTCTCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAAAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCaa                                                      <  2:600891/139‑1 (MQ=255)                                                                           tttCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCaa                                                      <  1:597053/139‑1 (MQ=255)                                                                              cAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAg                                                                                     >  1:613654/1‑105 (MQ=255)                                                                              cAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAg                                                                                     <  2:613654/105‑1 (MQ=255)                                                                                  gTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATcc                                                                                          <  1:178508/96‑1 (MQ=255)                                                                                  gTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATcc                                                                                          >  2:178508/1‑96 (MQ=255)                                                                                         cggcTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTg                                        >  2:473810/1‑139 (MQ=255)                                                                                             ttGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGcc                                    <  2:784068/139‑1 (MQ=255)                                                                                                  tatTGAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGCCAAATg                               <  2:67207/139‑1 (MQ=255)                                                                                                      gAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGCCAAATGCCAg                           >  2:92543/1‑139 (MQ=255)                                                                                                      gAACAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGCCAAATGCCAg                           >  1:795793/1‑139 (MQ=255)                                                                                                       aaCAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGCCAAATGCCAGc                          >  1:318506/1‑139 (MQ=255)                                                                                                       aaCAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGCCAAATGCCAGc                          >  2:436633/1‑139 (MQ=255)                                                                                                       aaCAGGAGGCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGCCAAATGCCAGc                          >  1:498860/1‑139 (MQ=255)                                                                                                          aggaggCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGcc                                    <  1:540859/126‑1 (MQ=255)                                                                                                          aggaggCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGcc                                    >  2:540859/1‑126 (MQ=255)                                                                                                          aggaggCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGCCAAATGCCAGCCTg                       <  2:225670/139‑1 (MQ=255)                                                                                                            gaggCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCg                                                                                   <  1:404334/77‑1 (MQ=255)                                                                                                            gaggCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCg                                                                                   >  2:404334/1‑77 (MQ=255)                                                                                                             aggCGTTATGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGCCAAATGCCAGCCTGCGt                    >  2:252180/1‑139 (MQ=255)                                                                                                                     tGCCAATTTATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGCCAAATGCCAGCCTGCGTGGCGAttt            <  1:436633/139‑1 (MQ=255)                                                                                                                             tATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAg                                                                                     <  2:643918/58‑1 (MQ=255)                                                                                                                             tATCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAg                                                                                     >  1:643918/1‑58 (MQ=255)                                                                                                                               tCAGATAAACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGCCAAATGCCAGCCTGCGTGGCGATTTTGGTCGTATc  >  2:792259/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                                    CTGCTGCTTTAGCGTTCGTTAGTGAAAAAGGTGAAACCTCGGTCGGTGAACTGCCTGCGCCGGGAAGCCTGTATTTTCAACGTTTTGGCGGCTTGAATATTGAACAGGAGGCGTTATGCCAATTTATCAGATAGACGGTCTGACTCCGGTTGTGCCAGAAGAGAGTTTTGTCCATCCGACAGCGGTATTGATCGGCGATGTTATTCTCGGCAAGGGCGTTTACGTTGGGCCAAATGCCAGCCTGCGTGGCGATTTTGGTCGTATC  >  NC_000913/1464356‑1464620

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