Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,485,171 | C→T | 21.3% | intergenic (+77/‑70) | purH → / → purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,485,171 | 0 | C | T | 21.3% | 66.0 / 14.0 | 38 | intergenic (+77/‑70) | purH/purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base C (11/18); new base T (3/5); total (14/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-01 |
TCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGC > NC_002947/5485027‑5485312 | tCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATTGTGGATTCACCGCCG‑CATTc < 7:185316/149‑1 (MQ=255) tCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTc < 5:303516/149‑1 (MQ=255) cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCgg < 2:91241/149‑1 (MQ=37) gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCgg < 5:136809/148‑1 (MQ=37) tgctgcCGACGAGGCCGGCATCGCGATGGCTTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGGCCGCGAATAAGGTGGTGGGTTGACCCGCG‑CATTAGCGGGc > 7:256922/1‑148 (MQ=255) ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGc < 3:22445/143‑1 (MQ=255) ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGc > 4:22445/1‑143 (MQ=255) cgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTc < 3:215101/149‑1 (MQ=255) tGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGgcgc < 8:177685/148‑1 (MQ=255) aCTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCCCGAATAAGATGGGGGGTTAACCGCCG‑CCTTAGCGGGCAACCCCCCGCCCCCAGTGGTAGGCGGCA‑GCCA‑‑Cgg > 2:298663/1‑149 (MQ=255) cACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCgg < 7:20116/100‑1 (MQ=17) cACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCgg > 8:20116/1‑100 (MQ=17) ccACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTcgc > 5:11485/1‑90 (MQ=255) ccACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTcgc < 6:11485/90‑1 (MQ=255) gcgTATCCCGTGGCAAGTCGAGATCCTGTGGGAGGGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGaaa < 7:219067/148‑1 (MQ=255) gcgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCGCATAACATGGTGGTTTAAACCCCGCCATAGCCGGGCAA‑CCCCCGCCCCCACAGTTGGGCGCCA‑GCCC‑‑CCCGGTCACATGAATACAGGGTTTTACaataaa > 8:219067/1‑143 (MQ=255) aGGCAAGTCGAGATCCTGTGGGCGCGGGCTTGCCCGCGAATACGAGGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcat < 8:256922/149‑1 (MQ=255) cAAGGCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCgg < 3:241020/149‑1 (MQ=255) gaTCCGGTGGGCGCGGGCTTGCCCGCGAATACTATGGTGTATTCACCGCCC‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTg < 2:10098/149‑1 (MQ=39) tCCTGTGGGCGCGGTGTTGCCCGATAATAGTATGGTGTATTCACCGCCT‑ATTTCGGGGGCATGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTgg < 1:323597/148‑1 (MQ=255) ttgGGCGCGGGTTTCCCCGCCAATAGGTTGTTGCTTTCACCGCCG‑CATTCCGGGGCATGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt < 3:221770/147‑1 (MQ=255) tgGGAGCGGGCTTGCCCGCTAATACGATGGTGTATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAACCCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt < 1:112623/148‑1 (MQ=38) ggAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCACCCGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa < 1:298663/149‑1 (MQ=255) aGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGAGGGCAAGCCCGCTCCCACAGTGTTAGGCGCAA‑GCCA‑‑CGGGATCACTTGAAACGAGGTTTTTACATGAAAGTTTTGATAATAGGCAGCGGTGGGCGTGAGCAc > 4:50684/1‑149 (MQ=38) gCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGtttt < 2:7792/106‑1 (MQ=255) gCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGtttt > 1:7792/1‑106 (MQ=255) ccGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa < 1:17295/134‑1 (MQ=25) ccGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa > 2:17295/1‑134 (MQ=25) ccGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAC‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCa > 4:97707/1‑121 (MQ=255) ccGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCa < 3:97707/121‑1 (MQ=255) gcgAATGCGGCGGTGAATCCACCATCG‑TATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTgg > 4:289308/1‑126 (MQ=25) gcgAATGCGGCGGTGAATCCACCATCG‑TATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTgg < 3:289308/126‑1 (MQ=25) ctaaTACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa < 3:50684/147‑1 (MQ=255) tACGATGGTGGATCCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc < 5:105081/149‑1 (MQ=255) ggATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCCCAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAc > 6:266307/1‑149 (MQ=255) aCCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGa > 4:136064/1‑68 (MQ=255) aCCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGa < 3:136064/68‑1 (MQ=255) aCCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcgc > 4:322490/1‑149 (MQ=255) ccg‑cATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCa > 4:132142/1‑135 (MQ=255) ccg‑cATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCa < 3:132142/135‑1 (MQ=255) | TCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGC > NC_002947/5485027‑5485312 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |