Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/a2/66b4d5680ee0859b26e0c89a5cbeff/fi_workdir/aih-tih-sc-0ce15f-R1_A23WJ53LT4_1.gtf
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/0d/02b371e63fd675b536d8a867c023b8/.command.sh
Downloading: s3://natera-rnd-pltf-dev-nextflow-scratch-01/work/0d/02b371e63fd675b536d8a867c023b8/.command.run
==> STAGING COMPLETE (3 inputs)
Using standard /usr/local/lib/perl5/site_perl/auto/share/dist/AGAT/agat_config.yaml file
------------------------------------------------------------------------------
| Another GFF Analysis Toolkit (AGAT) - Version: v1.2.0 |
| https://github.com/NBISweden/AGAT |
| National Bioinformatics Infrastructure Sweden (NBIS) - www.nbis.se |
------------------------------------------------------------------------------
------ Start parsing ------
-------------------------- parse options and metadata --------------------------
=> Accessing the feature_levels YAML file
Using standard /usr/local/lib/perl5/site_perl/auto/share/dist/AGAT/feature_levels.yaml file
=> Attribute used to group features when no Parent/ID relationship exists (i.e common tag):
* locus_tag
* gene_id
=> merge_loci option deactivated
=> Machine information:
This script is being run by perl v5.32.1
Bioperl location being used: /usr/local/lib/perl5/site_perl/Bio/
Operating system being used: linux
=> Accessing Ontology
No ontology accessible from the gff file header!
We use the SOFA ontology distributed with AGAT:
/usr/local/lib/perl5/site_perl/auto/share/dist/AGAT/so.obo
Read ontology /usr/local/lib/perl5/site_perl/auto/share/dist/AGAT/so.obo:
4 root terms, and 2596 total terms, and 1516 leaf terms
Filtering ontology:
We found 1861 terms that are sequence_feature or is_a child of it.
--------------------------------- parsing file ---------------------------------
=> Number of line in file: 1070
=> Number of comment lines: 0
=> Fasta included: No
=> Number of features lines: 1070
=> Number of feature type (3rd column): 2
* Level1: 0 =>
* level2: 0 =>
* level3: 2 => exon CDS
* unknown: 0 =>
=>Check because only level3 features:
* Number of feature with Parent attribute:0
* Number of feature with a common attribute:1070
=> Some common attributes and some Parent attributes missing.
/!\ For features where both are missing A single Level2 features (e.g. mRNA) and a single level1 (e.g. gene) will be created by AGAT, and all level3 feautres (e,g, CDS,exon) will be attached to them. This is probably not what you want...
see B. 2.2 and 3. at https://agat.readthedocs.io/en/latest/agat_how_does_it_work.html
/!\ For features where the common attribute or the parent attribute is missing, it would be fine as long as you do not expect isoforms in your annotation (Eukaryote). see B. 4. at https://agat.readthedocs.io/en/latest/agat_how_does_it_work.html
!! You might try to fix the issue by choosing a common tag attribute to use in order to group the features correctly (parameter --ct in agat_convert_sp_gxf2gxf.pl).
=> Version of the Bioperl GFF parser selected by AGAT: 2
------ End parsing (done in 0 second) ------
------ Start checks ------
---------------------------- Check1: feature types -----------------------------
----------------------------------- ontology -----------------------------------
All feature types in agreement with the Ontology.
------------------------------------- agat -------------------------------------
AGAT can deal with all the encountered feature types (3rd column)
------------------------------ done in 0 seconds -------------------------------
------------------------------ Check2: duplicates ------------------------------
None found
------------------------------ done in 0 seconds -------------------------------
-------------------------- Check3: sequential bucket ---------------------------
Nothing to check as sequential bucket!
------------------------------ done in 0 seconds -------------------------------
--------------------------- Check4: l2 linked to l3 ----------------------------
L1 and L2 created:
ABCA12--AC092844.1 HAVANA gene 1001 58400 . + . FI_gene_label "ABCA12^ENSG00000144452.15" ; ID "ABCA12--AC092844.1^ABCA12^ENSG00000144452.15" ; ccdsid "CCDS33372.1" ; exon_id "ENSE00001332842.5" ; exon_number 1 ; gene_id "ABCA12--AC092844.1^ABCA12^ENSG00000144452.15" ; gene_name ABCA12 ; gene_type protein_coding ; havana_gene "OTTHUMG00000154801.4" ; havana_transcript "OTTHUMT00000337111.2" ; hgnc_id "HGNC:14637" ; level 2 ; orig_coord_info "chr2,215138140,215138208,-" ; protein_id "ENSP00000272895.7" ; tag basic MANE_Select appris_principal_1 CCDS ; transcript_id "ABCA12--AC092844.1^ENST00000272895.12" ; transcript_name "ABCA12-201" ; transcript_support_level 1 ; transcript_type protein_coding
ABCA12--AC092844.1 HAVANA mRNA 1001 58400 . + . FI_gene_label "ABCA12^ENSG00000144452.15" ; ID "ABCA12--AC092844.1^ENST00000272895.12" ; Parent "ABCA12--AC092844.1^ABCA12^ENSG00000144452.15" ; ccdsid "CCDS33372.1" ; exon_id "ENSE00001332842.5" ; exon_number 1 ; gene_id "ABCA12--AC092844.1^ABCA12^ENSG00000144452.15" ; gene_name ABCA12 ; gene_type protein_coding ; havana_gene "OTTHUMG00000154801.4" ; havana_transcript "OTTHUMT00000337111.2" ; hgnc_id "HGNC:14637" ; level 2 ; orig_coord_info "chr2,215138140,215138208,-" ; protein_id "ENSP00000272895.7" ; tag basic MANE_Select appris_principal_1 CCDS ; transcript_id "ABCA12--AC092844.1^ENST00000272895.12" ; transcript_name "ABCA12-201" ; transcript_support_level 1 ; transcript_type protein_coding
L1 and L2 created:
POLR2J3--RASA4B HAVANA gene 13661 31288 . + . FI_gene_label "RASA4B^ENSG00000170667.15" ; ID "POLR2J3--RASA4B^RASA4B^ENSG00000170667.15" ; exon_id "ENSE00003705682.1" ; exon_number 3 ; gene_id "POLR2J3--RASA4B^RASA4B^ENSG00000170667.15" ; gene_name RASA4B ; gene_type protein_coding ; havana_gene "OTTHUMG00000022944.11" ; havana_transcript "OTTHUMT00000351773.1" ; hgnc_id "HGNC:35202" ; level 2 ; orig_coord_info "chr7,102510179,102510198,-" ; protein_id "ENSP00000303968.6" ; tag basic ; transcript_id "POLR2J3--RASA4B^ENST00000306682.6" ; transcript_name "RASA4B-201" ; transcript_support_level 1 ; transcript_type protein_coding
POLR2J3--RASA4B HAVANA mRNA 13661 31288 . + . FI_gene_label "RASA4B^ENSG00000170667.15" ; ID "POLR2J3--RASA4B^ENST00000306682.6" ; Parent "POLR2J3--RASA4B^RASA4B^ENSG00000170667.15" ; exon_id "ENSE00003705682.1" ; exon_number 3 ; gene_id "POLR2J3--RASA4B^RASA4B^ENSG00000170667.15" ; gene_name RASA4B ; gene_type protein_coding ; havana_gene "OTTHUMG00000022944.11" ; havana_transcript "OTTHUMT00000351773.1" ; hgnc_id "HGNC:35202" ; level 2 ; orig_coord_info "chr7,102510179,102510198,-" ; protein_id "ENSP00000303968.6" ; tag basic ; transcript_id "POLR2J3--RASA4B^ENST00000306682.6" ; transcript_name "RASA4B-201" ; transcript_support_level 1 ; transcript_type protein_coding
L1 and L2 created:
POLR2J3--RASA4B HAVANA gene 1001 10115 . + . FI_gene_label "POLR2J3^ENSG00000285437.1" ; ID "POLR2J3--RASA4B^POLR2J3^ENSG00000285437.1" ; exon_id "ENSE00003824619.1" ; exon_number 1 ; gene_id "POLR2J3--RASA4B^POLR2J3^ENSG00000285437.1" ; gene_name POLR2J3 ; gene_type protein_coding ; havana_gene "OTTHUMG00000192387.1" ; havana_transcript "OTTHUMT00000358591.2" ; hgnc_id "HGNC:33853" ; level 2 ; orig_coord_info "chr7,102572472,102572521,-" ; protein_id "ENSP00000422109.1" ; tag basic ; transcript_id "POLR2J3--RASA4B^ENST00000511313.5" ; transcript_name "POLR2J3-217" ; transcript_support_level 1 ; transcript_type protein_coding
POLR2J3--RASA4B HAVANA mRNA 1001 10115 . + . FI_gene_label "POLR2J3^ENSG00000285437.1" ; ID "POLR2J3--RASA4B^ENST00000511313.5" ; Parent "POLR2J3--RASA4B^POLR2J3^ENSG00000285437.1" ; exon_id "ENSE00003824619.1" ; exon_number 1 ; gene_id "POLR2J3--RASA4B^POLR2J3^ENSG00000285437.1" ; gene_name POLR2J3 ; gene_type protein_coding ; havana_gene "OTTHUMG00000192387.1" ; havana_transcript "OTTHUMT00000358591.2" ; hgnc_id "HGNC:33853" ; level 2 ; orig_coord_info "chr7,102572472,102572521,-" ; protein_id "ENSP00000422109.1" ; tag basic ; transcript_id "POLR2J3--RASA4B^ENST00000511313.5" ; transcript_name "POLR2J3-217" ; transcript_support_level 1 ; transcript_type protein_coding
L1 and L2 created:
STX16--NPEPL1 HAVANA gene 1021 7754 . + . FI_gene_label "STX16^ENSG00000124222.22" ; ID "STX16--NPEPL1^STX16^ENSG00000124222.22" ; exon_id "ENSE00003459479.1" ; exon_number 3 ; gene_id "STX16--NPEPL1^STX16^ENSG00000124222.22" ; gene_name STX16 ; gene_type protein_coding ; havana_gene "OTTHUMG00000033084.8" ; havana_transcript "OTTHUMT00000267911.2" ; hgnc_id "HGNC:11431" ; level 2 ; orig_coord_info "chr20,58667505,58667597,+" ; protein_id "ENSP00000312086.8" ; tag alternative_5_UTR mRNA_end_NF cds_end_NF ; transcript_id "STX16--NPEPL1^ENST00000312283.12" ; transcript_name "STX16-201" ; transcript_support_level 3 ; transcript_type protein_coding
STX16--NPEPL1 HAVANA mRNA 1021 7754 . + . FI_gene_label "STX16^ENSG00000124222.22" ; ID "STX16--NPEPL1^ENST00000312283.12" ; Parent "STX16--NPEPL1^STX16^ENSG00000124222.22" ; exon_id "ENSE00003459479.1" ; exon_number 3 ; gene_id "STX16--NPEPL1^STX16^ENSG00000124222.22" ; gene_name STX16 ; gene_type protein_coding ; havana_gene "OTTHUMG00000033084.8" ; havana_transcript "OTTHUMT00000267911.2" ; hgnc_id "HGNC:11431" ; level 2 ; orig_coord_info "chr20,58667505,58667597,+" ; protein_id "ENSP00000312086.8" ; tag alternative_5_UTR mRNA_end_NF cds_end_NF ; transcript_id "STX16--NPEPL1^ENST00000312283.12" ; transcript_name "STX16-201" ; transcript_support_level 3 ; transcript_type protein_coding
L1 and L2 created:
STX16--NPEPL1 HAVANA gene 18353 29399 . + . FI_gene_label "NPEPL1^ENSG00000215440.12" ; ID "STX16--NPEPL1^NPEPL1^ENSG00000215440.12" ; ccdsid "CCDS46621.1" ; exon_id "ENSE00001408052.7" ; exon_number 1 ; gene_id "STX16--NPEPL1^NPEPL1^ENSG00000215440.12" ; gene_name NPEPL1 ; gene_type protein_coding ; havana_gene "OTTHUMG00000033060.8" ; havana_transcript "OTTHUMT00000080402.7" ; hgnc_id "HGNC:16244" ; level 2 ; orig_coord_info "chr20,58692901,58693050,+" ; protein_id "ENSP00000348395.6" ; tag basic MANE_Select appris_principal_1 CCDS ; transcript_id "STX16--NPEPL1^ENST00000356091.11" ; transcript_name "NPEPL1-201" ; transcript_support_level 1 ; transcript_type protein_coding
STX16--NPEPL1 HAVANA mRNA 18353 29399 . + . FI_gene_label "NPEPL1^ENSG00000215440.12" ; ID "STX16--NPEPL1^ENST00000356091.11" ; Parent "STX16--NPEPL1^NPEPL1^ENSG00000215440.12" ; ccdsid "CCDS46621.1" ; exon_id "ENSE00001408052.7" ; exon_number 1 ; gene_id "STX16--NPEPL1^NPEPL1^ENSG00000215440.12" ; gene_name NPEPL1 ; gene_type protein_coding ; havana_gene "OTTHUMG00000033060.8" ; havana_transcript "OTTHUMT00000080402.7" ; hgnc_id "HGNC:16244" ; level 2 ; orig_coord_info "chr20,58692901,58693050,+" ; protein_id "ENSP00000348395.6" ; tag basic MANE_Select appris_principal_1 CCDS ; transcript_id "STX16--NPEPL1^ENST00000356091.11" ; transcript_name "NPEPL1-201" ; transcript_support_level 1 ; transcript_type protein_coding
L1 and L2 created:
TLK2--LINC01347 HAVANA gene 4784 35054 . + . FI_gene_label "TLK2^ENSG00000146872.18" ; ID "TLK2--LINC01347^TLK2^ENSG00000146872.18" ; ccdsid "CCDS62283.1" ; exon_id "ENSE00003568961.1" ; exon_number 2 ; gene_id "TLK2--LINC01347^TLK2^ENSG00000146872.18" ; gene_name TLK2 ; gene_type protein_coding ; havana_gene "OTTHUMG00000179176.3" ; havana_transcript "OTTHUMT00000445140.1" ; hgnc_id "HGNC:11842" ; level 2 ; orig_coord_info "chr17,62481126,62481206,+" ; protein_id "ENSP00000316512.9" ; tag basic CCDS ; transcript_id "TLK2--LINC01347^ENST00000326270.13" ; transcript_name "TLK2-201" ; transcript_support_level 1 ; transcript_type protein_coding
TLK2--LINC01347 HAVANA mRNA 4784 35054 . + . FI_gene_label "TLK2^ENSG00000146872.18" ; ID "TLK2--LINC01347^ENST00000326270.13" ; Parent "TLK2--LINC01347^TLK2^ENSG00000146872.18" ; ccdsid "CCDS62283.1" ; exon_id "ENSE00003568961.1" ; exon_number 2 ; gene_id "TLK2--LINC01347^TLK2^ENSG00000146872.18" ; gene_name TLK2 ; gene_type protein_coding ; havana_gene "OTTHUMG00000179176.3" ; havana_transcript "OTTHUMT00000445140.1" ; hgnc_id "HGNC:11842" ; level 2 ; orig_coord_info "chr17,62481126,62481206,+" ; protein_id "ENSP00000316512.9" ; tag basic CCDS ; transcript_id "TLK2--LINC01347^ENST00000326270.13" ; transcript_name "TLK2-201" ; transcript_support_level 1 ; transcript_type protein_coding
L1 and L2 created:
ABCA12--AC092844.1 HAVANA gene 61401 66632 . + . FI_gene_label "AC092844.1^ENSG00000286836.1" ; ID "ABCA12--AC092844.1^AC092844.1^ENSG00000286836.1" ; exon_id "ENSE00003853015.1" ; exon_number 1 ; gene_id "ABCA12--AC092844.1^AC092844.1^ENSG00000286836.1" ; gene_name "AC092844.1" ; gene_type lncRNA ; havana_gene "OTTHUMG00000195823.1" ; havana_transcript "OTTHUMT00000528925.1" ; level 2 ; orig_coord_info "chr2,215310679,215310947,-" ; tag basic TAGENE ; transcript_id "ABCA12--AC092844.1^ENST00000657931.1" ; transcript_name "AC092844.1-201" ; transcript_type lncRNA
ABCA12--AC092844.1 HAVANA RNA 61401 66632 . + . FI_gene_label "AC092844.1^ENSG00000286836.1" ; ID "ABCA12--AC092844.1^ENST00000657931.1" ; Parent "ABCA12--AC092844.1^AC092844.1^ENSG00000286836.1" ; exon_id "ENSE00003853015.1" ; exon_number 1 ; gene_id "ABCA12--AC092844.1^AC092844.1^ENSG00000286836.1" ; gene_name "AC092844.1" ; gene_type lncRNA ; havana_gene "OTTHUMG00000195823.1" ; havana_transcript "OTTHUMT00000528925.1" ; level 2 ; orig_coord_info "chr2,215310679,215310947,-" ; tag basic TAGENE ; transcript_id "ABCA12--AC092844.1^ENST00000657931.1" ; transcript_name "AC092844.1-201" ; transcript_type lncRNA
L1 and L2 created:
TLK2--LINC01347 HAVANA gene 41405 47393 . + . FI_gene_label "LINC01347^ENSG00000214837.8" ; ID "TLK2--LINC01347^LINC01347^ENSG00000214837.8" ; exon_id "ENSE00002389344.1" ; exon_number 1 ; gene_id "TLK2--LINC01347^LINC01347^ENSG00000214837.8" ; gene_name LINC01347 ; gene_type transcribed_unprocessed_pseudogene ; havana_gene "OTTHUMG00000039861.5" ; havana_transcript "OTTHUMT00000481524.1" ; hgnc_id "HGNC:50566" ; level 2 ; ont "PGO:0000005" "PGO:0000019" ; orig_coord_info "chr1,243097190,243097299,-" ; tag basic ; transcript_id "TLK2--LINC01347^ENST00000417964.2" ; transcript_name "LINC01347-201" ; transcript_support_level NA ; transcript_type transcribed_unprocessed_pseudogene
TLK2--LINC01347 HAVANA RNA 41405 47393 . + . FI_gene_label "LINC01347^ENSG00000214837.8" ; ID "TLK2--LINC01347^ENST00000417964.2" ; Parent "TLK2--LINC01347^LINC01347^ENSG00000214837.8" ; exon_id "ENSE00002389344.1" ; exon_number 1 ; gene_id "TLK2--LINC01347^LINC01347^ENSG00000214837.8" ; gene_name LINC01347 ; gene_type transcribed_unprocessed_pseudogene ; havana_gene "OTTHUMG00000039861.5" ; havana_transcript "OTTHUMT00000481524.1" ; hgnc_id "HGNC:50566" ; level 2 ; ont "PGO:0000005" "PGO:0000019" ; orig_coord_info "chr1,243097190,243097299,-" ; tag basic ; transcript_id "TLK2--LINC01347^ENST00000417964.2" ; transcript_name "LINC01347-201" ; transcript_support_level NA ; transcript_type transcribed_unprocessed_pseudogene
L1 and L2 created:
TLK2P1--AC110079.1 HAVANA gene 1001 4192 . + . FI_gene_label "TLK2P1^ENSG00000226049.3" ; ID "TLK2P1--AC110079.1^TLK2P1^ENSG00000226049.3" ; exon_id "ENSE00002174619.1" ; exon_number 1 ; gene_id "TLK2P1--AC110079.1^TLK2P1^ENSG00000226049.3" ; gene_name TLK2P1 ; gene_type processed_pseudogene ; havana_gene "OTTHUMG00000166422.1" ; havana_transcript "OTTHUMT00000389700.1" ; hgnc_id "HGNC:18048" ; level 1 ; ont "PGO:0000004" ; orig_coord_info "chr17,34036681,34039872,-" ; tag pseudo_consens basic ; transcript_id "TLK2P1--AC110079.1^ENST00000530992.1" ; transcript_name "TLK2P1-201" ; transcript_support_level NA ; transcript_type processed_pseudogene
TLK2P1--AC110079.1 HAVANA RNA 1001 4192 . + . FI_gene_label "TLK2P1^ENSG00000226049.3" ; ID "TLK2P1--AC110079.1^ENST00000530992.1" ; Parent "TLK2P1--AC110079.1^TLK2P1^ENSG00000226049.3" ; exon_id "ENSE00002174619.1" ; exon_number 1 ; gene_id "TLK2P1--AC110079.1^TLK2P1^ENSG00000226049.3" ; gene_name TLK2P1 ; gene_type processed_pseudogene ; havana_gene "OTTHUMG00000166422.1" ; havana_transcript "OTTHUMT00000389700.1" ; hgnc_id "HGNC:18048" ; level 1 ; ont "PGO:0000004" ; orig_coord_info "chr17,34036681,34039872,-" ; tag pseudo_consens basic ; transcript_id "TLK2P1--AC110079.1^ENST00000530992.1" ; transcript_name "TLK2P1-201" ; transcript_support_level NA ; transcript_type processed_pseudogene
L1 and L2 created:
TLK2P1--AC110079.1 HAVANA gene 7193 19707 . + . FI_gene_label "AC110079.1^ENSG00000260404.3" ; ID "TLK2P1--AC110079.1^AC110079.1^ENSG00000260404.3" ; exon_id "ENSE00002625896.1" ; exon_number 1 ; gene_id "TLK2P1--AC110079.1^AC110079.1^ENSG00000260404.3" ; gene_name "AC110079.1" ; gene_type transcribed_unprocessed_pseudogene ; havana_gene "OTTHUMG00000161164.4" ; havana_transcript "OTTHUMT00000364170.2" ; level 2 ; orig_coord_info "chr4,118591773,118591895,+" ; tag dotter_confirmed basic ; transcript_id "TLK2P1--AC110079.1^ENST00000567913.2" ; transcript_name "AC110079.1-201" ; transcript_support_level 5 ; transcript_type processed_transcript
TLK2P1--AC110079.1 HAVANA RNA 7193 19707 . + . FI_gene_label "AC110079.1^ENSG00000260404.3" ; ID "TLK2P1--AC110079.1^ENST00000567913.2" ; Parent "TLK2P1--AC110079.1^AC110079.1^ENSG00000260404.3" ; exon_id "ENSE00002625896.1" ; exon_number 1 ; gene_id "TLK2P1--AC110079.1^AC110079.1^ENSG00000260404.3" ; gene_name "AC110079.1" ; gene_type transcribed_unprocessed_pseudogene ; havana_gene "OTTHUMG00000161164.4" ; havana_transcript "OTTHUMT00000364170.2" ; level 2 ; orig_coord_info "chr4,118591773,118591895,+" ; tag dotter_confirmed basic ; transcript_id "TLK2P1--AC110079.1^ENST00000567913.2" ; transcript_name "AC110079.1-201" ; transcript_support_level 5 ; transcript_type processed_transcript
72 cases fixed where L3 features have parent feature(s) missing
------------------------------ done in 0 seconds -------------------------------
--------------------------- Check5: l1 linked to l2 ----------------------------
No problem found
------------------------------ done in 0 seconds -------------------------------
--------------------------- Check6: remove orphan l1 ---------------------------
We remove only those not supposed to be orphan
None found
------------------------------ done in 0 seconds -------------------------------
------------------------- Check7: all level3 locations -------------------------
------------------------------ done in 0 seconds -------------------------------
------------------------------ Check8: check cds -------------------------------
No problem found
------------------------------ done in 0 seconds -------------------------------
----------------------------- Check9: check exons ------------------------------
No exons created
No exons locations modified
No supernumerary exons removed
No level2 locations modified
------------------------------ done in 0 seconds -------------------------------
----------------------------- Check10: check utrs ------------------------------
134 UTRs created that were missing
No UTRs locations modified
No supernumerary UTRs removed
------------------------------ done in 0 seconds -------------------------------
------------------------ Check11: all level2 locations -------------------------
No problem found
------------------------------ done in 0 seconds -------------------------------
------------------------ Check12: all level1 locations -------------------------
We fixed 6 wrong level1 location cases
------------------------------ done in 0 seconds -------------------------------
---------------------- Check13: remove identical isoforms ----------------------
None found
------------------------------ done in 1 seconds -------------------------------
------ End checks (done in 1 second) ------
GFF3 file parsed