This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Incompatibility issue with fasta and GFF (fasta2GFF)

Good day to you all, I am about to run a script that requires a fasta and a gff file in the command. But the fatsa and gff I have downloaded from NCBI contain different loci names. for example fasta has the following type loci naming

Bf_V2_22 
Bf_V2_22

Gff has the following type

NW_003101570.1 
NW_003101570.1

I need both of them to have same names. What do I do? would I get the same cordinates in the gff if I convert fasta2gff? ( How do I convert fasta2gff). Or what is a better thing to do ? Thanks

fasta2gff mapping bedtools
• 2,399 views • link
written by jaqx008 jaqx008

Hi finswimmer. I downloaded them from ncbi as mentioned above.

• 1 views • link
written by jaqx008 jaqx008

NCBI has tones of data. So please be more specific.

• 1 views • link
written by finswimmer finswimmer

1 answer

You probably did not get the sequence from the same location.

$ zgrep ">" GCF_000003815.1_Version_2_genomic.fna.gz
>NW_003101570.1 Branchiostoma floridae genomic scaffold BRAFLscaffold_1, whole genome shotgun sequence
>NW_003101569.1 Branchiostoma floridae genomic scaffold BRAFLscaffold_2, whole genome shotgun sequence
>NW_003101568.1 Branchiostoma floridae genomic scaffold BRAFLscaffold_3, whole genome shotgun sequence
>NW_003101567.1 Branchiostoma floridae genomic scaffold BRAFLscaffold_4, whole genome shotgun sequence
>NW_003101566.1 Branchiostoma floridae genomic scaffold BRAFLscaffold_5, whole genome shotgun sequence
>NW_003101565.1 Branchiostoma floridae genomic scaffold BRAFLscaffold_6, whole genome shotgun sequence
>NW_003101564.1 Branchiostoma floridae genomic scaffold BRAFLscaffold_7, whole genome shotgun sequence
>NW_003101563.1 Branchiostoma floridae genomic scaffold BRAFLscaffold_8, whole genome shotgun sequence
>NW_003101562.1 Branchiostoma floridae genomic scaffold BRAFLscaffold_9, whole genome shotgun sequence
• 1 views • link
written by GenoMax GenoMax

Were did you get this from please?

• 1 views • link
written by jaqx008 jaqx008

I provided a link to the file above.

• 1 views • link
written by GenoMax GenoMax

I just found that the same script have to call a x.bed file. my bed file has the loci in the other format which is

Bf_V2_22 
Bf_V2_22

this is probably because the TE.bed was made from the fasta file with Bf_v type naming

• 1 views • link
written by jaqx008 jaqx008

Then you will either need to remake the bed file or translate those ID's so they match NCBI's identifiers.

• 1 views • link
written by GenoMax GenoMax

Do I have to do that manually? and Isnt there a possibility for mismatch of ID?

• 1 views • link
written by jaqx008 jaqx008

You will need to go back and re-make/process the data to get ID's that match.

• 1 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this