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.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
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.
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
-
by seqadmin
Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...-
Channel: Articles
04-04-2024, 04:25 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 08:47 AM
|
0 responses
12 views
0 likes
|
Last Post
by seqadmin
Yesterday, 08:47 AM
|
||
Started by seqadmin, 04-11-2024, 12:08 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
60 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
54 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
Comment