Command line: [exonerate -m p2g --showtargetgff -q SelP/SelP_U.fa -t SelP/SUBSEQ/KB229292.1_2_subseq.txt --exhaustive yes] Hostname: [luke] C4 Alignment: ------------ Query: SelP Target: KB229292.1:subseq(2319327,25215) Odobenus rosmarus divergens isolate Ivan unplaced genomic scaffold Scaffold415, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 1328 Query range: 0 -> 380 Target range: 25215 -> 17705 1 : MetTrpArgSerLeuGlyLeuAlaLeuAlaLeuCysLeuLeuProSerGlyGlyThrGlu : 20 ||||||||||||||||||||||||||||||||||||||||||||| ||||||...||| MetTrpArgSerLeuGlyLeuAlaLeuAlaLeuCysLeuLeuProTrpGlyGlyAlaGlu 25215 : ATGTGGAGAAGCCTGGGGCTTGCCCTGGCTCTCTGTCTTCTCccatggggaggagcagag : 25158 21 : SerGlnAspGlnSerSerLeuCysLysGlnProProAlaTrpSerIleArgAspGlnAsp : 40 |||||| |||||||||!!!||||||! !|||||||||||||||||| !!||||||:!! SerGlnGlyGlnSerSerPheCysLysProProProAlaTrpSerIleGlyAspGlnAsn 25157 : agccagggACAAAGCTCCTTTTGTAAGCCGCCTCCAGCCTGGAGCATAGGGGATCAAAAT : 25098 41 : ProMetLeuAsnSerAsnGlySerValThrValValAlaLeuLeuGlnAlaSerUnkTyr : 60 ||||||:!!||||||!:!|||||||||||||||||||||||||||||||||||| ||| ProMetValAsnSerSerGlySerValThrValValAlaLeuLeuGlnAlaSer***Tyr 25097 : CCAATGGTGAACTCAAGTGGTTCAGTGACAGTGGTGGCCCTTCTTCAAGCCAGCTGATAC : 25038 61 : LeuCysIleLeuGlnAlaSer{Ly} >>>> Target Intron 1 >>>> {s}Leu : 69 |||||||||||||||||||||{!:} 764 bp {!}||| LeuCysIleLeuGlnAlaSer{Ar}++ ++{g}Leu 25037 : CTGTGCATTCTGCAGGCATCC{AG}gt.........................ag{a}TTG : 24247 70 : GluAspLeuArgValLysLeuLysLysGluGlyTyrSerAsnIleSerTyrIleValVal : 89 ||||||||||||:::||||||:::||||||||||||! !||||||||||||:!!|||||| GluAspLeuArgIleLysLeuGluLysGluGlyTyrLeuAsnIleSerTyrValValVal 24246 : GAAGACCTGCGaataaaactggagaaagaaGGATATTTGAACATCTCTTATGTTGTTGTT : 24187 90 : AsnHisGlnGlyIleSerSerArgLeuLysTyrThrHisLeuLysAsnLysValSerGlu : 109 ||||||||||||:!!||||||:::||||||||| |||||||||||||||||||||||| AsnHisGlnGlyLeuSerSerGlnLeuLysTyrMetHisLeuLysAsnLysValSerGlu 24186 : AATCATCAAGGACTTTCTTctcaattaaaatatatgcatctTAAAAATAAGGTTTCAGAG : 24127 110 : HisIleProValTyrGlnGlnGluGluAsnGlnThrAspValTrpThrLeuLeuAsnGly : 129 :!!||| !!||||||||||||||||||||||||||||||||||||||||||||||||||| TyrIleAlaValTyrGlnGlnGluGluAsnGlnThrAspValTrpThrLeuLeuAsnGly 24126 : TATATTGCTGTTTATCAACaagaagaaaaccaaacagaTGTCTGGACGCTCTTAAATGGG : 24067 130 : SerLysAspAspPheLeuIleTyrAsp{Ar} >>>> Target Intron 2 >>>> : 139 !.!|||!!:||||||||||||||||||{||} 1464 bp LysLysGluAspPheLeuIleTyrAsp{Ar}++ ++ 24066 : AAAAAGGAGGACTTCCTCATATATGAC{AG}gt.........................ag : 22573 140 : {g}CysGlyArgLeuValTyrHisLeuGlyLeuProPheSerPheLeuThrPheProTyr : 158 {|}|||||||||||||||||||||||||||||||||!:!||||||||||||||||||||| {g}CysGlyArgLeuValTyrHisLeuGlyLeuProTyrSerPheLeuThrPheProTyr 22572 : {a}tgtggcCGTCTTGTATACCATCTTGGTTTGCCTTACTCCTTCCTAACTTTCCCGTAT : 22516 159 : ValGluGluAlaIleLysIleAlaTyrCysGluLysLysCysGlyAsnCysSerLeuThr : 178 |||||||||||||||||||||.!!||||||:!!:!!!:!||||||!:!|||||||||||| ValGluGluAlaIleLysIleThrTyrCysLysGluArgCysGlySerCysSerLeuThr 22515 : GTAGAAGAAGCCATTAAAATTACTTACTGTAAAGAGAGGTGTGGAAGCTGCTCTCTCacg : 22456 179 : >>>> Target Intron 3 >>>> ThrLeuLysAspGluAspPheCysLysArgV : 189 4109 bp ||||||:!!! !|||||||||||||||! !| ++ ++ThrLeuGluAlaGluAspPheCysLysMetV 22455 : gt.........................agACCCTGGAAGCTGAAGACTTCTGTAAAATGG : 18314 190 : alSerLeuAlaThrValAspLysThrValGluThrProSerProHisTyrHisHisGluH : 209 ||!!!||||||||||||! ! |||..!|||.!!|||...||||||:::||||||...| alCysLeuAlaThrValGlyAspThrThrGluAlaProGlnProHisHisHisHisHisH 18313 : TATGTTTGGCTACTGTGGGAGACACAACTGAGGCTCCCCAgccacatcaccaccaccacc : 18254 210 : isHisHisAsn<-><-><-><-><-><-><-><->HisGlyHisGlnHisLeuGlySerS : 221 ||||||||::: |||||||||||||||! !|||! !. isHisHisHisHisHisHisArgHisHisHisLysHisGlyHisGlnHisHisGlyIleA 18253 : accaccaccaccaccaccatcgtcATCATCACAAGCATGGGCATCAGCATCATGGGATTG : 18194 222 : erGluLeuSerGluAsnGlnGlnProGlyAlaProAsnAlaProThrHisProAlaProP : 241 .!! !:!!||||||||||||||||||! !||||||:!!.!! !! !||||||.!!||| spGlyValSerGluAsnGlnGlnProGluAlaProAspThrAlaGluHisProThrPro* 18193 : ATGGGGTTTCAGAGAATCAGCAACCAGAAGCCCCAGATACTGCTGAGCATCCCACTCCTT : 18134 242 : roGlyLeuHisHisHisHisLysHisLysGlyGlnHisArgGlnGlyHisProGluAsnA : 261 !|||||||||!.!|||||||||||||||..!!!.!!.|||:!!|||||||||!!:!.! **GlyLeuHisArgHisHisLysHisLysAsnHisGlnArgGluGlyHisProAspThr* 18133 : GAGGCCTCCATCGCCACCATAAGCACAAGAACCACCAAAGAGAGGGCCACCCAGACACCT : 18074 262 : rgAspMetProAla<->SerGluAspLeuGlnAspLeu<-><->GlnLysLysLeuCysA : 278 !!|||:!!|||||| ||||||..!! !||| !! ! ||||||!:!|||!! | **AspLeuProAlaGlySerGluSerSerGlnProSerValProGlnLysArgLeu***A 18073 : GAGATCTGCCAGCAGGAAGTGAAAGTTCACAACCTTCTGTTCCACAAAAGAGGCTCTGAC : 18014 279 : rgLysArgCysIleAsnGlnLeuLeuCysLysLeuProThrAspSerGluLeuAlaProA : 298 ||||| !!||||||||||||||||||||||||||||||! !||||||! !|||||||||! rgLysGlyCysIleAsnGlnLeuLeuCysLysLeuProLysAspSerGlyLeuAlaProA 18013 : GAAAGGGGTGTATAAATCAATTACTCTGCAAATTGCCTAAAGATTCAGGGTTGGCCCCTA : 17954 299 : rgSerUnkCysCysHisCysArgHisLeuIlePheGluLysThrGlySerAlaIleThrU : 318 .!||| |||!! |||!! ||||||||||||||||||||||||||||||||||||||| snSerCysCys***His***ArgHisLeuIlePheGluLysThrGlySerAlaIleThr* 17953 : ACAGCTGTTGCTGACACTGACGACATCTGATATTTGAAAAAACAGGGTCTGCAATCACCT : 17894 319 : nkGlnCysLysGluAsnLeuProSerLeuCysSerUnkGlnGlyLeuArgAlaGluGluA : 338 ||||||||||||!:!|||||||||||||||||| ||||||||| !!|||||||||| **GlnCysLysGluSerLeuProSerLeuCysSer***GlnGlyLeuTrpAlaGluGluA 17893 : GACAATGTAAAGAAAGCCTCCCCTCTCTATGCAGCTGACAGGGACTTTGGGCAGAGGAGA : 17834 339 : snIleThrGluSerCysGlnUnkArgLeuProProAlaAlaUnkGlnIleSerGlnGlnL : 358 ||:!!! !||||||!! ||| |||||||||||||||||| ||| !|||||||||| snValIleGluSer***Gln***ArgLeuProProAlaAla***GlnAlaSerGlnGlnL 17833 : ATGTCATTGAATCTTGACAGTGACGTTTACCTCCAGCTGCCTGACAAGCAAGTCAGCAGC : 17774 359 : euIleProThrGluAlaSerAlaSerUnkArgUnkLysAsnGlnAlaLysLysUnkGluU : 378 || |||||||||||||||.!!!:! !!! |||:!!:!!||||||! ! :!! euGlnProThrGluAlaSerThrAsn***Ser***LysHisLysAlaLysMet***Lys* 17773 : TTCAGCCCACAGAAGCCAGCACCAATTGAAGCTGAAAACATAAGGCCAAAATGTGAAAAT : 17714 379 : nkProSer : 380 |||||| **ProSer 17713 : GACCTTCA : 17706 vulgar: 0 380 . KB229292.1:subseq(2319327,25215) 25215 17705 - 1328 M 67 201 S 0 2 5 0 2 I 0 760 3 0 2 S 1 1 M 70 210 S 0 2 5 0 2 I 0 1460 3 0 2 S 1 1 M 39 117 5 0 2 I 0 4105 3 0 2 M 34 102 G 0 24 M 53 159 G 0 3 M 7 21 G 0 6 M 108 324 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2015-11-25 ##type DNA # # # seqname source feature start end score strand frame attributes # KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local gene 17706 25215 1328 - . gene_id 1 ; sequence ; gene_orientation + KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local cds 25013 25215 . - . KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local exon 25013 25215 . - . insertions 0 ; deletions 0 KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local splice5 25011 25012 . - . intron_id 1 ; splice_site "GT" KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local intron 24249 25012 . - . intron_id 1 KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local splice3 24249 24250 . - . intron_id 0 ; splice_site "ag" KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local cds 24036 24248 . - . KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local exon 24036 24248 . - . insertions 0 ; deletions 0 KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local splice5 24034 24035 . - . intron_id 2 ; splice_site "GT" KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local intron 22572 24035 . - . intron_id 2 KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local splice3 22572 22573 . - . intron_id 1 ; splice_site "ag" KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local cds 22454 22571 . - . KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local exon 22454 22571 . - . insertions 0 ; deletions 0 KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local splice5 22452 22453 . - . intron_id 3 ; splice_site "gt" KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local intron 18345 22453 . - . intron_id 3 KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local splice3 18345 18346 . - . intron_id 2 ; splice_site "AG" KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local cds 17706 18344 . - . KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local exon 17706 18344 . - . insertions 33 ; deletions 0 KB229292.1:subseq(2319327,25215) exonerate:protein2genome:local similarity 17706 25215 1328 - . alignment_id 1 ; Query ; Align 25216 1 201 ; Align 24248 69 210 ; Align 22571 140 117 ; Align 18345 179 102 ; Align 18219 213 159 ; Align 18057 266 21 ; Align 18030 273 324 # --- END OF GFF DUMP --- # -- completed exonerate analysis