Command line: [exonerate -m p2g --showtargetgff -q MSRB2.fa -t MSRB2.subseq] Hostname: [ws101069u] C4 Alignment: ------------ Query: SPP00000091_2.0 # Protein # Methionine-R-sufoxide reductase 2 (MSRB2) # Human Target: JSZB01014487.1:subseq(130000,120000) Manis javanica isolate MP_PG03-UM SCAFFOLD17423, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 648 Query range: 35 -> 182 Target range: 62373 -> 50873 36 : LeuGlyGluAlaGlySerLeuAlaThrCysGluLeuProLeuAlaLysSerGluTrpGln : 55 !!!..!.!..!!||||||! !.!!! ! !|||! ! !!|||.!!|||!!!||||||||| PheAsnAsnThrGlySerProThrLysHisGluProSerLeuThrLysArgGluTrpGln 62373 : TTTAACAATACAGGGTCCCCCACAAAGCATGAGCCTTCTCTCACCAAGAGAGAGTGGCAA : 62316 56 : LysLysLeuThrProGluGlnPheTyrValThrArgGluLysGlyThrGluPro >>>> : 74 |||:!!|||||| !!|||||||||||||||||||||! !||||||||||||||| LysGluLeuThrAlaGluGlnPheTyrValThrArgGlyLysGlyThrGluPro++ 62315 : AAGGAGCTGACTGCAGAACAATTCTATGTCACAAGAGGAAAGGGAACAGAACCGgt.... : 62257 75 : Target Intron 1 >>>> ProPheSerGlyIleTyrLeuAsnAsnLysGluAlaG : 86 4875 bp |||||||||||||||||||||||||||! !|||:!!| ++ProPheSerGlyIleTyrLeuAsnAsnThrGluSerG 62256 : .....................agCCTTTCAGTGGGATCTACCTGAATAACACGGAATCAG : 57348 87 : lyMetTyrHisCysValCysCysAspSerProLeuPhe{Se} >>>> Target Intr : 99 ||||||||||||||||||||||||||||||||||||||{||} 5236 bp lyMetTyrHisCysValCysCysAspSerProLeuPhe{Se}++ 57347 : GAATGTATCACTGTGTGTGTTGTGACAGCCCACTCTTC{AG}gt................ : 57305 100 : on 2 >>>> {r}SerGluLysLysTyrCysSerGlyThrGlyTrpProSerPheSerG : 115 {|}|||||||||||||||||||||||||||||||||||||||||||||| ++{r}SerGluLysLysTyrCysSerGlyThrGlyTrpProSerPheSerG 57304 : .........ag{T}TCTGAGAAAAAGTACTGCTCTGGCACTGGATGGCCTTCATTTTCTG : 52025 116 : luAlaHisGlyThrSerGlySerAspGluSerHisThrGlyIleLeuArgArgLeuAspT : 135 ||||||||||||||!!!|||||||||||||||:!!|||||||||||||||!.!:!!|||| luAlaHisGlyThrCysGlySerAspGluSerAsnThrGlyIleLeuArgHisValAspT 52024 : AGGCTCATGGTACGTGTGGCTCTGATGAAAGTAACACAGGGATCCTGAGACACGTGGACA : 51965 136 : hrSerLeuGlySerAlaArgThrGluValValCysLysGln >>>> Target Intro : 149 |||||||||||!!!!.!|||||||||||||||||||||||| 948 bp hrSerLeuGlyTrpValArgThrGluValValCysLysGln++ 51964 : CCTCGCTGGGATGGGTGCGCACAGAGGTTGTCTGCAAGCaggt................. : 51921 150 : n 3 >>>> CysGluAlaHisLeuGlyHisValPheProAspGlyProGlyProAsnGl : 165 ||||||||||||||||||||||||||||||||||||||||||||| !|| ++CysGluAlaHisLeuGlyHisValPheProAspGlyProGlyProAlaGl 51920 : ........agTGTGAAGCTCATCTTGGCCATGTGTTTCCTGATGGACCTGGGCCTGCTGG : 50927 166 : yGlnArgPheCysIleAsnSerValAlaLeuLysPheLysProArgLysHis : 182 |||||||||||||||||||||||||||||||||||||||||||||| !||| yGlnArgPheCysIleAsnSerValAlaLeuLysPheLysProArgGlyHis 50926 : TCAGAGGTTTTGCATCAACAGTGTGGCCCTGAAATTCAAGCCAAGGGGACAC : 50874 vulgar: SPP00000091_2.0 35 182 . JSZB01014487.1:subseq(130000,120000) 62373 50873 - 648 M 38 114 5 0 2 I 0 4871 3 0 2 M 25 75 S 0 2 5 0 2 I 0 5232 3 0 2 S 1 1 M 49 147 5 0 2 I 0 944 3 0 2 M 34 102 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-14 ##type DNA # # # seqname source feature start end score strand frame attributes # JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local gene 50874 62373 648 - . gene_id 1 ; sequence SPP00000091_2.0 ; gene_orientation + JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local cds 62260 62373 . - . JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local exon 62260 62373 . - . insertions 0 ; deletions 0 JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local splice5 62258 62259 . - . intron_id 1 ; splice_site "GT" JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local intron 57385 62259 . - . intron_id 1 JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local splice3 57385 57386 . - . intron_id 0 ; splice_site "AG" JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local cds 57308 57384 . - . JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local exon 57308 57384 . - . insertions 0 ; deletions 0 JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local splice5 57306 57307 . - . intron_id 2 ; splice_site "GT" JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local intron 52072 57307 . - . intron_id 2 JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local splice3 52072 52073 . - . intron_id 1 ; splice_site "AG" JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local cds 51924 52071 . - . JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local exon 51924 52071 . - . insertions 0 ; deletions 0 JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local splice5 51922 51923 . - . intron_id 3 ; splice_site "gt" JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local intron 50976 51923 . - . intron_id 3 JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local splice3 50976 50977 . - . intron_id 2 ; splice_site "AG" JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local cds 50874 50975 . - . JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local exon 50874 50975 . - . insertions 0 ; deletions 0 JSZB01014487.1:subseq(130000,120000) exonerate:protein2genome:local similarity 50874 62373 648 - . alignment_id 1 ; Query SPP00000091_2.0 ; Align 62374 36 114 ; Align 57385 74 75 ; Align 52071 100 147 ; Align 50976 149 102 # --- END OF GFF DUMP --- # -- completed exonerate analysis