Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,001,812	G→A	V490M (GTG→ATG)	xdhA →	putative xanthine dehydrogenase molybdenum‑binding subunit XdhA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,001,812	0	G	A	100.0%	176.3 / NA	48	V490M (GTG→ATG)	xdhA	putative xanthine dehydrogenase molybdenum‑binding subunit XdhA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/28);  total (20/28)

ACCATCAACGTGCAAAGCGGCGCGACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCGTGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTT  >  NC_000913/3001734‑3001889                                                                               |                                                                              aCCATCAACGTGCAAAGCGGCGCGACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATg                                                                             <  2:544989/81‑1 (MQ=255) aCCATCAACGTGCAAAGCGGCGCGACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATg                                                                             <  1:544988/81‑1 (MQ=255)    aTCAACGTGCAAAGCGGCGCGACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATgggg                                                                          >  1:545046/1‑81 (MQ=255)    aTCAACGTGCAAAGCGGCGCGACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATgggg                                                                          >  1:545045/1‑81 (MQ=255)    aTCAACGTGCAAAGCGGCGCGACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATgggg                                                                          >  2:545044/1‑81 (MQ=255)       aaCGTGCAAAGCGGCGCGACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAGCCATGGGGGtt                                                                       >  2:545042/1‑81 (MQ=255)       aaCGTGCAAAGCGGCGCGACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGtt                                                                       <  2:544990/81‑1 (MQ=255)             cAAAGCGGCGCGACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGt                                                                  >  2:545049/1‑80 (MQ=255)              aaaGCGGCGCGACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGt                                                                  >  2:545051/1‑79 (MQ=255)              aaaGCGGCGCGACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGt                                                                  >  1:545057/1‑79 (MQ=255)                  cggcgCGACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGa                                                            <  1:544991/81‑1 (MQ=255)                   ggcgCGACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGAc                                                           <  1:544992/81‑1 (MQ=255)                       cgACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTc                                                       <  1:544994/81‑1 (MQ=255)                       cgACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTc                                                       <  1:544993/81‑1 (MQ=255)                       cgACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTc                                                       >  1:545054/1‑81 (MQ=255)                        gACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTcg                                                      >  1:545071/1‑81 (MQ=255)                             aaaTCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTa                                                 <  2:544995/81‑1 (MQ=255)                              aaTCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTAt                                                <  2:544996/81‑1 (MQ=255)                                  ggTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTc                                             <  1:544997/80‑1 (MQ=255)                                  ggTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTc                                             >  1:545061/1‑80 (MQ=255)                                     cAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACt                                         >  2:545048/1‑81 (MQ=255)                                      aGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACt                                         >  1:545059/1‑80 (MQ=255)                                      aGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTc                                        <  2:544999/81‑1 (MQ=255)                                      aGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTc                                        <  1:544998/81‑1 (MQ=255)                                        ggTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTc                                        >  2:545075/1‑79 (MQ=255)                                         gTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAg                                     <  1:545000/81‑1 (MQ=255)                                           gCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGa                                    >  2:545058/1‑80 (MQ=255)                                               acacCGTCTTCTCGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCg                               <  2:545001/81‑1 (MQ=255)                                                         ctcGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACg                       <  1:545002/79‑1 (MQ=255)                                                         ctcGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACg                       <  1:545004/79‑1 (MQ=255)                                                         ctcGCAAATGGTGGCAGAAACCATGGGGATTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTAc                        <  2:545003/78‑1 (MQ=255)                                                          tcGCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCg                    >  1:545047/1‑81 (MQ=255)                                                            gCAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGcc                     <  2:545005/78‑1 (MQ=255)                                                             cAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTc                 <  2:545006/81‑1 (MQ=255)                                                             cAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTc                 <  2:545007/81‑1 (MQ=255)                                                             cAAATGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTc                 >  2:545065/1‑81 (MQ=255)                                                                aTGGTGGCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGAt              >  1:545076/1‑81 (MQ=255)                                                                 tggtggCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCTTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATc             <  2:545009/81‑1 (MQ=255)                                                                 tggtggCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATc             >  1:545063/1‑81 (MQ=255)                                                                 tggtggCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATc             <  1:545008/81‑1 (MQ=255)                                                                  ggtggCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATcc            >  1:545069/1‑81 (MQ=255)                                                                   gtggCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATccc           <  2:545011/81‑1 (MQ=255)                                                                   gtggCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATcc            <  2:545010/80‑1 (MQ=255)                                                                    tggCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCg          >  1:545050/1‑81 (MQ=255)                                                                      gCAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGgc        <  1:545012/81‑1 (MQ=255)                                                                       cAGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGgcg       <  2:545013/81‑1 (MQ=255)                                                                        aGAAACCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGgcgc      <  2:545014/81‑1 (MQ=255)                                                                            aCCATGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCAttt  <  1:545015/81‑1 (MQ=255)                                                                               |                                                                              ACCATCAACGTGCAAAGCGGCGCGACGGAAATCGGTCAGGGTGCCGACACCGTCTTCTCGCAAATGGTGGCAGAAACCGTGGGGGTTCCGGTCAGCGACGTTCGCGTTATTTCAACTCAAGATACCGACGTTACGCCGTTCGATCCCGGCGCATTT  >  NC_000913/3001734‑3001889

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