Command line: [exonerate --exhaustive yes -m p2g --showtargetgff -q treball/run/MSRB2/MSRB2.x.fa -t treball/run/MSRB2/MSRB2.extraction.KV884788.1.fa] Hostname: [sitdoc] C4 Alignment: ------------ Query: MSRB2 Protein # Methionine-R-sufoxide reductase 2 (MSRB2) # Zebrafish Target: KV884788.1:subseq(0,100000) Monopterus albus unplaced genomic scaffold scaffold100.1, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 697 Query range: 0 -> 179 Target range: 37817 -> 35897 1 : MetSerArgPheLeuValArgLeuSerThrValValSerLysGlyAlaThrGlyLysSer : 20 ||||||||||||:!!!.!||||||! !! !:!!!.!|||:!! !!||||||!.!!:!||| MetSerArgPheValAlaArgLeuPheIleLeuAlaSerGlnArgAlaThrAlaArgSer 37817 : atgtCACGATTTGTCGCGCGTTTGTTCATTCTCGCTTCTCAGCGCGCCACAGCCAGATCC : 37760 21 : ValLeuProGlnLysArgIlePheAlaGlyIleArgLeuIleSerSerSerThr{G} > : 39 |||:!! ! !!:!|||:!!||| ||||||! !:!!|||!!!!!! !{|} ValValLeuLeuArgArgValPhe------IleArgProLeuSerThrThrGln{G}++ 37759 : GTGGTGTTGTTGAGAAGGGTATTC------ATCCGTCCTCTATCCACAACTCAA{G}gt. : 37706 40 : >>> Target Intron 1 >>>> {ly}LeuGlnSerLeuThrArgTyrAspGlnSer : 49 519 bp {||}|||||||||||||||||||||..!:!!!:! ++{ly}LeuGlnSerLeuThrArgTyrSerGluAsn 37705 : ........................ag{GC}CTGCAGTCCCTCACACGATACAGTGAGAAT : 37160 50 : ThrAspTrpGlnArgLysLeuSerProGluGlnTyrValValThrArgGluLysGlyThr : 69 !|||||||||!:!||||||!!!|||||||||||||||||||||||||||:::|||||| GluAspTrpGlnLysLysLeuThrProGluGlnTyrValValThrArgGluArgGlyThr 37159 : GAAGACTGGCAAAAGAAACTCACCCCAGAACAGTATGTAGTcaccagagagagaggaact : 37100 70 : GluVal >>>> Target Intron 2 >>>> ProPheSerGlyIleTyrLeuAsnH : 80 |||||| 384 bp ||||||||||||||||||||||||| GluVal++ ++ProPheSerGlyIleTyrLeuAsnH 37099 : gaGGTGgt.........................agCCTTTTAGTGGGATCTACCTTAACC : 36683 81 : isAsnGluValGlyMetTyrHisCysValCysCysAspSerProLeuPhe{Se} >>>> : 97 || !!|||||||||||||||||||||||||||||||||:!!|||||||||{||} isTyrGluValGlyMetTyrHisCysValCysCysAspAlaProLeuPhe{Se}++ 36682 : ATTATGAAGTGGGGATGTACCACTGTGTCTGCTGTGATGCCCCACTCTTC{AG}gt.... : 36628 98 : Target Intron 3 >>>> {r}SerGluAlaLysTyrAspAlaGlyAlaGlyTrpP : 109 235 bp {|}||||||||||||||||||:!!|||.!!||||||| ++{r}SerGluAlaLysTyrAspSerGlyThrGlyTrpP 36627 : .....................ag{T}TCAGAGGCTAAATATGACTCTGGGACAGGCTGGC : 36361 110 : roSerPheCysGluAlaHisGlyThrTrpGluLysAspGluSerHisAlaAsnIleValA : 129 ||:!!||| |||||||||||||||||||||!:!||||||||||||.!!::!|||:!!| roAlaPheLysGluAlaHisGlyThrTrpGluArgAspGluSerHisThrSerIleIleA 36360 : CAGCATTCAAAGAGGCTCATGGGACATGGGAGAGGGATGAAAGCCACACTTCCATCATTC : 36301 130 : rgArgProAspAsnSerLeuGlySerThrGlyThrGluValIleCysLysHis >>>> : 147 ||||||||||||||||||||||||||.!!||||||||||||:!!||||||:!! rgArgProAspAsnSerLeuGlySerAlaGlyThrGluValLeuCysLysAsn++ 36300 : GTCGCCCTGACAACAGCCTGGGTAGTGCTGGGACAGAAGTCCTTTGTAAAAATgt..... : 36245 148 : Target Intron 4 >>>> CysAspAlaHisLeuGlyHisValPheGluAspGlyPr : 159 251 bp |||||||||||||||||||||||||||!!:|||||||| ++CysAspAlaHisLeuGlyHisValPheAspAspGlyPr 36244 : ....................agtgTGATGCTCACCTGGGTCATGTGTTTGATGATGGACC : 35960 160 : oAspProThrGlyGlnArgPheCysIleAsnSerValAlaLeuAsnPheLysProArgGl : 179 |||||||||||||||||||||||||||||||||||||||||||!!.|||||||||||||| oAspProThrGlyGlnArgPheCysIleAsnSerValAlaLeuLysPheLysProArgGl 35959 : AGACCCAACGGGTCAGCGATTCTGCATCAACAGTGTGGCACTCAAGTTTAAACCAAGAGA : 35900 180 : u : 179 | u 35899 : A : 35898 vulgar: MSRB2 0 179 . KV884788.1:subseq(0,100000) 37817 35897 - 697 M 28 84 G 2 0 M 8 24 S 0 1 5 0 2 I 0 515 3 0 2 S 1 2 M 32 96 5 0 2 I 0 380 3 0 2 M 25 75 S 0 2 5 0 2 I 0 231 3 0 2 S 1 1 M 49 147 5 0 2 I 0 247 3 0 2 M 33 99 # --- 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 # KV884788.1:subseq(0,100000) exonerate:protein2genome:local gene 35898 37817 697 - . gene_id 1 ; sequence MSRB2 ; gene_orientation + KV884788.1:subseq(0,100000) exonerate:protein2genome:local cds 37709 37817 . - . KV884788.1:subseq(0,100000) exonerate:protein2genome:local exon 37709 37817 . - . insertions 0 ; deletions 2 KV884788.1:subseq(0,100000) exonerate:protein2genome:local splice5 37707 37708 . - . intron_id 1 ; splice_site "GT" KV884788.1:subseq(0,100000) exonerate:protein2genome:local intron 37190 37708 . - . intron_id 1 KV884788.1:subseq(0,100000) exonerate:protein2genome:local splice3 37190 37191 . - . intron_id 0 ; splice_site "AG" KV884788.1:subseq(0,100000) exonerate:protein2genome:local cds 37092 37189 . - . KV884788.1:subseq(0,100000) exonerate:protein2genome:local exon 37092 37189 . - . insertions 0 ; deletions 0 KV884788.1:subseq(0,100000) exonerate:protein2genome:local splice5 37090 37091 . - . intron_id 2 ; splice_site "GT" KV884788.1:subseq(0,100000) exonerate:protein2genome:local intron 36708 37091 . - . intron_id 2 KV884788.1:subseq(0,100000) exonerate:protein2genome:local splice3 36708 36709 . - . intron_id 1 ; splice_site "AG" KV884788.1:subseq(0,100000) exonerate:protein2genome:local cds 36631 36707 . - . KV884788.1:subseq(0,100000) exonerate:protein2genome:local exon 36631 36707 . - . insertions 0 ; deletions 0 KV884788.1:subseq(0,100000) exonerate:protein2genome:local splice5 36629 36630 . - . intron_id 3 ; splice_site "GT" KV884788.1:subseq(0,100000) exonerate:protein2genome:local intron 36396 36630 . - . intron_id 3 KV884788.1:subseq(0,100000) exonerate:protein2genome:local splice3 36396 36397 . - . intron_id 2 ; splice_site "AG" KV884788.1:subseq(0,100000) exonerate:protein2genome:local cds 36248 36395 . - . KV884788.1:subseq(0,100000) exonerate:protein2genome:local exon 36248 36395 . - . insertions 0 ; deletions 0 KV884788.1:subseq(0,100000) exonerate:protein2genome:local splice5 36246 36247 . - . intron_id 4 ; splice_site "Gt" KV884788.1:subseq(0,100000) exonerate:protein2genome:local intron 35997 36247 . - . intron_id 4 KV884788.1:subseq(0,100000) exonerate:protein2genome:local splice3 35997 35998 . - . intron_id 3 ; splice_site "ag" KV884788.1:subseq(0,100000) exonerate:protein2genome:local cds 35898 35996 . - . KV884788.1:subseq(0,100000) exonerate:protein2genome:local exon 35898 35996 . - . insertions 0 ; deletions 0 KV884788.1:subseq(0,100000) exonerate:protein2genome:local similarity 35898 37817 697 - . alignment_id 1 ; Query MSRB2 ; Align 37818 1 84 ; Align 37734 31 24 ; Align 37188 40 96 ; Align 36708 72 75 ; Align 36395 98 147 ; Align 35997 147 99 # --- END OF GFF DUMP --- # -- completed exonerate analysis