Command line: [exonerate -m p2g --showtargetgff -q gpx3b_zf.aa.fa.canviat -t genomegpx3b_JXSJ01000009.1.fa] Hostname: [ws101097u] C4 Alignment: ------------ Query: SPP00000623_2.0 # Protein # Glutathione peroxidase 3b (GPx3b) # Zebrafish Target: JXSJ01000009.1:subseq(4832795,101397) Miichthys miiuy scaffold9, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 597 Query range: 0 -> 161 Target range: 51397 -> 49999 1 : MetAlaAlaAsnSerValTyrAspPheThrValGluThrLeuGluGluLysSerValSer : 20 |||||| !!!.|||:!!|||||||||:!!!.!|||||||||!!:! !:!! !||| !! MetAlaAsnLysSerIleTyrAspPheSerAlaGluThrLeuAspGlyGlnLeuValPro 51397 : ATGGCGAATAAATCAATCTATGATTTCTCTGCTGAGACCCTGGATGGACAGCTGGTTCCG : 51340 21 : LeuAspIlePheArgGlyLysValLeuLeuIleMetAsnValAlaThrPheUnkGlySer : 40 |||..!! !!:!|||||||||||||||||||||!!:||||||||||||||| |||||| LeuSerAsnTyrArgGlyLysValLeuLeuIleIleAsnValAlaThrPhe***GlySer 51339 : CTCAGTAACTACAGGGGTAAGGTGCTTCTCATCATCAATGTTGCCACCTTCTGAGGGTCA : 51280 41 : ThrGln >>>> Target Intron 1 >>>> TyrHisArgMetAsnAlaLeuMetA : 51 ||| ! 120 bp |||||||||:!!||||||||||||! ThrIle+- ++TyrHisArgLeuAsnAlaLeuMetG 51279 : ACAATAga.........................agTACCACAGACTGAATGCACTTATGG : 51127 52 : spMetTyrGlyGlyGlnArgPheThrValLeuGlyPheProCysAsnGlnPheGlyLeuG : 71 !:|||!:!|||! ! ||||||||||||||||||||| !!||||||||||||||||||| luMetPheGlyAspIleArgPheThrValLeuGlyPheSerCysAsnGlnPheGlyLeuG 51126 : AAATGTTTGGCGACATCAGATTTACTGTCCTTGGATTCTCCTGCAACCAATTCGGCCTTC : 51067 72 : lnSerPro{G} >>>> Target Intron 2 >>>> {ly}GluAsnHisGluThrL : 80 ||||||||{!} 715 bp { !}! !||||||||||||| lnSerPro{G}++ ++{lu}ValAsnHisGluThrL 51066 : AGTCACCT{G}gt.........................ag{AG}GTGAATCATGAAACCC : 50325 81 : euAsnValLeuGlnHisValArgProGlySerGlyPheLeuProLysPheProIlePheS : 100 |||||:!!|||:!!:!!||||||||||||.!!||||||:!!||||||||||||||||||: euAsnIleLeuLysTyrValArgProGlyGlyGlyPheValProLysPheProIlePheA 50324 : TTAATATTCTTAAGTATGTGAGACCTGGTGGAGGGTTTGTGCCAAAGTTTCCTATTTTCG : 50265 101 : erArgIleGluValAsnGlySerAspGluAspProLeuTyrAlaTyrLeuLys >>>> : 118 !!|||:!!||||||||||||! !:!!|||!!:||||||!:!.!!||||||||| laArgValGluValAsnGlyLeuAsnGluGluProLeuPheThrTyrLeuLys++ 50264 : CCAGAGTTGAAGTGAATGGATTAAATGAAGAGCCTCTTTTCACCTATCTAAAGgt..... : 50209 119 : Target Intron 3 >>>> GluSerLeuProPheValAsnProValIleGlyAspIl : 130 80 bp |||||||||||||||||||||||||||||||||||||| ++GluSerLeuProPheValAsnProValIleGlyAspIl 50208 : ....................agGAATCTTTGCCATTTGTGAATCCCGTAATTGGAGACAT : 50095 131 : eArgLysLeuTyrTrpSerProIleLysAlaAsnAspValArgTrpAsnPheGluLysPh : 150 |!:!|||!!!|||||||||||||||!!.!.!||||||:!!|||||||||||||||||||| eLysLysPheTyrTrpSerProIleAsnValAsnAspIleArgTrpAsnPheGluLysPh 50094 : AAAGAAATTCTACTGGTCTCCAATCAATGTCAATGACATCCGCTGGAATTTTGAGAAGTT : 50035 151 : eLeuIleThrAlaAspGlyArgProTyrLysArg : 161 |||||||||||||||||||! !|||!:!!:!||| eLeuIleThrAlaAspGlyMetProPheArgArg 50034 : TCTAATAACCGCTGATGGCATGCCCTTCAGGAGG : 50000 vulgar: SPP00000623_2.0 0 161 . JXSJ01000009.1:subseq(4832795,101397) 51397 49999 - 597 M 42 126 5 0 2 I 0 116 3 0 2 M 31 93 S 0 1 5 0 2 I 0 711 3 0 2 S 1 2 M 43 129 5 0 2 I 0 76 3 0 2 M 44 132 # --- 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 # JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local gene 50000 51397 597 - . gene_id 1 ; sequence SPP00000623_2.0 ; gene_orientation + JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local cds 51272 51397 . - . JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local exon 51272 51397 . - . insertions 0 ; deletions 0 JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local splice5 51270 51271 . - . intron_id 1 ; splice_site "GA" JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local intron 51152 51271 . - . intron_id 1 JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local splice3 51152 51153 . - . intron_id 0 ; splice_site "AG" JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local cds 51058 51151 . - . JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local exon 51058 51151 . - . insertions 0 ; deletions 0 JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local splice5 51056 51057 . - . intron_id 2 ; splice_site "GT" JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local intron 50343 51057 . - . intron_id 2 JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local splice3 50343 50344 . - . intron_id 1 ; splice_site "AG" JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local cds 50212 50342 . - . JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local exon 50212 50342 . - . insertions 0 ; deletions 0 JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local splice5 50210 50211 . - . intron_id 3 ; splice_site "gt" JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local intron 50132 50211 . - . intron_id 3 JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local splice3 50132 50133 . - . intron_id 2 ; splice_site "AG" JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local cds 50000 50131 . - . JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local exon 50000 50131 . - . insertions 0 ; deletions 0 JXSJ01000009.1:subseq(4832795,101397) exonerate:protein2genome:local similarity 50000 51397 597 - . alignment_id 1 ; Query SPP00000623_2.0 ; Align 51398 1 126 ; Align 51152 43 93 ; Align 50341 75 129 ; Align 50132 118 132 # --- END OF GFF DUMP --- # -- completed exonerate analysis