More posts like this
-
How do blastp step in MCScanX
written by carlosgonzalezcruz327 •
-
Ho can I add a standard desviation to my diversity values?
written by apl00028
-
I can not assembly this fasta file
written by apl00028
-
How can I see the low frequency SNP
written by apl00028
-
How can I get Allele Depth of my vcf file ?
written by apl00028
-
Blast wrong, RNA seq
written by apl00028
-
Incompatibility issue with fasta and GFF (fasta2GFF)
written by jaqx008
-
What program do I must to use to design primer for especific samples?
written by apl00028
-
Pipeline for genome annotation?
written by Buffo
-
.map file for Impute2
written by lmh327 •
As posed this question does not make a lot of sense. VCF files describe polymorphisms compared to a reference, where as GFF files are used for describing gene models/annotations. They are not directly inter-convertible.
What exactly are you trying to do?
You are right, I am sorry. I trying to analyze my polimorfism data using this commands in R:
For do that I need the vcf file, the reference sequence and the gff file of my vcf file, but I do not how to get that.
I got that code from here: https://cran.r-project.org/web/packages/vcfR/vignettes/intro_to_vcfR.html
You actually need a GFF file for the reference genome you used to create that VCF file. Generally recommended that you get GFF from the same source where you got the reference from, so there are no discrepancies with chromosomal identifiers etc.
which program should do I use for that?
What genome are you working with and where did you get the reference sequence from?
My reference sequence I got for an assembly of my data
You probably have some additional work to do then. Get the RefSeq assembly for this virus from NCBI. You will also find a GFF file in that same directory.
Compare your assembly with the reference and depending on changes you see you could modify the GFF file available for RefSeq genome and use that.
You need to tell us more about your data first. We need far more information.
What organism? Where did you get the VCF? How was the VCF derived? From what sequence?
The organism is Cucumber Mosaic Virus RNA3, I got this vcf file from bam file using these functions:
So this is a genome you de novo assembled?
Yes, using IDBA_UD program
Hi, are you the same person as the user apl0028?
No, I work with him. He can not reply to you because he spent all messages per day
I see. We will figure out if it is possible to remove this comment limit restriction.
By the way, if you've been up-voting his post and comments, please understand that such upvoting is not encouraged. Upvotes are for content with inherent objective value only, and almost none of apl0028's post/comments seem to have that value. In case such upvotes were given out, please go back and undo them.
Please, could you answer to me? How can I get this file?
genomax gave you an answer, see his comment up above.