Command line: [exonerate -m p2g --showtargetgff -q selenot1_zf.aa.fa.canviat -t genomeselenot1_JXSJ01000006.1.fa] Hostname: [ws116066U] C4 Alignment: ------------ Query: SPP00000649_2.0 # Protein # Selenoprotein T1 (SELENOT1) # Zebrafish Target: JXSJ01000006.1:subseq(6104489,103645) Miichthys miiuy scaffold6, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 754 Query range: 3 -> 178 Target range: 53645 -> 49999 4 : MetArgTrpLeuProPheSerAlaLeuLeuLeuTrpAlaLeuCysLeuHisSerAla<-> : 23 ||||||||||||! !|||||| !||| !||| !! !!!!! ||| !!!!||| MetArgTrpLeuArgPheSerLeuLeuAlaLeuGlyLeuPheSerLeuCysCysAlaThr 53645 : ATGAGGTGGTTGCGTTTTTCGCTCCTCGCGTTGGGGCTCTTCTCGCTATGCTGCGCCACG : 53588 24 : SerAlaAspAsnAsnGlyValLysLysMetLysMetGlnPheAlaThrGlyProLeuLeu : 42 |||!.!|||!:!!:!||||||||||||||||||||||||||||||||||||||||||||| SerGlyAspSerSerGlyValLysLysMetLysMetGlnPheAlaThrGlyProLeuLeu 53587 : TCTGGTGACAGCAGCGGcgtgaagaagatgaagatgcaGTTTGCGACGGGACCCCTCCTC : 53528 43 : LysPheGlnIle{Cy} >>>> Target Intron 1 >>>> {s}ValSerUnkGly : 51 ||||||||||||{||} 916 bp {|}:!!||| ||| LysPheGlnIle{Cy}++ ++{s}IleSer***Gly 53527 : AAGTTTCAAATT{TG}gt.........................ag{C}ATTTCCTGAGGG : 52585 52 : TyrLysArgValPheGluGluTyrThrGlnAlaLeuTyrGlnArgTyrProAspIleArg : 71 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| TyrLysArgValPheGluGluTyrThrGlnAlaLeuTyrGlnArgTyrProAspIleArg 52584 : TACAAGCGGGTGTTTGAGGAGTACACGCAGGCCTTGTACCAGCGGTACCCGGACATCCGC : 52525 72 : IleGluGlyGluAsnTyrLeuProLeuProLeuTyr{Ar} >>>> Target Intron : 84 ||||||||||||||||||||||||:!!|||||||||{||} 1049 bp IleGluGlyGluAsnTyrLeuProIleProLeuTyr{Ar}++ 52524 : ATTGAAGGGGAGAACTATCTTCCTATACCCCTCTAT{CG}gt.................. : 52482 85 : 2 >>>> {g}HisIleAlaSerPheLeuSerMetPheLysLeuLeuLeuIleGlyVal : 100 {|}||||||||||||||||||||||||||||||||||||:!!||||||||| ++{g}HisIleAlaSerPheLeuSerMetPheLysLeuLeuValIleGlyVal 52481 : .......ag{a}CACATTGCTTCCTTCCTGTCTATGTTCAAACTGCTGGTGATTGGGGTG : 51389 101 : IleIleLeuGlyLysAspProPheAlaLeuCysGlyMetGlnAlaProGlyIleTrpVal : 120 ||||||:!!|||||||||||||||||||||! !||||||||||||||||||||||||! ! IleIleIleGlyLysAspProPheAlaLeuPheGlyMetGlnAlaProGlyIleTrpGlu 51388 : ATTATTATCGGCAAGGACCCATTCGCCCTCTTTGGTATGCAAGCCCCGGGCATCTGGGAG : 51329 121 : TrpSerGlnGluAsnLys >>>> Target Intron 3 >>>> IleTyrAlaCysM : 131 |||.!!|||! !|||||| 1041 bp ||||||||||||| TrpGlyGlnGlyAsnLys++ ++IleTyrAlaCysM 51328 : TGGGGCCAGGGAAATAAGgt.........................agATATATGCCTGCA : 50255 132 : etMetValPhePhePheSerAsnMetIleGluAsnGlnCysMetSerThrGlyAlaPheG : 151 ||||||||||||||.!!|||||||||||||||||||||! ||||||||||||||||||| etMetValPhePheLeuSerAsnMetIleGluAsnGlnLeuMetSerThrGlyAlaPheG 50254 : TGATGGTGTTCTTTCTCAGCAATATGATTGAAAACCAGTTAATGTCTACAGGGGCTTTTG : 50195 152 : luIleThrLeuAsn{A} >>>> Target Intron 4 >>>> {sp}ValProValT : 160 ||||||||||||||{|} 112 bp {||}|||||||||| luIleThrLeuAsn{A}++ ++{sp}ValProValT 50194 : AAATCACATTAAAC{g}gt.........................ag{AT}GTGCCAGTGT : 50056 161 : rpSerLysLeuGluSerGlyHisLeuProSerMetGlnGlnLeuValGlnIleLeu : 178 |||||||||||||||||||||||||||||||||||||||||||||||||||||||| rpSerLysLeuGluSerGlyHisLeuProSerMetGlnGlnLeuValGlnIleLeu 50055 : GGTCAAAACTGGAGTCAGGTCACCTCCCCTCCATGCAGCAGCTTGTGCAAATCCTG : 50000 vulgar: SPP00000649_2.0 3 178 . JXSJ01000006.1:subseq(6104489,103645) 53645 49999 - 754 M 19 57 G 0 3 M 24 72 S 0 2 5 0 2 I 0 912 3 0 2 S 1 1 M 36 108 S 0 2 5 0 2 I 0 1045 3 0 2 S 1 1 M 42 126 5 0 2 I 0 1037 3 0 2 M 29 87 S 0 1 5 0 2 I 0 108 3 0 2 S 1 2 M 22 66 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-17 ##type DNA # # # seqname source feature start end score strand frame attributes # JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local gene 50000 53645 754 - . gene_id 1 ; sequence SPP00000649_2.0 ; gene_orientation + JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local cds 53512 53645 . - . JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local exon 53512 53645 . - . insertions 3 ; deletions 0 JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local splice5 53510 53511 . - . intron_id 1 ; splice_site "GT" JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local intron 52596 53511 . - . intron_id 1 JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local splice3 52596 52597 . - . intron_id 0 ; splice_site "AG" JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local cds 52485 52595 . - . JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local exon 52485 52595 . - . insertions 0 ; deletions 0 JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local splice5 52483 52484 . - . intron_id 2 ; splice_site "GT" JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local intron 51436 52484 . - . intron_id 2 JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local splice3 51436 51437 . - . intron_id 1 ; splice_site "ag" JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local cds 51309 51435 . - . JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local exon 51309 51435 . - . insertions 0 ; deletions 0 JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local splice5 51307 51308 . - . intron_id 3 ; splice_site "GT" JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local intron 50268 51308 . - . intron_id 3 JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local splice3 50268 50269 . - . intron_id 2 ; splice_site "AG" JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local cds 50180 50267 . - . JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local exon 50180 50267 . - . insertions 0 ; deletions 0 JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local splice5 50178 50179 . - . intron_id 4 ; splice_site "gt" JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local intron 50068 50179 . - . intron_id 4 JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local splice3 50068 50069 . - . intron_id 3 ; splice_site "AG" JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local cds 50000 50067 . - . JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local exon 50000 50067 . - . insertions 0 ; deletions 0 JXSJ01000006.1:subseq(6104489,103645) exonerate:protein2genome:local similarity 50000 53645 754 - . alignment_id 1 ; Query SPP00000649_2.0 ; Align 53646 4 57 ; Align 53586 23 72 ; Align 52595 48 108 ; Align 51435 85 126 ; Align 50268 127 87 ; Align 50066 157 66 # --- END OF GFF DUMP --- # -- completed exonerate analysis