Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,503,790 | A→G | 13.6% | intergenic (+39/‑59) | xanP → / → tusD | xanthine permease/sulfur transfer protein complex subunit TusD |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,503,790 | 0 | A | G | 13.6% | 99.5 / 9.0 | 44 | intergenic (+39/‑59) | xanP/tusD | xanthine permease/sulfur transfer protein complex subunit TusD |
Reads supporting (aligned to +/- strand): ref base A (16/22); new base G (4/2); total (20/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.87e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.10e-02 |
GCACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCGG > NC_002947/4503643‑4503937 | tctcaTGCCGGCGGCGATTCTCAAGGTGCTGGATGCGGGGGTGGCTTCCGGGGGGTTTTGTGCGCTGGTGCTGATCTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGcac < 4:60014/146‑1 (MQ=255) acaTGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGcac > 6:398643/1‑147 (MQ=255) ccGGCGGCGAGTCGCAATGTGCTGGAGGCGGGTGTTGCTACCGGGGTGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGATGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAgcgg < 4:238959/148‑1 (MQ=255) gcggcgATTCGCAATGTGCTGGAGGCGGGTGGGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGTGCAAgcgggcgc < 7:388297/149‑1 (MQ=255) gCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCttttt < 5:176211/148‑1 (MQ=255) cAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCGTGGACAAGCGGGGGGGGGGttttttt > 5:39544/1‑148 (MQ=255) actgtgCTGGAGTCGGGGGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGc < 1:371299/147‑1 (MQ=255) cTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTAt > 7:283372/1‑148 (MQ=255) tGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGATTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCa < 4:6435/149‑1 (MQ=255) ggtgtggggTGTGGCTACGGGGGGTTTTTTTGCGGTGTTGCTGAACTTTTTTTTGCCGGTGGTAAAGTATCTGGTTTTAACTTTTAGACTGAATTATATTAGCCCGCGCCCGCTTGAACACGGGGGGGCGGGCTTTTTTGCTTTATCAt < 8:132670/143‑1 (MQ=255) aGTCGGGTTTGTCTACCGGGGGGATTTGTGCGCTTGTGCTGAAGTGGTTCTTGCCGGAGGGCAAGGAGCAGTCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGATTGCACAAGCTGTCGCTGGCTTTTTTGCTTTATCATg < 8:232520/148‑1 (MQ=255) ggTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGTGCGCGGGCTTTTTTTCTTTATAATTTCAcc > 4:110593/1‑148 (MQ=255) cTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCttt < 3:331696/149‑1 (MQ=255) ccGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTg < 8:337797/147‑1 (MQ=255) ccGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCa > 3:298687/1‑149 (MQ=255) tgtgCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGcc < 5:398643/148‑1 (MQ=255) gTGCTTAACTGGTTCTTGACGGATGGCAAGGATCAGGCTTGAACCTTGATAATGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTc < 3:110593/148‑1 (MQ=255) gAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATc > 5:394709/1‑149 (MQ=255) cTGGTTCTTGCCGGAGGGCAAGGAGCAGTCTTGAACGTTGATACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTg < 1:24395/119‑1 (MQ=255) cTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCCCTTGCACAAGCCGGCGCCGGCTTTTTTGCTTTATCAAGGCACCATAACTTTg > 2:24395/1‑119 (MQ=255) gTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGgtg > 2:98023/1‑149 (MQ=255) ttCTTGCCGGAGGGCAAGGAGCAGGCTTTAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATcg‑‑cg < 5:50963/145‑1 (MQ=255) cTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTtc > 7:301508/1‑149 (MQ=255) ttGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGtt > 4:304149/1‑147 (MQ=255) ttGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGtt < 3:304149/147‑1 (MQ=255) ggAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCg < 8:283372/149‑1 (MQ=255) gAGTGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCGCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCgg < 7:300890/149‑1 (MQ=255) gAGGGCAAGGAGCAGGCTTGAACCTTGAGCCTGAATGACAGTAGCCCGCGCCCGCGTGCACCACCGCGCGCGGGCGTTTTTGCTTTAACCTGGACCCATTACTTTGCAAGAGTTAGCCAAGAAATTAGCGATCG‑‑AGTTGGTAGCCCCgg > 6:102319/1‑149 (MQ=255) ggAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTAGCTATCG‑‑CGGTGTTCTCCCCGGCCCATgtg > 2:121891/1‑147 (MQ=255) aGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCT‑CCCGGCCCATGCGCCctct < 4:298687/149‑1 (MQ=255) ttGAACCTTGAGACTGAATGAGAGTAGCCCGCCCCCCCATGCACAAGCGGGCGCGGGCTTTTTTGCATTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATAG‑‑AGGTGTTCTCCCCGGCCCATGCGCCCACATcg > 6:180514/1‑149 (MQ=255) gAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTcg < 1:98023/147‑1 (MQ=255) gAACCTTGAGACGTAATGAGCGTCGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTcgc < 5:102319/148‑1 (MQ=255) ttGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGAcgcg > 1:303408/1‑149 (MQ=255) gagaCTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGAcgcg < 5:180514/147‑1 (MQ=255) gagTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGcc < 6:394709/149‑1 (MQ=255) agTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTTTCATGGCACCATTCCTTTGCATGAGTTTGCCATGAAATTAGCTATAG‑‑CGGTGTTGTCCCCCGCCCAAGCGCCCGCTTAGCGACGCGCCCGGCGCGTAGCCg > 8:120261/1‑149 (MQ=255) cccgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCT‑‑CGGTGTTCTCCCCGGCCCATGCCCCCTCTTCGCCACGCGCCCTGCGCTTCGCCGAGCCg > 6:127225/1‑149 (MQ=255) ccgcgcccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCgg > 2:76431/1‑149 (MQ=255) cgcgcccgcTTGTGCAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGt > 5:90807/1‑149 (MQ=255) cgcgcccgcTTGTGCAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGt > 2:259999/1‑149 (MQ=255) gcgcccgcTTGTGCAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTg > 4:117391/1‑149 (MQ=255) gcgcccgcTTGTGCAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCAGAGGCGGTg > 7:263324/1‑149 (MQ=255) cccgcTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTg < 7:120261/149‑1 (MQ=255) cgcTTGTGCAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGc < 6:90807/149‑1 (MQ=255) cgcTTGCACAAGCGGGCGCGGGGTTTTTTTCTTTATCAAGGCACCATTCCTTTTCATGAGTTATCCATGAAAATCGCTATCGCCCTGTGTTCTCCCCGGCCCATGTGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTg < 1:121891/149‑1 (MQ=255) cacaAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCcggcgg < 2:303408/149‑1 (MQ=255) | GCACATGCCGGCGGCGATTCGCAATGTGCTGGAGTCGGGTGTGGCTACCGGGGGGATTTGTGCGCTGGTGCTGAACTGGTTCTTGCCGGAGGGCAAGGAGCAGGCTTGAACCTTGAGACTGAATGAGAGTAGCCCGCGCCCGCTTGCACAAGCGGGCGCGGGCTTTTTTGCTTTATCATGGCACCATTCCTTTGCATGAGTTAGCCATGAAATTCGCTATCG‑‑CGGTGTTCTCCCCGGCCCATGCGCCCTCTTCGCGACGCGCCCTGCGCTTCGCCGAGGCGGTGCTGGCCGGCGG > NC_002947/4503643‑4503937 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |