Question: Adding the first line in a gff3 file
gravatar for relindisg
8 months ago by
relindisg0 wrote:

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 • 276 views
ADD COMMENTlink written 8 months ago by relindisg0
echo  "##gff-version 3" > new.gff && cat old.gff >> new.gff
ADD REPLYlink written 8 months ago by Pierre Lindenbaum128k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1867 users visited in the last hour