Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	543,951	A→G	D32D (GAT→GAC)	allE ←	(S)‑ureidoglycine aminohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	543,951	0	A	G	100.0%	85.5 / NA	25	D32D (GAT→GAC)	allE	(S)‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/9);  total (16/9)

AAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTATTTCTGCCAGGCAAGTTGAT  >  NC_000913/543820‑544079                                                                                                                                    |                                                                                                                                 aaGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGt                                                                                                                          >  1:413331/1‑140 (MQ=255)           ggttgCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCg                                                                                                                >  2:592922/1‑140 (MQ=255)            gttgCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGa                                                                                                               >  2:732448/1‑140 (MQ=255)                 cACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTa                                                                                                          >  1:693270/1‑140 (MQ=255)                   ccGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTAcc                                                                                                        <  2:693270/140‑1 (MQ=255)                                gtgtGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATc                                                                                           >  2:21739/1‑140 (MQ=255)                                   tgACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCa                                                                                        <  1:21739/140‑1 (MQ=255)                                         gATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTa                                                                                  >  2:747876/1‑140 (MQ=255)                                         gATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTa                                                                                  >  2:630494/1‑140 (MQ=255)                                           tAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGc                                                                                >  2:236557/1‑140 (MQ=255)                                             aTCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCa                                                                                        <  1:369018/130‑1 (MQ=255)                                             aTCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCa                                                                                        >  2:369018/1‑130 (MQ=255)                                                          gCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTcgcg                                                                 >  2:317736/1‑140 (MQ=255)                                                          gCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTcgcg                                                                 <  1:236557/140‑1 (MQ=255)                                                                           ggAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTAt                                                <  2:399041/140‑1 (MQ=255)                                                                            gagaGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTAtt                                               >  1:727937/1‑140 (MQ=255)                                                                                      tCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGt                                                   >  2:483533/1‑126 (MQ=255)                                                                                      tCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGt                                                   <  1:483533/126‑1 (MQ=255)                                                                                           tCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATc                                       <  1:657393/133‑1 (MQ=255)                                                                                           tCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATc                                       >  2:657393/1‑133 (MQ=255)                                                                                           tCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTa                                >  2:447008/1‑140 (MQ=255)                                                                                           tCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTa                                >  2:627059/1‑140 (MQ=255)                                                                                                           ccGGAATAATATTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTATTTCTg                <  2:579110/140‑1 (MQ=255)                                                                                                                  aataTTTTTTACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTATTTCTGCCAGGCa         <  1:675129/140‑1 (MQ=255)                                                                                                                         tttACCAGACCGTCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTATTTCTGCCAGGCAAGTTGAt  >  1:478083/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                 AAGCCCTGTTGGTTGCCACCGTTTTGATGCAGTGTGACCAGATAATCAACAAAAGAGGCACCCAGCTTTGGCGTGGAGAGGATTGTCGCGTCACAATTTTCAAAGCCCGGAATAATATTTTTTACCAGACCATCCGGGGTTAACAGTGCGAAATTACCGTGTTTAACAATCGCACGGTTAGCCAGTAAATCTTCGCGGTAACCGGTGACGTTATTTAAATATCCCATTTATGACTCCTTATTTCTGCCAGGCAAGTTGAT  >  NC_000913/543820‑544079

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