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

ATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGCCGATGCCATCAACCTGCT  >  NC_000913/3002391‑3002672                                                                                                                                            |                                                                                                                                               aTGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACTCGGTGAAGATGGCTACCGATGTTGcaat                                                                                                                                       <  5:57023/149‑1 (MQ=255)         tCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATacac                                                                                                                               >  7:69543/1‑149 (MQ=255)             ccACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTACCGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGcc                                                                                                                           >  6:166503/1‑149 (MQ=255)               acaACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaatcaa                                                                                                                                    >  8:85889/1‑138 (MQ=255)               acaACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaatcaa                                                                                                                                    <  7:85889/138‑1 (MQ=255)                  aCTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGa                                                                                                                      >  2:214899/1‑149 (MQ=255)                    tGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaa                                                                                                                                        >  5:104260/1‑129 (MQ=255)                    tGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaa                                                                                                                                        <  6:104260/129‑1 (MQ=255)                      gAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGcc                                                                                                                  <  2:15486/149‑1 (MQ=255)                                  gtcgAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTtata                                                                                                      <  8:69543/149‑1 (MQ=255)                                  gtcgAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTtata                                                                                                      >  6:66433/1‑149 (MQ=255)                                   tcgAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTatat                                                                                                     >  5:109227/1‑149 (MQ=255)                                   tcgAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTatat                                                                                                     <  1:41889/149‑1 (MQ=255)                                    cgAAATCAATGAGCCGCAATCCGAATACGGACATAAGTCACTGTTTGAGCCACCAATAATTCCTGTTGCCGCAGATATTCTTATCGCGGTTAAGATGGCTACCGATGTTGCAATCAATACAATGCCGCTGACGCCAAAAAGGTAatat                                                                                                     >  7:170293/1‑148 (MQ=255)                                     gAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGCTGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATg                                                                                                    >  7:68408/1‑148 (MQ=255)                                       aaTCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATa                                                                                                                                  <  3:98481/116‑1 (MQ=255)                                       aaTCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATa                                                                                                                                  >  4:98481/1‑116 (MQ=255)                                           aaTGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGtt                                                                                             >  7:112342/1‑149 (MQ=255)                                            aTGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTc                                                                                            <  2:204314/149‑1 (MQ=255)                                                gCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGaa                                                                                                  <  6:68942/139‑1 (MQ=255)                                                gCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGaa                                                                                                  >  5:68942/1‑139 (MQ=255)                                                gCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATc                                                                                        >  4:140343/1‑149 (MQ=255)                                                gCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGGGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATc                                                                                        <  8:106389/149‑1 (MQ=255)                                                 ccGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCt                                                                                       <  6:109227/149‑1 (MQ=255)                                                  cGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTg                                                                                      >  5:130400/1‑149 (MQ=255)                                                   gCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTACCGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTgg                                                                                     <  5:166503/149‑1 (MQ=255)                                                      aTCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAg                                                                                  <  3:140343/149‑1 (MQ=255)                                                              aCGGACATAAGCCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTtata                                                                                                      <  4:36047/121‑1 (MQ=255)                                                              aCGGACATAAGCCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTtata                                                                                                      >  3:36047/1‑121 (MQ=255)                                                               cGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGAttt                                                                         <  8:112342/149‑1 (MQ=255)                                                                 gACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGa                                                                       <  3:76900/149‑1 (MQ=255)                                                                      aaGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGAt                                                                   >  7:54523/1‑148 (MQ=255)                                                                               ggTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAAcatcat                                                           <  1:214899/147‑1 (MQ=255)                                                                               ggTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGt                                                         >  6:65317/1‑149 (MQ=255)                                                                               ggTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGt                                                         <  5:65317/149‑1 (MQ=255)                                                                                ttGAGCCACCAATAATTCCTGTGGCCGCTGCTATTCGTTACTAGGTTAAGATTGATTGCGATGTTGCAATCAATACACTGACGCTGACTCCAAAAAGGTTATATCTAGAGTTCCATCTGGCATTTTTGATTTGAGGATAACATAATGtt                                                        <  8:170293/148‑1 (MQ=255)                                                                                       accaATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGAtttt                                                 >  6:106422/1‑149 (MQ=255)                                                                                        ccaATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAAc                                                                >  6:42162/1‑133 (MQ=255)                                                                                        ccaATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAAc                                                                <  5:42162/133‑1 (MQ=255)                                                                                        ccaATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTg                                                >  2:34436/1‑149 (MQ=255)                                                                                         caATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGc                                               <  5:34124/149‑1 (MQ=255)                                                                                          aataatTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGCCGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGct                                              >  4:16083/1‑149 (MQ=255)                                                                                          aataatTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGct                                              >  3:101452/1‑149 (MQ=255)                                                                                                   tGTTGCCGCTGATATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTACATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCa                                      >  3:31812/1‑148 (MQ=255)                                                                                                        ccGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCgcgc                                <  1:34436/149‑1 (MQ=255)                                                                                                        ccGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCgcgc                                <  5:66433/149‑1 (MQ=255)                                                                                                         cGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCgca                               <  4:101452/149‑1 (MQ=255)                                                                                                          gctgctATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCgcag                              <  8:54523/149‑1 (MQ=255)                                                                                                                 ttCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCt                       >  7:131431/1‑149 (MQ=255)                                                                                                                    gTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGc                    <  5:106422/149‑1 (MQ=255)                                                                                                                    gTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGc                    <  6:130400/149‑1 (MQ=255)                                                                                                                         gcgGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGCCGATg               <  2:172319/149‑1 (MQ=255)                                                                                                                                  aTGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCtt                      >  4:169216/1‑133 (MQ=255)                                                                                                                                  aTGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCtt                      <  3:169216/133‑1 (MQ=255)                                                                                                                                      cTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTtat                                                                                                       >  2:217366/1‑48 (MQ=255)                                                                                                                                      cTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTtat                                                                                                       <  1:217366/48‑1 (MQ=255)                                                                                                                                      cTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGCCGATGCCATCAACctgct  <  4:31812/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                               ATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGCCGATGCCATCAACCTGCT  >  NC_000913/3002391‑3002672

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