Command line: [exonerate -m p2g --showtargetgff -q SelKHsapiensprot.fa -t sojae_tall.fa --exhaustive yes]
Hostname: [ws101137.aulescm.upf.edu]

C4 Alignment:
------------
         Query: SPP00000016_1.0 # Protein # Selenoprotein K (SelK) # Homo sapiens # Complete
        Target: gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618) Phytophthora sojae strain P6497 Psojae_Cont630, whole genome shotgun sequence
         Model: protein2genome:local
     Raw score: 105
   Query range: 0 -> 88
  Target range: 999 -> 19599

     1 : MetValTyrIleSerAsnGlyGlnValLeuAspSerArgSerGlnSerProTrpArgLeu :    20
         |||..!|||:!:!!!..!|||.!.|||:!!! !!.!|||      |||||||||||||||
         MetThrTyrValThrGlyGlyAspValValAlaLysArg------SerProTrpArgLeu
  1000 : ATGACGTACGTGACCGGCGGCGACGTGGTGGCCAAGCGC------TCGCCGTGGCGGCTC :  1051

    21 : SerLeuIleThrAspPhePheTrpGlyIleAlaGluPheValValLeuPhePhe  >>>> :    39
         |||:!!:!:||||||   |||||||||:!!  !.!.||||||! !|||!!.!:!      
         SerIleValThrAspAlaPheTrpGlyValIleAsnPheValGlyLeuLeuTyr++    
  1052 : TCCATCGTGACGGACGCGTTCTGGGGCGTCATCAACTTTGTCGGACTCTTGTACgt.... :  1110

    40 :  Target Intron 1 >>>>  LysThrLeuLeuGlnGlnAsp{V}  >>>> Target :    46
              14286 bp          :!!  !|||  !|||!:!   { }            1
                              ++GluLeuLeuAlaGlnArgPro{T}-+           
  1111 : .....................agGAGCTGCTCGCTCAACGGCCA{T}ct........... : 15418

    47 :  Intron 2 >>>>  {al}LysLysArgArgSerTyrGlyAsnSerSerAspSerArgT :    60
         24 bp           { !}   ::!|||!:! !!:!!|||.....! !!.!. !!|||:
                       ++{yr}AlaArgArgGlnArgHisGlyArgGlnArgGlnArgArgH
 15419 : ..............ag{AC}GCTCGGAGACAGCGTCATGGCCGGCAACGACAACGGCGCC : 15581

    61 : yrAspAspGlyArgGlyProProGlyAsnProPro{Ar}  >>>> Target Intron  :    72
         !! !!! !! !|||! !|||! !|||.!.! !! !{||}           3932 bp   
         isHisGlyAspArgGluProHisGlyGluArgArg{Ar}++                   
 15582 : ACCATGGTGACCGAGAGCCTCACGGAGAGCGCCGC{CG}gt................... : 15621

    73 : 3 >>>>  {g}ArgMetGlyArgIleAsnHisLeuArgGlyProSerProProProMet :    88
                 {|}  !:!!  !! !:!!|||||||||!:!  !||| !!|||||||||:!:
               -+{g}PheLeu***ThrValAsnHisLeuGlnProProProProProProLeu
 15622 : ......cg{A}TTCTTGTAAACAGTCAATCATTTGCAACCACCTCCACCACCTCCACTA : 19599

vulgar: SPP00000016_1.0 0 88 . gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618) 999 19599 + 105 M 13 39 G 2 0 M 23 69 5 0 2 I 0 14282 3 0 2 M 7 21 S 0 1 5 0 2 I 0 120 3 0 2 S 1 2 M 25 75 S 0 2 5 0 2 I 0 3928 3 0 2 S 1 1 M 16 48
# --- START OF GFF DUMP ---
#
#
##gff-version 2
##source-version exonerate:protein2genome:local 2.2.0
##date 2010-03-04
##type DNA
#
#
# seqname source feature start end score strand frame attributes
#
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	gene	1000	19599	105	+	.	gene_id 1 ; sequence SPP00000016_1.0 ; gene_orientation +
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	cds	1000	1107	.	+	.	
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	exon	1000	1107	.	+	.	insertions 0 ; deletions 2
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	splice5	1108	1109	.	+	.	intron_id 1 ; splice_site "GT"
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	intron	1108	15393	.	+	.	intron_id 1
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	splice3	15392	15393	.	+	.	intron_id 0 ; splice_site "AG"
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	cds	15394	15415	.	+	.	
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	exon	15394	15415	.	+	.	insertions 0 ; deletions 0
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	splice5	15416	15417	.	+	.	intron_id 2 ; splice_site "CT"
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	intron	15416	15539	.	+	.	intron_id 2
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	splice3	15538	15539	.	+	.	intron_id 1 ; splice_site "AG"
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	cds	15540	15618	.	+	.	
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	exon	15540	15618	.	+	.	insertions 0 ; deletions 0
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	splice5	15619	15620	.	+	.	intron_id 3 ; splice_site "GT"
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	intron	15619	19550	.	+	.	intron_id 3
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	splice3	19549	19550	.	+	.	intron_id 2 ; splice_site "CG"
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	cds	19551	19599	.	+	.	
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	exon	19551	19599	.	+	.	insertions 0 ; deletions 0
gi|113926118|gb|AAQY01000630.1|:subseq(30511,32618)	exonerate:protein2genome:local	similarity	1000	19599	105	+	.	alignment_id 1 ; Query SPP00000016_1.0 ; Align 1000 1 39 ; Align 1039 16 69 ; Align 15394 39 21 ; Align 15542 47 75 ; Align 19552 73 48
# --- END OF GFF DUMP ---
#
-- completed exonerate analysis