Command line: [exonerate --exhaustive yes -m p2g --showtargetgff -q treball/run/DIO3b/DIO3b.x.fa -t treball/run/DIO3b/DIO3b.extraction.KV884701.1.fa] Hostname: [sitdoc] C4 Alignment: ------------ Query: DIO3b # Protein # Iodothyronine deiodinase 3b (DIO3b) # Zebrafish Target: KV884701.1:subseq(3433220,100000) Monopterus albus unplaced genomic scaffold scaffold13.1, whole genome shotgun sequence Model: protein2genome:local Raw score: 846 Query range: 0 -> 253 Target range: 49999 -> 50740 1 : ArgAlaLeuLysAsnAlaLeuValCysLeuLeuIleLeuProArgPheLeuValAlaAla : 20 !:!|||:!!|||||||||:!!|||||||||:!!:!:||||||||||||||||||!.!||| LysAlaIleLysAsnAlaValValCysLeuValLeuLeuProArgPheLeuValValAla 50000 : AAAGCTATTAAAAATGCAGTAGTCTGCCTGGTCCTGCTGCCCCGTTTTCTGGTGGTAGCA : 50057 21 : PheMetLeuTrpCysLeuAspPheLeuCysValArgLysArgValLeuValHisLeuGln : 40 !!!!:|||||| |||||||||||||||:!!||||||||||||!!!:!! ValIleLeuTrpLeuLeuAspPheLeuCysIleArgLysArgValPheLeu--------- 50058 : GTCATATTGTGGCTACTTGACTTTTTATGCATAAGGAAAAGGGTTTTCCTA--------- : 50111 41 : GluArgAlaGluGluGluGluGluAspAlaGluGluGluGluGluProLeuCysIleSer : 60 ||| !:!!! !:!!|||! |||||| !!: ... |||||||||||| ---ArgMetLysAlaGlnGluGlyAspAlaIleAspProSer------LeuCysIleSer 50112 : ---AGGATGAAGGCGCAGGAGGGCGATGCCATCGATCCTTCC------CTCTGCATATCG : 50159 61 : AspSerSerArgMetPheSerTrpGluSerLeuLysAlaValPheHisGlyHisLysLeu : 80 ||||||!:!|||:!:|||||| !||||||||||||||||||!::||||||||||||||| AspSerAsnArgLeuPheSerLeuGluSerLeuLysAlaValTrpHisGlyHisLysLeu 50160 : GACTCCAATCGTCTCTTCAGCCTGGAGTCCCTCAAAGCTGTGTGGCACGGCCATAAACTG : 50219 81 : AspTyrMetLysSerAlaArgLeuGlyHisAlaAlaProAspSerGluValPheProLeu : 100 |||!:!:!:|||:!!|||!.!||||||||| !||||||:!!!!!|||||| !!! !||| AspPheLeuLysAlaAlaHisLeuGlyHisArgAlaProAsnThrGluValValGlnLeu 50220 : GACTTTCTTAAAGCAGCGCATCTCGGACACAGAGCGCCCAACACCGAAGTTGTTCAGCTG : 50279 101 : AlaGluProArgArgGlyArgValLeuGluPheAlaArgGlyHisArgProLeuValLeu : 120 ! !!!: !||||||.!!|||:!!|||!!:!:!|||!:!! ! ! |||||||||:!!||| GluAspLysArgArgSerArgIleLeuAspTyrAlaLysAspLysArgProLeuIleLeu 50280 : GAGGATAAAAGGCGCAGCCGTATCCTTGATTACGCAAAGGACAAGAGACCGCTTATCCTC : 50339 121 : SerPheGlySerCysSerUnkProProPheMetArgArgLeuLysAlaPheArgArgLeu : 140 !:!||||||||||||!!! |||||||||||| !|||||||||||||||!:! !!:!! AsnPheGlySerCysThr***ProProPheMetAlaArgLeuLysAlaPheGlnGlyVal 50340 : AACTTTGGCAGCTGCACCTGACCACCGTTCATGGCTCGTCTCAAGGCTTTCCAGGGAGTT : 50399 141 : ValLeuArgTyrAlaAspValAlaAspAlaLeuLeuIleTyrIleGluGluAlaHisPro : 160 |||! !!:! !!||||||:!!||||||:!!:!!:!!:!:||||||||||||||||||||| ValGlnGlnAsnAlaAspIleAlaAspSerValValValTyrIleGluGluAlaHisPro 50400 : GTGCAGCAGAATGCAGACATCGCAGACTCTGTAGTAGTGTACATCGAGGAAGCGCACCCT : 50459 161 : SerAspGlyTrpArgSerSerAspAlaProHisGlnIleArgArgHisArgSerLeuGlu : 180 ||||||||||||! !|||!!!|||||||||:!!||||||! !!:!||||||!!!|||||| SerAspGlyTrpMetSerThrAspAlaProTyrGlnIleProLysHisArgCysLeuGlu 50460 : TCCGACGGCTGGATGAGCACTGACGCGCCATATCAAATCCCCAAACACCGGTGTCTGGAG : 50519 181 : GluArgLeuSerAlaAlaArgLeuMetGluArgGluAlaProGlyCysAlaValValAla : 200 !!:||||||!:!||||||!:!||||||.!.! !|||!.!||||||||| !||||||||| AspArgLeuAsnAlaAlaGlnLeuMetHisLeuGluValProGlyCysLeuValValAla 50520 : GACAGGCTGAACGCTGCGCAGCTGATGCACCTGGAGGTACCCGGCTGCCTGGTCGTGGCG : 50579 201 : AspGlyMetGluAsnAlaAlaAsnSerAlaTyrGlyAlaTyrPheAspArgLeuTyrIle : 220 |||.!!|||||||||:!!:!!|||:!!||||||||||||||||||||||||||||||||| AspSerMetGluAsnSerSerAsnAlaAlaTyrGlyAlaTyrPheAspArgLeuTyrIle 50580 : GACAGCATGGAAAACTCCTCCAACGCTGCGTACGGAGCTTATTTCGACAGGCTTTATATA : 50639 221 : ValGlnAspGlyArgValValTyrGlnGlyGlyArgGlyProGluGlyPheArgIleSer : 240 :!!|||||||||!:!:!!||||||||||||||||||||||||||||||!:!||||||||| LeuGlnAspGlyLysIleValTyrGlnGlyGlyArgGlyProGluGlyTyrArgIleSer 50640 : CTGCAGGATGGAAAGATAGTTTACCAGGGTGGCAGAGGACCTGAGGGGTATCGGATCTCA : 50699 241 : GluLeuArgHisTrpLeuAspArgTyrArgGluArgLeu : 253 ||||||||| !!|||||||||!:!||||||:!!|||||| GluLeuArgAspTrpLeuAspGlnTyrArgGlnArgLeu 50700 : GAGCTTAGAGACTGGCTGGATCAATACAGACAAAGGCTG : 50740 vulgar: DIO3b 0 253 . KV884701.1:subseq(3433220,100000) 49999 50740 + 846 M 37 111 G 4 0 M 13 39 G 2 0 M 197 591 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-12-01 ##type DNA # # # seqname source feature start end score strand frame attributes # KV884701.1:subseq(3433220,100000) exonerate:protein2genome:local gene 50000 50740 846 + . gene_id 1 ; sequence DIO3b ; gene_orientation . KV884701.1:subseq(3433220,100000) exonerate:protein2genome:local cds 50000 50740 . + . KV884701.1:subseq(3433220,100000) exonerate:protein2genome:local exon 50000 50740 . + . insertions 0 ; deletions 6 KV884701.1:subseq(3433220,100000) exonerate:protein2genome:local similarity 50000 50740 846 + . alignment_id 1 ; Query DIO3b ; Align 50000 1 111 ; Align 50111 42 39 ; Align 50150 57 591 # --- END OF GFF DUMP --- # -- completed exonerate analysis