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%	174.3 / NA	55	A187A (GCG→GCA)	hyuA	D‑stereospecific phenylhydantoinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/26);  total (29/26)

AGTCGAGGAAGGACTGAGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAA  >  NC_000913/3010453‑3010722                                                                                                                                        |                                                                                                                                       aGTCGAGGAAGGACTGAGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCt                                                                                                                                     >  2:1101903/1‑139 (MQ=255)   tCGAGGAAGGACTGAGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTAt                                                                                                                                   >  1:791982/1‑139 (MQ=255)       ggaaggaCTGAGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGcc                                                                                                                               >  1:598355/1‑139 (MQ=255)             aCTGAGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAg                                                                                                                         >  1:111168/1‑139 (MQ=255)                 aGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTCTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAgcggg                                                                                                                     >  2:1561457/1‑139 (MQ=255)                      ttttAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGaa                                                                                                                >  1:784685/1‑139 (MQ=255)                      ttttAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGaa                                                                                                                >  1:1236706/1‑139 (MQ=255)                           aaCTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTAt                                                                                                           <  1:1101903/139‑1 (MQ=255)                           aaCTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTAt                                                                                                           >  1:811700/1‑139 (MQ=255)                           aaCTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTAt                                                                                                           <  2:791982/139‑1 (MQ=255)                                tATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATcgccg                                                                                                      >  2:486095/1‑139 (MQ=255)                                tATTTAACCTATCAATACAAACTCAAAGATGACGAGGTTTTGCAAGCATTACGCCGCTTGCATGAATCAGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATcgccg                                                                                                      >  2:100333/1‑139 (MQ=255)                                  tttAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCa                                                                                                    <  2:1236706/139‑1 (MQ=255)                                      aCCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGt                                                                                                >  1:917964/1‑139 (MQ=255)                                            cAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccg                                                                                          >  1:1598740/1‑139 (MQ=255)                                            cAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCACTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccg                                                                                          >  2:1442319/1‑139 (MQ=255)                                            cAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCAGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccg                                                                                          >  1:872245/1‑139 (MQ=255)                                                 cAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcgccg                                                                                     <  2:598355/139‑1 (MQ=255)                                                 cAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcgccg                                                                                     <  2:934706/139‑1 (MQ=255)                                                 cAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcgccg                                                                                     >  1:480607/1‑139 (MQ=255)                                                   aaCTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcgccgcg                                                                                   >  1:1102668/1‑139 (MQ=255)                                                    aCTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACcgcgccgcgc                                                                                  >  1:1430454/1‑139 (MQ=255)                                                     cTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAgcgggcg                                                                                                                   <  2:1657417/105‑1 (MQ=255)                                                     cTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAgcgggcg                                                                                                                   >  1:1657417/1‑105 (MQ=255)                                                       cAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTAt                                                                               <  1:388293/139‑1 (MQ=255)                                                        aaCGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATc                                                                              <  2:784685/139‑1 (MQ=255)                                                         aCGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCa                                                                             >  1:590268/1‑139 (MQ=255)                                                         aCGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCa                                                                             <  2:480607/139‑1 (MQ=255)                                                          cGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCAt                                                                            >  1:1515843/1‑139 (MQ=255)                                                                       ttGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGcc                                                               >  2:1207867/1‑139 (MQ=255)                                                                                       gcTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGc                                               >  2:911071/1‑139 (MQ=255)                                                                                        cTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcg                                                                                        >  2:154417/1‑97 (MQ=255)                                                                                        cTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccgcg                                                                                        <  1:154417/97‑1 (MQ=255)                                                                                         ttGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccg                                                                                          >  2:867872/1‑94 (MQ=255)                                                                                         ttGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAAccg                                                                                          <  1:867872/94‑1 (MQ=255)                                                                                                aaTCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAg                                                                    >  1:472138/1‑109 (MQ=255)                                                                                                aaTCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAg                                                                    <  2:472138/109‑1 (MQ=255)                                                                                                     ggCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGccg                                 <  2:211917/139‑2 (MQ=255)                                                                                                     ggCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGccc                                 <  1:854778/139‑1 (MQ=255)                                                                                                     ggCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGccc                                 <  1:911071/139‑1 (MQ=255)                                                                                                     ggCGCGCTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCCAAGCGGAAGCCATCGccc                                 <  2:917964/139‑1 (MQ=255)                                                                                                      gcgcgcTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCg                                <  2:1102668/139‑1 (MQ=255)                                                                                                      gcgcgcTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCg                                <  1:486095/139‑1 (MQ=255)                                                                                                      gcgcgcTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCg                                <  2:242784/139‑1 (MQ=255)                                                                                                      gcgcgcTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCg                                <  2:1598740/139‑1 (MQ=255)                                                                                                           cTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGc                               >  2:596909/1‑135 (MQ=255)                                                                                                           cTGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGc                               <  1:596909/135‑1 (MQ=255)                                                                                                            tGACCACCGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGAt                          <  2:590268/139‑1 (MQ=255)                                                                                                            tGACCACCGTGCACCAGGAAAATGATGCAGTTATCGCCAGCAAGCGGGCGGAATTTATCGCCGTAGGGTTAACCGCTCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGTATGAt                          <  2:1231520/139‑1 (MQ=255)                                                                                                                caccGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCTGCATGATTAAc                      <  1:1561457/139‑1 (MQ=255)                                                                                                                caccGTGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAAc                      <  1:1207867/139‑1 (MQ=255)                                                                                                                     tGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTAt                                                                                                           >  2:179296/1‑49 (MQ=255)                                                                                                                     tGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTAt                                                                                                           <  1:179296/49‑1 (MQ=255)                                                                                                                     tGCACCCGGAAAATGATGCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGc                 >  2:1503671/1‑139 (MQ=255)                                                                                                                                    atgCAGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTGCCGGTaa  >  2:1496946/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                       AGTCGAGGAAGGACTGAGCAGTTTTAAACTCTATTTAACCTATCAATACAAACTCAACGATGACGAGGTTTTGCAGGCATTACGCCGCTTGCATGAATCCGGCGCGCTGACCACCGTGCACCCGGAAAATGATGCGGCTATCGCCAGCAAGCGGGCGGAATTTATCGCCGCAGGGTTAACCGCGCCGCGCTATCATGCCTTGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAA  >  NC_000913/3010453‑3010722

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