Command line: [exonerate -m p2g --showtargetgff -q ./SPP00000314_2.0.fa -t ./SPP00000314_2.0/fastasubseq/KV390107.1_fastasubseq.fa] Hostname: [ws101073u] C4 Alignment: ------------ Query: SPP00000314_2.0 # Protein # Glutathione peroxidase 5 (GPx5) # Marmoset Target: KV390107.1:subseq(503168,100218) Cebus capucinus imitator isolate Cc_AM_T3 unplaced genomic scaffold Scaffold681, whole genome shotgun sequence Model: protein2genome:local Raw score: 999 Query range: 0 -> 194 Target range: 45657 -> 50218 1 : PheGlnMetAspCysHisArgAspGlnLysAspThrIleTyrAspTyrGluAlaIleMet : 20 ||||||||||||||||||||||||||||||! !||||||||||||||||||||||||||| PheGlnMetAspCysHisArgAspGlnLysGlyThrIleTyrAspTyrGluAlaIleMet 45658 : ttccagatGGATTGTCACAGAGATCAAAAAGGCACCATCTATGACTATGAGGCCATCATG : 45715 21 : LeuAsnGluLysGluTyrValProPheLysArgTyrValGlyLysHisValLeuPheVal : 40 ||||||:!!||||||||||||||||||||||||||||||||||||||||||||||||||| LeuAsnLysLysGluTyrValProPheLysArgTyrValGlyLysHisValLeuPheVal 45716 : CTTAATAAGAAGGAATACGTTCCCTTCAAGCGGTATGTGGGCAAGCACGTCCTCTTTGTC : 45775 41 : AsnValAlaThrTyrCysGlyLeuThrAlaGlnTyrPro{G} >>>> Target Intr : 54 |||||||||||||||||||||||||||||||||||||||{|} 2023 bp AsnValAlaThrTyrCysGlyLeuThrAlaGlnTyrPro{G}++ 45776 : AACGTGGCTACCTACTGTGGTCTGACAGCACAATATCCT{G}gt................ : 45820 55 : on 1 >>>> {lu}LeuAsnAlaLeuGlnGluGluLeuLysProPheGlyLeuValVal : 69 {||}||||||||||||||||||||||||||||||||||||||||||||| ++{lu}LeuAsnAlaLeuGlnGluGluLeuLysProPheGlyLeuValVal 45821 : .........ag{AA}CTGAATGCACTCCAGGAGGAGCTGAAGCCCTTTGGTCTAGTTGTG : 47885 70 : LeuGlyPheProCysAsnGlnPheGlyLysGlnGluProGlyAspAsnLysGluIleLeu : 89 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| LeuGlyPheProCysAsnGlnPheGlyLysGlnGluProGlyAspAsnLysGluIleLeu 47886 : CTGGGCTTTCCCTGCAACCAGTTTGGAAAGCAAGAACCAGGAGATAACAAAGAGattctt : 47945 90 : ProGlyLeu{Ly} >>>> Target Intron 2 >>>> {s}TyrValArgProGly : 98 |||||||||{||} 418 bp {|}||||||||||||||| ProGlyLeu{Ly}++ ++{s}TyrValArgProGly 47946 : cctgggctc{aa}gt.........................ag{G}TATGTCCGTCCAGGG : 48390 99 : GlyGlyPheValProAsnPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGlu : 118 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| GlyGlyPheValProAsnPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGlu 48391 : GGAGGATTTGTACCTAATTTCCAGCTTTTTGAGAAAGGAGATGTGAATggtgaaaaagaa : 48450 119 : GlnLysValPheSerPheLeuLys >>>> Target Intron 3 >>>> HisAlaC : 129 |||||||||||||||||||||||| 1538 bp |||:!!| GlnLysValPheSerPheLeuLys++ ++HisSerC 48451 : cagaaagtctTCAGCTTCTTGAAGgt.........................agCACTCCT : 50021 130 : ysProHisProSerGluIleLeuGlyThrPheLysSerIleSerTrpAspProValLysV : 149 |||||||||||||||||||||||||||||||||||||||||||||||:!!|||||||||| ysProHisProSerGluIleLeuGlyThrPheLysSerIleSerTrpAsnProValLysV 50022 : GTCCTCATCCCTCTGAGATTTTGGGCACATTCAAATCTATATCCTGGAACCCTGTCAAAG : 50081 150 : alHisAspValArgTrpAsnPheGluLysPheLeuValGlyProAspGlyIleProValM : 169 ||||||||:!!|||||||||||||||||||||||||||||||||:!!||||||||||||| alHisAspIleArgTrpAsnPheGluLysPheLeuValGlyProAsnGlyIleProValM 50082 : TCCATGACATCCGTTGGAACTTTGAAAAGTTCCTGGTGGGGCCCAATGGAATCCCTGTCA : 50141 170 : etArgTrpSerHisGlnAlaThrValSerIleValLysAlaAspIleValAlaAsnLeuL : 189 |||||||||||||||||||||||||||||:!!||||||||||||||||||||| !!|||| etArgTrpSerHisGlnAlaThrValSerLeuValLysAlaAspIleValAlaTyrLeuL 50142 : TGCGGTGGTCCCACCAGGCCACAGTCAGCTTAGTCAAGGCAGACATCGTGGCCTACTTGA : 50201 190 : ysGlnPheLysThrGln : 194 ||||||||||||||:!! ysGlnPheLysThrLys 50202 : AGCAATTCAAAACCAAA : 50218 vulgar: SPP00000314_2.0 0 194 . KV390107.1:subseq(503168,100218) 45657 50218 + 999 M 53 159 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 68 204 # --- 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(503168,100218) exonerate:protein2genome:local gene 45658 50218 999 + . gene_id 1 ; sequence SPP00000314_2.0 ; gene_orientation + KV390107.1:subseq(503168,100218) exonerate:protein2genome:local cds 45658 45817 . + . KV390107.1:subseq(503168,100218) exonerate:protein2genome:local exon 45658 45817 . + . insertions 0 ; deletions 0 KV390107.1:subseq(503168,100218) exonerate:protein2genome:local splice5 45818 45819 . + . intron_id 1 ; splice_site "GT" KV390107.1:subseq(503168,100218) exonerate:protein2genome:local intron 45818 47840 . + . intron_id 1 KV390107.1:subseq(503168,100218) exonerate:protein2genome:local splice3 47839 47840 . + . intron_id 0 ; splice_site "AG" KV390107.1:subseq(503168,100218) exonerate:protein2genome:local cds 47841 47958 . + . KV390107.1:subseq(503168,100218) exonerate:protein2genome:local exon 47841 47958 . + . insertions 0 ; deletions 0 KV390107.1:subseq(503168,100218) exonerate:protein2genome:local splice5 47959 47960 . + . intron_id 2 ; splice_site "gt" KV390107.1:subseq(503168,100218) exonerate:protein2genome:local intron 47959 48376 . + . intron_id 2 KV390107.1:subseq(503168,100218) exonerate:protein2genome:local splice3 48375 48376 . + . intron_id 1 ; splice_site "AG" KV390107.1:subseq(503168,100218) exonerate:protein2genome:local cds 48377 48476 . + . KV390107.1:subseq(503168,100218) exonerate:protein2genome:local exon 48377 48476 . + . insertions 0 ; deletions 0 KV390107.1:subseq(503168,100218) exonerate:protein2genome:local splice5 48477 48478 . + . intron_id 3 ; splice_site "GT" KV390107.1:subseq(503168,100218) exonerate:protein2genome:local intron 48477 50014 . + . intron_id 3 KV390107.1:subseq(503168,100218) exonerate:protein2genome:local splice3 50013 50014 . + . intron_id 2 ; splice_site "AG" KV390107.1:subseq(503168,100218) exonerate:protein2genome:local cds 50015 50218 . + . KV390107.1:subseq(503168,100218) exonerate:protein2genome:local exon 50015 50218 . + . insertions 0 ; deletions 0 KV390107.1:subseq(503168,100218) exonerate:protein2genome:local similarity 45658 50218 999 + . alignment_id 1 ; Query SPP00000314_2.0 ; Align 45658 1 159 ; Align 47843 55 114 ; Align 48378 94 99 ; Align 50015 127 204 # --- END OF GFF DUMP --- # C4 Alignment: ------------ Query: SPP00000314_2.0 # Protein # Glutathione peroxidase 5 (GPx5) # Marmoset Target: KV390107.1:subseq(503168,100218) Cebus capucinus imitator isolate Cc_AM_T3 unplaced genomic scaffold Scaffold681, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 733 Query range: 0 -> 189 Target range: 26795 -> 20285 1 : PheGlnMetAspCysHisArgAspGlnLysAspThrIleTyrAspTyrGluAlaIleMet : 20 ||||||:!!|||||| !!!:!:!! !! !! !|||||||||!!:|||! !|||:!!! PheGlnValAspCysAspLysAsnValThrGlyThrIleTyrGluTyrGlyAlaLeuThr 26795 : TTCCAGGTGGATTGCGACAAAAACGTGACAGGCACCATCTATGAGTATGGAGCCCTCACC : 26738 21 : LeuAsnGluLysGluTyrValProPheLysArgTyrValGlyLysHisValLeuPheVal : 40 |||!:!! :!!||||||:!!! !||||||!:!!:!!.!||||||! !:!!||||||||| LeuSerGlyGluGluTyrIleGlnPheLysGlnPheAlaGlyLysProIleLeuPheVal 26737 : CTCAGCGGCGAGGAGTACATCCAGTTCAAGCAGTTTGCAGGCAAGCCCATCCTGTTTGTC : 26678 41 : AsnValAlaThrTyrCysGlyLeuThrAlaGlnTyrPro{G} >>>> Target Intr : 54 ||||||||||||||||||.!!||||||||||||||||||{|} 4236 bp AsnValAlaThrTyrCysSerLeuThrAlaGlnTyrPro{G}++ 26677 : AATGTAGCCACCTACTGCAGCTTGACAGCTCAGTATCCT{G}gt................ : 26633 55 : on 1 >>>> {lu}LeuAsnAlaLeuGlnGluGluLeuLysProPheGlyLeuValVal : 69 {||}||||||||||||||||||||||||||| !.!.|||:!!:!!||| ++{lu}LeuAsnAlaLeuGlnGluGluLeuLysSerLeuGlyValIleVal 26632 : .........ag{AA}CTGAATGCGCTCCAGGAGGAGCTGAAGAGCCTGGGAGTCATTGTG : 22355 70 : LeuGlyPheProCysAsnGlnPheGlyLysGlnGluProGlyAspAsnLysGluIleLeu : 89 ||||||||||||||||||||||||||||||||||||||||||! |||..!||||||||| LeuGlyPheProCysAsnGlnPheGlyLysGlnGluProGlyAlaAsnSerGluIleLeu 22354 : CTGGGCTTTCCCTGCAACCAGTTTGGAAAACAAGAACCAGGAGCAAACTCAGAGATACTT : 22295 90 : ProGlyLeu{Ly} >>>> Target Intron 2 >>>> {s}TyrValArgProGly : 98 ! !||||||{||} 500 bp {|}||||||||||||||| LeuGlyLeu{Ly}++ ++{s}TyrValArgProGly 22294 : CTTGGTCTC{AA}gt.........................ag{G}TATGTGCGTCCAGGT : 21768 99 : GlyGlyPheValProAsnPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGlu : 118 .!!||||||||||||||||||||||||||||||||||||||||||||||||||||||||| SerGlyPheValProAsnPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGlu 21767 : AGCGGCTTTGTACCCAATTTCCAGCTCTTTGAGAAAGGGGATGTGaatggagaaaaagaa : 21708 119 : GlnLysValPheSerPheLeuLys >>>> Target Intron 3 >>>> HisAlaC : 129 ||||||||||||!!!||||||||| 1207 bp :!!|||| GlnLysValPheThrPheLeuLys++ ++AsnAlaC 21707 : cagaaggtCTTTACTTTCCTGAAGgt.........................agAACGCCT : 20468 130 : ysProHisProSerGluIleLeuGlyThrPheLysSerIleSerTrpAspProValLysV : 149 |||||! !!|||!!::!!||||||:!!! !!...!:!!! !|||!!:|||:!!|||| ysProProThrSerAspLeuLeuGlySerSerAsnGlnLeuPheTrpGluProMetLysV 20467 : GCCCTCCGACCTCTGATCTTTTGGGCTCATCAAACCAGCTCTTCTGGGAGCCCATGAAGG : 20408 150 : alHisAspValArgTrpAsnPheGluLysPheLeuValGlyProAspGlyIleProValM : 169 ||||||||:!!||||||||||||||||||||||||||||||||| !!|||:!!||||||| alHisAspIleArgTrpAsnPheGluLysPheLeuValGlyProHisGlyValProValM 20407 : TCCATGACATCCGCTGGAACTTTGAAAAATTTCTGGTGGGGCCCCATGGAGTCCCTGTCA : 20348 170 : etArgTrpSerHisGlnAlaThrValSerIleValLysAlaAspIleValAlaAsnLeuL : 189 ||||||||! !|||... ! !!||||||! !||||||:!!||||||:!!!.! !!|||| etArgTrpPheHisSerIleProValSerThrValLysSerAspIleLeuValTyrLeuL 20347 : TGCGCTGGTTCCACAGTATTCCAGTCAGCACAGTCAAGTCAGACATCCTGGTGTACCTAA : 20288 190 : ys : 189 || ys 20287 : AG : 20286 vulgar: SPP00000314_2.0 0 189 . KV390107.1:subseq(503168,100218) 26795 20285 - 733 M 53 159 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(503168,100218) exonerate:protein2genome:local gene 20286 26795 733 - . gene_id 1 ; sequence SPP00000314_2.0 ; gene_orientation + KV390107.1:subseq(503168,100218) exonerate:protein2genome:local cds 26636 26795 . - . KV390107.1:subseq(503168,100218) exonerate:protein2genome:local exon 26636 26795 . - . insertions 0 ; deletions 0 KV390107.1:subseq(503168,100218) exonerate:protein2genome:local splice5 26634 26635 . - . intron_id 1 ; splice_site "GT" KV390107.1:subseq(503168,100218) exonerate:protein2genome:local intron 22400 26635 . - . intron_id 1 KV390107.1:subseq(503168,100218) exonerate:protein2genome:local splice3 22400 22401 . - . intron_id 0 ; splice_site "AG" KV390107.1:subseq(503168,100218) exonerate:protein2genome:local cds 22282 22399 . - . KV390107.1:subseq(503168,100218) exonerate:protein2genome:local exon 22282 22399 . - . insertions 0 ; deletions 0 KV390107.1:subseq(503168,100218) exonerate:protein2genome:local splice5 22280 22281 . - . intron_id 2 ; splice_site "GT" KV390107.1:subseq(503168,100218) exonerate:protein2genome:local intron 21782 22281 . - . intron_id 2 KV390107.1:subseq(503168,100218) exonerate:protein2genome:local splice3 21782 21783 . - . intron_id 1 ; splice_site "AG" KV390107.1:subseq(503168,100218) exonerate:protein2genome:local cds 21682 21781 . - . KV390107.1:subseq(503168,100218) exonerate:protein2genome:local exon 21682 21781 . - . insertions 0 ; deletions 0 KV390107.1:subseq(503168,100218) exonerate:protein2genome:local splice5 21680 21681 . - . intron_id 3 ; splice_site "GT" KV390107.1:subseq(503168,100218) exonerate:protein2genome:local intron 20475 21681 . - . intron_id 3 KV390107.1:subseq(503168,100218) exonerate:protein2genome:local splice3 20475 20476 . - . intron_id 2 ; splice_site "AG" KV390107.1:subseq(503168,100218) exonerate:protein2genome:local cds 20286 20474 . - . KV390107.1:subseq(503168,100218) exonerate:protein2genome:local exon 20286 20474 . - . insertions 0 ; deletions 0 KV390107.1:subseq(503168,100218) exonerate:protein2genome:local similarity 20286 26795 733 - . alignment_id 1 ; Query SPP00000314_2.0 ; Align 26796 1 159 ; Align 22398 55 114 ; Align 21781 94 99 ; Align 20475 127 189 # --- END OF GFF DUMP --- # -- completed exonerate analysis