Command line: [exonerate -m p2g --showtargetgff -q SPS2/SPS2_U.fa -t SPS2/SUBSEQ/KB228942.1_0_subseq.txt --exhaustive yes] Hostname: [luke] C4 Alignment: ------------ Query: SPS2 Target: KB228942.1:subseq(3077522,51376) Odobenus rosmarus divergens isolate Ivan unplaced genomic scaffold Scaffold65, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 2187 Query range: 25 -> 483 Target range: 51376 -> 49999 26 : LeuThrLeuGlyArgSerPheSerAsnTyrArgProPheGluProGlnAlaLeuGlyLeu : 45 |||:!!|||||||||.!!||||||!:!||||||||||||||||||||||||||||||||| LeuSerLeuGlyArgGlyPheSerSerTyrArgProPheGluProGlnAlaLeuGlyLeu 51376 : TTGTCTCTGGGCCGGGGCTTCTCGAGCTACCGGCCCTTCGAGCCCCAGGCGCTGGGCCTC : 51319 46 : SerProSerTrpArgLeuThrGlyPheSerGlyMetLysGlyUnkGlyCysLysValPro : 65 |||||||||||||||||||||||||||||||||||||||||| ||||||||||||||| SerProSerTrpArgLeuThrGlyPheSerGlyMetLysGly***GlyCysLysValPro 51318 : AGCCCGAGCTGGCGGCTGACCGGCTTCTCCGGCATGAAGGGCTGAGGCTGCAAGGTCCCC : 51259 66 : GlnGluAlaLeuLeuLysLeuLeuAlaGlyLeuThrArgProAspValArgProProLeu : 85 ||||||.!!||||||||||||||||||||||||||||||||||||||||||||||||||| GlnGluThrLeuLeuLysLeuLeuAlaGlyLeuThrArgProAspValArgProProLeu 51258 : CAGGAGACGCTGCTCAAACTCCTGGCGGGACTGACGCGACCAGACGTGCGGCCCccgctg : 51199 86 : GlyArgGlyLeuValGlyGlyGlnGluGluAlaSerGlnGluAlaGlyLeuProAlaGly : 105 ||| |||||||||||||||!!.|||||||||:!!||||||!.!||||||||||||||| GlyLeuGlyLeuValGlyGlyHisGluGluAlaAlaGlnGluGlyGlyLeuProAlaGly 51198 : ggcctgggcctggtcGGAGGCCATGAAGAGGCGGCCCAGGAAGGCGGCCTACCGGCCGGG : 51139 106 : AlaGlyProSerProThrPheProAlaLeuGlyIleGlyMetAspSerCysValIlePro : 125 |||! !|||!:!|||||||||||||||||||||||||||||||||||||||||||||||| AlaGluProAsnProThrPheProAlaLeuGlyIleGlyMetAspSerCysValIlePro 51138 : GCCGAACCCAACCCCACCTTTCCAGCCCTGGGCATTGGGATGGACTCCTGCGTCATACCC : 51079 126 : LeuArgHisGlyGlyLeuSerLeuValGlnThrThrAspPhePheTyrProLeuValGlu : 145 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| LeuArgHisGlyGlyLeuSerLeuValGlnThrThrAspPhePheTyrProLeuValGlu 51078 : CTGAGGCACGGGGGTCTGTCACTGGTGCAGACTACAGACTTCTTTTACCCCTTGGTGGAA : 51019 146 : AspProTyrMetMetGlyArgIleAlaCysAlaAsnValLeuSerAspLeuTyrAlaMet : 165 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| AspProTyrMetMetGlyArgIleAlaCysAlaAsnValLeuSerAspLeuTyrAlaMet 51018 : GATCCCTACATGATGGGGCGCATAGCGTGTGCCAACGTGCTGAGTGACCTCTACGCCATG : 50959 166 : GlyIleThrGluCysAspAsnMetLeuMetLeuLeuSerValSerGlnSerMetSerGlu : 185 ||||||||||||||||||||||||||||||||||||||||||||||||||||||!!!||| GlyIleThrGluCysAspAsnMetLeuMetLeuLeuSerValSerGlnSerMetThrGlu 50958 : GGCATTACTGAATGTGACAACATGTTGATGTTACTCAGCGTCAGCCAGAGTATGACTGAG : 50899 186 : GluGluArgGluLysValThrProLeuMetValLysGlyPheArgAspAlaAlaGluGlu : 205 |||||||||||||||:!!||||||||||||:!!||||||||||||||||||||||||||| GluGluArgGluLysIleThrProLeuMetIleLysGlyPheArgAspAlaAlaGluGlu 50898 : GAGGAGCGAGAAAAGATAACACCACTCATGATCAAAGGCTTTCGTGATGCTGCTGAGGAA : 50839 206 : GlyGlyThrAlaValThrGlyGlyGlnThrValValAsnProTrpIleIleIleGlyGly : 225 |||||||||||||||||||||||||||||||||||||||||||||||||||:!!|||||| GlyGlyThrAlaValThrGlyGlyGlnThrValValAsnProTrpIleIleValGlyGly 50838 : GGAGGGACTGCAGTGACTGGTGGACAAACGGTGGTAAACCCTTGGATTATCGTCGGTGGG : 50779 226 : ValAlaThrValValCysGlnProAsnGluPheIleMetProAspSerAlaValValGly : 245 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| ValAlaThrValValCysGlnProAsnGluPheIleMetProAspSerAlaValValGly 50778 : GTTGCCACTGTGGTATGTCAGCCAAATGAATTCATCATGCCCGACAGTGCAGTTGTTGGG : 50719 246 : AspValLeuValLeuThrLysProLeuGlyThrGlnValAlaValAsnAlaHisGlnTrp : 265 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| AspValLeuValLeuThrLysProLeuGlyThrGlnValAlaValAsnAlaHisGlnTrp 50718 : GATGTGCTCGTGTTAACCAAACCCTTAGGAACCCAAGTTGCTGTCAATGCCCACCAATGG : 50659 266 : LeuAspAsnProGluArgTrpAsnLysValLysMetValValSerArgGluGluValGlu : 285 |||||||||||||||||||||||||||:::...||||||||||||||||||||||||||| LeuAspAsnProGluArgTrpAsnLysIleAsnMetValValSerArgGluGluValGlu 50658 : CTGGATAATCCGGAgagatggaataaaataaatatggtgGTCTCCAGAGAAGAAGTAGAG : 50599 286 : LeuAlaTyrGlnGluAlaMetPheAsnMetAlaThrLeuAsnArgThrAlaAlaGlyLeu : 305 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| LeuAlaTyrGlnGluAlaMetPheAsnMetAlaThrLeuAsnArgThrAlaAlaGlyLeu 50598 : CTGGCTTATCAGGAAGCCATGTTCAACATGGCTACCCTAAACAGAACTGCTGCTGGATTA : 50539 306 : MetHisThrPheAsnAlaHisAlaAlaThrAspIleThrGlyPheGlyIleLeuGlyHis : 325 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| MetHisThrPheAsnAlaHisAlaAlaThrAspIleThrGlyPheGlyIleLeuGlyHis 50538 : ATGCACACGTTTAATGCACACGCAGCCACGGATATCACAGGCTTTGGCATCCTAGGACAC : 50479 326 : SerGlnAsnLeuAlaLysGlnGlnArgAsnGluValSerPheValIleHisAsnLeuPro : 345 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| SerGlnAsnLeuAlaLysGlnGlnArgAsnGluValSerPheValIleHisAsnLeuPro 50478 : TCTCAGAACCTTGCAAAACAACAAAGAAACGAGGTGTCCTTTGTCATTCATAATCTGCCA : 50419 346 : IleIleAlaLysMetAlaAlaValSerLysAlaSerGlyArgPheGlyLeuLeuGlnGly : 365 |||||||||||||||||||||:!!|||||||||||||||||||||||||||||||||||| IleIleAlaLysMetAlaAlaIleSerLysAlaSerGlyArgPheGlyLeuLeuGlnGly 50418 : ATTATTGCCAAGATGGCTGCCATCAGCAAGGCCAGTGGGCGCTTTGGGCTCCTTCAAGGA : 50359 366 : ThrSerAlaGluThrSerGlyGlyLeuLeuIleCysLeuProArgGluGlnAlaAlaArg : 385 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| ThrSerAlaGluThrSerGlyGlyLeuLeuIleCysLeuProArgGluGlnAlaAlaArg 50358 : ACCTCAGCAGAAACTTCTGGGGGATTACTGATTTGTCTGCCAAGAGAACAGGCAGCTCGC : 50299 386 : PheCysSerGluIleLysSerSerLysTyrGlyGluGlyHisGlnAlaTrpIleValGly : 405 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| PheCysSerGluIleLysSerSerLysTyrGlyGluGlyHisGlnAlaTrpIleValGly 50298 : TTTTGTTCTGAAATCAAATCTTCCAAGTATGGAGAGGGTCACCAAGCATGGATCGTCGGT : 50239 406 : IleValGluLysGlyAsnArgThrAlaArgIleIleAspLysProArgValIleGluVal : 425 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| IleValGluLysGlyAsnArgThrAlaArgIleIleAspLysProArgValIleGluVal 50238 : ATTGTGGAAAAGGGAAACCGGACAGCCCGGATCATTGACAAGCCTCGAGTTATTGAGGTC : 50179 426 : LeuProArgGlyAlaThrAlaAlaValLeuAlaProAspSerSerAsnAlaSerSerGlu : 445 |||||||||||||||||||||.!!!.!||||||||||||!:!|||!:!|||||||||||| LeuProArgGlyAlaThrAlaThrAlaLeuAlaProAspAsnSerSerAlaSerSerGlu 50178 : CTACCTCGTGGGGCCACTGCTACTGCTCTTGCTCCTGACAATTCCAGTGCCTCCTCTGAG : 50119 446 : ProSerSerUnkAspGluArgThrGluValValTrpThrLeuGluProLeuSerThrIle : 465 ||||||||| |||! !! !.!!|||||||||||||||||||||||||||! !|||||| ProSerSer***AspGlyIleAlaGluValValTrpThrLeuGluProLeuTyrThrIle 50118 : CCTAGTTCGTGAGATGGAATAGCAGAAGTTGTTTGGACCTTAGAGCCTCTGTACACAATC : 50059 466 : ThrAspGlySerGlnGluLeuIleVal<->ArgAsnPheGlnArgArgLeuProAla : 483 ! !:!!! !||||||......|||::: ||||||||||||||||||||||||||| IleAsnAspSerGlnAsnPheIleLeuArgArgAsnPheGlnArgArgLeuProAla 50058 : ATTAACGATTCTCAGAactttattttaagaagaaatttccAAAGAAGGCTGCCTGCA : 50000 vulgar: 25 483 . KB228942.1:subseq(3077522,51376) 51376 49999 - 2187 M 449 1347 G 0 3 M 9 27 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2015-11-26 ##type DNA # # # seqname source feature start end score strand frame attributes # KB228942.1:subseq(3077522,51376) exonerate:protein2genome:local gene 50000 51376 2187 - . gene_id 1 ; sequence ; gene_orientation . KB228942.1:subseq(3077522,51376) exonerate:protein2genome:local cds 50000 51376 . - . KB228942.1:subseq(3077522,51376) exonerate:protein2genome:local exon 50000 51376 . - . insertions 3 ; deletions 0 KB228942.1:subseq(3077522,51376) exonerate:protein2genome:local similarity 50000 51376 2187 - . alignment_id 1 ; Query ; Align 51377 26 1347 ; Align 50027 475 27 # --- END OF GFF DUMP --- # -- completed exonerate analysis