Command line: [exonerate -m p2g --showtargetgff -q SPP00000009_1.0.fa -t SPP00000009_1.0.GL397312.1/GL397312.1.subseq] Hostname: [upf-OptiPlex-755] C4 Alignment: ------------ Query: SPP00000009_1.0 # Protein # Glutathione peroxidase 6 (GPx6) # Homo sapiens # Complete Target: GL397312.1:subseq(15866587,103211) dna:supercontig supercontig:Nleu1.0:GL397312.1:1:16802348:1 REF:[revcomp] Model: protein2genome:local Raw score: 302 Query range: 35 -> 216 Target range: 53211 -> 49999 36 : ValThrGlyThrIleTyrGluTyrGlyAlaLeuThrLeuAsnGlyGluGluTyrIleGln : 55 :!!.!! !:!!.!!|||!!: !.!!|||:!!!:!|||:!!||||||:!! :!!.!. IleAlaLysSerPheTyrAspLeuSerAlaIleSerLeuAspGlyGluLys---ValAsp 53211 : ATTGCCAAGTCCTTCTATGACCTCAGTGCAATCAGCCTGGATGGGGAGAAG---GTAGAT : 53157 56 : PheLysGlnPheAlaGlyLysHisValLeuPheValAsnValAlaAlaTyrUnkGlyLeu : 75 |||!!. !||| !|||!:! !||||||.!!! !|||||||||:!! ! ||| ! PheAsnThrPheArgGlyArgAlaValLeuIleGluAsnValAlaSerLeu***GlyThr 53156 : TTCAATACGTTCCGGGGCAGGGCCGTGCTGATTGAGAATGTGGCTTCACTCTGAGGCACA : 53097 76 : AlaAlaGln<->TyrProGluLeuAsnAlaLeuGlnGluGluLeuLysAsnPheGlyVal : 94 .!!.!!!:! !:! !!:!!||||||! !|||||| ...!!! !.. !!:!! ThrThrArgAspPheThrGlnLeuAsnGluLeuGlnCysArgPheProArg---ArgLeu 53096 : ACCACCCGGGACTTCACCCAGCTCAACGAGCTGCAATGCCGCTTTCCCAGG---CGCCTG : 53040 95 : IleValLeuAlaPheProCysAsnGlnPheGlyLysGln >>>> Target Intron : 108 :!:||||||!.!||||||||||||||||||||| ! ||| 2660 bp ValValLeuGlyPheProCysAsnGlnPheGlyHisGln++ 53039 : GTGGTCCTTGGCTTCCCCTGCAACCAATTTGGACATCAGgt................... : 52996 109 : 1 >>>> GluProGlyThrAsnSerGluIleLeuLeuGlyLeuLysTyrValCysProG : 125 ||| ! !! !|||..!||||||||| !.!!|||||||||||| !!|||| ++GluAsnCysGlnAsnGluGluIleLeuAsnSerLeuLysTyrValArgProG 52995 : ......agGAGAACTGTCAGAATGAGGAGATCCTGAACAGTCTCAAGTATGTCCGTCCTG : 50287 126 : lySerGlyPheValProSerPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysG : 145 ||.!!|||!:! !|||!!!||| ||| !!:!!||| !!!!:|||||||||:!!!!.| lyGlyGlyTyrGlnProThrPheThrLeuValGlnLysCysGluValAsnGlyGlnAsnG 50286 : GGGGTGGATACCAGCCCACCTTCACCCTTGTCCAAAAATGTGAGGTGAATGGGCAGAACG : 50227 146 : luGlnLysValPheThrPheLeuLysAsnSerCysPro<->ProThrSerAsp<-><->< : 162 ||!!. ||||||.!!!:!||||||:!!!.! !||| ||| !..!||| luHisProValPheAlaTyrLeuLysAspLysLeuProTyrProTyrAspAspProPheS 50226 : AGCATCCTGTCTTCGCCTACCTGAAGGACAAGCTCCCCTACCCTTATGATGACCCATTTT : 50167 163 : ->LeuLeuGlySerSerSerGlnLeuPheTrpGluProMetLysValHisAspIleArgT : 181 |||:!! !..! !!!.!! :!!.!!|||...|||:!!::: ! |||:!: !| erLeuMetThrAspProLysLeuIleIleTrpSerProValArgArgSerAspValAlaT 50166 : CCCTCATGACCGATCCCAAGCTCATCATTTGGAGCCCTGTGCGCCGCTCAGATGTGGCCT : 50107 182 : rpAsnPheGluLysPheLeuValGlyProAspGlyValProValMetHisTrpPheHisG : 201 ||||||||||||||||||||:!!||||||!!:|||! !||| !! !.!!:: !!.! rpAsnPheGluLysPheLeuIleGlyProGluGlyGluProPheArgArgTyrSerArgT 50106 : GGAACTTTGAGAAGTTCCTCATAGGGCCGGAAGGAGAGCCCTTCCGACGCTACAGCCGCA : 50047 202 : lnAlaProValSerThrValLysSerAspIleLeuGluTyrLeuLys : 216 !|||..!! !!.!:!!:!! !!|||||| !... !|||||| hrPheProThrIleAsnIleGluProAspIleLysArgLeuLeuLys 50046 : CCTTCCCAACCATCAACATTGAGCCTGACATCAAGCGCCTCCTTAAA : 50000 vulgar: SPP00000009_1.0 35 216 . GL397312.1:subseq(15866587,103211) 53211 49999 - 302 M 17 51 G 1 0 M 25 75 G 0 3 M 13 39 G 1 0 M 15 45 5 0 2 I 0 2656 3 0 2 M 50 150 G 0 3 M 4 12 G 0 9 M 55 165 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2012-11-29 ##type DNA # # # seqname source feature start end score strand frame attributes # GL397312.1:subseq(15866587,103211) exonerate:protein2genome:local gene 50000 53211 302 - . gene_id 1 ; sequence SPP00000009_1.0 ; gene_orientation + GL397312.1:subseq(15866587,103211) exonerate:protein2genome:local cds 52999 53211 . - . GL397312.1:subseq(15866587,103211) exonerate:protein2genome:local exon 52999 53211 . - . insertions 3 ; deletions 2 GL397312.1:subseq(15866587,103211) exonerate:protein2genome:local splice5 52997 52998 . - . intron_id 1 ; splice_site "GT" GL397312.1:subseq(15866587,103211) exonerate:protein2genome:local intron 50339 52998 . - . intron_id 1 GL397312.1:subseq(15866587,103211) exonerate:protein2genome:local splice3 50339 50340 . - . intron_id 0 ; splice_site "AG" GL397312.1:subseq(15866587,103211) exonerate:protein2genome:local cds 50000 50338 . - . GL397312.1:subseq(15866587,103211) exonerate:protein2genome:local exon 50000 50338 . - . insertions 12 ; deletions 0 GL397312.1:subseq(15866587,103211) exonerate:protein2genome:local similarity 50000 53211 302 - . alignment_id 1 ; Query SPP00000009_1.0 ; Align 53212 36 51 ; Align 53161 54 75 ; Align 53083 79 39 ; Align 53044 93 45 ; Align 50339 108 150 ; Align 50186 158 12 ; Align 50165 162 165 # --- END OF GFF DUMP --- # -- completed exonerate analysis