Command line: [exonerate -m p2g --showtargetgff -q SELENOH.fa.txt -t SELENOH.subseq.txt] Hostname: [ws116044u] C4 Alignment: ------------ Query: SPP00000070_2.0 # Protein # Selenoprotein H (SELENOH) # Human Target: JSZB01009980.1:subseq(26642,101000) Manis javanica isolate MP_PG03-UM SCAFFOLD11697, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 517 Query range: 0 -> 121 Target range: 50559 -> 49978 1 : MetAlaProArgGlyArgLysArgLysAlaGluAlaAlaValValAlaValAlaGluLys : 20 |||||| !!||||||||||||||||||||||||||||||!.!||||||||| |||||| MetAlaSerArgGlyArgLysArgLysAlaGluAlaAlaAlaValAlaVal---GluLys 50559 : ATGGCTTCCCGCGGGAGGAAGCGCAAGGCGGAGGCAGCGGCGGTGGCGGTC---GAGAAG : 50505 21 : ArgGluLysLeuAlaAsnGlyGlyGluGlyMetGluGluAlaThrValValIleGluHis : 40 |||||||||||||||!:!||| !!||||||:!!|||||||||.!!||||||||||||||| ArgGluLysLeuAlaSerGlyArgGluGlyValGluGluAlaAlaValValIleGluHis 50504 : CGAGAGAAATTGGCTAGCGGCCGAGAGGGAGTGGAGGAGGCGGCCGTCGTTATAGAGCAT : 50445 41 : {Cy} >>>> Target Intron 1 >>>> {s}ThrSerUnkArgValTyrGlyArg : 49 {||} 83 bp {|}|||||| ||||||||||||||| {Cy}++ ++{s}ThrSer***ArgValTyrGlyArg 50444 : {TG}gt.........................ag{C}ACCAGCTGACGCGTCTACGGGCGC : 50335 50 : AsnAlaAlaAlaLeuSerGlnAlaLeuArgLeuGluAlaProGluLeuProValLysVal : 69 ||||||||||||||||||||||||||||||||||||||||||:!!|||||||||!:!||| AsnAlaAlaAlaLeuSerGlnAlaLeuArgLeuGluAlaProGlnLeuProValArgVal 50334 : AACGCCGCGGCCCTGAGCCAGGCGCTGCGCCTGGAGGCCCCGCAGCTGCCAGTGAGGGTG : 50275 70 : AsnProThrLysProArgArgGlySerPheGluValThrLeuLeuArgProAspGlySer : 89 ||||||.!!||||||||||||||||||||||||||||||||||||||||||||||||||| AsnProAlaLysProArgArgGlySerPheGluValThrLeuLeuArgProAspGlySer 50274 : AACCCTGCCAAACCCCGGAGGGGCAGCTTCGAGGTGACGCTGCTGCGCCCGGACGGCAGC : 50215 90 : {S} >>>> Target Intron 2 >>>> {er}AlaGluLeuTrpThrGlyIleLys : 98 {|} 138 bp {||}!.!|||||||||||||||:!!||| {S}++ ++{er}ValGluLeuTrpThrGlyValLys 50214 : {A}gt.........................ag{GT}GTGGAGCTCTGGACGGGGGTTAAG : 50050 99 : LysGlyProProArgLysLeuLysPheProGluProGlnGluValValGluGluLeuLys : 118 ||||||||||||||||||||||||||||||||||||:!!||||||||||||||||||||| LysGlyProProArgLysLeuLysPheProGluProGluGluValValGluGluLeuLys 50049 : AAGGGGCCCCCACGCAAACTCAAGTTCCCTGAGCCTGAAGAGGTGGTGGAGGAGCTGAAG : 49990 119 : LysTyrLeu : 121 ||||||||| LysTyrLeu 49989 : AAGTACCTG : 49979 vulgar: SPP00000070_2.0 0 121 . JSZB01009980.1:subseq(26642,101000) 50559 49978 - 517 M 17 51 G 1 0 M 22 66 S 0 2 5 0 2 I 0 79 3 0 2 S 1 1 M 48 144 S 0 1 5 0 2 I 0 134 3 0 2 S 1 2 M 31 93 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-22 ##type DNA # # # seqname source feature start end score strand frame attributes # JSZB01009980.1:subseq(26642,101000) exonerate:protein2genome:local gene 49979 50559 517 - . gene_id 1 ; sequence SPP00000070_2.0 ; gene_orientation + JSZB01009980.1:subseq(26642,101000) exonerate:protein2genome:local cds 50441 50559 . - . JSZB01009980.1:subseq(26642,101000) exonerate:protein2genome:local exon 50441 50559 . - . insertions 0 ; deletions 1 JSZB01009980.1:subseq(26642,101000) exonerate:protein2genome:local splice5 50439 50440 . - . intron_id 1 ; splice_site "GT" JSZB01009980.1:subseq(26642,101000) exonerate:protein2genome:local intron 50358 50440 . - . intron_id 1 JSZB01009980.1:subseq(26642,101000) exonerate:protein2genome:local splice3 50358 50359 . - . intron_id 0 ; splice_site "AG" JSZB01009980.1:subseq(26642,101000) exonerate:protein2genome:local cds 50212 50357 . - . JSZB01009980.1:subseq(26642,101000) exonerate:protein2genome:local exon 50212 50357 . - . insertions 0 ; deletions 0 JSZB01009980.1:subseq(26642,101000) exonerate:protein2genome:local splice5 50210 50211 . - . intron_id 2 ; splice_site "GT" JSZB01009980.1:subseq(26642,101000) exonerate:protein2genome:local intron 50074 50211 . - . intron_id 2 JSZB01009980.1:subseq(26642,101000) exonerate:protein2genome:local splice3 50074 50075 . - . intron_id 1 ; splice_site "AG" JSZB01009980.1:subseq(26642,101000) exonerate:protein2genome:local cds 49979 50073 . - . JSZB01009980.1:subseq(26642,101000) exonerate:protein2genome:local exon 49979 50073 . - . insertions 0 ; deletions 0 JSZB01009980.1:subseq(26642,101000) exonerate:protein2genome:local similarity 49979 50559 517 - . alignment_id 1 ; Query SPP00000070_2.0 ; Align 50560 1 51 ; Align 50509 19 66 ; Align 50357 42 144 ; Align 50072 91 93 # --- END OF GFF DUMP --- # -- completed exonerate analysis