Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,462,292	A→G	E56E (GAA→GAG)	paaK →	phenylacetate‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,462,292	0	A	G	100.0%	201.3 / NA	54	E56E (GAA→GAG)	paaK	phenylacetate‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/31);  total (23/31)

CGTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGA  >  NC_000913/1462162‑1462425                                                                                                                                   |                                                                                                                                      cgTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGAc                                                                                                                              >  1:1027465/1‑140 (MQ=255) cgTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGAc                                                                                                                              <  1:1082810/140‑1 (MQ=255)     tGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGc                                                                                                                          >  2:762339/1‑140 (MQ=255)            gaGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAAtt                                                                                                                   >  1:837465/1‑140 (MQ=255)              gTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTc                                                                                                                 <  2:895078/140‑1 (MQ=255)               ttACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTcc                                                                                                                <  2:346203/140‑1 (MQ=255)                   aGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGt                                                                                                            >  1:765143/1‑140 (MQ=255)                    ggCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTa                                                                                                           <  1:762339/140‑1 (MQ=255)                       gTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTacca                                                                                                        <  2:1129227/140‑1 (MQ=255)                        ttGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTaccac                                                                                                       <  2:796975/140‑1 (MQ=255)                        ttGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTaccac                                                                                                       <  2:1072758/140‑1 (MQ=255)                        ttGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTaccac                                                                                                       <  2:915820/140‑1 (MQ=255)                        ttGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTaccac                                                                                                       <  1:764432/140‑1 (MQ=255)                        ttGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTaccac                                                                                                       <  1:692189/140‑1 (MQ=255)                         tGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTaccacc                                                                                                      >  2:912813/1‑140 (MQ=255)                         tGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTaccacc                                                                                                      >  1:80003/1‑140 (MQ=255)                            aaaCACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCaaa                                                                                                   <  2:288481/140‑1 (MQ=255)                             aaCACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTACCCTGTACCACCAAAc                                                                                                  <  1:330858/140‑1 (MQ=255)                                acaGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAgg                                                                                               <  1:1045368/140‑1 (MQ=255)                                   gCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATc                                                                                            <  1:704874/140‑1 (MQ=255)                                   gCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATc                                                                                            <  1:417081/140‑1 (MQ=255)                                    cGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCt                                                                                           >  1:589406/1‑140 (MQ=255)                                        tGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTgcgc                                                                                       >  1:933058/1‑140 (MQ=255)                                         gAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGcgcg                                                                                      <  2:906159/140‑1 (MQ=255)                                                 cGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTAt                                                                              >  2:241828/1‑140 (MQ=255)                                                 cGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCCTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTAt                                                                              <  2:1027465/140‑1 (MQ=255)                                                  gCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATc                                                                             <  1:999784/140‑1 (MQ=255)                                                    tCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATccc                                                                           >  2:40532/1‑140 (MQ=255)                                                                      aTGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTg                                                         >  1:554354/1‑140 (MQ=255)                                                                      aTGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTg                                                         >  2:895242/1‑140 (MQ=255)                                                                      aTGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTg                                                         <  2:589406/140‑1 (MQ=255)                                                                           ccGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAAt                                                    <  2:80003/140‑1 (MQ=255)                                                                           ccGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAAt                                                    >  2:679500/1‑140 (MQ=255)                                                                               tGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGaa                                                <  2:765143/140‑1 (MQ=255)                                                                                        gcAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAgt                                           <  2:1138691/136‑1 (MQ=255)                                                                                        gcAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAgt                                           >  1:1138691/1‑136 (MQ=255)                                                                                           aaTTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTgcgc                                    <  2:511979/140‑1 (MQ=255)                                                                                                          gCGTACATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCg                     <  1:912813/140‑1 (MQ=255)                                                                                                               cATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCgg                    <  2:443243/136‑1 (MQ=255)                                                                                                               cATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCgg                    >  1:443243/1‑136 (MQ=255)                                                                                                               cATCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAAc                <  1:40532/140‑1 (MQ=255)                                                                                                                 tCCTGATGATTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACca              >  1:337373/1‑140 (MQ=255)                                                                                                                    tgatgaTTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGTAACCACAg           >  1:628885/1‑140 (MQ=255)                                                                                                                    tgatgaTTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAg           >  1:354911/1‑140 (MQ=255)                                                                                                                    tgatgaTTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAg           >  2:1147818/1‑140 (MQ=255)                                                                                                                      atgatTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCggg                   <  2:31187/130‑1 (MQ=255)                                                                                                                      atgatTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCggg                   >  1:31187/1‑130 (MQ=255)                                                                                                                      atgatTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACca              >  2:799534/1‑135 (MQ=255)                                                                                                                      atgatTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACca              <  1:799534/135‑1 (MQ=255)                                                                                                                      atgatTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGt         <  1:1147818/140‑1 (MQ=255)                                                                                                                       tgatTTCAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTa        >  1:1079167/1‑140 (MQ=255)                                                                                                                            tCAGGTAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCg   >  1:1061402/1‑140 (MQ=255)                                                                                                                             cAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGa  <  2:312649/140‑1 (MQ=255)                                                                                                                             cAGGGAGCTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGa  <  1:895242/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                      CGTCTGTTGATGAGTTACAGGCGTTGCAAACACAGCGTCTGAAATGGACGCTCAAACACGCGTATGAAAATGTGCCGATGTATCGGCGCAAATTCGACGCAGCAGGCGTACATCCTGATGATTTCAGGGAACTTTCAGACCTGCGTAAATTTCCCTGTACCACCAAACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGA  >  NC_000913/1462162‑1462425

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