Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,485,162 | T→C | 27.5% | intergenic (+68/‑79) | 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,162 | 0 | T | C | 27.5% | 50.8 / 22.3 | 40 | intergenic (+68/‑79) | purH/purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base T (12/17); new base C (5/6); total (17/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01 |
TGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAC > NC_002947/5485016‑5485306 | tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCa > 7:177685/1‑149 (MQ=255) gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccat > 1:68692/1‑147 (MQ=38) tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgccg > 6:303516/1‑149 (MQ=255) tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGAATACGATGGTGGATTAAccgccg > 1:292825/1‑149 (MQ=255) gAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcg > 1:91241/1‑144 (MQ=38) 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 > 4:22445/1‑143 (MQ=255) ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGc < 3:22445/143‑1 (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 > 2:17295/1‑134 (MQ=25) ccGCGAATGCGGCGGTGAATCCACCATCG‑TATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa < 1:17295/134‑1 (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 < 3:289308/126‑1 (MQ=25) gcgAATGCGGCGGTGAATCCACCATCG‑TATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTgg > 4:289308/1‑126 (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) | TGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAC > NC_002947/5485016‑5485306 |
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: ‑ |