Command line: [exonerate --exhaustive yes -m p2g --showtargetgff -q /export/home/u114420/treball/fastas/MsrA2.x.fa -t MsrA2.extraction.KV884890.1.fa] Hostname: [sitdoc] C4 Alignment: ------------ Query: MsrA2 Protein # Methionine sulfoxide reductase A (MsrA) # Zebrafish Target: KV884890.1:subseq(802559,100000) Monopterus albus unplaced genomic scaffold scaffold202.1, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 961 Query range: 0 -> 232 Target range: 82460 -> 49838 1 : MetValAlaArgThrSerValArgLeuIleTrpArgGlnPheIleGlnSerArgMetGly : 20 ||||||..!! !:!!!!!|||||||||:!:||||||!::.!!:!!:!!|||||||||||| MetValCysLeuSerArgValArgLeuValTrpArgArgLeuValLysSerArgMetGly 82460 : ATGGTCTGCCTCTCCAGAGTCCGTCTAGTGTGGCGTCGTCTTGTCAAAAGCAGGATGGGA : 82403 21 : GluMetSerSerLysValGlnMetIleSerProGluGluAlaLeuProGlyArgGluGln : 40 !!:|||:!!||||||||||||||| !!!|||||| !||||||||||||||| .!. AspMetAlaSerLysValGlnMetProThrProGluThrAlaLeuProGlyArgThrAsp 82402 : GACATGGCCTCTAAAGTCCAGATGCCGACTCCAGAGACAGCGCTGCCCGGCCGCACTGAC : 82343 41 : SerIleLysValSer{A} >>>> Target Intron 1 >>>> {la}LysHisAsp : 49 |||||||||||||||{|} 18582 bp {||}|||||| !! SerIleLysValSer{A}++ ++{la}LysHisHis 82342 : AGCATCAAAGTCTCG{G}gt.........................ag{CC}AAACATCAT : 63734 50 : ValAsnGlyAsnArgThrValProProPheProGluGlyLeuGlnMetValLeuPhe{G} : 69 ||||||||||||||||||||||||||||||||||||||| !:!!||||||:!!|||{|} ValAsnGlyAsnArgThrValProProPheProGluGlyThrGluMetValValPhe{G} 63733 : GTGAATGGCAACAGGACTGTTCCACCTTTTCCAGAGGGGACAGAGATGGTCGTGTTT{G} : 63674 70 : >>>> Target Intron 2 >>>> {ly}MetGlyCysPheTrpGlyAlaGluArg : 78 249 bp {||}||||||||||||||||||||||||||| ++ ++{ly}MetGlyCysPheTrpGlyAlaGluArg 63673 : gt.........................ag{GG}ATGGGCTGTTTCTGGGGAGCGGAGAGG : 63398 79 : LysPheTrpArgGlnLysGlyValTyrSerThrGlnValGlyTyrSerGlyGlyTyrThr : 98 ||||||||||||||||||||||||||||||||||||||||||||||||||||||:!!||| LysPheTrpArgGlnLysGlyValTyrSerThrGlnValGlyTyrSerGlyGlyHisThr 63397 : AAGTTTTGGAGACAAAAAGGAGTTTATTCCACCCAAGTAGGCTACTCTGGAGGACACACA : 63338 99 : ProAsnProThrTyrGluGluValCysThr{G} >>>> Target Intron 3 >>>> : 109 |||||||||:!!|||:!!||||||||||||{|} 4168 bp ProAsnProSerTyrLysGluValCysThr{G}++ 63337 : CCTAACCCCTCATACAAGGAAGTCTGCACA{G}gt......................... : 63302 110 : {ly}LysThrGlyHisThrGluValValArgValValPheGluProGlnLysIleLys : 127 {||}!:!|||||||||.!!|||||||||||||||||||||.!.|||.!.:!!||| ++{ly}ArgThrGlyHisAlaGluValValArgValValPheHisProAspGlnIleCys 63301 : ag{GT}AGGACAGGCCACGCTGAGGTGGTGAGAGTTGTTTTCCATCCCGACCAGATATGC : 59083 128 : PheSerGluLeuLeuLysValPheTrpGluSerHisAsnProThrGln{G} >>>> Ta : 144 |||:!!...|||||||||||||||||||||||||||:!!|||||||||{|} PheAlaSerLeuLeuLysValPheTrpGluSerHisAspProThrGln{G}++ 59082 : TTCGCAAGCCTCCTCAAAGTTTTCTGGGAAAGCCATGACCCGACTCAA{G}gt....... : 59029 145 : rget Intron 4 >>>> {ly}MetArgGlnGlyAsnAspValGlyThrThrTyrArg : 156 694 bp {||}||||||||||||||||||:!!||||||..!|||||| ++{ly}MetArgGlnGlyAsnAspIleGlyThrValTyrArg 59028 : ..................ag{ga}aTGCGGCAGGGGAACGATATTGGGACAGTATACCGC : 58302 157 : SerSerIleTyrThrAsnThrGlnGluGlnLeuGluGlnAlaLeuGlnSerArgGluGlu : 176 |||:!!||||||||| !!|||:!!:!!|||||||||:!!|||||| |||!:!!!::!! SerAlaIleTyrThrTyrThrLysGlnGlnLeuGluGluAlaLeuAlaSerLysAspGln 58301 : TCTGCCATCTACACCTACACGAAGCAGCAGCTCGAGGAAGCTTTGGCCTCCAAAGATCAA : 58242 177 : TyrGlnLys >>>> Target Intron 5 >>>> ValLeuThrGluGluGlyPheG : 187 ||||||||| 8233 bp |||||||||||||||||||||| TyrGlnLys++ ++ValLeuThrGluGluGlyPheG 58241 : TACCAGAAGgt.........................agGTACTAACAGAAGAAGGTTTTG : 49976 188 : lyAlaIleThrThrGluIleThrMetAlaLysGluPheTyrTyrAlaGluAspTyrHisG : 207 || !!:!!||||||||||||:!! |||!!. !||||||||||||||||||||||||| lyProValThrThrGluIleSerGluAlaAsnProPheTyrTyrAlaGluAspTyrHisG 49975 : GTCCAGTTACAACAGAAATATCTGAAGCCAATCCATTCTACTATGCTGAAGACTACCACC : 49916 208 : lnGlnTyrLeuSerLysAsnProAspGlyTyrCysGlyLeuGlyGlyThrGlyValSerC : 227 |||||||||||||||||:!!|||||||||||||||||||||||||||||||||||||||| lnGlnTyrLeuSerLysHisProAspGlyTyrCysGlyLeuGlyGlyThrGlyValSerC 49915 : AGCAGTACCTGAGCAAACACCCTGATGGATACTGTGGCCTGGGAGGGACAGGAGTCTCTT : 49856 228 : ysProIleGlyLeuLys : 232 |||||||||||:!!||| ysProIleGlyIleLys 49855 : GTCCAATAGGAATCAAG : 49839 vulgar: MsrA2 0 232 . KV884890.1:subseq(802559,100000) 82460 49838 - 961 M 45 135 S 0 1 5 0 2 I 0 18578 3 0 2 S 1 2 M 22 66 S 0 1 5 0 2 I 0 245 3 0 2 S 1 2 M 39 117 S 0 1 5 0 2 I 0 4164 3 0 2 S 1 2 M 34 102 S 0 1 5 0 2 I 0 690 3 0 2 S 1 2 M 35 105 5 0 2 I 0 8229 3 0 2 M 53 159 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-23 ##type DNA # # # seqname source feature start end score strand frame attributes # KV884890.1:subseq(802559,100000) exonerate:protein2genome:local gene 49839 82460 961 - . gene_id 1 ; sequence MsrA2 ; gene_orientation + KV884890.1:subseq(802559,100000) exonerate:protein2genome:local cds 82325 82460 . - . KV884890.1:subseq(802559,100000) exonerate:protein2genome:local exon 82325 82460 . - . insertions 0 ; deletions 0 KV884890.1:subseq(802559,100000) exonerate:protein2genome:local splice5 82323 82324 . - . intron_id 1 ; splice_site "GT" KV884890.1:subseq(802559,100000) exonerate:protein2genome:local intron 63743 82324 . - . intron_id 1 KV884890.1:subseq(802559,100000) exonerate:protein2genome:local splice3 63743 63744 . - . intron_id 0 ; splice_site "AG" KV884890.1:subseq(802559,100000) exonerate:protein2genome:local cds 63674 63742 . - . KV884890.1:subseq(802559,100000) exonerate:protein2genome:local exon 63674 63742 . - . insertions 0 ; deletions 0 KV884890.1:subseq(802559,100000) exonerate:protein2genome:local splice5 63672 63673 . - . intron_id 2 ; splice_site "GT" KV884890.1:subseq(802559,100000) exonerate:protein2genome:local intron 63425 63673 . - . intron_id 2 KV884890.1:subseq(802559,100000) exonerate:protein2genome:local splice3 63425 63426 . - . intron_id 1 ; splice_site "AG" KV884890.1:subseq(802559,100000) exonerate:protein2genome:local cds 63305 63424 . - . KV884890.1:subseq(802559,100000) exonerate:protein2genome:local exon 63305 63424 . - . insertions 0 ; deletions 0 KV884890.1:subseq(802559,100000) exonerate:protein2genome:local splice5 63303 63304 . - . intron_id 3 ; splice_site "GT" KV884890.1:subseq(802559,100000) exonerate:protein2genome:local intron 59137 63304 . - . intron_id 3 KV884890.1:subseq(802559,100000) exonerate:protein2genome:local splice3 59137 59138 . - . intron_id 2 ; splice_site "AG" KV884890.1:subseq(802559,100000) exonerate:protein2genome:local cds 59032 59136 . - . KV884890.1:subseq(802559,100000) exonerate:protein2genome:local exon 59032 59136 . - . insertions 0 ; deletions 0 KV884890.1:subseq(802559,100000) exonerate:protein2genome:local splice5 59030 59031 . - . intron_id 4 ; splice_site "GT" KV884890.1:subseq(802559,100000) exonerate:protein2genome:local intron 58338 59031 . - . intron_id 4 KV884890.1:subseq(802559,100000) exonerate:protein2genome:local splice3 58338 58339 . - . intron_id 3 ; splice_site "ag" KV884890.1:subseq(802559,100000) exonerate:protein2genome:local cds 58231 58337 . - . KV884890.1:subseq(802559,100000) exonerate:protein2genome:local exon 58231 58337 . - . insertions 0 ; deletions 0 KV884890.1:subseq(802559,100000) exonerate:protein2genome:local splice5 58229 58230 . - . intron_id 5 ; splice_site "GT" KV884890.1:subseq(802559,100000) exonerate:protein2genome:local intron 49998 58230 . - . intron_id 5 KV884890.1:subseq(802559,100000) exonerate:protein2genome:local splice3 49998 49999 . - . intron_id 4 ; splice_site "AG" KV884890.1:subseq(802559,100000) exonerate:protein2genome:local cds 49839 49997 . - . KV884890.1:subseq(802559,100000) exonerate:protein2genome:local exon 49839 49997 . - . insertions 0 ; deletions 0 KV884890.1:subseq(802559,100000) exonerate:protein2genome:local similarity 49839 82460 961 - . alignment_id 1 ; Query MsrA2 ; Align 82461 1 135 ; Align 63741 47 66 ; Align 63423 70 117 ; Align 59135 110 102 ; Align 58336 145 105 ; Align 49998 180 159 # --- END OF GFF DUMP --- # -- completed exonerate analysis