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_0009133,003,4731.G100.0% 48.4 / NA 15A287G (GCG→GGG) xdhBputative xanthine dehydrogenase FAD‑binding subunit XdhB
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (10/5);  total (10/5)

CAAACGCTGGAAGCTATCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTA  >  NC_000913/3003337‑3003602
                                                                                                                                         |                                                                                                                                 
cAAACGCTGGAAGCTATCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCggg                                                                                                                                <  1:14176/140‑1 (MQ=255)
                         cTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGc                                                                                                        <  2:158921/140‑1 (MQ=255)
                          tGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCg                                                                                                       <  1:390784/140‑1 (MQ=255)
                                                               gCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGGAAATTGCAATGAAACACAGCGAAACAATTACCATCGAATGCACCCTTAACGGGGTGcc                                                                  >  1:23636/1‑140 (MQ=255)
                                                                          gTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTc                                                       >  1:524665/1‑140 (MQ=255)
                                                                             tCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCAcg                                                    >  1:279647/1‑140 (MQ=255)
                                                                                          aTCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCAc                                       >  2:68055/1‑140 (MQ=255)
                                                                                                  gATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGctct                               <  2:524665/140‑1 (MQ=255)
                                                                                                    tGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCAcg                                      >  1:364599/1‑131 (MQ=255)
                                                                                                    tGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCAcg                                      <  2:364599/131‑1 (MQ=255)
                                                                                                    tGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCg                             >  1:325654/1‑140 (MQ=255)
                                                                                                      aCCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGa                           >  2:526132/1‑140 (MQ=255)
                                                                                                          aaaaaGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTa                       >  2:361055/1‑140 (MQ=255)
                                                                                                                     aGCGAAGCCGTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACa            >  1:431385/1‑140 (MQ=255)
                                                                                                                              gTCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATACCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCt   <  1:241024/140‑1 (MQ=255)
                                                                                                                               tCGCCGCGGCGGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTa  >  1:336122/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                 
CAAACGCTGGAAGCTATCAGCGAATCTGTCCTGCAAGATGTCGCCCCGCGTTCTTCATGGCGGGCCAGTAAAGAGTTTCGTCTGCATCTCATCCAGACGATGACCAAAAAAGTGATTAGCGAAGCCGTCGCCGCGGCGGGGGGAAAATTGCAATGAATCACAGCGAAACAATTACCATCGAATGCACCATTAACGGGATGCCTTTTCAGCTTCACGCCGCACCAGGCACGCCGCTCTCGGAATTACTCCGCGAACAAGGACTGCTA  >  NC_000913/3003337‑3003602

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