Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,485,153 | A→G | 16.4% | 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 | 16.4% | 93.1 / 9.8 | 45 | intergenic (+59/‑88) | purH/purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base A (17/19); new base G (4/3); total (21/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.98e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01 |
TCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCC > NC_002947/5485008‑5485300 | tCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGGGTTGCCC‑CGGATACGGtgg > 7:191746/1‑149 (MQ=255) gTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGAATACGATGGTGGAtt > 6:163329/1‑148 (MQ=255) tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGAATACGATGGTGGATTc < 6:83307/148‑1 (MQ=255) tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGAATACGATGGTGGATTc > 5:83307/1‑148 (MQ=255) gATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGGTTTGCCGCGCATAAGAAGGGGGTTTccc > 3:261169/1‑147 (MQ=255) aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGCATACGGTGGGGGATTGACCGCCGCCTTcg > 1:80926/1‑149 (MQ=255) tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGC‑GGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc < 7:286367/149‑1 (MQ=14) tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGAATACGATGGTGGATTCACCACCGCATTcgcg < 1:315726/149‑1 (MQ=255) gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg < 4:62226/148‑1 (MQ=37) tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGc > 8:200811/1‑149 (MQ=25) ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGc < 7:200811/147‑1 (MQ=25) ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGcc < 4:85979/148‑1 (MQ=255) cgacgaGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCg < 2:47541/149‑1 (MQ=255) cgacgaGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCCCGAATAAGAATGTGGATTCACCCCCGCAATAGCCGGCAAGCCCg > 2:269051/1‑149 (MQ=255) gCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTTCCCGCGAATACGATGGTGGATTAACCGCCGCAATAGCGGGCAAGCCCCCTCCcac > 6:155116/1‑149 (MQ=255) ccGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTTCCCGCGAATACGATTGGGGATTCACCGCCGCATTAGCGGGCCAGCCCGCGCCcacc > 1:276693/1‑148 (MQ=255) aTCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCGTGGCCCCCAATACGAAGGTGGGATAACAGCCGCCATAACGGGCAAGCCCCCTCCCCCAat‑gg > 4:211974/1‑147 (MQ=255) tCGCTATGGTCTTCACTGGAAT‑GCGCCACTTCCGCCACAAATTACGCGCATCCCGAGCCAAGTCGAGATC‑CTGTGGTCGC‑GGGCTTGCCCGCTAATACGATGGTGGATTCCCCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGT‑GTTc < 8:191746/149‑1 (MQ=255) cgATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg > 1:328217/1‑123 (MQ=255) cgATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg < 2:328217/123‑1 (MQ=255) aCTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACg < 5:337879/148‑1 (MQ=255) aCTTCCGCCACTAATTCCGCGGATCCCGAGGCATGTCTAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCTAATACTATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATc < 4:261169/149‑1 (MQ=38) ccGCCACTAATTACGCGGATCCCGAGGCACGTCGAGATC‑CTGGGGGAGC‑GGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAg < 5:155116/148‑1 (MQ=255) cgcgGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCTAATACGATGGGGGATTCACCGCCGAATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa < 1:269051/149‑1 (MQ=255) tCCCGAGCCAAGTCGATATC‑CTGGGGGAGC‑GGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg < 2:80926/149‑1 (MQ=255) aGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcat < 5:68943/149‑1 (MQ=255) ccAAGTCGATATC‑CTGTGGGCGC‑GGGCTTCCCCGCTAATAGGTTGGTGTATTCCCCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcatc < 3:154642/147‑1 (MQ=255) aaGTCGATATC‑CTGTGGGAGC‑GGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTCGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa < 1:267336/142‑1 (MQ=255) ccagaTC‑CTGTGGGTGG‑GGGCTGGCCCGCTAATACGAGGGTGTATTCACCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCg < 3:211974/147‑1 (MQ=255) tC‑CTGTGGGAGC‑GGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGc < 1:40601/149‑1 (MQ=255) tatgtgAGG‑GGGCTTGCCGGCTAATACGTTGGTGTATTCACCGCGGCATTCGGGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCTAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt < 5:3889/144‑1 (MQ=255) gC‑GGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACg < 5:163329/149‑1 (MQ=255) c‑GGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGaa > 8:165833/1‑95 (MQ=255) c‑GGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGaa < 7:165833/95‑1 (MQ=255) gggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGcc > 6:290946/1‑149 (MQ=255) ggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGccc > 3:124426/1‑149 (MQ=255) ggCTTGCCCGCGAATACGATGGTGGATTCACCGACGCATTCGCGGGCAAGCCCACTCCCACATT‑GTACGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATCAAAGTTTTGATNATCGGCAGCGGTGGCAGTGAGCACGca > 4:114544/1‑147 (MQ=255) cTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATGCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg < 7:64518/149‑1 (MQ=255) ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctgg > 8:177285/1‑149 (MQ=255) ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 6:94631/1‑149 (MQ=37) ccGCGAATACCATGGTGGATTCACCCCCGCATTCGCGGGCAAGCCCGCTCCCACACT‑GTTCGGCGCAACCCACGGGAACACTTGAATCCAGGTTTTGACAGGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcaag > 7:236610/1‑146 (MQ=255) cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg < 2:276693/149‑1 (MQ=255) aaTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAg > 1:18316/1‑149 (MQ=255) gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc > 3:175247/1‑149 (MQ=255) gATGGTGGATTCACAGCCGCATTCGCGGTCAAGCCAGCTCCCACAGT‑GTTCGGCCCAAGCCCAGGGATCAATTGAATCAAGTTTTTGACATGAAAGTTTTGATCATCGGCAGCGCTGGGCGTGAGCACGCCCTGCACTGGAAACTCTcc > 6:263075/1‑149 (MQ=255) | TCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCAT‑GCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGC‑GGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGT‑GTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCC > NC_002947/5485008‑5485300 |
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: ‑ |