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_0009133,002,5300GA100.0% 87.9 / NA 26G729D (GGT→GAT) xdhAxanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/13);  total (13/13)

CGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGC  >  NC_000913/3002386‑3002654
                                                                                                                                                |                                                                                                                            
cGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGtt                                                                                                                          <  1:53983/149‑1 (MQ=255)
                   acaACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGc                                                                                                       >  5:184930/1‑149 (MQ=255)
                       cTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGAc                                                                                                   <  2:153580/149‑1 (MQ=255)
                                                   aGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCAt                                                                       <  5:61860/149‑1 (MQ=255)
                                                       gCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTgg                                                                   >  5:2582/1‑149 (MQ=255)
                                                                 tACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGgattgat                                                         <  6:2582/149‑1 (MQ=255)
                                                                       cATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACTCTGCCGCTGACGCCAAAACGGTtat                                                                                     >  6:129247/1‑115 (MQ=255)
                                                                       cATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACTCTGCCGCTGACGCCAAAACGGTtat                                                                                     <  5:129247/115‑1 (MQ=255)
                                                                                   ggTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGt                                       >  3:46323/1‑149 (MQ=255)
                                                                                        gCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTg                                                                                                     >  3:98379/1‑82 (MQ=255)
                                                                                        gCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTg                                                                                                     <  4:98379/82‑1 (MQ=255)
                                                                                           accaATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGgattgatt                                                        >  1:148546/1‑124 (MQ=255)
                                                                                           accaATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGgattgatt                                                        <  2:148546/124‑1 (MQ=255)
                                                                                            ccaATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTc                                                                          >  6:59822/1‑105 (MQ=255)
                                                                                               ataatTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaa                                                                                                                      <  2:8432/58‑1 (MQ=255)
                                                                                               ataatTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaa                                                                                                                      >  1:8432/1‑58 (MQ=255)
                                                                                               ataatTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGa                                                                                                    <  7:103255/76‑1 (MQ=255)
                                                                                               ataatTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGa                                                                                                    >  8:103255/1‑76 (MQ=255)
                                                                                                        gTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGgattgat                                                         <  2:86423/110‑1 (MQ=255)
                                                                                                        gTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGgattgat                                                         >  1:86423/1‑110 (MQ=255)
                                                                                                         ttGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAATGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGAt                                                 <  1:94507/117‑1 (MQ=255)
                                                                                                         ttGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAAAGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGAt                                                 >  2:94507/1‑117 (MQ=255)
                                                                                                                    aTTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAAccc      >  2:121727/1‑149 (MQ=255)
                                                                                                                        gTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGc  <  4:46323/149‑1 (MQ=255)
                                                                                                                                          cTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTg                                                                    >  3:55003/1‑65 (MQ=255)
                                                                                                                                          cTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTg                                                                    <  4:55003/65‑1 (MQ=255)
                                                                                                                                                |                                                                                                                            
CGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGC  >  NC_000913/3002386‑3002654

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