SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Any Available uniform cufflinks masking file Bardj RNA Sequencing 1 08-18-2013 10:13 AM
Mapping vs De Novo & Masking or not Eric Fournier General 0 07-05-2011 07:49 AM
Blast query masking behavior fungs Bioinformatics 0 04-09-2010 07:03 AM
Masking reads before aligning is imp. thinkRNA Bioinformatics 0 03-22-2010 10:52 PM
BWA and masking sequences... spindrift Bioinformatics 4 02-17-2010 10:20 AM

Reply
 
Thread Tools
Old 06-01-2011, 01:05 AM   #1
Amy S
Junior Member
 
Location: UCSF

Join Date: Oct 2009
Posts: 2
Default Repeat masking

Hi everyone,

We have whole genome data from a number of tumour-normal pairs, and are looking at our score distributions for somatic changes. We find that some of the highest scores look like errors due to read pile up in repetitive regions when viewed in IGV. Thus we are beginning to think about masking repetitive parts of the genome using repeat masker or something similar.

Does anyone have any tips on how to go about this? I've had a brief look, and to me it seems I could mask up to 50% of the genome using the least stringent criteria. Does anyone have any experience to get me started? Which programs/repeat libraries are the best to use.

Thank you for your help,
Amy.
Amy S is offline   Reply With Quote
Old 06-01-2011, 02:14 AM   #2
jochensupper
Junior Member
 
Location: Germany

Join Date: Nov 2009
Posts: 7
Default

For homology based repeat masking I would use RepeatMasker (www.repeatmasker.org/) together with Repbase (www.girinst.org/repbase/). Its easy to use and usually provides very good results.

If you want to uncover de novo repeats there are many programs you could choose from, e.g., LTRfinder, PILER, RepeatScout.
jochensupper is offline   Reply With Quote
Old 06-02-2011, 03:32 PM   #3
mattanswers
Member
 
Location: Boston

Join Date: Oct 2009
Posts: 65
Default

When you do your alignment, you could align for only unique sequences (sequences that match to only 1 place on the genome). This can be done in bowtie by inserting -m 1.

Also, have a look at The Uniqueome: A mappability resource for short-tag sequencing,http://bioinformatics.oxfordjournals...s.btq640.short and here for a turorial:http://grimmond.imb.uq.edu.au/unique...ary_File_2.pdf
mattanswers is offline   Reply With Quote
Reply

Tags
repeat, repeated sequences, repeatmasker, repeats

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 06:45 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO