Command line: [/mnt/NFS_UPF/soft/exonerate/exonerate-2.2.0-x86_64/bin/exonerate -m p2g --showtargetgff --exhaustive yes -q query3 -t q_subseq] Hostname: [am33111-116073u] C4 Alignment: ------------ Query: SPP00000100_1.0 # Protein # Glutathione peroxidase 3 (GPx3) # Mus musculus # Complete Target: PVKD010003256.1:subseq(35591,123257) Cricetomys gambianus isolate US082 CriGam_scaffold_6508, whole genome shotgun sequence Model: protein2genome:local Raw score: 1104 Query range: 0 -> 223 Target range: 93060 -> 100227 1 : MetAlaArgIleLeuArgAlaSerCysLeuLeuSerLeuLeuLeuAlaGlyPheValP : 20 |||||||||:!!|||||||||||||||||||||||||||||||||||||||||||||| MetAlaArgLeuLeuArgAlaSerCysLeuLeuSerLeuLeuLeuAlaGlyPheValP 93061 : ATGGCCCGGCTCCTCCGGGCATCCTGCCTCCTGTCCCTGCTCCTGGCCGGGTTCGTTC : 93118 21 : roProGlyArgGlyGlnGluLysSerLys >>>> Target Intron 1 >>>> : 30 ||||||||||||||||||||||||||||| 4452 bp roProGlyArgGlyGlnGluLysSerLys++ ++ 93119 : CGCCCGGCCGGGGGCAAGAGAAGTCCAAGgt.........................ag : 97598 31 : ThrAspCysHisGlyGlyMetSerGlyThrIleTyrGluTyrGlyAlaLeuThrIleA : 49 .!!||||||||||||||||||||||||||||||||||||||||||||||||||||||| AlaAspCysHisGlyGlyMetSerGlyThrIleTyrGluTyrGlyAlaLeuThrIleA 97599 : GCAGACTGCCACGGAGGCATGAGTGGTACCATCTATGAGTACGGGGCCCTCACCATTG : 97657 50 : spGlyGluGluTyrIleProPheLysGlnTyrAlaGlyLysTyrIleLeuPheValAs : 68 ||||||||||||||||||||||||||||||||||||||||||||:!!||||||||||| spGlyGluGluTyrIleProPheLysGlnTyrAlaGlyLysTyrValLeuPheValAs 97658 : ACGGAGAGGAATACATTCCTTTTAAGCAGTATGCAGGCAAATATGTCCTCTTTGTCAA : 97714 69 : nValAlaSerTyr<-><-><->GlyAspGlnTyrLeu{G} >>>> Target Intr : 78 ||||||||||||| !||||||||||||{|} 1345 bp nValAlaSerTyr***GlyLeuThrAspGlnTyrLeu{G}++ 97715 : CGTAGCTAGCTACTGAGGTCTGACAGACCAATACCTT{G}gt................ : 97756 79 : on 2 >>>> {lu}LeuAsnAlaLeuGlnGluGluLeuGlyProPheGlyLeuValI : 93 {||}||||||||||||||||||||||||||||||||||||||||||| ++{lu}LeuAsnAlaLeuGlnGluGluLeuGlyProPheGlyLeuValI 97757 : .........ag{AA}CTGAATGCACTACAAGAAGAACTTGGGCCATTCGGCCTGGTCA : 99143 94 : leLeuGlyPheProSerAsnGlnPheGlyLysGlnGluProGlyGluAsnSerGluIl : 112 ||||||||||||||!!!||||||||||||||||||||||||||||||||||||||||| leLeuGlyPheProCysAsnGlnPheGlyLysGlnGluProGlyGluAsnSerGluIl 99144 : TCCTGGGCTTCCCCTGCAACCAATTTGGCAAACAGGAGCCAGGCGAGAACTCGGAGAT : 99200 113 : eLeuProSerLeu{Ly} >>>> Target Intron 3 >>>> {s}TyrValArg : 120 |||||||||||||{||} 362 bp {|}||||||||| eLeuProSerLeu{Ly}++ ++{s}TyrValArg 99201 : CCTACCCAGTCTC{AA}gt.........................ag{G}TATGTCCGA : 99586 121 : ProGlyGlyGlyPheValProAsnPheGlnLeuPheGluLysGlyAspValAsnGlyG : 140 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||| ProGlyGlyGlyPheValProAsnPheGlnLeuPheGluLysGlyAspValAsnGlyG 99587 : CCAGGTGGGGGTTTTGTCCCTAACTTCCAGCTCTTTGAGAAAGGAGACGTGAATGGGG : 99646 141 : luLysGluGlnLysPheTyrThrPheLeuLys >>>> Target Intron 4 >>> : 151 |||||||||||||||||||||||||||||||| 330 bp luLysGluGlnLysPheTyrThrPheLeuLys++ 99647 : AGAAAGAGCAGAAGTTCTACACTTTCTTGAAGgt........................ : 99681 152 : > AsnSerCysProProThrAlaGluLeuLeuGlySerProGlyArgLeuPheTrpG : 169 |||||||||||||||||||||||||||||||||||||||! !||||||||||||| ++AsnSerCysProProThrAlaGluLeuLeuGlySerProAspArgLeuPheTrpG 99682 : .agAACTCCTGCCCTCCCACTGCGGAGCTCCTGGGCTCACCTGACCGCCTCTTTTGGG : 100063 170 : luProMetLysIleHisAspIleArgTrpAsnPheGluLysPheLeuValGlyProAs : 188 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||| luProMetLysIleHisAspIleArgTrpAsnPheGluLysPheLeuValGlyProAs 100064 : AACCCATGAAGATCCATGACATCCGCTGGAACTTTGAGAAGTTTCTGGTGGGGCCAGA : 100120 189 : pGlyIleProValMetArgTrpTyrHisArgThrThrValSerAsnValLysMetAsp : 207 ||||||||||:!!||||||||||||||||||||||||||||||||||||||||||||| pGlyIleProIleMetArgTrpTyrHisArgThrThrValSerAsnValLysMetAsp 100121 : TGGCATACCCATTATGCGCTGGTACCACCGGACCACAGTCAGCAACGTCAAGATGGAC : 100177 208 : IleLeuSerTyrMetArgArgGlnAlaAlaLeuSerAlaArgGlyLys : 223 |||||||||||||||||||||||||||||||||.!!|||||||||||| IleLeuSerTyrMetArgArgGlnAlaAlaLeuGlyAlaArgGlyLys 100178 : ATCCTGTCCTACATGAGGCGGCAGGCTGCCCTGGGGGCCAGAGGAAAG : 100227 vulgar: SPP00000100_1.0 0 223 . PVKD010003256.1:subseq(35591,123257) 93060 100227 + 1104 M 29 87 5 0 2 I 0 4448 3 0 2 M 43 129 G 0 9 M 5 15 S 0 1 5 0 2 I 0 1341 3 0 2 S 1 2 M 38 114 S 0 2 5 0 2 I 0 358 3 0 2 S 1 1 M 33 99 5 0 2 I 0 326 3 0 2 M 73 219 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2019-11-20 ##type DNA # # # seqname source feature start end score strand frame attributes # PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local gene 93061 100227 1104 + . gene_id 1 ; sequence SPP00000100_1.0 ; gene_orientation + PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local cds 93061 93147 . + . PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local exon 93061 93147 . + . insertions 0 ; deletions 0 PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local splice5 93148 93149 . + . intron_id 1 ; splice_site "GT" PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local intron 93148 97599 . + . intron_id 1 PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local splice3 97598 97599 . + . intron_id 0 ; splice_site "AG" PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local cds 97600 97753 . + . PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local exon 97600 97753 . + . insertions 9 ; deletions 0 PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local splice5 97754 97755 . + . intron_id 2 ; splice_site "GT" PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local intron 97754 99098 . + . intron_id 2 PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local splice3 99097 99098 . + . intron_id 1 ; splice_site "AG" PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local cds 99099 99216 . + . PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local exon 99099 99216 . + . insertions 0 ; deletions 0 PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local splice5 99217 99218 . + . intron_id 3 ; splice_site "GT" PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local intron 99217 99578 . + . intron_id 3 PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local splice3 99577 99578 . + . intron_id 2 ; splice_site "AG" PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local cds 99579 99678 . + . PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local exon 99579 99678 . + . insertions 0 ; deletions 0 PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local splice5 99679 99680 . + . intron_id 4 ; splice_site "GT" PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local intron 99679 100008 . + . intron_id 4 PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local splice3 100007 100008 . + . intron_id 3 ; splice_site "AG" PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local cds 100009 100227 . + . PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local exon 100009 100227 . + . insertions 0 ; deletions 0 PVKD010003256.1:subseq(35591,123257) exonerate:protein2genome:local similarity 93061 100227 1104 + . alignment_id 1 ; Query SPP00000100_1.0 ; Align 93061 1 87 ; Align 97600 30 129 ; Align 97738 73 15 ; Align 99101 79 114 ; Align 99580 118 99 ; Align 100009 151 219 # --- END OF GFF DUMP --- # -- completed exonerate analysis