Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,003,097	C→T	H162Y (CAT→TAT)	xdhB →	putative xanthine dehydrogenase FAD‑binding subunit XdhB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,003,097	0	C	T	100.0%	116.2 / NA	33	H162Y (CAT→TAT)	xdhB	putative xanthine dehydrogenase FAD‑binding subunit XdhB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/16);  total (17/16)

AGATCCACTCCCCGCGCGGTGTTCGTTTCGTCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGCATGACGAAATCCTCGTCGCCTTTCATTTTCCGCCACAGCCGAAAGAACACGCGGGCAGCGCGCATT  >  NC_000913/3003020‑3003163                                                                              |                                                                   agaTCCACTCCCCGCGCGGTGTTCGTTTCGTCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTAt                                                                  >  2:491563/1‑80 (MQ=255)   aTCCACTCCCCGCGCGGTGTTCGTTTCGTCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGAc                                                               >  2:491577/1‑81 (MQ=255)    tCCACTCCCCGCGCGGTGTTCGTTTCGTCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACg                                                              <  1:491523/81‑1 (MQ=255)      cACTCCCCGCGCGGTGTTCGTTTCGTCCCGATTAATGGCTTTCATACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGa                                                             <  2:491523/80‑1 (MQ=255)      cACTCCCCGCGCGGTGTTCGTTTCGTCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGa                                                             >  2:491581/1‑80 (MQ=255)       aCTCCCCGCGCGGTGTTCGTTTCGTCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGaaa                                                           >  2:491561/1‑81 (MQ=255)           cccGCGCGGTGTTCGTTTCGTCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCt                                                       >  1:491569/1‑81 (MQ=255)              gcgcgGTGTTCGTTTCGTCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCtcgt                                                    <  1:491525/81‑1 (MQ=255)               cgcgGTGTTCGTTTCGTCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCtcgtc                                                   <  2:491525/81‑1 (MQ=255)                gcgGTGTTCGTTTCGTCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCtcgtcg                                                  >  2:491569/1‑81 (MQ=255)                  ggTGTTCGTTTCGTCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGcc                                                <  2:491526/81‑1 (MQ=255)                    tgtTCGTTTCGTCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCtt                                              <  2:491527/81‑1 (MQ=255)                      ttcgtttcgtCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCttt                                             >  1:491575/1‑80 (MQ=255)                      ttcgtttcgtCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCttt                                             >  1:491576/1‑80 (MQ=255)                      ttcgtttcgtCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCttt                                             >  2:491571/1‑80 (MQ=255)                       tcgtttcgtCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCa                                           >  1:491556/1‑81 (MQ=255)                        cgtttcgtCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCAt                                          >  2:491573/1‑81 (MQ=255)                         gtttcgtCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCAtt                                         >  2:491579/1‑81 (MQ=255)                            tcgtCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCAtttt                                       <  2:491528/80‑1 (MQ=255)                              gtCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCATTTTCCg                                    >  2:491565/1‑81 (MQ=255)                                 ccGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCATTTTCCGcc                                  <  1:491530/80‑1 (MQ=255)                                    aTTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCATTTTCCGCCACAg                              <  2:491530/81‑1 (MQ=255)                                    aTTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCATTTTCCGCCACAg                              >  1:491590/1‑81 (MQ=255)                                     ttAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCATTTTCCGCCACAGc                             <  2:491531/81‑1 (MQ=255)                                       aaTGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCATTTTCCGCCACAGcc                            <  1:491533/80‑1 (MQ=255)                                         tGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCATTTTCCGCCACAGCCGaa                         <  1:491534/81‑1 (MQ=255)                                          ggCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCATTTTCCGCCACAGCCGa                          >  2:491580/1‑79 (MQ=255)                                            cTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCATTTTCCGCCACAGCCGAAAg                       >  2:491583/1‑80 (MQ=255)                                               tCACACCGGGCCGGGCAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCATTTTCCGCCACAGCCGAAAGaa                     <  1:491535/79‑1 (MQ=255)                                                          cgggcAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCATTTTCCGCCACAGCCGAAAGAACACGCGGGCAgcg        <  2:491535/81‑1 (MQ=255)                                                          cgggcAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCATTTTCCGCCACAGCCGAAAGAACACGCGGGCAgcg        >  2:491576/1‑81 (MQ=255)                                                             gcAAAGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCATTTTCCGCCACAGCCGAAAGAACACGCGGGCAgcgcgc     <  2:491536/81‑1 (MQ=255)                                                                aaGTGTCTCTTGAGTATGACGAAATCCTCGTCGCCTTTCATTTTCCGCCACAGCCGAAAGAACACGCGGGCAGCGCGCAtt  <  2:491537/81‑1 (MQ=255)                                                                              |                                                                   AGATCCACTCCCCGCGCGGTGTTCGTTTCGTCCCGATTAATGGCTTTCACACCGGGCCGGGCAAAGTGTCTCTTGAGCATGACGAAATCCTCGTCGCCTTTCATTTTCCGCCACAGCCGAAAGAACACGCGGGCAGCGCGCATT  >  NC_000913/3003020‑3003163

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