How to create bed file from list of variants with indels and bp values
2 answers
Hope this is useful, if you have the original VCF file:
vcf2bed --snvs < variants.vcf > snvs.bed
vcf2bed --insertions < variants.vcf > insertions.bed
vcf2bed --deletions < variants.vcf > deletions.bed
Ref.: https://bedops.readthedocs.io/en/latest/content/reference/file-management/conversion/vcf2bed.html
• 0 views
•
link
written
by
Alex Reynolds
Log in to answer this question.
More posts like this
-
Inherited pipeline uses RNA-seq variant calls for WASP filtering—shouldn’t it be external SNPs?
written by jonas.andersson •
-
Where do I gather indel and frameshift information in SnpEff output?
written by mgranada3 •
-
How to draw genome atlas on Linux?
written by marongiu.luigi
-
Identify Alternative Splicing Event using exon and junction count table
written by rhonddaskl •
-
KEGG analysis in Hordeum
written by martta95 •
-
How to export the result contents (analyzed part only) of clr function to the CSV file ? (The micro…
written by ohtang7 •
-
Driver Phylogenetic Trees of TRACERx project
written by Vasei •
-
Coverage plot/histogram of a complete genome
written by Michal Nevo •
-
Tutorial: Extracting genomic feature sequences from GTF/GFF files with AGAT
written by Juke34
-
vcftools --TajimaD does not start at 0 and misses variants
written by grey •
A bed file is a white space delimited file with three required fields - CHR, START and END. You have two fields from which the three can be derived. What is your problem? I don't see any
+, only-s signifying an insertion of deletion like VCF pre 4.2 convention.Should I just subtract or add bases for the Pos_End? Right, they have not given + or -. I wanted to search these variants in my VCF file.
Look up some of these on dbSNP or gnomAD to understand how base positioning works - you'll need some math to calculate the End Pos, and some more math to account for the
-Indels.For example,
2:47803501delC(C>-) is in reality2:47803500AC>A(AC>A). (https://gnomad.broadinstitute.org/variant/2-47803500-AC-A?dataset=gnomad_r3)You can of course skip this step if standardization in not something you need. You could do just the End Pos calculation using Start Pos and the number of REF bases (I'll leave calculating the exact formula to you).