Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,995,567	A→G	T635A (ACC→GCC)	xdhA →	xanthine dehydrogenase molybdenum‑binding subunit XdhA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,995,567	0	A	G	100.0%	98.9 / NA	30	T635A (ACC→GCC)	xdhA	xanthine dehydrogenase molybdenum‑binding subunit XdhA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/13);  total (17/13)

AAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATG  >  NZ_CP009273/2995481‑2995650                                                                                       |                                                                                    aaaGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCa                                                                                  <  1:110704/90‑1 (MQ=255)     cTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCaa                                                                              <  2:151136/90‑1 (MQ=255)     cTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCaa                                                                              >  1:709852/1‑90 (MQ=255)       ccATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAAcc                                                                            >  2:117289/1‑90 (MQ=255)               accaccACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTc                                                                    <  1:325145/90‑1 (MQ=255)                ccaccaCTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCa                                                                   <  2:639089/90‑1 (MQ=255)                 caccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCaa                                                                  <  1:266194/90‑1 (MQ=255)                  accacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAAc                                                                 >  1:163081/1‑90 (MQ=255)                   ccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACg                                                                >  1:543522/1‑90 (MQ=255)                   ccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACg                                                                <  2:543522/90‑1 (MQ=255)                        aaCCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCAt                                                           <  1:589006/90‑1 (MQ=255)                             aCCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTc                                                      >  1:460049/1‑90 (MQ=255)                                           gTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTa                                        >  2:251137/1‑90 (MQ=255)                                           gTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTa                                        >  2:610481/1‑90 (MQ=255)                                            tACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTaa                                       >  1:213710/1‑90 (MQ=255)                                                 ttgttgATCTGACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCAc                                  <  2:527043/90‑1 (MQ=255)                                                           gACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGa                        >  2:190096/1‑90 (MQ=255)                                                           gACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGa                        <  2:639095/90‑1 (MQ=255)                                                           gACGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGa                        >  1:135937/1‑90 (MQ=255)                                                            aCGGTCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGaa                       <  1:303817/90‑1 (MQ=255)                                                                tCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTc                   >  2:74826/1‑90 (MQ=255)                                                                tCGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACCGCTGGCAGAAGGTc                   >  2:543239/1‑90 (MQ=255)                                                                 cGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                              >  1:617916/1‑62 (MQ=255)                                                                 cGATATTGCGCTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCata                                              <  2:617916/62‑1 (MQ=255)                                                                           cTGTGCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggc        <  1:74826/90‑1 (MQ=255)                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg      >  1:701275/1‑90 (MQ=255)                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg      >  1:489665/1‑90 (MQ=255)                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg      >  1:406119/1‑90 (MQ=255)                                                                             gtgCAAAGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg      >  1:374195/1‑90 (MQ=255)                                                                                 aaaGTCGCCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCggaatg  <  1:195228/90‑1 (MQ=255)                                                                                       |                                                                                    AAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATG  >  NZ_CP009273/2995481‑2995650

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