Command line: [exonerate --exhaustive yes -m p2g --showtargetgff -q /export/home/u114420/treball/fastas/GP1a.x.fa -t GP1a.extraction.KV884777.1.fa] Hostname: [sitdoc] C4 Alignment: ------------ Query: GP1a # Protein # Glutathione peroxidase 1a (GPx1a) # Zebrafish Target: KV884777.1:subseq(695402,100000) Monopterus albus unplaced genomic scaffold scaffold89.1, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 816 Query range: 4 -> 187 Target range: 51143 -> 49661 5 : LysSerPheTyrAspIleThrAlaLysThrLeuThrGlyGluGluPheLysPheSerSer : 24 |||!.!|||:!!!!:.!!||||||||| !|||||||||||| !|||!!.||||||||| LysLysPheHisGluPheThrAlaLysLeuLeuThrGlyGluThrPheAsnPheSerSer 51143 : AAGAAATTTCACGAGTTCACCGCCAAACTCTTAACGGGAGAAACGTTTAATTTTTCCTCC : 51086 25 : LeuGlnGlyLysValValLeuIleGluAsnValAlaSerLeuUnkGlyThrThrSerArg : 44 |||||||||||||||||||||||||||||||||||||||||| ||||||||| !||| LeuGlnGlyLysValValLeuIleGluAsnValAlaSerLeu***GlyThrThrValArg 51085 : CTGCAGGGCAAAGTTGTCCTCATTGAGAATGTGGCGTCTCTTTGAGGTACGACCGTCAGG : 51026 45 : AspTyrThrGlnMetAsnGluLeuHisGluArgPheAlaGluLysGlyLeuValValLeu : 64 |||||||||||||||||||||||||||||||||!:!|||...||||||||||||:!!||| AspTyrThrGlnMetAsnGluLeuHisGluArgTyrAlaSerLysGlyLeuValIleLeu 51025 : GATTACACCCAGATGAACGAGCTCCACGAGCGGTACGCCAGCAAGGGGCTCGTAATCCTG : 50966 65 : GlyValProCysAsnGlnPheGlyTyrGln >>>> Target Intron 1 >>>> G : 75 ||||||||||||||||||||||||:!!||| 933 bp | GlyValProCysAsnGlnPheGlyHisGln++ ++G 50965 : GGAGTGCCCTGCAACCAGTTCGGCCATCAGgt.........................agG : 50000 76 : luAsnCysThrAsnGluGluIleLeuLeuSerLeuLysTyrValArgProGlyAsnGlyT : 95 ||||||||! !||||||||||||||||||||||||||||||||||||||||||||||||! luAsnCysLysAsnGluGluIleLeuLeuSerLeuLysTyrValArgProGlyAsnGlyP 49999 : AGAACTGCAAGAATGAAGAAATTCTTCTGTCTCTCAAGTACGTCCGTCCTGGAAATGGCT : 49940 96 : yrGluProLysPheGlnLeuLeuGluLysValAspValAsnGlyLysAsnAlaHisProL : 115 :!|||||||||||||||||||||||||||||||||||||||||||||:!!|||||||||| heGluProLysPheGlnLeuLeuGluLysValAspValAsnGlyLysAspAlaHisProL 49939 : TTGAGCCCAAGTTTCAGCTCCTGGAGAAGGTAGATGTGAATGGGAAGGATGCCCATCCCT : 49880 116 : euPheThrPheLeuLysGluLysLeuProPheProSerAspGluProMetProPheMetS : 135 |||||..!|||||||||||||||||||||! !|||||||||||||||.!. !!.!.|||! euPheValPheLeuLysGluLysLeuProSerProSerAspGluProPheAlaLeuMetT 49879 : TGTTTGTGTTCCTGAAAGAGAAGCTCCCGTCCCCCAGTGACGAGCCATTCGCCCTGATGA : 49820 136 : erAspProLysPheIleIleTrpSerProValCysArgAsnAspIleAlaTrpAsnPheG : 155 !!||||||||||||||||||||||||||||||||||||||||||:!::!!|||||||||| hrAspProLysPheIleIleTrpSerProValCysArgAsnAspValSerTrpAsnPheG 49819 : CTGACCCGAAGTTTATCATCTGGAGCCCGGTCTGCAGGAATGACGTGTCCTGGAACTTTG : 49760 156 : luLysPheLeuIleGlySerAspGlyValProPheLysArgTyrSerArgArgPheLeuT : 175 ||||||||||||||||| !!||||||! !||||||||||||||||||||||||||||||| luLysPheLeuIleGlyProAspGlyGluProPheLysArgTyrSerArgArgPheLeuT 49759 : AGAAGTTCCTCATTGGGCCTGATGGAGAGCCGTTCAAGCGTTACAGCAGGAGGTTCCTCA : 49700 176 : hrSerGlyIleAspGlyAspIleLysLysLeuLeuSer : 187 |||||! !:!!!!:|||||||||||||||||||||||| hrSerAspLeuGluGlyAspIleLysLysLeuLeuSer 49699 : CCAGCGACCTTGAGGGAGACATCAAGAAGCTCCTCAGC : 49662 vulgar: GP1a 4 187 . KV884777.1:subseq(695402,100000) 51143 49661 - 816 M 70 210 5 0 2 I 0 929 3 0 2 M 113 339 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-22 ##type DNA # # # seqname source feature start end score strand frame attributes # KV884777.1:subseq(695402,100000) exonerate:protein2genome:local gene 49662 51143 816 - . gene_id 1 ; sequence GP1a ; gene_orientation + KV884777.1:subseq(695402,100000) exonerate:protein2genome:local cds 50934 51143 . - . KV884777.1:subseq(695402,100000) exonerate:protein2genome:local exon 50934 51143 . - . insertions 0 ; deletions 0 KV884777.1:subseq(695402,100000) exonerate:protein2genome:local splice5 50932 50933 . - . intron_id 1 ; splice_site "GT" KV884777.1:subseq(695402,100000) exonerate:protein2genome:local intron 50001 50933 . - . intron_id 1 KV884777.1:subseq(695402,100000) exonerate:protein2genome:local splice3 50001 50002 . - . intron_id 0 ; splice_site "AG" KV884777.1:subseq(695402,100000) exonerate:protein2genome:local cds 49662 50000 . - . KV884777.1:subseq(695402,100000) exonerate:protein2genome:local exon 49662 50000 . - . insertions 0 ; deletions 0 KV884777.1:subseq(695402,100000) exonerate:protein2genome:local similarity 49662 51143 816 - . alignment_id 1 ; Query GP1a ; Align 51144 5 210 ; Align 50001 75 339 # --- END OF GFF DUMP --- # -- completed exonerate analysis