Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,023,896	G→T	G861C (GGC→TGC)	xdhD →	fused putative xanthine/hypoxanthine oxidase: molybdopterin‑binding subunit and Fe‑S binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,023,896	0	G	T	100.0%	82.9 / NA	25	G861C (GGC→TGC)	xdhD	fused putative xanthine/hypoxanthine oxidase: molybdopterin‑binding subunit and Fe‑S binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (10/15);  total (10/15)

ACTTCGCTGAAGTTGCCGTCAACACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGA  >  NC_000913/3023759‑3024026                                                                                                                                          |                                                                                                                                   aCTTCGCTGAAGTTGCCGTCAACACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGc                                                                                                                                  <  1:82397/140‑1 (MQ=255)               gCCGTCAACACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGa                                                                                                                    <  2:132576/140‑1 (MQ=255)                ccGTCAACACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGaa                                                                                                                   <  2:239831/140‑1 (MQ=255)                ccGTCAACACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGaa                                                                                                                   <  1:530379/140‑1 (MQ=255)                  gTCAACACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAaga                                                                                                                 <  1:272336/140‑1 (MQ=255)                            gTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTAt                                                                                                       <  1:21993/140‑1 (MQ=255)                                 ggTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGc                                                                                                  <  1:475195/140‑1 (MQ=255)                                            ccTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCAcc                                                                                       >  2:642111/1‑140 (MQ=255)                                                gATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGtt                                                                                   >  1:265607/1‑140 (MQ=255)                                                     aTTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAAcgc                                                                              <  1:92160/140‑1 (MQ=255)                                                               ttGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACg                                                                    <  1:426968/140‑1 (MQ=255)                                                               ttGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACg                                                                    <  1:642111/140‑1 (MQ=255)                                                                          cGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGgcg                                                         >  1:235555/1‑140 (MQ=255)                                                                            gTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCa                                                       <  2:265607/140‑1 (MQ=255)                                                                               cacCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCg                                                    >  2:83303/1‑140 (MQ=255)                                                                                  cGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGc                                                                   <  2:249320/122‑1 (MQ=255)                                                                                  cGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGc                                                                   >  1:249320/1‑122 (MQ=255)                                                                                         tCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTg                                          >  2:182692/1‑140 (MQ=255)                                                                                         tCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCAAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTg                                          >  2:470912/1‑140 (MQ=255)                                                                                                  aGCGTTGGGACAAATCTACGGTGCCTCCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCcgc                                 <  1:582243/140‑1 (MQ=255)                                                                                                   gCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCcgcg                                >  1:693549/1‑140 (MQ=255)                                                                                                               aTCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCgcgc                    <  2:322732/140‑1 (MQ=255)                                                                                                               aTCTACGGTGCCACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCgcgc                    <  2:693549/140‑1 (MQ=255)                                                                                                                         ccACCCTGCGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTg          >  2:354684/1‑140 (MQ=255)                                                                                                                                 cGTGCTATCTGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGcga  >  1:132425/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                   ACTTCGCTGAAGTTGCCGTCAACACGCGTACGGGTGAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGA  >  NC_000913/3023759‑3024026

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