Command line: [exonerate -m p2g --showtargetgff -q 1B.query.dir/SP.SelH.1.fa.x.fa -t 3B.fastasub.dir/1.SP.SelH.1.fa.x.fa]
Hostname: [ws116063]

C4 Alignment:
------------
         Query: SPP00000070_2.0 # Protein # Selenoprotein H (SelH) # Human (Homo sapiens)
        Target: gi|398921812|gb|JH767791.1|:subseq(1259632,100522) Ceratotherium simum simum unplaced genomic scaffold scaffold00069, whole genome shotgun sequence
         Model: protein2genome:local
     Raw score: 525
   Query range: 0 -> 122
  Target range: 49933 -> 50522

     1 : MetAlaProArgGlyArgLysArgLysAlaGluAlaAlaValValAlaValAlaGluLys :    20
         |||||| !!|||||||||||||||||||||||||||||||||..!|||   |||||||||
         MetAlaSerArgGlyArgLysArgLysAlaGluAlaAlaValThrAla---AlaGluLys
 49934 : ATGGCTTCCCGCGGGAGGAAGCGGAAGGCGGAGGCGGCGGTCACGGCG---GCCGAGAAG : 49988

    21 : ArgGluLysLeuAlaAsnGlyGlyGluGlyMetGluGluAlaThrValValIleGluHis :    40
         |||||||||||||||..!||||||:!!|||:!!|||||||||||||||||||||||||||
         ArgGluLysLeuAlaGlyGlyGlyLysGlyValGluGluAlaThrValValIleGluHis
 49989 : CGGGAGAAGCTGGCCGGCGGCGGGAAGGGAGTGGAGGAGGCGACGGTCGTTATCGAGCAT : 50048

    41 : {Cy}  >>>> Target Intron 1 >>>>  {s}ThrSerUnkArgValTyrGlyArg :    49
         {||}            87 bp            {|}..!|||   |||||||||||||||
         {Cy}++                         ++{s}ValSer***ArgValTyrGlyArg
 50049 : {TG}gt.........................ag{C}GTAAGCTGACGCGTCTACGGGCGC : 50162

    50 : AsnAlaAlaAlaLeuSerGlnAlaLeuArgLeuGluAlaProGluLeuProValLysVal :    69
         ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
         AsnAlaAlaAlaLeuSerGlnAlaLeuArgLeuGluAlaProGluLeuProValLysVal
 50163 : AACGCCGCGGCCCTGAGCCAGGCGCTGCGCCTGGAGGCCCCCGAGCTGCCAGTGAAGGTG : 50222

    70 : AsnProThrLysProArgArgGlySerPheGluValThrLeuLeuArgProAspGlySer :    89
         ||||||.!!|||||||||||||||||||||||||||||||||||||||||||||||||||
         AsnProAlaLysProArgArgGlySerPheGluValThrLeuLeuArgProAspGlySer
 50223 : AACCCTGCCAAGCCCCGGAGGGGCAGCTTCGAGGTGACGCTGCTGCGCCCCGACGGCAGC : 50282

    90 : {S}  >>>> Target Intron 2 >>>>  {er}AlaGluLeuTrpThrGlyIleLys :    98
         {|}            139 bp           {||}!.!|||||||||||||||||||||
         {S}++                         ++{er}ValGluLeuTrpThrGlyIleLys
 50283 : {A}gt.........................ag{GT}GTGGAGCTCTGGACTGGGATCAAG : 50448

    99 : LysGlyProProArgLysLeuLysPheProGluProGlnGluValValGluGluLeuLys :   118
         ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
         LysGlyProProArgLysLeuLysPheProGluProGlnGluValValGluGluLeuLys
 50449 : AAGGGGCCCCCACGCAAACTCAAGTTCCCTGAGCCTCAAGAGGTGGTGGAGGAGCTGAAG : 50508

   119 : LysTyrLeuSer :   122
         ||||||||||||
         LysTyrLeuSer
 50509 : AAGTACCTTTCG : 50522

vulgar: SPP00000070_2.0 0 122 . gi|398921812|gb|JH767791.1|:subseq(1259632,100522) 49933 50522 + 525 M 16 48 G 1 0 M 23 69 S 0 2 5 0 2 I 0 83 3 0 2 S 1 1 M 48 144 S 0 1 5 0 2 I 0 135 3 0 2 S 1 2 M 32 96
# --- START OF GFF DUMP ---
#
#
##gff-version 2
##source-version exonerate:protein2genome:local 2.2.0
##date 2014-11-25
##type DNA
#
#
# seqname source feature start end score strand frame attributes
#
gi|398921812|gb|JH767791.1|:subseq(1259632,100522)	exonerate:protein2genome:local	gene	49934	50522	525	+	.	gene_id 1 ; sequence SPP00000070_2.0 ; gene_orientation +
gi|398921812|gb|JH767791.1|:subseq(1259632,100522)	exonerate:protein2genome:local	cds	49934	50052	.	+	.	
gi|398921812|gb|JH767791.1|:subseq(1259632,100522)	exonerate:protein2genome:local	exon	49934	50052	.	+	.	insertions 0 ; deletions 1
gi|398921812|gb|JH767791.1|:subseq(1259632,100522)	exonerate:protein2genome:local	splice5	50053	50054	.	+	.	intron_id 1 ; splice_site "GT"
gi|398921812|gb|JH767791.1|:subseq(1259632,100522)	exonerate:protein2genome:local	intron	50053	50139	.	+	.	intron_id 1
gi|398921812|gb|JH767791.1|:subseq(1259632,100522)	exonerate:protein2genome:local	splice3	50138	50139	.	+	.	intron_id 0 ; splice_site "AG"
gi|398921812|gb|JH767791.1|:subseq(1259632,100522)	exonerate:protein2genome:local	cds	50140	50285	.	+	.	
gi|398921812|gb|JH767791.1|:subseq(1259632,100522)	exonerate:protein2genome:local	exon	50140	50285	.	+	.	insertions 0 ; deletions 0
gi|398921812|gb|JH767791.1|:subseq(1259632,100522)	exonerate:protein2genome:local	splice5	50286	50287	.	+	.	intron_id 2 ; splice_site "GT"
gi|398921812|gb|JH767791.1|:subseq(1259632,100522)	exonerate:protein2genome:local	intron	50286	50424	.	+	.	intron_id 2
gi|398921812|gb|JH767791.1|:subseq(1259632,100522)	exonerate:protein2genome:local	splice3	50423	50424	.	+	.	intron_id 1 ; splice_site "AG"
gi|398921812|gb|JH767791.1|:subseq(1259632,100522)	exonerate:protein2genome:local	cds	50425	50522	.	+	.	
gi|398921812|gb|JH767791.1|:subseq(1259632,100522)	exonerate:protein2genome:local	exon	50425	50522	.	+	.	insertions 0 ; deletions 0
gi|398921812|gb|JH767791.1|:subseq(1259632,100522)	exonerate:protein2genome:local	similarity	49934	50522	525	+	.	alignment_id 1 ; Query SPP00000070_2.0 ; Align 49934 1 48 ; Align 49982 18 69 ; Align 50141 42 144 ; Align 50427 91 96
# --- END OF GFF DUMP ---
#
-- completed exonerate analysis
