Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,462,463	G→A	W113* (TGG→TGA)	paaK →	phenylacetyl‑CoA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,462,463	0	G	A	100.0%	126.5 / NA	41	W113* (TGG→TGA)	paaK	phenylacetyl‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/18);  total (23/18)

AACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGGGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGCCTGTTTACTGGTGGGCTGGGTGCGCACTATGGTGCCGAACGTCTGGGCGCTACGG  >  NC_000913/1462327‑1462596                                                                                                                                         |                                                                                                                                      aaCAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGc                                                                                                                                     >  2:1581674/1‑139 (MQ=255)            cgcgACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAg                                                                                                                          <  1:1581674/139‑1 (MQ=255)                             tGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTtgcgtgcg                                                                                                         >  1:1193503/1‑139 (MQ=255)                                 acCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAg                                                                                                     <  2:602264/139‑1 (MQ=255)                                     ttGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAggggg                                                                                                 >  2:9196/1‑139 (MQ=255)                                           tGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGc                                                                                                                                     <  2:1462504/97‑1 (MQ=255)                                             ccAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAAGGATATTGATAACTGAGCCAATATTGTAGCGCGTTATTTGCGTGCGGCAGGGGGCTCGCCCa                                                                                         >  1:1023904/1‑139 (MQ=255)                                             ccAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCa                                                                                         >  1:1475144/1‑139 (MQ=255)                                                      cAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATAtt                                                                                                                              >  2:1539901/1‑93 (MQ=255)                                                      cAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATAtt                                                                                                                              <  1:1539901/93‑1 (MQ=255)                                                            gtgCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGc                                                                                                                                     <  1:101139/80‑1 (MQ=255)                                                             tgCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACTCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCAt                                                                           <  1:953556/137‑1 (MQ=255)                                                             tgCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCa                                                                                         <  2:444161/123‑1 (MQ=255)                                                             tgCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCa                                                                                         >  1:444161/1‑123 (MQ=255)                                                             tgCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCAt                                                                           >  2:953556/1‑137 (MQ=255)                                                                 cATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGTGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGcc                                                                     >  2:1188496/1‑139 (MQ=255)                                                                    tCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTAt                                                                  >  1:1532397/1‑139 (MQ=255)                                                                           tcttcGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACg                                                           >  1:1087970/1‑139 (MQ=255)                                                                              tcGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGa                                                                                     <  2:144703/110‑1 (MQ=255)                                                                              tcGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGa                                                                                     >  1:144703/1‑110 (MQ=255)                                                                              tcGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGcc                                                        >  2:729856/1‑139 (MQ=255)                                                                                     ccACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGCCTGtt                                                    >  1:1586771/1‑136 (MQ=255)                                                                                     ccACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGCCTGtt                                                    <  2:1586771/136‑1 (MQ=255)                                                                                            tAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCAt                                                                           >  1:574936/1‑106 (MQ=255)                                                                                            tAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCAt                                                                           <  2:574936/106‑1 (MQ=255)                                                                                             aaaCCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGCCTGTTTACTGGTGGGc                                         <  1:720831/139‑1 (MQ=255)                                                                                              aaCCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGTGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGCCTGTTTACTGGTGGGCt                                        <  1:1188496/139‑1 (MQ=255)                                                                                                ccGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCaa                                                                                                                                  <  2:1152413/47‑1 (MQ=255)                                                                                                ccGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGAGCCaa                                                                                                                                  >  1:1152413/1‑47 (MQ=255)                                                                                                                 taCGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGtt                                                                        <  1:419961/88‑1 (MQ=255)                                                                                                                 taCGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGtt                                                                        >  2:419961/1‑88 (MQ=255)                                                                                                                 taCGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGCCTGTTTACTGGTGGGCTGGGTGCGCACTATGGTGcc                     >  2:1300801/1‑139 (MQ=255)                                                                                                                 taCGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGCCTGTTTACTGGTGGGCTGGGTGCGCACTATGGTGcc                     >  2:670773/1‑139 (MQ=255)                                                                                                                 taCGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGCCTGTTTACTGGTGGGCTGGGTGCGCACTATGGTGcc                     >  1:388652/1‑139 (MQ=255)                                                                                                                   cGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGCCTGTTTACTGGTGGGCTGGGt                                    <  2:1094880/122‑1 (MQ=255)                                                                                                                   cGCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGCCTGTTTACTGGTGGGCTGGGt                                    >  1:1094880/1‑122 (MQ=255)                                                                                                                    gCAAAACGATATTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGCCTGTTTACTGGTGGGCTGGGTGCGCACTATGGTGCCGaa                  >  1:1199673/1‑139 (MQ=255)                                                                                                                            atatTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGCCTGTTTACTGGTGGGCTGGGTGCGCACTATGGTGCCGAACGTCTggg          <  2:1193503/139‑1 (MQ=255)                                                                                                                             tatTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGtt                                                                        >  2:1450673/1‑76 (MQ=255)                                                                                                                             tatTGATAACTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGtt                                                                        <  1:1450673/76‑1 (MQ=255)                                                                                                                                    aaCTGAGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGCCTGTTTACTGGTGGGCTGGGTGCGCACTATGGTGCCGAACGTCTGGGCGCTACgg  <  2:1532397/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                      AACAGGATCTGCGCGACAACTATCCCTTTGACACCTTTGCCGTGCCAATGGAACAAGTGGTGCGCATTCATGCTTCTTCGGGAACCACAGGTAAACCGACAGTTGTCGGCTATACGCAAAACGATATTGATAACTGGGCCAATATTGTAGCGCGTTCTTTGCGTGCGGCAGGGGGCTCGCCCAAAGACAAAATTCATGTTGCCTATGGTTACGGCCTGTTTACTGGTGGGCTGGGTGCGCACTATGGTGCCGAACGTCTGGGCGCTACGG  >  NC_000913/1462327‑1462596

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