Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,002,530	G→A	G729D (GGT→GAT)	xdhA →	xanthine dehydrogenase, molybdenum binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,002,530	0	G	A	100.0%	110.6 / NA	32	G729D (GGT→GAT)	xdhA	xanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/17);  total (15/17)

TGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCA  >  NC_000913/3002392‑3002643                                                                                                                                           |                                                                                                                  tGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaatc                                                                                                         <  8:71227/149‑1 (MQ=255)                    ggAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGc                                                                                      >  5:34374/1‑149 (MQ=255)                                     aaaTCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaatcaat                                                                                                      >  3:2074/1‑116 (MQ=255)                                     aaaTCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaatcaat                                                                                                      <  4:2074/116‑1 (MQ=255)                                           aTGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTc                                                               >  1:56710/1‑149 (MQ=255)                                              aGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCAt                                                            >  2:102594/1‑149 (MQ=255)                                                  gCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGt                                                                             >  5:12786/1‑128 (MQ=255)                                                  gCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGt                                                                             <  6:12786/128‑1 (MQ=255)                                                      tCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAAc                                                                                >  8:118332/1‑121 (MQ=255)                                                      tCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAAc                                                                                <  7:118332/121‑1 (MQ=255)                                                        cGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGgat                                                  >  6:17623/1‑149 (MQ=255)                                                              cGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGgattgatt                                             <  6:34374/148‑1 (MQ=255)                                                              cGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGAttt                                            <  5:17623/149‑1 (MQ=255)                                                                 aCATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAg                                         <  1:102594/149‑1 (MQ=255)                                                                                tGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGttt                          <  2:56710/149‑1 (MQ=255)                                                                                   gCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATg                                                                       >  3:99667/1‑101 (MQ=255)                                                                                   gCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATg                                                                       <  4:99667/101‑1 (MQ=255)                                                                                      accaATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGt                                                                             <  6:114272/92‑1 (MQ=255)                                                                                      accaATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGt                                                                             >  5:114272/1‑92 (MQ=255)                                                                                        caATAATTGCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACTGATGTTGCAATCATTACACTGCCGCTGACTCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGc                  <  7:28406/149‑1 (MQ=255)                                                                                                               aTTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCgca  <  1:151977/142‑1 (MQ=255)                                                                                                               aTTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCgca  >  2:151977/1‑142 (MQ=255)                                                                                                                ttCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAAcatc                                <  2:91866/111‑1 (MQ=255)                                                                                                                ttCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAAcatc                                >  1:91866/1‑111 (MQ=255)                                                                                                                ttCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGctt                >  4:145787/1‑127 (MQ=255)                                                                                                                ttCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGctt                <  3:145787/127‑1 (MQ=255)                                                                                                                      aCGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTatat                                                                        <  5:13696/65‑1 (MQ=255)                                                                                                                      aCGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTatat                                                                        >  6:13696/1‑65 (MQ=255)                                                                                                                          ggTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTa            >  6:109737/1‑121 (MQ=255)                                                                                                                          ggTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTa            <  5:109737/121‑1 (MQ=255)                                                                                                                                      tACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGG‑TATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCgca  <  4:56180/118‑1 (MQ=255)                                                                                                                                      tACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGG‑TATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCgca  >  3:56180/1‑118 (MQ=255)                                                                                                                                           |                                                                                                                  TGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCA  >  NC_000913/3002392‑3002643

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