I want to manipulate (acutally add information to) the vcf file and its header information. I tried using sed in linux but just couldn't understand it.
It original vcf file looks like:
CHROM POS REF ALT
1 69 A G
1 75 C T
1 85 G A
I want to add the "ID" field after the 2nd column and replace the value for this field with "." Additionally I want to add other fields "QUAL", "FILTER" and "INFO" after the last column. All the columns should be tab separated and finally should look like.
CHROM POS ID REF ALT QUAL FILTER INFO
1 69 . A C . . .
1 75 . G A . . .
1 87 . C A . . .
1 88 . T C . . .
Thanks much in advance !
•
written
by
kirannbishwa01
1 answer
You're looking to modify a stream of text. Substitution using sed works like this:
sed 's/old_text_or_regex/new_text_or_regex/modifiers/' input_file >output_file
But this scenario is better addressed using awk, where you can process individual columns and just add column at the appropriate position.
Once again, I'm providing no code because that would cripple learning.
Log in to answer this question.
-
Linux command to edit CHROM column in VCF file
written by anilkumarboddapati •
-
Reomve the first occurance of overalapping variants from VCF file
written by Frieda •
-
Converting GP field in a VCF to GT
written by Garan
-
Extract AD thatt ' support the ALTERNATIVE ALLELE ' from VCF file
written by Maryem •
-
Extract the DP from VCF file along with the chromosome postion and alteration
written by manojkumarbioinfo
-
How to remove chr from VCF file
written by MAPK
-
How to seperate vcf annoation info field
written by raya.girish •
-
Comparing samtools,gatk and varscan against gold standard 1000G snps and need to decide the best me…
written by bioinforesearchquestions
-
Calculate Allele balance by sample
written by kirannbishwa01
-
A tool to annotate one VCF file with INFO records of another VCF taking SNP into account?
written by tsukanoffkirill •
