Command line: [exonerate -m p2g --showtargetgff -q ./SPP00000009_1.0.fa -t ./SPP00000009_1.0/fastasubseq/KV390107.1_fastasubseq.fa] Hostname: [ws101072u] C4 Alignment: ------------ Query: SPP00000009_1.0 # Protein # Glutathione peroxidase 6 (GPx6) # Homo sapiens # Complete Target: KV390107.1:subseq(473451,100206) Cebus capucinus imitator isolate Cc_AM_T3 unplaced genomic scaffold Scaffold681, whole genome shotgun sequence Model: protein2genome:local Raw score: 747 Query range: 21 -> 220 Target range: 75356 -> 79932 22 : ThrLeuLysProGlnAsnArgLysValAspCysAsnLysGlyValThrGlyThrIleTyr : 41 :::|||||| ::: :::::!||||||:!!!:!! ! !! !|||||||||||| SerLeuLysLysLysIlePheGlnMetAspCysHisArgAspGlnLysGlyThrIleTyr 75357 : tcacttaaaaaaaaaatcttccagatGGATTGTCACAGAGATCAAAAAGGCACCATCTAT : 75414 42 : GluTyrGlyAlaLeuThrLeuAsnGlyGluGluTyrIleGlnPheLysGlnPheAlaGly : 61 !!:|||! !|||:!!! !|||||| !:!!||||||:!!! ||||||!:!!:!!.!||| AspTyrGluAlaIleMetLeuAsnLysLysGluTyrValProPheLysArgTyrValGly 75415 : GACTATGAGGCCATCATGCTTAATAAGAAGGAATACGTTCCCTTCAAGCGGTATGTGGGC : 75474 62 : LysHisValLeuPheValAsnValAlaAlaTyrUnkGlyLeuAlaAlaGlnTyrPro{G} : 81 |||||||||||||||||||||||||||.!!||| ||||||.!!||||||||||||{|} LysHisValLeuPheValAsnValAlaThrTyrCysGlyLeuThrAlaGlnTyrPro{G} 75475 : AAGCACGTCCTCTTTGTCAACGTGGCTACCTACTGTGGTCTGACAGCACAATATCCT{G} : 75534 82 : >>>> Target Intron 1 >>>> {lu}LeuAsnAlaLeuGlnGluGluLeuLys : 90 2023 bp {||}||||||||||||||||||||||||||| ++ ++{lu}LeuAsnAlaLeuGlnGluGluLeuLys 75535 : gt.........................ag{AA}CTGAATGCACTCCAGGAGGAGCTGAAG : 77584 91 : AsnPheGlyValIleValLeuAlaPheProCysAsnGlnPheGlyLysGlnGluProGly : 110 !||||||:!!:!!||||||!.!|||||||||||||||||||||||||||||||||||| ProPheGlyLeuValValLeuGlyPheProCysAsnGlnPheGlyLysGlnGluProGly 77585 : CCCTTTGGTCTAGTTGTGCTGGGCTTTCCCTGCAACCAGTTTGGAAAGCAAGAACCAGGA : 77644 111 : ThrAsnSerGluIleLeuLeuGlyLeu{Ly} >>>> Target Intron 2 >>>> : 120 !|||!.!||||||||| ||||||{||} 418 bp AspAsnLysGluIleLeuProGlyLeu{Ly}++ ++ 77645 : GATAACAAAGAGattcttcctgggctc{aa}gt.........................ag : 78092 121 : {s}TyrValCysProGlySerGlyPheValProSerPheGlnLeuPheGluLysGlyAsp : 139 {|}|||||| !!||||||.!!||||||||||||!:!|||||||||||||||||||||||| {s}TyrValArgProGlyGlyGlyPheValProAsnPheGlnLeuPheGluLysGlyAsp 78093 : {G}TATGTCCGTCCAGGGGGAGGATTTGTACCTAATTTCCAGCTTTTTGAGAAAGGAGAT : 78149 140 : ValAsnGlyGluLysGluGlnLysValPheThrPheLeuLys >>>> Target Intr : 154 ||||||||||||||||||||||||||||||!:!||||||||| 1538 bp ValAsnGlyGluLysGluGlnLysValPheSerPheLeuLys++ 78150 : GTGAATggtgaaaaagaacagaaagtctTCAGCTTCTTGAAGgt................ : 78196 155 : on 3 >>>> AsnSerCysProProThrSerAspLeuLeuGlySerSerSerGlnLeuP : 170 :!!|||||||||! ! !!|||!!::!!||||||!!!! !!.!...:!!! ++HisSerCysProHisProSerGluIleLeuGlyThrPheLysSerIleS 78197 : .........agCACTCCTGTCCTCATCCCTCTGAGATTTTGGGCACATTCAAATCTATAT : 79780 171 : heTrpGluProMetLysValHisAspIleArgTrpAsnPheGluLysPheLeuValGlyP : 190 !|||.!.|||:!:|||||||||||||||||||||||||||||||||||||||||||||| erTrpAsnProValLysValHisAspIleArgTrpAsnPheGluLysPheLeuValGlyP 79781 : CCTGGAACCCTGTCAAAGTCCATGACATCCGTTGGAACTTTGAAAAGTTCCTGGTGGGGC : 79840 191 : roAspGlyValProValMetHisTrpPheHisGlnAlaProValSerThrValLysSerA : 210 ||:!!|||:!!|||||||||!..|||! !||||||||| !!|||||| !||||||:!!| roAsnGlyIleProValMetArgTrpSerHisGlnAlaThrValSerLeuValLysAlaA 79841 : CCAATGGAATCCCTGTCATGCGGTGGTCCCACCAGGCCACAGTCAGCTTAGTCAAGGCAG : 79900 211 : spIleLeuGluTyrLeuLysGlnPheAsnThr : 220 |||||:!!! |||||||||||||||!!.||| spIleValAlaTyrLeuLysGlnPheLysThr 79901 : ACATCGTGGCCTACTTGAAGCAATTCAAAACC : 79932 vulgar: SPP00000009_1.0 21 220 . KV390107.1:subseq(473451,100206) 75356 79932 + 747 M 59 177 S 0 1 5 0 2 I 0 2019 3 0 2 S 1 2 M 38 114 S 0 2 5 0 2 I 0 414 3 0 2 S 1 1 M 33 99 5 0 2 I 0 1534 3 0 2 M 67 201 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-16 ##type DNA # # # seqname source feature start end score strand frame attributes # KV390107.1:subseq(473451,100206) exonerate:protein2genome:local gene 75357 79932 747 + . gene_id 1 ; sequence SPP00000009_1.0 ; gene_orientation + KV390107.1:subseq(473451,100206) exonerate:protein2genome:local cds 75357 75534 . + . KV390107.1:subseq(473451,100206) exonerate:protein2genome:local exon 75357 75534 . + . insertions 0 ; deletions 0 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local splice5 75535 75536 . + . intron_id 1 ; splice_site "GT" KV390107.1:subseq(473451,100206) exonerate:protein2genome:local intron 75535 77557 . + . intron_id 1 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local splice3 77556 77557 . + . intron_id 0 ; splice_site "AG" KV390107.1:subseq(473451,100206) exonerate:protein2genome:local cds 77558 77675 . + . KV390107.1:subseq(473451,100206) exonerate:protein2genome:local exon 77558 77675 . + . insertions 0 ; deletions 0 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local splice5 77676 77677 . + . intron_id 2 ; splice_site "gt" KV390107.1:subseq(473451,100206) exonerate:protein2genome:local intron 77676 78093 . + . intron_id 2 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local splice3 78092 78093 . + . intron_id 1 ; splice_site "AG" KV390107.1:subseq(473451,100206) exonerate:protein2genome:local cds 78094 78193 . + . KV390107.1:subseq(473451,100206) exonerate:protein2genome:local exon 78094 78193 . + . insertions 0 ; deletions 0 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local splice5 78194 78195 . + . intron_id 3 ; splice_site "GT" KV390107.1:subseq(473451,100206) exonerate:protein2genome:local intron 78194 79731 . + . intron_id 3 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local splice3 79730 79731 . + . intron_id 2 ; splice_site "AG" KV390107.1:subseq(473451,100206) exonerate:protein2genome:local cds 79732 79932 . + . KV390107.1:subseq(473451,100206) exonerate:protein2genome:local exon 79732 79932 . + . insertions 0 ; deletions 0 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local similarity 75357 79932 747 + . alignment_id 1 ; Query SPP00000009_1.0 ; Align 75357 22 177 ; Align 77560 82 114 ; Align 78095 121 99 ; Align 79732 154 201 # --- END OF GFF DUMP --- # C4 Alignment: ------------ Query: SPP00000009_1.0 # Protein # Glutathione peroxidase 6 (GPx6) # Homo sapiens # Complete Target: KV390107.1:subseq(473451,100206) Cebus capucinus imitator isolate Cc_AM_T3 unplaced genomic scaffold Scaffold681, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 990 Query range: 0 -> 216 Target range: 60071 -> 50002 1 : MetPheGlnGlnPheGlnAlaSerCysLeuValLeuPhePheLeuValGlyPheAlaGln : 20 |||.!!!:!|||||||||||||||||||||||||||||||||||||||.!!||||||||| MetIleArgGlnPheGlnAlaSerCysLeuValLeuPhePheLeuValSerPheAlaGln 60071 : ATGATCCGGCAGTTCCAAGCGTCCTGTCTTGTCCTGTTTTTCCTGGTTAGCTTTGCTCAG : 60014 21 : GlnThrLeuLysProGlnAsnArgLys >>>> Target Intron 1 >>>> ValA : 31 |||||||||||||||:!!:!!|||||| 3478 bp |||| GlnThrLeuLysProGluAspArgLys++ ++ValA 60013 : CAGACCCTAAAGCCTGAAGATAGGAAGgt.........................agGTGG : 56503 32 : spCysAsnLysGlyValThrGlyThrIleTyrGluTyrGlyAlaLeuThrLeuAsnGlyG : 51 |||||:!!|||..!|||||||||||||||||||||||||||||||||||||||!:!|||| spCysAspLysAsnValThrGlyThrIleTyrGluTyrGlyAlaLeuThrLeuSerGlyG 56502 : ATTGCGACAAAAACGTGACAGGCACCATCTATGAGTATGGAGCCCTCACCCTCAGCGGCG : 56443 52 : luGluTyrIleGlnPheLysGlnPheAlaGlyLysHisValLeuPheValAsnValAlaA : 71 |||||||||||||||||||||||||||||||||||! !:!!||||||||||||||||||. luGluTyrIleGlnPheLysGlnPheAlaGlyLysProIleLeuPheValAsnValAlaT 56442 : AGGAGTACATCCAGTTCAAGCAGTTTGCAGGCAAGCCCATCCTGTTTGTCAATGTAGCCA : 56383 72 : laTyrUnkGlyLeuAlaAlaGlnTyrPro{G} >>>> Target Intron 2 >>>> : 81 !!||| .!!|||.!!||||||||||||{|} 4236 bp hrTyrCysSerLeuThrAlaGlnTyrPro{G}++ + 56382 : CCTACTGCAGCTTGACAGCTCAGTATCCT{G}gt.........................a : 52118 82 : {lu}LeuAsnAlaLeuGlnGluGluLeuLysAsnPheGlyValIleValLeuAlaPheP : 100 {||}|||||||||||||||||||||||||||!:!.!.|||||||||||||||!.!|||| +{lu}LeuAsnAlaLeuGlnGluGluLeuLysSerLeuGlyValIleValLeuGlyPheP 52117 : g{AA}CTGAATGCGCTCCAGGAGGAGCTGAAGAGCCTGGGAGTCATTGTGCTGGGCTTTC : 52060 101 : roCysAsnGlnPheGlyLysGlnGluProGlyThrAsnSerGluIleLeuLeuGlyLeu{ : 120 ||||||||||||||||||||||||||||||||.!!||||||||||||||||||||||||{ roCysAsnGlnPheGlyLysGlnGluProGlyAlaAsnSerGluIleLeuLeuGlyLeu{ 52059 : CCTGCAACCAGTTTGGAAAACAAGAACCAGGAGCAAACTCAGAGATACTTCTTGGTCTC{ : 52000 121 : Ly} >>>> Target Intron 3 >>>> {s}TyrValCysProGlySerGlyPheV : 129 ||} 500 bp {|}|||||| !!|||||||||||||||| Ly}++ ++{s}TyrValArgProGlySerGlyPheV 51999 : AA}gt.........................ag{G}TATGTGCGTCCAGGTAGCGGCTTTG : 51473 130 : alProSerPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGluGlnLysValP : 149 |||||!:!|||||||||||||||||||||||||||||||||||||||||||||||||||| alProAsnPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGluGlnLysValP 51472 : TACCCAATTTCCAGCTCTTTGAGAAAGGGGATGTGaatggagaaaaagaacagaaggtCT : 51413 150 : heThrPheLeuLys >>>> Target Intron 4 >>>> AsnSerCysProProTh : 159 |||||||||||||| 1207 bp |||:!!||||||||||| heThrPheLeuLys++ ++AsnAlaCysProProTh 51412 : TTACTTTCCTGAAGgt.........................agAACGCCTGCCCTCCGAC : 50176 160 : rSerAspLeuLeuGlySerSerSerGlnLeuPheTrpGluProMetLysValHisAspIl : 179 ||||||||||||||||||||||!:!||||||||||||||||||||||||||||||||||| rSerAspLeuLeuGlySerSerAsnGlnLeuPheTrpGluProMetLysValHisAspIl 50175 : CTCTGATCTTTTGGGCTCATCAAACCAGCTCTTCTGGGAGCCCATGAAGGTCCATGACAT : 50116 180 : eArgTrpAsnPheGluLysPheLeuValGlyProAspGlyValProValMetHisTrpPh : 199 |||||||||||||||||||||||||||||||||| !!|||||||||||||||!.!||||| eArgTrpAsnPheGluLysPheLeuValGlyProHisGlyValProValMetArgTrpPh 50115 : CCGCTGGAACTTTGAAAAATTTCTGGTGGGGCCCCATGGAGTCCCTGTCATGCGCTGGTT : 50056 200 : eHisGlnAlaProValSerThrValLysSerAspIleLeuGluTyrLeuLys : 216 ||||... !||||||||||||||||||||||||||||||! !||||||||| eHisSerIleProValSerThrValLysSerAspIleLeuValTyrLeuLys 50055 : CCACAGTATTCCAGTCAGCACAGTCAAGTCAGACATCCTGGTGTACCTAAAG : 50003 vulgar: SPP00000009_1.0 0 216 . KV390107.1:subseq(473451,100206) 60071 50002 - 990 M 29 87 5 0 2 I 0 3474 3 0 2 M 51 153 S 0 1 5 0 2 I 0 4232 3 0 2 S 1 2 M 38 114 S 0 2 5 0 2 I 0 496 3 0 2 S 1 1 M 33 99 5 0 2 I 0 1203 3 0 2 M 63 189 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-16 ##type DNA # # # seqname source feature start end score strand frame attributes # KV390107.1:subseq(473451,100206) exonerate:protein2genome:local gene 50003 60071 990 - . gene_id 1 ; sequence SPP00000009_1.0 ; gene_orientation + KV390107.1:subseq(473451,100206) exonerate:protein2genome:local cds 59985 60071 . - . KV390107.1:subseq(473451,100206) exonerate:protein2genome:local exon 59985 60071 . - . insertions 0 ; deletions 0 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local splice5 59983 59984 . - . intron_id 1 ; splice_site "GT" KV390107.1:subseq(473451,100206) exonerate:protein2genome:local intron 56507 59984 . - . intron_id 1 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local splice3 56507 56508 . - . intron_id 0 ; splice_site "AG" KV390107.1:subseq(473451,100206) exonerate:protein2genome:local cds 56353 56506 . - . KV390107.1:subseq(473451,100206) exonerate:protein2genome:local exon 56353 56506 . - . insertions 0 ; deletions 0 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local splice5 56351 56352 . - . intron_id 2 ; splice_site "GT" KV390107.1:subseq(473451,100206) exonerate:protein2genome:local intron 52117 56352 . - . intron_id 2 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local splice3 52117 52118 . - . intron_id 1 ; splice_site "AG" KV390107.1:subseq(473451,100206) exonerate:protein2genome:local cds 51999 52116 . - . KV390107.1:subseq(473451,100206) exonerate:protein2genome:local exon 51999 52116 . - . insertions 0 ; deletions 0 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local splice5 51997 51998 . - . intron_id 3 ; splice_site "GT" KV390107.1:subseq(473451,100206) exonerate:protein2genome:local intron 51499 51998 . - . intron_id 3 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local splice3 51499 51500 . - . intron_id 2 ; splice_site "AG" KV390107.1:subseq(473451,100206) exonerate:protein2genome:local cds 51399 51498 . - . KV390107.1:subseq(473451,100206) exonerate:protein2genome:local exon 51399 51498 . - . insertions 0 ; deletions 0 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local splice5 51397 51398 . - . intron_id 4 ; splice_site "GT" KV390107.1:subseq(473451,100206) exonerate:protein2genome:local intron 50192 51398 . - . intron_id 4 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local splice3 50192 50193 . - . intron_id 3 ; splice_site "AG" KV390107.1:subseq(473451,100206) exonerate:protein2genome:local cds 50003 50191 . - . KV390107.1:subseq(473451,100206) exonerate:protein2genome:local exon 50003 50191 . - . insertions 0 ; deletions 0 KV390107.1:subseq(473451,100206) exonerate:protein2genome:local similarity 50003 60071 990 - . alignment_id 1 ; Query SPP00000009_1.0 ; Align 60072 1 87 ; Align 56507 30 153 ; Align 52115 82 114 ; Align 51498 121 99 ; Align 50192 154 189 # --- END OF GFF DUMP --- # -- completed exonerate analysis