Parse Non Hashed VCF info and Store It in csv File using BASH

Can anyone help me with such a problem?

I want to parse a VCF file from 1000 human genome projects to write the non hashed line (##) escaping this meta info and write the information into comma-separated file >> 1000_genome,CSV

parse vcf bash

• 1,086 views • link

updated 4.8 years ago by Pierre Lindenbaum 17K • | written 4.8 years ago by Ramy Ali 0

man grep

man tr

this meta info and write the information into comma-separated file

there is a good chance that this is a bad idea.

• 0 views • link

written 4.8 years ago by Pierre Lindenbaum 17K

0 answers

No answers yet.

More posts like this

problem with bcftools syntax
written 3.1 years ago by Barista 1

Hi all! I am having difficulty with creating a bcftools command. I have a .vcf.gz file downloaded from the 1000G site and a csv file …
Why is NA18498 sample missing from 1000 genome GRCh38 vcf file?
written 3.6 years ago by octpus616 14

Hi I am trying to compare the 1000 genome phase3 GRCh19 vcf file and the GRCh38 （low coverage） vcf file. I noticed that one sample, …
Parsing a GTF file with BCBio-gff: AttributeError
written 3.9 years ago by Flavia 0

I'm trying to parse a gtf file using this code: ```py from BCBio import GFF gtf_rec = [] in_file = 'cuffcmp.combined.gtf' out_file = 'extract.gtf' with …
Error: PLINK does not support more than 2^31 - 3 variants.
written 4.9 years ago by Adarsh Kuamr 6

Hi there, I was converting my vcf file into bfiles in plink, and I got an error 'Error: PLINK does not support more than 2^31 …
Mutect2 vcf.stats file missing
written 5.6 years ago by rugarem 0

Hi everyone, Just wanted to ask if anyone can help me locate the `.vcf.stats` file supposedly made by Mutect2 along with the .vcf. I've called …
How to use UCSC net file to get genome rearrangement info between two species
written 10.5 years ago by Rongbin 1

Hi all, I have a fundamental problem. Actually I just want the information of rearrangement between human and mouse genome. It seems UCSC provide this …
Using Plink to test 1000 human Genome population data distance, most recent common ancestor and tim…
written 10.6 years ago by lyle 0

Hi everyone, I'm trying to use vcf files from the 1000 Human Genomes Project to establish distance between populations, most recent common ancestor and a …
Using VCFTools (v0.1.13) with ExAC (release 0.3.1)
written 10.6 years ago by areyoujokingme 7

I'm trying to extract allele frequency information for all variants in the latest ExAC data release. I'm currently using vcftools version 0.1.13, and am coming …
Extract each annotation in 'info' column in snpEff output
written 11.1 years ago by jan 19

Hi, I have used snpEff to annotate my vcf file. I previously use Annovar to annotate my vcf file and the fields are nicely separated …
Tool: Variantannotation Package (Bioconductor)
written 14.2 years ago by Malachi Griffith 2.0K

Using the Bioconductor [VariantAnnotation][1] Package to [Annotate Genetic Variants][2] > The VariantAnnotation package has facilities for reading in all or portions of Variant Call Format …