Yum, tasty mutations...

mutation t@sting

 documentation

Prediction

disease causing

 Model: simple_aae, prob: 0.999999999233087 (classification due to ClinVar, real probability is shown anyway)      (explain)
Summary
  • amino acid sequence changed
  • known disease mutation at this position (HGMD CM021494)
  • known disease mutation: rs9148 (pathogenic)
  • protein features (might be) affected
  • splice site changes
 hyperlink
analysed issue analysis result
name of alteration no title
alteration (phys. location) chr15:73027508C>AN/A show variant in all transcripts   IGV
HGNC symbol BBS4
Ensembl transcript ID ENST00000539603
Genbank transcript ID N/A
UniProt peptide Q96RK4
alteration type single base exchange
alteration region CDS
DNA changes c.1055C>A
cDNA.1170C>A
g.48982C>A
AA changes A352E Score: 107 explain score(s)
position(s) of altered AA
if AA alteration in CDS		 352
frameshift no
known variant Reference ID: rs28938468
Allele 'A' was neither found in ExAC nor 1000G.
known disease mutation: rs9148 (pathogenic for Bardet-Biedl syndrome 4) dbSNP  NCBI variation viewer
known disease mutation at this position, please check HGMD for details (HGMD ID CM021494)

known disease mutation at this position, please check HGMD for details (HGMD ID CM021494)
known disease mutation at this position, please check HGMD for details (HGMD ID CM021494)
regulatory features DNase1, Open Chromatin, DNase1 Hypersensitive Site
Gene Associated, Regulatory Feature, Gene associated regulatory feature
H3K36me3, Histone, Histone 3 Lysine 36 Tri-Methylation
phyloP / phastCons
PhyloPPhastCons
(flanking)3.1371
5.9541
(flanking)-2.6980.001
explain score(s) and/or inspect your position(s) in in UCSC Genome Browser
splice sites
effectgDNA positionscorewt detection sequence exon-intron border
Donor increased48978wt: 0.21 / mu: 0.99wt: ACGCAGAAGCAGTCC
mu: ACGCAGAAGAAGTCC		 GCAG|aagc
distance from splice site 16
Kozak consensus sequence altered? N/A
conservation
protein level for non-synonymous changes
speciesmatchgeneaaalignment
Human      352DIENAKRAYAEAVHLDKCNPLVNL
mutated  not conserved    352DIENAKRAYAEEVHLDKCNP
Ptroglodytes  all identical  ENSPTRG00000007255  372DIENAKRAYAEAVHLDKCNPLVN
Mmulatta  all identical  ENSMMUG00000013685  364DIENAKRAYAEAVHLDKCNPLVN
Fcatus  no alignment  ENSFCAG00000012249  n/a
Mmusculus  all identical  ENSMUSG00000025235  364DIENARRAYVEAVRLDKCNPLVN
Ggallus  all identical  ENSGALG00000001798  342DIENAKRSYEQAVALDKCNPLVN
Trubripes  no homologue    
Drerio  all identical  ENSDARG00000063522  364DAENARRSYEQAVQIDESSPLVN
Dmelanogaster  all identical  FBgn0033578  377DMENAFVALERASSMATGQQ
Celegans  all conserved  F58A4.14  395DDKNARLAYKKSIELKK-NPSTI
Xtropicalis  all identical  ENSXETG00000005200  364DIENAKSSYQQAASLDQ
protein features
start (aa)end (aa)featuredetails 
338519REGIONRequired for localization to centrosomes.lost
339371REPEATTPR 9.lost
373408REPEATTPR 10.might get lost (downstream of altered splice site)
length of protein normal
AA sequence altered yes
position of stopcodon in wt / mu CDS 1524 / 1524
position (AA) of stopcodon in wt / mu AA sequence 508 / 508
position of stopcodon in wt / mu cDNA 1639 / 1639
poly(A) signal N/A
conservation
nucleotide level for all changes - no scoring up to now		 N/A
position of start ATG in wt / mu cDNA 116 / 116
chromosome 15
strand 1
last intron/exon boundary 1530
theoretical NMD boundary in CDS 1364
length of CDS 1524
coding sequence (CDS) position 1055
cDNA position
(for ins/del: last normal base / first normal base)		 1170
gDNA position
(for ins/del: last normal base / first normal base)		 48982
chromosomal position
(for ins/del: last normal base / first normal base)		 73027508
original gDNA sequence snippet CAAGAGAGCCTACGCAGAAGCAGTCCACCTGGATAAGTATG
altered gDNA sequence snippet CAAGAGAGCCTACGCAGAAGAAGTCCACCTGGATAAGTATG
original cDNA sequence snippet CAAGAGAGCCTACGCAGAAGCAGTCCACCTGGATAAGTGTA
altered cDNA sequence snippet CAAGAGAGCCTACGCAGAAGAAGTCCACCTGGATAAGTGTA
wildtype AA sequence MIWNTRNSTY IFTAPEFPIL EKQNWLIHLH YIRKDYEACK AVIKEQLQET QGLCEYAIYV
QALIFRLEGN IQESLELFQT CAVLSPQSAD NLKQVARSLF LLGKHKAAIE VYNEAAKLNQ
KDWEISHNLG VCYIYLKQFN KAQDQLHNAL NLNRHDLTYI MLGKIHLLEG DLDKAIEVYK
KAVEFSPENT ELLTTLGLLY LQLGIYQKAF EHLGNALTYD PTNYKAILAA GSMMQTHGDF
DVALTKYRVV ACAVPESPPL WNNIGMCFFG KKKYVAAISC LKRANYLAPF DWKILYNLGL
VHLTMQQYAS AFHFLSAAIN FQPKMGELYM LLAVALTNLE DIENAKRAYA EAVHLDKCNP
LVNLNYAVLL YNQGEKKNAL AQYQEMEKKV SLLKDNSSLE FDSEMVEMAQ KLGAALQVGE
ALVWTKPVKD PKSKHQTTST SKPASFQQPL GSNQALGQAM SSAAAYRTLP SGAGGTSQFT
KPPSLPLEPE PAVESSPTET SEQIREK*
mutated AA sequence MIWNTRNSTY IFTAPEFPIL EKQNWLIHLH YIRKDYEACK AVIKEQLQET QGLCEYAIYV
QALIFRLEGN IQESLELFQT CAVLSPQSAD NLKQVARSLF LLGKHKAAIE VYNEAAKLNQ
KDWEISHNLG VCYIYLKQFN KAQDQLHNAL NLNRHDLTYI MLGKIHLLEG DLDKAIEVYK
KAVEFSPENT ELLTTLGLLY LQLGIYQKAF EHLGNALTYD PTNYKAILAA GSMMQTHGDF
DVALTKYRVV ACAVPESPPL WNNIGMCFFG KKKYVAAISC LKRANYLAPF DWKILYNLGL
VHLTMQQYAS AFHFLSAAIN FQPKMGELYM LLAVALTNLE DIENAKRAYA EEVHLDKCNP
LVNLNYAVLL YNQGEKKNAL AQYQEMEKKV SLLKDNSSLE FDSEMVEMAQ KLGAALQVGE
ALVWTKPVKD PKSKHQTTST SKPASFQQPL GSNQALGQAM SSAAAYRTLP SGAGGTSQFT
KPPSLPLEPE PAVESSPTET SEQIREK*
speed 0.41 s
All positions are in basepairs (bp) if not explicitly stated differently.
AA/aa: amino acid; CDS: coding sequence; mu: mutated; NMD: nonsense-mediated mRNA decay; nt: nucleotide; wt: wildtype; TGP: 1000 Genomes Project