Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,010,588	G→A	A187A (GCG→GCA)	hyuA →	D‑stereospecific phenylhydantoinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,010,588	0	G	A	100.0%	214.0 / NA	67	A187A (GCG→GCA)	hyuA	D‑stereospecific phenylhydantoinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (36/31);  total (36/31)

GAAGGACTGAGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAA  >  NC_000913/3010460‑3010722                                                                                                                                 |                                                                                                                                       gaaggaCTGAGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCa                                                                                                                              >  2:1114087/1‑139 (MQ=255)   aggaCTGAGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGc                                                                                                                            >  2:1742232/1‑139 (MQ=255)     gaCTGAGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCaa                                                                                                                          <  1:1742232/139‑1 (MQ=255)      aCTGAGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAg                                                                                                                         >  2:480510/1‑139 (MQ=255)          aGCAGTTTTACACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTa                                                                                                                                    >  1:2889/1‑124 (MQ=255)          aGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTa                                                                                                                                    <  2:2889/124‑1 (MQ=255)          aGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAgcggg                                                                                                                     <  1:1114087/139‑1 (MQ=255)          aGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAgcggg                                                                                                                     <  1:254452/139‑1 (MQ=255)               ttttAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGaa                                                                                                                >  1:611684/1‑139 (MQ=255)                   aaaCTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTa                                                                                                            <  1:480510/139‑1 (MQ=255)                         tATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATcgccg                                                                                                      >  2:1336632/1‑139 (MQ=255)                         tATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATcgccg                                                                                                      <  2:245894/139‑1 (MQ=255)                          aTTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATcgccgc                                                                                                     <  1:1246052/139‑1 (MQ=255)                                   aTCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAc                                                                                            >  2:30992/1‑139 (MQ=255)                                   aTCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAc                                                                                            >  1:180835/1‑139 (MQ=255)                                     cAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccg                                                                                          >  2:763204/1‑139 (MQ=255)                                     cAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccg                                                                                          >  2:113259/1‑139 (MQ=255)                                     cAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccg                                                                                          >  1:691075/1‑139 (MQ=255)                                     cAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCACTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccg                                                                                          >  1:1022937/1‑139 (MQ=255)                                       aTACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcg                                                                                        <  1:1336632/139‑1 (MQ=255)                                       aTACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCACTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcg                                                                                        <  2:1022937/139‑1 (MQ=255)                                          cAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcgccg                                                                                     <  2:283363/139‑1 (MQ=255)                                           aaaCTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcgccgc                                                                                    >  1:583403/1‑139 (MQ=255)                                            aaCTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcgccgcg                                                                                   <  2:611684/139‑1 (MQ=255)                                             aCTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcgccgcg                                                                                   <  2:822452/138‑1 (MQ=255)                                             aCTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcgccgcg                                                                                   >  1:822452/1‑138 (MQ=255)                                             aCTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACcgcgccgcgc                                                                                  >  2:420124/1‑139 (MQ=255)                                             aCTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACcgcgccgcgc                                                                                  >  1:830140/1‑139 (MQ=255)                                              cTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcgccg                                                                                     >  1:1045648/1‑135 (MQ=255)                                              cTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcgccg                                                                                     <  2:1045648/135‑1 (MQ=255)                                                 aaCGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATc                                                                              <  2:583403/139‑1 (MQ=255)                                                      tGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGcc                                                                         >  2:443208/1‑139 (MQ=255)                                                         cGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAgcggg                                                                                                                     <  2:939542/92‑1 (MQ=255)                                                         cGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAgcggg                                                                                                                     >  1:939542/1‑92 (MQ=255)                                                           aGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAg                                                                    <  2:180835/139‑1 (MQ=255)                                                             gTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTc                                                                  <  2:1203578/139‑1 (MQ=255)                                                                           aCGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAgcgggcg                                                                                                                   >  1:1763521/1‑76 (MQ=255)                                                                           aCGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAgcgggcg                                                                                                                   <  2:1763521/76‑1 (MQ=255)                                                                           aCGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGc                                                    <  1:146523/139‑1 (MQ=255)                                                                            cgccgcTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCg                                                   <  2:830140/139‑1 (MQ=255)                                                                                gcTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATcgcc                                                                                                       >  1:186272/1‑83 (MQ=255)                                                                                gcTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATcgcc                                                                                                       <  2:186272/83‑1 (MQ=255)                                                                                 cTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCg                                              >  2:1428785/1‑139 (MQ=255)                                                                                    gCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGaa                                           >  1:1631730/1‑139 (MQ=255)                                                                                           tCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGc                                                    <  1:1232714/123‑1 (MQ=255)                                                                                           tCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGc                                                    >  2:1232714/1‑123 (MQ=255)                                                                                              ggCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGccc                                 <  1:113259/139‑1 (MQ=255)                                                                                               gcgcgcTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCg                                >  2:612255/1‑139 (MQ=255)                                                                                                cgcgcTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGc                               >  1:1303165/1‑139 (MQ=255)                                                                                                    cTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCg                                >  2:1427602/1‑134 (MQ=255)                                                                                                    cTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCg                                >  1:57585/1‑134 (MQ=255)                                                                                                    cTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCg                                <  1:1427602/134‑1 (MQ=255)                                                                                                    cTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCg                                <  2:57585/134‑1 (MQ=255)                                                                                                    cTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGa                           <  1:30992/139‑1 (MQ=255)                                                                                                     tGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcgccg                                                                                     >  1:130048/1‑80 (MQ=255)                                                                                                     tGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcgccg                                                                                     <  2:130048/80‑1 (MQ=255)                                                                                                                 aCCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGcc                                                                         <  2:1632582/80‑1 (MQ=255)                                                                                                                 aCCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGcc                                                                         >  1:1632582/1‑80 (MQ=255)                                                                                                                 aCCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCgg                                             <  2:1780790/108‑1 (MQ=255)                                                                                                                 aCCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCgg                                             >  1:1780790/1‑108 (MQ=255)                                                                                                                   ccGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCACaaa            >  2:811978/1‑139 (MQ=255)                                                                                                                      gAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATcgccgc                                                                                                     <  1:1701325/47‑1 (MQ=255)                                                                                                                      gAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATcgccgc                                                                                                     >  2:1701325/1‑47 (MQ=255)                                                                                                                      gAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTg         >  1:98670/1‑139 (MQ=255)                                                                                                                             atgCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCg                                                   <  1:1027287/90‑1 (MQ=255)                                                                                                                             atgCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCg                                                   >  2:1027287/1‑90 (MQ=255)                                                                                                                             atgCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTGCCGGTaa  >  1:1653215/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                       GAAGGACTGAGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAA  >  NC_000913/3010460‑3010722

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