Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
I wrote a little Perl script to generate an mlRho profile file from samtools mpileup. You can find pileup2pro.pl at https://github.com/douglasgscofield/...master/scripts. Bernhard also links to it from the mlRho site.Last edited by dgscofield; 08-07-2012, 08:42 AM.
-
(I would like to obtain the number of A,C,G and T at each position of each contig from a SAM/BAM file.)
I think GATK tools does that. I am also using it to calculate the breakdown of each nucleotide at every position matching reference.
Leave a comment:
-
The author Bernhard Haubold created a script to go from SAMTOOLS pileup to number of A,T,C,G available on asking.
Thanks for your help mgogol anyway
Leave a comment:
-
One way I know of is to run samtools mpileup, but you'll have to write another script to interpret and format the output.
Leave a comment:
-
From BAM/SAM to number of A,C,G,T at each location
Hi,
I would like to obtain the number of A,C,G and T at each position of each contig from a SAM/BAM file.
This is tu run the program mlRho from Bernhard Haubold, Peter Pfaffelhuber and Michael Lynch (http://guanine.evolbio.mpg.de/cgi-bi...o/mlRho.cgi.pl).
The program can also use ACE format but I don't know how to obtain ACE from BAM.
I don't know if there already is some tool available to do that, otherwise I will have to code it myself.
Many thanks in advance for your helpTags: None
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, Today, 08:47 AM
|
0 responses
12 views
0 likes
|
Last Post
by seqadmin
Today, 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
59 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
|
Leave a comment: