Questions tagged [variants]

The tag has no usage guidance.

Filter by
Sorted by
Tagged with
1 vote
0 answers
19 views

How to make a POST "vep/:species/region" request that is able to accomodate multiple alternate alleles for a single variant?

There is an R code example on Ensembl's website on how to execute a POST "vep/:species/region" request for variants that only have one alternate allele. The part that holds the variant data ...
Derk's user avatar
  • 11
2 votes
1 answer
48 views

Standardizing Variant IDs best practise

I'm currently pondering on the best approach to handle and standardize variant IDs within our department to mitigate the challenges associated with rsids and the potential ambiguity of just using ...
user27815's user avatar
  • 121
0 votes
0 answers
24 views

spliceAI not giving output value while running using vep (Variant Ensemble Predictor)

While generating the predictions using spliceAI from vep, I am not getting any values. I used this command to generate the output, ...
Sruthi's user avatar
  • 1
0 votes
0 answers
20 views

What statistical analysis should I perform on DNA variants from unequal number of case (34) and control (11) samples?

I have performed NGS on 34 patient samples and 11 controls. Now I have variants from both groups but there is obvious difference in number of variants. How can I compare variants of both groups ...
Saad Farooq's user avatar
4 votes
1 answer
221 views

How to get a GISAID account? I registered months ago, still no reply!

Inspired by amateur variant hunters, I would like to join the Pango lineage proposal community and help contribute to variant surveillance. However, I cannot seem to get access to GISAID, the platform ...
AppleBees's user avatar
4 votes
1 answer
57 views

Identifying somatic mutations in cell lines

I would like to identify the somatic mutations present in a cell line and characterise the genes that are potentially affected by those mutations. For example, are there oncogenes mutated in a ...
Macintosh's user avatar
  • 160
0 votes
0 answers
76 views

Combine gvcf for whole population or for different phenotypic groups?

I created 400 or so gvcf files (g.vcf.gz) using the gatk HaplotypeCaller. These files show populations for three distinct traits (HR, MR, HS). Now, I want to use gatk GenomicDBImport to combine gvcfs ...
Learner's user avatar
  • 51
2 votes
1 answer
693 views

Compare VCF files

I have two strains, X and Y, and I would like to compare the variants found in both. I have my own parent strain, Z, that was sequenced and assembled to be used as the reference genome. I first ...
rimo's user avatar
  • 923
1 vote
1 answer
60 views

Name this alternative INDEL format

I have recently been working with variants which are represented in a non-standard format. This alternative format uses a period instead of prefixing variants with a reference base. What is this ...
Moss Richardson's user avatar
0 votes
1 answer
70 views

syndip dataset for benchmark variant

I have a question about syndip dataset : https://github.com/lh3/CHM-eval . I'm struggling to find the syndip vcf. In the release ( https://github.com/lh3/CHM-eval/releases ), we have a file named : <...
Maryem Life's user avatar
0 votes
0 answers
158 views

Extract SNPs from multiple sequence alignment

I’m wondering if folks have recommendations for tools/scripts to extract SNP sites from a multiple sequence alignment of consensus bacterial genomes? Specifically, I am interested in a multiple ...
ksw's user avatar
  • 111
3 votes
1 answer
54 views

pathogeniticy prediction tool

I have a question related to the evaluation of variant pathogenicity in the context of human genetics. There are many predictors available, however, is there a specific predictor that takes into ...
Tomas Kucera's user avatar
1 vote
0 answers
39 views

HGVS python package : How to parse complex insertions?

In HGVS package, how to create more complex variants ? I just wonder how can I parse more complex insertions such as those : ...
jossefaz's user avatar
  • 111
3 votes
3 answers
243 views

Correspondence of SARS-CoV-2 annotations (Nextstrain clades - Pango lineages)

Is there any annotation file providing correspondence between Nextclade and Pangolin variant nomenclature/annotations, to annotate some SARS-CoV-2 genomes from Gisaid with both these? For the moment I ...
aechchiki's user avatar
  • 2,676
1 vote
1 answer
246 views

Dictionary and index of vcf for base recalibration step

I need suggestions in creating index and dictionary with vcf files. For the base recalibration step, I downloaded Homo_sapiens_assembly38.known_indels.vcf.gz from the given link: https://console.cloud....
Atom's user avatar
  • 11
2 votes
1 answer
52 views

Counting the co-occurrence of variants A and B in an aligned sequencing read

I need some help getting started on this project. To simplify we want to be able to quantify the occurrence of 3 variants on each sequencing read in an alignment file as a proxy measurement for ...
Bioreeb's user avatar
  • 21
0 votes
2 answers
627 views

Remove repetitive region vcf file using repeatmasker bed file [duplicate]

I have a 1000 genomes vcf file for chromosome 14. And I want to remove variants in repetitive regions flagged by repeatMasker. I have a bed file for those repeats. I downloaded that file from UCSC. ...
John's user avatar
  • 115
5 votes
0 answers
52 views

Parse RNA variant effect annotations ("r." format)

I've got annotations for splicing variants in a format like this (this is one variant): Variant: NM_004092.3:c.88+5G>A Effect: Retention; r.87_88ins1_88+10:p.(Ala31Glufs*23) I want to extract ...
Ron's user avatar
  • 51
1 vote
1 answer
1k views

Find common SNPs in multiple VCF files

I have 3 VCF files. A1.vcf A2.vcf A3.vcf I want to get the common SNPs that are present in all these three files. And output must be in vcf format. output: ...
Mendel's user avatar
  • 75
7 votes
1 answer
49 views

How to represent random sequence elements in SBOL?

In our lab, we work with synthetic biology components with partially random sequences (similar to work in directed evolution). So, for example, we have a plasmid design with several components, ...
David Ross's user avatar
1 vote
1 answer
279 views

genotyping or variant calling in R

Are there any good options for calling variants or even just a pileup from R? R is not ideal for this, but I'd like to integrate with other functions. I found deepSNV::bam2R which roughly does what I ...
burger's user avatar
  • 2,169
0 votes
1 answer
40 views

Inheritance annotation for nsSNPs

I am relatively new to SNP analysis. Is there a database to find annotation about the Mendelian inheritance of SNPs? I have a small list of nsSNPs (non-synonymous SNPs) and I need to find how they are ...
Hamza Umut Karakurt's user avatar
1 vote
0 answers
211 views

About getting rs id from chromosome and position

Hi I have a dataset with chromosome and position. VARIANT_ID chromosome position 17_26797415_147499 17 26797415 17_26797556_147500 17 26797556 How can I ...
susuauidikd's user avatar
2 votes
2 answers
452 views

Produce .vcf file of ALL mutations in .bam file

I got tasked with a problem, that i thought would be quite simple to solve, but turned out to be quite tricky. Our lab is running targeted mutagenesis experiments in yeast using crispr base editors. ...
Angelo Limeta's user avatar
1 vote
1 answer
68 views

How to compact variant data to their genes without bias?

I have a dataset of genes I am trying to collect data on from public databases, to use as features in machine learning. I am trying to take some features from UCSC genome browser (e.g. number of CpG ...
DN1's user avatar
  • 45
2 votes
1 answer
358 views

SNPEFF genome build with specific codon usage

Failing to build my own Arabidopsis thaliana reference with Mt using 'codon.Mitochondrial' and Pt using 'codon.Bacterial_and_Plant_Plastid' but the same commands works well with both lines with codon ...
splaisan's user avatar
  • 141
0 votes
2 answers
72 views

Q: What analyses can I perform with a completely phased genome assembly?

Phasing a genome is the process of determining which variants lie on which copy of each chromosome. For example, position 5,430,500 might have an A on the paternal chromosome, but a T base on the ...
conchoecia's user avatar
  • 3,141
1 vote
0 answers
65 views

Filtering out 8-oxoG-derived artefacts in variant calls from cfDNA

I have cfDNA samples which were sequenced (at about 300-1000X) on a 30-gene panel via a capture-seq protocol for the purposes of calling variants on circulating tumour DNA. However, the resulting ...
user36196's user avatar
  • 291