Adding the first line in a gff3 file
0
0
Entering edit mode
4.6 years ago
relindisg • 0

I am working on WGBS for strawberry (Fragaria vesca) genome. I want to annotate CpG,CHG and CHH islands in the genome. The recent version of the genome has a transcript gff3 file and a genes gff3 file. I am using the methylkit to annotate differentially methylated bases /regions and methylkit requires a bed file for annotation.I have converted the transcript gff3 file to bed.However, I am not able to convert the genes.gff3 file to a bed because the file lacks the first line.When i cheched the file using gff3 validator, it shows that the file lacks the first line denoted as " ##gff-version 3". My questions are:

  1. How can I insert the line" ##gff-version 3" in the gff3 file.
  2. Can anyone help calrify if any of the above mentioned gff3 files contain information about CpG, CHH and CHG features? if not.
  3. Where can I get a gff3 file or bed file that can be used for annotating CpG,CHG and CHH.

I appreciate your help.

sequencing • 1.1k views
ADD COMMENT
1
Entering edit mode
echo  "##gff-version 3" > new.gff && cat old.gff >> new.gff
ADD REPLY

Login before adding your answer.

Traffic: 2950 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6