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.0 / NA	31	G729D (GGT→GAT)	xdhA	xanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/14);  total (17/14)

CGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGCCGATGCCATCA  >  NC_000913/3002386‑3002665                                                                                                                                                 |                                                                                                                                        cGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGtt                                                                                                                                     >  5:80843/1‑149 (MQ=255)  gACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTg                                                                                                                                    >  4:127920/1‑149 (MQ=255)                   cacaACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGcc                                                                                                                    >  1:95835/1‑148 (MQ=255)                   cacaACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGcc                                                                                                                    <  2:95835/148‑1 (MQ=255)                    acaACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGc                                                                                                                  >  1:268696/1‑149 (MQ=255)                           gAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGc                                                                                                            <  2:268696/148‑1 (MQ=255)                                    ttCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGtt                                                                                                  <  3:127920/149‑1 (MQ=255)                                       gtcgAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTtata                                                                                               >  8:181201/1‑149 (MQ=255)                                        tcgAAATCAATGAGCCGCAATCCGCATACGGCCATAAGTCACTGGGTGCGCCACCAATAATTCCTGTTGCCGCTGCTATTAGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTatat                                                                                              >  5:126675/1‑149 (MQ=255)                                           aaaTCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATg                                                                                             >  5:185233/1‑147 (MQ=255)                                             aTCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGtt                                                                                                  <  5:59141/140‑1 (MQ=255)                                             aTCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGtt                                                                                                  >  6:59141/1‑140 (MQ=255)                                                      ccGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTc                                                                                     >  5:181813/1‑144 (MQ=255)                                                      ccGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTc                                                                                     <  6:181813/144‑1 (MQ=255)                                                                     ggACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTg                                                                 <  1:109222/149‑1 (MQ=255)                                                                       aCATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAg                                                               <  6:126675/149‑1 (MQ=255)                                                                           aaGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaa                                                                                                                                 >  7:144393/1‑79 (MQ=255)                                                                           aaGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaa                                                                                                                                 <  8:144393/79‑1 (MQ=255)                                                                                          ccaccaATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGAtt                                            >  2:38985/1‑149 (MQ=255)                                                                                              caATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGc                                        >  4:204259/1‑149 (MQ=255)                                                                                                   atTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaatcaat                                                                                                                            >  5:229137/1‑60 (MQ=255)                                                                                                   atTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaatcaat                                                                                                                            <  6:229137/60‑1 (MQ=255)                                                                                                            gCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGAtt                                            >  7:119079/1‑131 (MQ=255)                                                                                                            gCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGAtt                                            <  8:119079/131‑1 (MQ=255)                                                                                                               gctgctATTCGTAACGCGGTGAAGATGGCTACCTATGTTGCAATCAATAAACTGCCGCTGACGACAAACAGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCgcag                       >  8:120881/1‑149 (MQ=255)                                                                                                                   ctATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGGTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAAc                   <  6:80843/149‑1 (MQ=255)                                                                                                                     aTTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAAccc                 >  2:249939/1‑149 (MQ=255)                                                                                                                         gTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGc             >  2:122977/1‑149 (MQ=255)                                                                                                                             cgcgGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGCCGAt         <  7:181201/149‑1 (MQ=255)                                                                                                                                 gTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGCCGATGCCa     <  1:122977/149‑1 (MQ=255)                                                                                                                                     aGATGGCTACCGATTTTGCAATCAATACACTGCCGCTGACGCCATAACGTTTATATGAAGAGTTCCATCTTGCAGTATTTATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAAACCTTGCCGATGCCATCa  <  6:185233/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                        CGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGCCGATGCCATCA  >  NC_000913/3002386‑3002665

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: ‑