gep.wustl.edu

If all four restriction digests matched and there are no issues with the assembly in this region (e.g. high quality discrepancies, unaligned high quality, etc), then I would not be too concerned about the inconsistent mate pairs. In particular, if the inconsistent mate pairs are in the correct relat...

Many of the genuine high quality discrepancies I have seen in the D. ananassae projects are placed in regions that matched known transposons in D. ananassae . The consensus sequences in the repetitive regions are marked with blue "repeat" tags. Consequently, some of the high quality discrepancies co...

The calc_total_fragSizes.pl script parses the file which contains the observed fragment sizes in the real digests (i.e. fragSizes.txt ) and calculates the sum of the fragment sizes for each of the real digests. We expect the real digest sizes to be approximately the same for all four restriction dig...

You are correct that the GEP annotation strategy would not allow you to annotate species-specific gene (e.g. genes that are only found in D. biarmipes). The problem is that we cannot confidently construct the gene models using only RNA-seq data. For example, while the oases and cufflinks transcripts...

The GEP web framework has been reset for the Spring 2013 semester. Summary of key changes 1. New D. ananassae annotation projects from the 3L reference region available 2. Updated GEP Web Framework and annotation materials (Note: new NCBI BLAST output format) 3. D. ananassae sequence improvement pro...

The change in the genomic sequence is unlikely to be an error in the consensus because there multiple RNA-seq reads aligned to this region and they agree with the consensus. A FlyBase TBLASTN search of the first CDS of the C isoform (17_1548_0) against the genomic assemblies in different Drosophila ...

In general, we expect to find some short spurious matches in the dot plot. Some of these matches can occur by chance and others can be attributed to repetitive or low complexity sequences. In fact, this is one of the reasons why we use the BLAST E-value threshold and the low complexity filter to try...

Missing exons in CG33978 Given that exons 3_1556_1 and 2_1566_1 overlap with each other in D. melanogaster , I would try to map the larger exon (3_1556_1) first. From your previous analysis of the adjacent exons, we know that we need this exon to have a phase 1 acceptor and a phase 0 donor and that...

In both cases, the exons are likely to be too small or too weakly conserved to be detected using BLAST searches. Given the lack of expression data, conservation and computational predictions, we will use conservation of exon length and biological constraints to help us pick the best candidate. Missi...

Yes, in frame stop codons in the gene model is typically caused by incompatible splice donor and acceptor sites. When the sum of the donor and acceptor phases is not either 0 or 3, it will introduce a frame shift in the translation. In most cases, translating in the incorrect reading frame would res...

> Will I miss important info if I only use the blastx from the genome browser or should I still do a NCBI blastx as a double check of the coordinates? Yes, because the blastx alignment track on the genome browser is created by searching the contig sequence against the set of full-length D. melanogas...

Because D. biarmipes is evolutionarily further away from D. melanogaster than D. erecta , we anticipate that the D. biarmipes projects will be more difficult to annotate than D. erecta . The difficulty levels of the D. biarmipes projects range from 4 to 6, with 4 being the easiest. D. biarmipes proj...

As you have mentioned, the gene record for the D. melanogaster sls gene has changed in June 2012. flybase_sls_update.png However, the most likely reason that you cannot detect the first exon (43_374_0) with blastx is that the low complexity filter and compositional adjustments were on when you perfo...

> Do I need to recheck the gene model? No, you can ignore the warning from the Gene Model Checker in this case. Having a stop codon at the beginning of the last coding exon is unusual and it merits additional explanations in the GEP Annotation Report. However, the explanation in the Annotation Repor...

Hello Everyone, The Project Management System has been reset for Fall 2012. A new set of 70 projects derived from the Muller F element of D. biarmipes is now available to be claimed on the Project Management System. You can view the projects using the GEP UCSC Genome Browser Mirror by selecting “D. ...