Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,485,153 | A→G | 27.6% | intergenic (+59/‑88) | 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,153 | 0 | A | G | 27.6% | 55.0 / 22.6 | 41 | intergenic (+59/‑88) | purH/purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base A (10/19); new base G (6/5); total (16/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.95e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-01 |
TTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTC > NC_002947/5485007‑5485297 | ttCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtg < 4:210743/149‑1 (MQ=255) tGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGa < 5:147724/149‑1 (MQ=255) 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) | TTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCG‑CATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAA‑GCCA‑‑CGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTC > NC_002947/5485007‑5485297 |
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: ‑ |