Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,003,479	(G)6→7	coding (866/879 nt)	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,473	1	.	G	100.0%	89.4 / NA	27	A287G (GCG→GGG)	xdhB	putative xanthine dehydrogenase FAD‑binding subunit XdhB
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (17/10);  total (17/10)

GCAAACGCTGGAAGCTATCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGC‑GGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTAAG  >  NC_000913/3003336‑3003604                                                                                                                                           |                                                                                                                                    gCAAACGCTGGAAGCTATCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGGCAGTAAAGAGCTTCGACTGCATTTTATCCAGACGGTGAGCAAAAAAGTGGTTAGCGAAGCCCCCACCGcggc‑gg                                                                                                                                   >  2:546767/1‑140 (MQ=255)  cAAACGCTGGAAGCTATCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGC‑ggg                                                                                                                                  <  2:5855/140‑1 (MQ=255)     aCGCTGGAAGCTATCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGC‑gggggg                                                                                                                               <  1:614619/140‑1 (MQ=255)           gAAGCTATCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGGCAGGACAGAGGTTCGTCTGGGTCTCAACCAGAAGAACCCCCAAAAAGAGAATAACGGAGACGGCGCCGCGGC‑ggggggggaaaa                                                                                                                         >  2:550128/1‑134 (MQ=255)                          tCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGGAAAATTCAATGAAACACAg                                                                                                           >  2:804474/1‑140 (MQ=255)                                           cGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGGTTCGTCTGCATCTCATCCAGACGATGACCCAAAAAGTGGTTAGCGAAACCGTCGCCACGGCGGGGGGGAAAAACTCAATTAACTACAGCGAAAAAACTaaccacc                                                                                          >  2:364616/1‑134 (MQ=255)                                                 gcgTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGGAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCa                                                                                    >  1:806814/1‑140 (MQ=255)                                                  cgTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGGAAAATGCAAAGAATCACAGCGGAACAAATACCATCCaaaggac                                                                                   >  2:547520/1‑135 (MQ=255)                                                                 gCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGGAAAATGCAATGAAACACAGCGAAAAAATTACCATCGAATGCACCCTTAAAGGGGTGcc                                                                    >  2:796758/1‑140 (MQ=255)                                                                           agTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCtt                                                          >  2:462562/1‑140 (MQ=255)                                                                             tttCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCa                                                        >  1:627439/1‑140 (MQ=255)                                                                                cGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCAcgc                                                     >  2:169137/1‑140 (MQ=255)                                                                                   cTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCAcgccgc                                                  <  1:804474/140‑1 (MQ=255)                                                                                    tGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCa                                                 >  2:450920/1‑140 (MQ=255)                                                                                             tCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAGTGCACCATGAACTGGAGGCCTTTTCTGCTTAACGCCGCACCAGGCAcg                                        >  2:246907/1‑140 (MQ=255)                                                                                               cAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCAcgcc                                      <  2:806814/140‑1 (MQ=255)                                                                                                  aCGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGct                                   >  1:458215/1‑140 (MQ=255)                                                                                                    gaggACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAACCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGctct                                 <  1:207353/137‑1 (MQ=255)                                                                                                    gATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGctct                                 <  1:547520/140‑1 (MQ=255)                                                                                                    gATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGctct                                 <  1:169137/140‑1 (MQ=255)                                                                                                      tGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCg                               >  2:595134/1‑140 (MQ=255)                                                                                                      tGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCg                               >  2:70672/1‑140 (MQ=255)                                                                                                      tGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCg                               >  2:735566/1‑140 (MQ=255)                                                                                                        aCCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCCAATGCACCATTAGCGGGAGGCGTTTGCAGTTTTACGCCTCATCTGGCACGCCGCTTTCGGa                             >  2:212613/1‑140 (MQ=255)                                                                                                               aaGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTc                      <  1:450920/140‑1 (MQ=255)                                                                                                                 gTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAACCAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCcg                    <  1:735566/140‑1 (MQ=255)                                                                                                                      tAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAAc               <  1:70672/140‑1 (MQ=255)                                                                                                                      tAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAAc               <  1:483450/140‑1 (MQ=255)                                                                                                                         caaaGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACAATTAAGGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAg            <  1:550128/138‑1 (MQ=255)                                                                                                                           aaGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGa          >  2:115230/1‑140 (MQ=255)                                                                                                                                 tCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCAGCACCAGGTACGACGCACTCGGAATTACTCGGCGAACAAGGACTGCTa    >  2:219290/1‑140 (MQ=255)                                                                                                                                 tCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACAGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTa    >  2:790026/1‑140 (MQ=255)                                                                                                                                   gccgcGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTAAg  >  1:793372/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                    GCAAACGCTGGAAGCTATCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGC‑GGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTAAG  >  NC_000913/3003336‑3003604

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑