Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,024,034 G→A V907I (GTC→ATC)  xdhD → fused putative xanthine/hypoxanthine oxidase: molybdopterin‑binding subunit and Fe‑S binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,024,0340GA94.4% 41.8 / NA 18V907I (GTC→ATC) xdhDfused putative xanthine/hypoxanthine oxidase: molybdopterin‑binding subunit and Fe‑S binding subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (7/10);  minor base C (1/0);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTC  >  NC_000913/3023794‑3024278
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
gAAAGCCGCCTTGATAAATTCTTCGCGTTGCGTGCCTGCGGTACACCGGTCAATCCAGAGTTAGTGGTGGGACAAATCTACGGTGCCCCCCGTCGTGCTATCTGCCGCAGTAGGAGCGAAGAGATCATGTATGACTCCGAAGGTTATCCGTTAACTCGTGATTTACGCAGGGTCTGCGCACCTAAAATGGGTGACATTCCTCGTGATTTCCTCGCTGTGCTGGTGCCGAGCGACGATAAAATCGGCCCGtt                                                                                                                                                                                                                                            <  1:151809/251‑1 (MQ=255)
                   tCTACGCCTTGCTGGACTGCGGGACACCGGTCAATCCAGAGTTATCGTTGGGACAAATCTACGGTGCCACCCTGCGTGATATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTCCCGCGCTGTGCTGGGGCCGAGCGACGATAAAATCGGCACGTTCGGGGCGAAATCGAtctc                                                                                                                                                                                                                          <  1:14707/250‑1 (MQ=255)
                                                               gCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAATCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCgg                                                                                                                                                                                             >  1:146168/1‑235 (MQ=255)
                                                                                            ccgatttATCGGCAACAGTATGAGCGAAGAGATCATGTATGACGCCGAAGGTCACCCGTGAACGCGTGATGTACGCAGTTACGGCGCACCGAAAATGGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGGGCCGAGCGACGATAAAATCGGCCCGTTCGGGGCTAAATCGATCGCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGAGTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTAcgc                                                                                                                                                 <  1:147091/244‑1 (MQ=255)
                                                                                                             gTATGAGCGAAGAGATCATTTATGACGGCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGGGCCGAGCGACGATAAAATCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCACGCGGCATCTCGTTACGCGAATGGCATTTCACACCg                                                                                                                               <  1:176065/251‑1 (MQ=255)
                                                                                                                                    gACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAATCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCg                                                                                                                                                               <  1:56716/196‑1 (MQ=255)
                                                                                                                                                                                     cGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAATCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCTGCGGCCTGTTTCTCTGTTTg                                                       <  1:167095/251‑1 (MQ=255)
                                                                                                                                                                                      gAAAATTGGTTACATGCCGCGTGATTTCCTCGCTGTGCGGGTGCCGCGCGACGAGAAAATCGGCCCGGTCGGGGGTAAAGCGTTTTCTGAAATCGGTGGAACAGGCGCTGCTCCGTCGATTGCTACCGCAATTCACGATGCATGCGGTATCTGGTTAc                                                                                                                                                   <  1:37695/158‑1 (MQ=255)
                                                                                                                                                                                                                 ccGCGCTGTGCTGGTGCCGAGCGACGATAAAATCGGCCCGTTCGGGGCGAAATCGATCTCGGCAATCGGTGTAACTGGCGCAGCTCCGGCGATTGCTCCCGCAATTCAc                                                                                                                                                                         >  1:201809/1‑109 (MQ=255)
                                                                                                                                                                                                                    cgcTGTGCTGGTGCCGCGCGCCGATAAACTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGCTTGCTACCGCACTTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATg                         >  1:70419/1‑250 (MQ=255)
                                                                                                                                                                                                                            tGGTGCCGAGCGACGATAAAATCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATg                                                                                                                                            <  1:121342/127‑1 (MQ=255)
                                                                                                                                                                                                                                gCCGAGCGACGATAAAATCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTtc            <  1:30436/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ccGAGCGACGATAAAATCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTtct           >  1:13690/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cGAGCGACGATAAAATCGGCCCGTTCGGGGCGAAATCGATCGCGGAAATCGGTGTAAAGGGCTCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGTCATATCACTCCGGAGAAAAGACTCGCTGCGCTGGAAAAAATAGAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGttt                                                        <  1:32670/205‑1 (MQ=255)
                                                                                                                                                                                                                                     gCGACGATAAAATCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCa         >  1:213891/1‑249 (MQ=255)
                                                                                                                                                                                                                                          naTAAAATCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTc  >  1:135456/2‑251 (MQ=255)
                                                                                                                                                                                                                                          gaTAAAATCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCCGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCCTCTGGTTACGCGAATGGCATTTCACACCGGTGAAAATACTCACTGCGCTGGACAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTc  >  1:112418/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gaTAAAATCGGCCCGTTCGGGGAGAAATGGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAACATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTc  >  1:181598/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
GAAATCCGCCTGGATAAATTCTACGCCTTGCTGGACTGCGGTACACCGGTCAATCCAGAGTTAGCGTTGGGACAAATCTACGGTGCCACCCTGCGTGCTATCGGCCACAGTATGAGCGAAGAGATCATTTATGACGCCGAAGGTCACCCGTTAACGCGTGATTTACGCAGTTACGGCGCACCGAAAATTGGTGACATTCCGCGTGATTTCCGCGCTGTGCTGGTGCCGAGCGACGATAAAGTCGGCCCGTTCGGGGCGAAATCGATCTCGGAAATCGGTGTAAATGGCGCAGCTCCGGCGATTGCTACCGCAATTCACGATGCATGCGGCATCTGGTTACGCGAATGGCATTTCACACCGGAGAAAATACTCACTGCGCTGGAAAAAATATAAATAATAAATGGCCGCTTCGGCGGCCTGTTTCTCTGTTTGAACGATGACTAAATAGTGACTGCGAAAATGAGATAAGCCGTTCTCAAGGCGTC  >  NC_000913/3023794‑3024278

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