Re: [Geuvadis_rna_analysis] Geuvadis RNAseq analysis group TC - slides

Dear all, In attachment the result of mapping for 32-36bp length reads. Most of the hits are to supercontigs. Here is the possible explanation: In alignments a lot of reads map to GL000220.1. My colleague Mat Davis found the following explanation online from Ian Dunham (Ensembl) regarding someone else's question: "In human the ribosomal RNA gene clusters containing the 28S, 18S and 5.8S RNAs are located on the short arms of the acrocentric chromosomes (13, 14, 15, 21 and 22) in long repeating arrays of 40 odd kb for a single unit. You have many copies in the arrays and copy number is variable between individuals. These regions are difficult to assemble into what are now called supercontigs because of this repeating structure and other satellite repeats that lie on these chromosome arms, so although there have been some BACs sequenced representing the clusters, the contigs are not integrated into the chromosome supercontigs. Instead they lie on the set of unplaced supercontigs. Supercontig GL000220.1 for instance is one of these. There are others. In the sequence databases there are several BACS that cover these clusters." So it seems the mapping to rRNAs and to supercontigs is related. Regards, Natalja -- Natalja Kurbatova, PhD Functional Genomics Group, EBI