Question: Negative value in "phase" line of a gff3 file.What does it mean?
gravatar for aramis.1994
2.5 years ago by
aramis.199410 wrote:

Hi you all,

I am trying to perform featureCount over miRNA reads from CHO cells for differential expression analysis purposes. I have created my own annotation file by mapping hairpin miRNA sequences downloaded from mirBase ( against ENSEMBL reference genome for CHO cell. For that purpose, I have used gmap -D:

gmap -D ~/miRNA/crigri_gmap -d crigri_gmap -f 2 -n 0 -t 16 --gff3-cds=genomic hairpin_crigri_dna_mod.fa > trial_1.gff3

I retrieve a .gff file with 9 columns containing the 9 expected fields (seqid,source,type,start,end,score,strand,phase and attributes), but for some scaffolds, I get a "-1" value in the "phase" field. According to ENSEMBLE ( 'One of '0', '1' or '2'. '0' indicates that the first base of the feature is the first base of a codon, '1' that the second base is the first base of a codon, and so on.' So I don't know how to interpret this value.

This is how the head of my .gff3 file looks:

scaffold_6      crigri_gmap     gene    62246878        62246977        .       +       .       ID=cgr_let_7a_MI0020368.path1;Name=cgr_let_7a_MI0020368

scaffold_6      crigri_gmap     mRNA    62246878        62246977        .       +       .       ID=cgr_let_7a_MI0020368.mrna1;Name=cgr_let_7a_MI0020368;Parent=cgr_let_7a_MI0020368.path1;coverage=100.0;identity=100.0;matches=100;mismatches=0;indels=0;unknowns=0

scaffold_6      crigri_gmap     exon    1       62246878        100     +       .       ID=cgr_let_7a_MI0020368.mrna1.exon1;Name=cgr_let_7a_MI0020368;Parent=cgr_let_7a_MI0020368.mrna1;Target=cgr_let_7a_MI0020368 1 1 +

scaffold_6      crigri_gmap     CDS     62246878        62246975        100     +       -1      ID=cgr_let_7a_MI0020368.mrna1.cds1;Name=cgr_let_7a_MI0020368;Parent=cgr_let_7a_MI0020368.mrna1;Target=cgr_let_7a_MI0020368 1 98 +

Any help with this will be much appreciated. Thank you very much!!

.gff file rna-seq • 1.0k views
ADD COMMENTlink modified 6 months ago by Biostar ♦♦ 20 • written 2.5 years ago by aramis.199410

This is strange, no idea why there is negative value... But to fix the phases you can use from AGAT

ADD REPLYlink written 6 months ago by Juke344.4k
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: 1170 users visited in the last hour