Command line: [exonerate --exhaustive yes -m p2g --showtargetgff -q ../GPx6-human.fa -t GPx6-NIVO01038789.1.subseq] Hostname: [sitdoc] C4 Alignment: ------------ Query: SPP00000009_1.0 # Protein # Glutathione peroxidase 6 (GPx6) # Homo sapiens # Complete Target: NIVO01038789.1:subseq(406729,133097) Ammotragus lervia isolate NAZ scaffold232_len768847_cov70, whole genome shotgun sequence Model: protein2genome:local Raw score: 818 Query range: 0 -> 220 Target range: 75121 -> 83097 1 : MetPheGlnGlnPheGlnAlaSerCysLeuValLeuPhePheLeuValGlyPheAlaGln : 20 ||| ! !|||!!.::!|||||||||||| |||...||||||..!||||||!.!||| MetThrValGlnLeuArgAlaSerCysLeuPheLeuLeuPheLeuAlaGlyPheValGln 75122 : ATGACTGTACAGTTAAGGGCCtcttgtctttttctccttttcctagcTGGCTTTGTGCAA : 75179 21 : GlnThrLeuLysProGlnAsn >>>> Target Intron 1 >>>> ArgLys >> : 30 !!.!! !!..! !:!!!!. 230 bp !:!! ThrAsnTyrSerLeuGluLys++ ++ValGlu++ 75180 : ACAAATTACAGTCTAGAGAAGgt.........................agGTAGAGgt.. : 75441 31 : >> Target Intron 2 >>>> ValAspCysAsnLysGlyValThrGlyThrIleTy : 41 2570 bp :!!|||||| !!|||! !|||! !||||||||||| ++MetAspCysTyrLysAspValLysGlyThrIleTy 75442 : .......................agATGGATTGCTACAAAGATGTGAAAGGCACCATCTA : 78042 42 : rGluTyrGlyAlaLeuThrLeuAsnGlyGluGluTyrIleGlnPheLysGlnPheAlaGl : 61 |!!:|||! !|||!!!||||||||||||:!!|||:!!|||||||||||||||!:!||||| rAspTyrAspAlaPheThrLeuAsnGlyLysGluHisIleGlnPheLysGlnTyrAlaGl 78043 : TGATTATGACGCTTTTACTCTTAATGGAAAGGAACACATTCAGTTCAAGCAGTATGCGGG : 78102 62 : yLysHisValLeuPheValAsnValAlaAlaTyrUnkGlyLeuAlaAlaGlnTyrPro{G : 81 ||||||||||||||||||||||||||||.!!||| ||||||.!!||||||||||||{| yLysHisValLeuPheValAsnValAlaThrTyrCysGlyLeuThrAlaGlnTyrPro{G 78103 : CAAACATGTCCTTTTTGTCAATGTGGCCACCTATTGTGGTCTGACAGCTCAATATCCT{G : 78162 82 : } >>>> Target Intron 3 >>>> {lu}LeuAsnAlaLeuGlnGluGluLeuLy : 90 } 2920 bp {||}|||||||||||||||||||||||||| }++ ++{lu}LeuAsnAlaLeuGlnGluGluLeuLy 78163 : }gt.........................ag{AA}CTGAATGCACTACAGGAGGAGCTGAA : 81109 91 : sAsnPheGlyValIleValLeuAlaPheProCysAsnGlnPheGlyLysGlnGluProGl : 110 | !||||||:!!:!!||||||!.!||||||||||||||||||||||||||||||||||| sProPheGlyLeuValValLeuGlyPheProCysAsnGlnPheGlyLysGlnGluProGl 81110 : GCCCTTTGGCTTAGTTGTGTTGGGATTTCCCTGTAACCAATTCGGAAAGCAAGAACCAGG : 81169 111 : yThrAsnSerGluIleLeuLeuGlyLeu{Ly} >>>> Target Intron 4 >>>> : 120 | !|||||||||||||||! !||||||{||} 481 bp yGluAsnSerGluIleLeuProGlyLeu{Ly}++ + 81170 : AGAAAACTCAGAAATCCTTCCAGGACTG{AA}gt.........................a : 81680 121 : {s}TyrValCysProGlySerGlyPheValProSerPheGlnLeuPheGluLysGlyAs : 139 {|}|||||| !!||||||.!!|||!:!||||||!:!||||||||||||||||||||||| +{s}TyrValArgProGlyGlyGlyTyrValProAsnPheGlnLeuPheGluLysGlyAs 81681 : g{G}TATGTCCGTCCAGGAGGAGGATACGTACCTAATTTCCAGCTGTTTGAGAAAGGGGA : 81737 140 : pValAsnGlyGluLysGluGlnLysValPheThrPheLeuLys >>>> Target Int : 154 |||||||||||||! !||||||||||||||||||||||||||| 1115 b pValAsnGlyGluThrGluGlnLysValPheThrPheLeuLys++ 81738 : TGTGAATGGTGAAACAGAGCAGAAAGTCTTCACCTTCTTAAAGgt............... : 81784 155 : ron 5 >>>> AsnSerCysProProThrSerAspLeuLeuGlySerSerSerGlnLeu : 169 p .!.|||||||||! ! !!|||!!:!!!:!!||||||! !!.!!!.:!! ++GlnSerCysProHisProSerGluPheMetGlySerIleLysHisIle 81785 : ..........agCAATCCTGTCCTCACCCTTCTGAGTTTATGGGCTCAATCAAACATATA : 82942 170 : PheTrpGluProMetLysValHisAspIleArgTrpAsnPheGluLysPheLeuValGly : 189 ! !|||||||||!!:! !|||!.!|||||||||||||||||||||||||||||||||||| SerTrpGluProIleMetValArgAspIleArgTrpAsnPheGluLysPheLeuValGly 82943 : TCCTGGGAACCCATCATGGTCCGTGACATCCGCTGGAATTTTGAAAAGTTCCTAGTGGGA : 83002 190 : ProAspGlyValProValMetHisTrpPheHisGlnAlaProValSerThrValLysSer : 209 |||||||||||||||||||||!.!|||||||||!:!.!!||||||||||||||||||||| ProAspGlyValProValMetArgTrpPheHisArgThrProValSerThrValLysSer 83003 : CCTGATGGCGTCCCTGTCATGCGCTGGTTTCATCGGACTCCAGTCAGTACAGTCAAGTCA : 83062 210 : AspIleLeuGluTyrLeuLysGlnPheAsnThr : 220 |||||||||! !|||:!!|||||||||.!.||| AspIleLeuAlaTyrMetLysGlnPheGlnThr 83063 : GATATTCTGGCATACATGAAACAGTTCCAAACC : 83097 vulgar: SPP00000009_1.0 0 220 . NIVO01038789.1:subseq(406729,133097) 75121 83097 + 818 M 27 81 5 0 2 I 0 226 3 0 2 M 2 6 5 0 2 I 0 2566 3 0 2 M 51 153 S 0 1 5 0 2 I 0 2916 3 0 2 S 1 2 M 38 114 S 0 2 5 0 2 I 0 477 3 0 2 S 1 1 M 33 99 5 0 2 I 0 1111 3 0 2 M 67 201 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-17 ##type DNA # # # seqname source feature start end score strand frame attributes # NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local gene 75122 83097 818 + . gene_id 1 ; sequence SPP00000009_1.0 ; gene_orientation + NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local cds 75122 75202 . + . NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local exon 75122 75202 . + . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice5 75203 75204 . + . intron_id 1 ; splice_site "GT" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local intron 75203 75432 . + . intron_id 1 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice3 75431 75432 . + . intron_id 0 ; splice_site "aG" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local cds 75433 75438 . + . NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local exon 75433 75438 . + . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice5 75439 75440 . + . intron_id 2 ; splice_site "GT" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local intron 75439 78008 . + . intron_id 2 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice3 78007 78008 . + . intron_id 1 ; splice_site "ag" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local cds 78009 78162 . + . NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local exon 78009 78162 . + . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice5 78163 78164 . + . intron_id 3 ; splice_site "GT" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local intron 78163 81082 . + . intron_id 3 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice3 81081 81082 . + . intron_id 2 ; splice_site "AG" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local cds 81083 81200 . + . NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local exon 81083 81200 . + . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice5 81201 81202 . + . intron_id 4 ; splice_site "GT" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local intron 81201 81681 . + . intron_id 4 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice3 81680 81681 . + . intron_id 3 ; splice_site "AG" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local cds 81682 81781 . + . NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local exon 81682 81781 . + . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice5 81782 81783 . + . intron_id 5 ; splice_site "GT" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local intron 81782 82896 . + . intron_id 5 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice3 82895 82896 . + . intron_id 4 ; splice_site "AG" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local cds 82897 83097 . + . NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local exon 82897 83097 . + . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local similarity 75122 83097 818 + . alignment_id 1 ; Query SPP00000009_1.0 ; Align 75122 1 81 ; Align 75433 28 6 ; Align 78009 30 153 ; Align 81085 82 114 ; Align 81683 121 99 ; Align 82897 154 201 # --- END OF GFF DUMP --- # C4 Alignment: ------------ Query: SPP00000009_1.0 # Protein # Glutathione peroxidase 6 (GPx6) # Homo sapiens # Complete Target: NIVO01038789.1:subseq(406729,133097) Ammotragus lervia isolate NAZ scaffold232_len768847_cov70, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 968 Query range: 0 -> 220 Target range: 56726 -> 49999 1 : MetPheGlnGlnPheGlnAlaSerCysLeuValLeuPhePheLeuValGlyPheAlaGln : 20 |||.!!! !|||.!! !|||||||||||| !!! !.!.|||||||||.!!!!.|||||| MetIleLeuGlnLeuTrpAlaSerCysLeuPheProLeuPheLeuValSerLeuAlaGln 56726 : ATGATCCTGCAGCTCTGGGCCTCCTGTCTTTTCCCACTGTTCCTGGTCAGCTTAGCTCAG : 56669 21 : GlnThrLeuLysProGlnAsnArgLys >>>> Target Intron 1 >>>> ValA : 31 ! !|||! !|||! !|||.!.! !||| 2330 bp |||| LeuThrProLysGlnGlnGlnMetLys++ ++ValA 56668 : CTGACTCCGAAGCAGCAACAGATGAAGgt.........................agGTGG : 54306 32 : spCysAsnLysGlyValThrGlyThrIleTyrGluTyrGlyAlaLeuThrLeuAsnGlyG : 51 ||||| !!|||||||||||||||||||||||||||||||||||||||||||||||||||| spCysTyrLysGlyValThrGlyThrIleTyrGluTyrGlyAlaLeuThrLeuAsnGlyG 54305 : ATTGCTATAAGGGGGTGACAGGCACCATCTATGAGTATGGAGCCCTCACCCTCAATGGTG : 54246 52 : luGluTyrIleGlnPheLysGlnPheAlaGlyLysHisValLeuPheValAsnValAlaA : 71 |||||||||||||||||||||||!:!|||||||||||||||||||||||||||||||||. luGluTyrIleGlnPheLysGlnTyrAlaGlyLysHisValLeuPheValAsnValAlaT 54245 : AGGAGTATATCCAGTTCAAGCAGTATGCAGGCAAGCATGTCCTGTTTGTCAATGTGGCCA : 54186 72 : laTyrUnkGlyLeuAlaAlaGlnTyrPro{G} >>>> Target Intron 2 >>>> : 81 !!||| ||||||.!!||||||||||||{|} 1830 bp hrTyr***GlyLeuThrAlaGlnTyrPro{G}++ + 54185 : CCTATTGAGGCTTGACAGCTCAGTATCCA{G}gt.........................a : 52327 82 : {lu}LeuAsnAlaLeuGlnGluGluLeuLysAsnPheGlyValIleValLeuAlaPheP : 100 {||}||||||||||||||||||||||||||| !|||||||||:!!||||||!.!|||| +{lu}LeuAsnAlaLeuGlnGluGluLeuLysProPheGlyValValValLeuGlyPheP 52326 : g{AA}CTGAATGCACTACAGGAGGAGCTGAAGCCTTTCGGTGTAGTTGTGTTGGGCTTTC : 52269 101 : roCysAsnGlnPheGlyLysGlnGluProGlyThrAsnSerGluIleLeuLeuGlyLeu{ : 120 |||||||||||||||||||||||||||||!.!! !|||||||||||||||:!!||||||{ roCysAsnGlnPheGlyLysGlnGluProAlaLysAsnSerGluIleLeuMetGlyLeu{ 52268 : CCTGCAACCAATTTGGAAAACAAGAGCCAGCAAAGAACTCAGAGATCCTTATGGGGCTC{ : 52209 121 : Ly} >>>> Target Intron 3 >>>> {s}TyrValCysProGlySerGlyPheV : 129 ||} 1156 bp {|}|||||| !!||||||.!!||||||| Ly}++ ++{s}TyrValArgProGlyGlyGlyPheV 52208 : AA}gt.........................ag{G}TATGTACGTCCAGGTGGTGGCTTTG : 51026 130 : alProSerPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGluGlnLysValP : 149 |||||!:!|||||||||||||||||||||||||||||||||||||||||||||||||||| alProAsnPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGluGlnLysValP 51025 : TTCCCAATTTTCAGCTCTTTGAGAAAGGGGATGTGaatggagaaaaagaacagaaggtCT : 50966 150 : heThrPheLeuLys >>>> Target Intron 4 >>>> AsnSerCysProProTh : 159 |||||||||||||| 751 bp ||||||||||||||||| heThrPheLeuLys++ ++AsnSerCysProProTh 50965 : TTACCTTCCTGAAGgt.........................agAACTCCTGCCCTCCAAC : 50185 160 : rSerAspLeuLeuGlySerSerSerGlnLeuPheTrpGluProMetLysValHisAspIl : 179 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| rSerAspLeuLeuGlySerSerSerGlnLeuPheTrpGluProMetLysValHisAspIl 50184 : CTCTGATCTTTTGGGCTCATCAAGTCAGCTCTTCTGGGAGCCCATGAAGGTCCATGACAT : 50125 180 : eArgTrpAsnPheGluLysPheLeuValGlyProAspGlyValProValMetHisTrpPh : 199 ||||||||||||||||||||||||||||||||||||||||||||||||||||!.!|||!: eArgTrpAsnPheGluLysPheLeuValGlyProAspGlyValProValMetArgTrpTy 50124 : TCGCTGGAACTTTGAGAAGTTTCTAGTGGGGCCTGATGGTGTCCCTGTCATGCGTTGGTA : 50065 200 : eHisGlnAlaProValSerThrValLysSerAspIleLeuGluTyrLeuLysGlnPheAs : 219 !|||!:!||| !!|||||||||||||||||||||!!:|||||||||||||||||||||!! rHisArgAlaSerValSerThrValLysSerAspMetLeuGluTyrLeuLysGlnPheLy 50064 : CCACCGGGCTTCAGTCAGCACAGTTAAATCAGACATGCTGGAGTACCTGAAGCAGTTCAA : 50005 220 : nThr : 220 .!:! sSer 50004 : AAGC : 50000 vulgar: SPP00000009_1.0 0 220 . NIVO01038789.1:subseq(406729,133097) 56726 49999 - 968 M 29 87 5 0 2 I 0 2326 3 0 2 M 51 153 S 0 1 5 0 2 I 0 1826 3 0 2 S 1 2 M 38 114 S 0 2 5 0 2 I 0 1152 3 0 2 S 1 1 M 33 99 5 0 2 I 0 747 3 0 2 M 67 201 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-17 ##type DNA # # # seqname source feature start end score strand frame attributes # NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local gene 50000 56726 968 - . gene_id 1 ; sequence SPP00000009_1.0 ; gene_orientation + NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local cds 56640 56726 . - . NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local exon 56640 56726 . - . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice5 56638 56639 . - . intron_id 1 ; splice_site "GT" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local intron 54310 56639 . - . intron_id 1 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice3 54310 54311 . - . intron_id 0 ; splice_site "ag" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local cds 54156 54309 . - . NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local exon 54156 54309 . - . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice5 54154 54155 . - . intron_id 2 ; splice_site "GT" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local intron 52326 54155 . - . intron_id 2 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice3 52326 52327 . - . intron_id 1 ; splice_site "AG" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local cds 52208 52325 . - . NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local exon 52208 52325 . - . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice5 52206 52207 . - . intron_id 3 ; splice_site "GT" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local intron 51052 52207 . - . intron_id 3 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice3 51052 51053 . - . intron_id 2 ; splice_site "AG" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local cds 50952 51051 . - . NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local exon 50952 51051 . - . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice5 50950 50951 . - . intron_id 4 ; splice_site "GT" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local intron 50201 50951 . - . intron_id 4 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local splice3 50201 50202 . - . intron_id 3 ; splice_site "AG" NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local cds 50000 50200 . - . NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local exon 50000 50200 . - . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406729,133097) exonerate:protein2genome:local similarity 50000 56726 968 - . alignment_id 1 ; Query SPP00000009_1.0 ; Align 56727 1 87 ; Align 54310 30 153 ; Align 52324 82 114 ; Align 51051 121 99 ; Align 50201 154 201 # --- END OF GFF DUMP --- # -- completed exonerate analysis