Gene or accession: CG12276 Gene annotation error Gene CG12276 corresponds to AF193554 (FBgn0029512) Comments: CG12276|CT17770| corresponds to Aos1 (FBgn0029512). ClustalW alignments of the Aos1 and CG12276|CT17770| CDSs are being sent to flybase-updates in a separate message. Subject: Aos1-CG12276|CT17770 ClustalW alignment Aos1 corresponds to CG12276. Aos1 and CG12276|CT17770 ClustalW alignment. Pairwise Scores: CLUSTAL W (1.8) Multiple Sequence Alignments Sequence format is Pearson Sequence 1: gi|6694276|gb|AAF25198.1|AF193 337 aa Sequence 2: gi|7299631|gb|AAF54815.1| 337 aa Start of Pairwise alignments Aligning... Sequences ( 1:2 ) Aligned. Score: 99 Guide tree file created: /net/nfs0/vol1/production/w3nobody/tmp/1047959028373.dnd Start of Multiple Alignment There are 1 groups Aligning... Group 1: Sequences: 2 Score:5488 Alignment Score 2029 CLUSTAL-Alignment file created /net/nfs0/vol1/production/w3nobody/tmp/1047959028373.aln CLUSTAL W (1.8) multiple sequence alignment gi|6694276|gb|AAF25198.1|AF193 MVVDMDTSETAVELTEAENELYDRQIRLWGLESQKRLRTAKILIAGLCGL 50 gi|7299631|gb|AAF54815.1| MVVDMDTSETAVELTEAENELYDRQIRLWGLESQKRLRTAKILIAGLCGL 50 gi|6694276|gb|AAF25198.1|AF193 GAEITKNIILSGVNSVKLLDDKDVTEEDFCSQFLVPRESLNTNRAEASLT 100 gi|7299631|gb|AAF54815.1| GAEITKNIILSGVNSVKLLDDKDVTEEDFCSQFLVPRESLNTNRAEASLT 100 gi|6694276|gb|AAF25198.1|AF193 RARALNPMVDISADREPLKEKTSEFFGQFDVVVVNGATNEELLRIDTICR 150 gi|7299631|gb|AAF54815.1| RARALNPMVDISADREPLKEKTSEFFGQFDVVVVNGATNEELLRIDTICR 150 gi|6694276|gb|AAF25198.1|AF193 DLGVKFIATDVWGTFGFYFASLQKHSYVEDVIKHKVVANSEKKKKYETVS 200 gi|7299631|gb|AAF54815.1| DLGVKFIATDVWGTFGFYFASLQKHSYVEDVIKHKVVANSEKKKKYETVS 200 gi|6694276|gb|AAF25198.1|AF193 IPTQRDVDYPGYSAWLDFDVTEPSYLRKLKRNGPGVLLLSVLQKFRTTHK 250 gi|7299631|gb|AAF54815.1| IPTQRDVDYPGYSAWLDFDVTEPSYLRKLKRNGPGVLLLSVLQKFRTTHK 250 gi|6694276|gb|AAF25198.1|AF193 RDPSYKTREADLELLRGIRDELLPNSILGNEALGLIFAQISPAVAVVGGV 300 gi|7299631|gb|AAF54815.1| RDPSYKTREADLELLRGIRDELLPNSILGDEALGLIFAQISPAVAVVGGV 300 gi|6694276|gb|AAF25198.1|AF193 VAQEVIKVVTKLEAPHRNLFVFDPETCAGYVEAIGAK 337 gi|7299631|gb|AAF54815.1| VAQEVIKVVTKLEAPHRNLFVFDPETCAGYVEAIGAK 337 Subject: Aos1 and SAE1 (plus CG12276) Michael and the Cambridge curators, By alignment, Aos1 and SAE1 seem to be the same gene. The attached file contains the BLASTN and ClustalW alignments for the nucleotide and protein sequences for Aos1 and SAE1. I have just sent in a FB-updates/error report indicating that Aos1 corresponds to CG12276|CT17770. Leyla ================================================================================ SAE1 AF218863 AAF31703.1 Aos1 AF193554 AAF25198.1 SAE1 and Aos1 BLASTN alignment Sequence 1 gi 6934293 Length 1169 (1 .. 1169) (SAE1) Sequence 2 gi 6694275Length 1169 (1 .. 1169) (Aos1) Score = 2173 bits (1130), Expect = 0.0 Identities = 1134/1136 (99%), Positives = 1134/1136 (99%) Query: 4 ctgaaagttgttgcaagaaatttgtattaagttgtgaaaaatacattttaaacttgataa 63 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 9 ctgaaagttgttgcaagaaatttgtattaagttgtgaaaaatacattttaaacttgataa 68 Query: 64 ataaataatggtcgtggatatggacaccagcgagacagccgtggagctcactgaggcgga 123 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 69 ataaataatggtcgtggatatggacaccagcgagacagccgtggagctcactgaggcgga 128 ubiquitin-like protein activati> 1 M V V D M D T S E T A V E L T E A E Query: 124 gaacgagctgtacgacagacaaatccgactctggggtctggaatctcagaaacgccttcg 183 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 129 gaacgagctgtacgacagacaaatccgactctggggtctggaatctcagaaacgccttcg 188 ubiquitin-like protein activati> 19 N E L Y D R Q I R L W G L E S Q K R L R Query: 184 cacggccaagattctcatcgctggactgtgcggcctgggcgccgagataaccaagaacat 243 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 189 cacggccaagattctcatcgctggactgtgcggcctgggcgccgagataaccaagaacat 248 ubiquitin-like protein activati> 39 T A K I L I A G L C G L G A E I T K N I Query: 244 catcctgtccggagtgaactccgtgaagctgctggatgacaaggacgtaaccgaggagga 303 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 249 catcctgtccggagtgaactccgtgaagctgctggatgacaaggacgtaaccgaggagga 308 ubiquitin-like protein activati> 59 I L S G V N S V K L L D D K D V T E E D Query: 304 cttttgttcacaattccttgtcccccgtgaatcgctgaacaccaaccgggccgaagcatc 363 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 309 cttttgttcacaattccttgtcccccgtgaatcgctgaacaccaaccgggccgaagcatc 368 ubiquitin-like protein activati> 79 F C S Q F L V P R E S L N T N R A E A S Query: 364 attgacacgggcacgtgctctcaatcccatggtggacatctccgccgaccgcgagccctt 423 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 369 attgacacgggcacgtgctctcaatcccatggtggacatctccgccgaccgcgagccctt 428 ubiquitin-like protein activati> 99 L T R A R A L N P M V D I S A D R E P L Query: 424 gaaggagaagacctctgagttcttcggtcagttcgacgttgtggtggtcaatggcgcgac 483 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 429 gaaggagaagacctctgagttcttcggtcagttcgacgttgtggtggtcaatggcgcgac 488 ubiquitin-like protein activati> 119 K E K T S E F F G Q F D V V V V N G A T Query: 484 caacgaggaactgttgcgcattgacaccatttgccgggacctgggcgtcaagttcatagc 543 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 489 caacgaggaactgttgcgcattgacaccatttgccgggacctgggcgtcaagttcatagc 548 ubiquitin-like protein activati> 139 N E E L L R I D T I C R D L G V K F I A Query: 544 caccgatgtgtggggcaccttcgggttctactttgccagtctgcagaagcacagctacgt 603 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 549 caccgatgtgtggggcaccttcgggttctactttgccagtctgcagaagcacagctacgt 608 ubiquitin-like protein activati> 159 T D V W G T F G F Y F A S L Q K H S Y V Query: 604 cgaggatgttattaatcacaaggtcgtagccaattcggagaagaaaaagaagtatgaaac 663 ||||||||||||||| |||||||||||||||||||||||||||||||||||||||||||| Sbjct: 609 cgaggatgttattaagcacaaggtcgtagccaattcggagaagaaaaagaagtatgaaac 668 ubiquitin-like protein activati> 179 E D V I K H K V V A N S E K K K K Y E T Query: 664 cgtgtccattccaacgcagcgtgatgtggattatcccggctactctgcctggcttgattt 723 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 669 cgtgtccattccaacgcagcgtgatgtggattatcccggctactctgcctggcttgattt 728 ubiquitin-like protein activati> 199 V S I P T Q R D V D Y P G Y S A W L D F Query: 724 cgatgttaccgagcctagttatttgcgaaaattgaagcgcaatggcccgggtgttttgct 783 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 729 cgatgttaccgagcctagttatttgcgaaaattgaagcgcaatggcccgggtgttttgct 788 ubiquitin-like protein activati> 219 D V T E P S Y L R K L K R N G P G V L L Query: 784 attaagtgtcctgcaaaagttccgcacaacccacaagcgggatcccagctacaagacccg 843 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 789 attaagtgtcctgcaaaagttccgcacaacccacaagcgggatcccagctacaagacccg 848 ubiquitin-like protein activati> 239 L S V L Q K F R T T H K R D P S Y K T R Query: 844 agaagctgacctggaattgcttcgcggaattcgcgacgaactgctgcccaactctatact 903 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 849 agaagctgacctggaattgcttcgcggaattcgcgacgaactgctgcccaactctatact 908 ubiquitin-like protein activati> 259 E A D L E L L R G I R D E L L P N S I L Query: 904 gggcgacgaagctctaggtctgatttttgcacagatctcgccagcggtcgctgttgtcgg 963 |||| ||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 909 gggcaacgaagctctaggtctgatttttgcacagatctcgccagcggtcgctgttgtcgg 968 ubiquitin-like protein activati> 279 G N E A L G L I F A Q I S P A V A V V G Query: 964 cggcgttgtggcacaagaggtcatcaaggtggttactaaattggaggcgccacaccgtaa 1023 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 969 cggcgttgtggcacaagaggtcatcaaggtggttactaaattggaggcgccacaccgtaa 1028 ubiquitin-like protein activati> 299 G V V A Q E V I K V V T K L E A P H R N Query: 1024 tctgttcgtttttgacccagagacttgtgccggatacgttgaggccatcggagcgaagtg 1083 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 1029 tctgttcgtttttgacccagagacttgtgccggatacgttgaggccatcggagcgaagtg 1088 ubiquitin-like protein activati> 319 L F V F D P E T C A G Y V E A I G A K ^^ Query: 1084 aatttaccacaaagggcttatttataaaaatctgaatatcaaataacaccagttgt 1139 |||||||||||||||||||||||||||||||||||||||||||||||||||||||| Sbjct: 1089 aatttaccacaaagggcttatttataaaaatctgaatatcaaataacaccagttgt 1144 ubiquitin-like protein activati> 338 ^ CPU time: 0.11 user secs. 0.02 sys. secs 0.13 total secs. Gapped Lambda K H 1.33 0.621 1.12 Matrix: blastn matrix:1 \-2 Gap Penalties: Existence: 5, Extension: 2 Number of Hits to DB: 1 Number of Sequences: 0 Number of extensions: 1 Number of successful extensions: 1 Number of sequences better than 10.0: 1 length of query: 1169 length of database: 706,884,452 effective HSP length: 23 effective length of query: 1146 effective length of database: 692,976,559 effective search space: 794151136614 effective search space used: 794151136614 T: 0 A: 0 X1: 6 (11.5 bits) X2: 26 (50.0 bits) S1: 12 (23.8 bits) S2: 19 (37.2 bits) _______________________________________________________________________ CLUSTAL W (1.8) Multiple Sequence Alignments Sequence format is Pearson Sequence 1: gi|6694276|gb|AAF25198.1|AF193 337 aa Sequence 2: gi|6934294|gb|AAF31703.1|AF218 337 aa Start of Pairwise alignments Aligning... Sequences ( 1:2 ) Aligned. Score: 99 Guide tree file created: /net/nfs0/vol1/production/w3nobody/tmp/23494959027282.dnd Start of Multiple Alignment There are 1 groups Aligning... Group 1: Sequences: 2 Score:5480 Alignment Score 2024 CLUSTAL-Alignment file created /net/nfs0/vol1/production/w3nobody/tmp/23494959027282.aln Your guide tree: 23494959027282.dnd (gi|6694276|gb|AAF25198.1| AF193:0.00297 ,gi|6934294|gb|AAF31703.1| AF218:0.00297 ); Your Multiple Sequence Alignment: 23494959027282.aln CLUSTAL W (1.8) multiple sequence alignment gi|6694276|gb|AAF25198.1|AF193 MVVDMDTSETAVELTEAENELYDRQIRLWGLESQKRLRTAKILIAGLCGL 50 gi|6934294|gb|AAF31703.1|AF218 MVVDMDTSETAVELTEAENELYDRQIRLWGLESQKRLRTAKILIAGLCGL 50 gi|6694276|gb|AAF25198.1|AF193 GAEITKNIILSGVNSVKLLDDKDVTEEDFCSQFLVPRESLNTNRAEASLT 100 gi|6934294|gb|AAF31703.1|AF218 GAEITKNIILSGVNSVKLLDDKDVTEEDFCSQFLVPRESLNTNRAEASLT 100 gi|6694276|gb|AAF25198.1|AF193 RARALNPMVDISADREPLKEKTSEFFGQFDVVVVNGATNEELLRIDTICR 150 gi|6934294|gb|AAF31703.1|AF218 RARALNPMVDISADREPLKEKTSEFFGQFDVVVVNGATNEELLRIDTICR 150 gi|6694276|gb|AAF25198.1|AF193 DLGVKFIATDVWGTFGFYFASLQKHSYVEDVIKHKVVANSEKKKKYETVS 200 gi|6934294|gb|AAF31703.1|AF218 DLGVKFIATDVWGTFGFYFASLQKHSYVEDVINHKVVANSEKKKKYETVS 200 gi|6694276|gb|AAF25198.1|AF193 IPTQRDVDYPGYSAWLDFDVTEPSYLRKLKRNGPGVLLLSVLQKFRTTHK 250 gi|6934294|gb|AAF31703.1|AF218 IPTQRDVDYPGYSAWLDFDVTEPSYLRKLKRNGPGVLLLSVLQKFRTTHK 250 gi|6694276|gb|AAF25198.1|AF193 RDPSYKTREADLELLRGIRDELLPNSILGNEALGLIFAQISPAVAVVGGV 300 gi|6934294|gb|AAF31703.1|AF218 RDPSYKTREADLELLRGIRDELLPNSILGDEALGLIFAQISPAVAVVGGV 300 gi|6694276|gb|AAF25198.1|AF193 VAQEVIKVVTKLEAPHRNLFVFDPETCAGYVEAIGAK 337 gi|6934294|gb|AAF31703.1|AF218 VAQEVIKVVTKLEAPHRNLFVFDPETCAGYVEAIGAK 337