美文网首页ggplot集锦
Gff3_to_gtf && gff3_to_bed

Gff3_to_gtf && gff3_to_bed

作者: 余绕 | 来源:发表于2021-10-27 11:54 被阅读0次

    1.gff3_to_bed

    Approadch 1:

    利用TransDecoder的gff3_file_to_bed.pl脚本完成(conda 安装TransDecoder)

    gff3_file_to_bed.pl /public/home/qtxu/rice_anotation/Rice_annotation/Tige7.0.gff3 >Tigr_MSU_7_raw.bed
    
    #format the bed files!
    
    sed 's/;LOC\S*\t/\t/g' Tigr_MSU_7_raw.bed |sed 's/ID=//g' >Tigr_MSU_7.bed
    
    
    Approadch 2:

    下载 Gff3ToBed(https://github.com/YoannPa/Gff3ToBed

    sudo chmod +x Gff3ToBed
    
    # Example: Extract all data to a Bed file
    ./Gff3ToBed.sh gencode.v26.annotation.gff3 gff3 all both
    # 
    
    
    Approadch 3:
    #下载# [BEDOPS v2.4.40](https://bedops.readthedocs.io/en/latest/content/reference/file-management/conversion/gff2bed.html#downloads)
    
    #We can convert it to sorted BED data in the following manner:
    
    $ gff2bed < foo.gff3
    chr1    1049    1500    exon00002       .       -       USA     exon    0       ID=exon00002;Ontology_term="GO:0046703";Ontology_term="GO:0046704"
    chr1    1299    1300    exon00001       .       +       Canada  exon    .       ID=exon00001;score=1;zeroLengthInsertion=True
    chr1    2999    3902    exon00003       .       ?       Canada  exon    2       ID=exon00003;score=4;Name=foo
    chr1    4999    5500    exon00004       .       .       .       exon    .       ID=exon00004;Gap=M8 D3 M6 I1 M6
    chr1    6999    9000    exon00005       10      +       .       exon    1       ID=exon00005;Dbxref="NCBI_gi:10727410"
    
    

    2.gff3_to_GTF

    gff 2 gtf

    
    gffread  gencode.v19.annotation.gff3  -T  -o   gencode.v19.gtf
    
    

    相关文章

      网友评论

        本文标题:Gff3_to_gtf && gff3_to_bed

        本文链接:https://www.haomeiwen.com/subject/fkoualtx.html