Command line: [exonerate -q dselM_prot.fasta -t D_sim_secDNA.fasta -m protein2genome -W 3 --showtargetgff yes] C4 Alignment display: Model: protein2genome Raw score: 1027 Aligned positions 0->249 of query Aligned positions 138->997 of target Query: dSelM Target: gnl|Dsim-washu-w501-asm|Contig1042.3 1 : MetProProLysArgAsnLysLysAlaGlu >>>> Target Intron 1 >>>> AlaPr : 11 ||||||||||||||||||||||||.!!||| 59 bp ||||| MetProProLysArgAsnLysLysThrGlu AlaPr 139 : ATGCCACCAAAACGGAACAAAAAAACAGAGgt.........................agGCGCC : 230 12 : oIleAlaGluArgAspAlaGlyGluGluLeuAspProAsnAlaProValLeuTyrValGluHis : 32 |:!!|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| oValAlaGluArgAspAlaGlyGluGluLeuAspProAsnAlaProValLeuTyrValGluHis 231 : AGTCGCTGAGCGGGATGCAGGCGAAGAGCTGGATCCCAATGCGCCGGTGCTCTACGTAGAGCAC : 293 33 : CysArgSerUnkArgValPheArgArgArgAlaGluGluLeuHisSerAlaLeuArgGluArgG : 54 ||||||||| |||||||||||||||||||||||||||||||||||||||||||||||||||| CysArgSer***ArgValPheArgArgArgAlaGluGluLeuHisSerAlaLeuArgGluArgG 294 : TGCCGATCCTGACGAGTCTTTCGTCGTCGAGCGGAGGAGCTGCACTCCGCTCTCCGTGAGCGTG : 359 55 : lyLeuGlnGlnLeuGlnLeuGlnLeuAsnAlaLeuGlyAlaProArgArgGlyAlaPheGluLe : 75 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| lyLeuGlnGlnLeuGlnLeuGlnLeuAsnAlaLeuGlyAlaProArgArgGlyAlaPheGluLe 360 : GTCTCCAACAGCTTCAGCTCCAGCTGAATGCACTTGGAGCACCGCGGCGAGGAGCCTTCGAGTT : 422 76 : uSerLeuSerAlaGlyGlyMetGlyLysGlnGluGlnValAlaLeuTrpSerGlyLeuLysArg : 96 ||||||||||||||||||||||||||||||||||||||||! !||||||||||||||||||||| uSerLeuSerAlaGlyGlyMetGlyLysGlnGluGlnValGluLeuTrpSerGlyLeuLysArg 423 : GAGTTTGAGCGCTGGCGGAATGGGGAAGCAGGAGCAGGTGGAACTTTGGTCGGGTCTAAAGCGG : 485 97 : GlyProProArgAlaArgLysPheProThrValGluGluValTyrAspArgIleValGlyIleL : 118 |||||||||||||||! !|||||||||||||||||||||||||||||||||||||||||||||| GlyProProArgAlaLeuLysPheProThrValGluGluValTyrAspArgIleValGlyIleL 486 : GGTCCACCACGTGCTCTCAAGTTTCCCACTGTGGAGGAGGTATACGATCGGATTGTAGGGATTC : 551 119 : euGlyAspGlnGlnGluSerLysGluGlnIleAsnThrGlnLysLeuSerLysIleAspLeuPr : 139 ||..!|||||||||||||||||||||! !! !||||||||||||! !||||||||||||||| ! euAsnAspGlnGlnGluSerLysGluLeuThrAsnThrGlnLysSerSerLysIleAspLeuTh 552 : TAAATGACCAGCAGGAGTCCAAGGAACTGACCAACACGCAGAAGTCATCGAAAATTGATTTAAC : 614 140 : oGlySerGluAlaLeuAlaSerProLysLysSerGluSerThrGluGluAlaGlnGluAsnGlu : 160 !! !||||||||| !.!!||||||||||||||||||||||||||||||||||||! !!.!||| rGluSerGluAlaLysThrSerProLysLysSerGluSerThrGluGluAlaGlnGlyThrGlu 615 : GGAGTCCGAGGCAAAAACATCGCCAAAGAAGTCGGAGTCCACAGAGGAAGCACAGGGAACTGAG : 677 161 : Ala >>>> Target Intron 2 >>>> ProThrSerThrSerThrSerArgLysSerLy : 172 !.! 53 bp !!! !:!! !!|||||||||||||||||||| Val ThrLysAlaProSerThrSerArgLysSerLy 678 : GTAgt.........................aaACAAAGGCACCTTCAACTTCGCGTAAATCCAA : 766 173 : sLysGluGlnLysSerGluGluGluProThrGlnValAspSerLysGluAlaLysGlnSerLys : 193 |||||||||||||||||||||||||! !|||||||||||||||||||||||||||||||||||| sLysGluGlnLysSerGluGluGluGlnThrGlnValAspSerLysGluAlaLysGlnSerLys 767 : AAAGGAACAGAAATCCGAGGAAGAGCAAACTCAAGTCGATTCAAAGGAAGCGAAACAATCCAAA : 829 194 : GluLeuValLysThrLysArgGlnProLysAlaGlnLysLysGlnAlaLysAlaSerGluSerG : 215 ||||||!.!|||||||||||||||||||||||||||||||||! !||||||||||||||||||| GluLeuAlaLysThrLysArgGlnProLysAlaGlnLysLysProAlaLysAlaSerGluSerG 830 : GAGCTGGCCAAAACGAAGAGACAACCGAAAGCCCAAAAGAAGCCGGCAAAGGCATCCGAATCCC : 895 216 : lnGluGluValAlaGluAspLysProProSerSerGlnLysArgLysArgThrThrArgSerSe : 236 ||:!!||||||.!!|||!!:|||||| !||||||||||||||||||||||||||||||||||| lnLysGluValThrGluGluLysProValSerSerGlnLysArgLysArgThrThrArgSerSe 896 : AAAAAGAAGTCACCGAAGAGAAGCCTGTAAGCAGCCAGAAACGGAAGCGGACCACCAGATCCTC : 958 237 : rThrAspGluAlaThrAlaGlyAlaLysArgArgArg : 249 |||||||||||||||||||.!!||||||||||||||| rThrAspGluAlaThrAlaSerAlaLysArgArgArg 959 : CACAGACGAGGCTACTGCCAGTGCCAAGCGGAGGAGG : 997 vulgar: dSelM 0 249 . gnl|Dsim-washu-w501-asm|Contig1042.3 138 997 + 1027 M 10 30 5 0 2 I 0 55 3 0 2 M 152 456 5 0 2 I 0 49 3 0 2 M 87 261