Command line: [exonerate --exhaustive yes -m p2g --showtargetgff -q ../SelS-human.fa -t SelS-NIVO01025374.1.subseq] Hostname: [sitdoc] C4 Alignment: ------------ Query: SPP00000024_1.0 # Protein # Selenoprotein S (SelS) # Homo sapiens # Complete Target: NIVO01025374.1:subseq(1283219,100530) Ammotragus lervia isolate NAZ scaffold355_len2295547_cov70, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 820 Query range: 0 -> 189 Target range: 50530 -> 49960 1 : MetGluArgGlnGluGluSerLeuSerAlaArgProAlaLeuGluThrGluGlyLeuArg : 20 |||||||||.!.! !!!:..!||||||.!!||||||.!!|||:!!||||||||||||!.! MetGluArgAspGlyAspGlnLeuSerThrArgProThrLeuLysThrGluGlyLeuHis 50530 : ATGGAGCGGGACGGGGATCAGCTGTCCACACGGCCCACCCTCAAGACCGAGGGGCTGCAT : 50473 21 : PheLeuHisThrThrValGlySerLeuLeuAlaThrTyrGlyTrpTyrIleValPheSer : 40 |||||||||..!|||||||||||||||! !|||||||||||||||||||||:!!|||||| PheLeuHisValThrValGlySerLeuProAlaThrTyrGlyTrpTyrIleIlePheSer 50472 : TTCCTTCACGTCACAGTGGGCTCCCTGCCGGCCACTTACGGCTGGTACATCATCTTCAGC : 50413 41 : CysIleLeuLeuTyrValValPheGlnLysLeuSerAlaArgLeuArgAlaLeuArgGln : 60 ||||||||||||||||||||||||||||||||||||.!!||||||||||||||||||||| CysIleLeuLeuTyrValValPheGlnLysLeuSerThrArgLeuArgAlaLeuArgGln 50412 : TGCATCCTTCTCTACGTGGTGTTTCAGAAGCTTTCCACCCGGTTGAGGGCCTTGAGGCAG : 50353 61 : ArgGlnLeuAspArgAlaAlaAlaAlaValGluProAspValValValLysArgGlnGlu : 80 |||!!.||||||!:!||||||||||||:!!|||||||||:!!|||||||||!:!|||||| ArgHisLeuAspGlnAlaAlaAlaAlaLeuGluProAspIleValValLysGlnGlnGlu 50352 : AGGCACTTAGATCAAGCTGCCGCTGCTCTGGAACCTGATATTGTTGTTAAACAACAGGAG : 50293 81 : AlaLeuAlaAlaAlaArgLeuLysMetGlnGluGluLeuAsnAlaGlnValGluLysHis : 100 |||||||||||||||!.!|||||||||||||||||||||||||||||||||||||||||| AlaLeuAlaAlaAlaHisLeuLysMetGlnGluGluLeuAsnAlaGlnValGluLysHis 50292 : GCTTTAGCAGCTGCTCATTTGAAAATGCAAGAAGAATTAAATGCACAAGTTGAaaaacat : 50233 101 : LysGluLysLeuLysGlnLeuGluGluGluLysArgArgGlnLysIleGluMetTrpAsp : 120 ||||||||||||:::||||||||| !!||||||||||||||||||||||||||||||||| LysGluLysLeuArgGlnLeuGlu***GluLysArgArgGlnLysIleGluMetTrpAsp 50232 : aaagagaaattaagacagCTTGAATAAGAAAAAAGGAGACAGAAGATTGAAATGTGGGAT : 50173 121 : SerMetGlnGluGlyLysSerTyrLysGlyAsnAlaLysLysProGlnGluGluAspSer : 140 |||||||||||||||||||||||||||||||||.!!!:!:!!|||||||||||||||||| SerMetGlnGluGlyLysSerTyrLysGlyAsnThrArgGluProGlnGluGluAspSer 50172 : AGCATGCAAGAAGGAAAAAGTTATAAAGGAAATACAAGAGAGCCTCAGGAAGAAGACAGT : 50113 141 : ProGlyProSerThrSerSerValLeu<->LysArgLysSerAspArgLysProLeuArg : 159 ||||||||||||||||||||||||:!! |||||||||||||||||||||||||||||| ProGlyProSerThrSerSerValIleProLysArgLysSerAspArgLysProLeuArg 50112 : CCTGGGCCTTCTACTTCATCAGTCATCCCGAAACGAAAATCCGACAGAAAGCCTTTGCGA : 50053 160 : GlyGlyGlyTyrAsnProLeuSerGlyGluGlyGlyGlyAlaCysSerTrpArgProGly : 179 |||||||||! !|||||||||||||||||||||||||||.!!|||||||||||||||||| GlyGlyGlySerAsnProLeuSerGlyGluGlyGlyGlyThrCysSerTrpArgProGly 50052 : GGAGGAGGTTCCAATCCTTTGTCTGGTGAAGGTGGTGGAACCTGCTCCTGGAGACCTGGA : 49993 180 : ArgArgGlyProSerSerGlyGlyUnkGly : 189 ||| !!|||||||||||||||||| ||| Arg***GlyProSerSerGlyGly***Gly 49992 : CGCTGAGGACCCTCatctggtggatgaggc : 49961 vulgar: SPP00000024_1.0 0 189 . NIVO01025374.1:subseq(1283219,100530) 50530 49960 - 820 M 149 447 G 0 3 M 40 120 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-17 ##type DNA # # # seqname source feature start end score strand frame attributes # NIVO01025374.1:subseq(1283219,100530) exonerate:protein2genome:local gene 49961 50530 820 - . gene_id 1 ; sequence SPP00000024_1.0 ; gene_orientation . NIVO01025374.1:subseq(1283219,100530) exonerate:protein2genome:local cds 49961 50530 . - . NIVO01025374.1:subseq(1283219,100530) exonerate:protein2genome:local exon 49961 50530 . - . insertions 3 ; deletions 0 NIVO01025374.1:subseq(1283219,100530) exonerate:protein2genome:local similarity 49961 50530 820 - . alignment_id 1 ; Query SPP00000024_1.0 ; Align 50531 1 447 ; Align 50081 150 120 # --- END OF GFF DUMP --- # -- completed exonerate analysis