SEQanswers

Go Back   SEQanswers > General



Similar Threads
Thread Thread Starter Forum Replies Last Post
tabix and 1000 genomes data Alessandra Bioinformatics 8 09-05-2013 12:13 AM
Interpreting 1000 Genomes data ashkot Bioinformatics 5 01-05-2012 12:27 AM
Annotating 1000 Genomes data ashkot Bioinformatics 4 12-22-2011 10:53 AM
1000 Genomes Data/ Exon targetted Firebird Bioinformatics 27 02-17-2011 12:08 PM
need 1000 genomes data for just one gene michelle.lupton Bioinformatics 11 08-09-2010 01:00 PM

Reply
 
Thread Tools
Old 03-16-2011, 08:12 AM   #1
RichardRocca
Junior Member
 
Location: USA

Join Date: Mar 2011
Posts: 2
Default 1000 Genomes Data

I am trying to count short tandem repeats (STR) from 1000 Genomes Project data using the low coverage BAM files, but am getting the same repeats even on known fast STRs. Should I be looking at other files within the project to do this? I am using SAMTools.
RichardRocca is offline   Reply With Quote
Old 03-16-2011, 12:11 PM   #2
krobison
Senior Member
 
Location: Boston area

Join Date: Nov 2007
Posts: 747
Default

Depending on the aligner & parameters used, these may not have aligned well. In the worst case, no single read can contain the entire repeat array AND unique sequence on either side to enable conclusive mapping -- even if the repeat is not strictly too long, very long repeats will be unlikely to have randomly fragmented in such a manner that they can be uniquely mapped.

STRs are definitely findable in the Watson & Venter BAM files & you might test out your method there. These were done with long read (454 & Sanger respectively) technologies that can read through all but the most ginormous repeats.

My intern last summer developed a tool for this purpose; anyone interested in beta testing it can PM me about access. I'm trying to get some additional data generated to test it prior to releasing & publishing it -- the problem is a dearth of data in which the actual repeat length is known as tested by conventional methods. I know we looked at 1K genome data & thought she was able to pull variation out, but I'd need to check.
krobison is offline   Reply With Quote
Reply

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 02:25 AM.


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