Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
I noticed that the referenced Galaxy code does not check for the '>' and '<' characters. Can anybody please provide some help how to process those? I.e., I don't really understand what "reference skip" means and would think that such a symbol means that the read does not cover the the respective position. Any help would be appreciated, thanks!
-
Originally posted by nekrut View Post
Best wishes,
Graham
Leave a comment:
-
-
Originally posted by Pepe View PostYou cannot open a file in Pileup format with Bio:B::sam
You should convert your SAM file to BAM format (samtools view) and sort it (samtools sort), maybe index it too? I'm not sure (samtools index).
Then follow the instructions in the CPAN website, which are very, veryhelpful. Basically:
use Bio:B::Sam;
my $sam = Bio:B::Sam->new(-bam => "my_sorted_bam_file.bam",
-fasta => "my_ref.fasta");
my @targets = $sam->seq_ids;
foreach my $chr (@targets){
$sam->pileup($chr,$my_subroutine);
}
where $my_subroutine is a subroutine that you'll need to get to do what you want: snp calling, coverage calculation...
Thanks for your answer. I'd realised that Bio:B:Sam couldn't parse the pileup, so this is why I am asking if anyone has, or knows about, anything that could parse it, before I spent time writing a Samtools parser myself.
I have a 'Parser.pm' module that parses various pileup formats (maq, bowtie, novoalign, etc), but I haven't added one for Samtools format yet.
Cheers,
Graham
Leave a comment:
-
You cannot open a file in Pileup format with Bio:B::sam
You should convert your SAM file to BAM format (samtools view) and sort it (samtools sort), maybe index it too? I'm not sure (samtools index).
Then follow the instructions in the CPAN website, which are very, veryhelpful. Basically:
use Bio:B::Sam;
my $sam = Bio:B::Sam->new(-bam => "my_sorted_bam_file.bam",
-fasta => "my_ref.fasta");
my @targets = $sam->seq_ids;
foreach my $chr (@targets){
$sam->pileup($chr,$my_subroutine);
}
where $my_subroutine is a subroutine that you'll need to get to do what you want: snp calling, coverage calculation...
Leave a comment:
-
Samtools Pileup Parser
Hi,
I'm trying to parse the Samtools Pileup format in Perl.
I've created a paired-end assembly using BWA (aln, sampe) and then used Samtools to create a pileup (using import - sort - index - pileup).
What I want to do is to open up the pileup (which is obviously in the Samtools Pileup format) and iterate over each position of the pileup, examining differences in base-calls.
Does anyone know of anything available, or have and sample code they could post?
Many thanks.
Latest Articles
Collapse
-
by seqadmin
The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...-
Channel: Articles
07-08-2024, 03:19 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 06:46 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
Yesterday, 06:46 AM
|
||
Started by seqadmin, 07-24-2024, 11:09 AM
|
0 responses
25 views
0 likes
|
Last Post
by seqadmin
07-24-2024, 11:09 AM
|
||
Started by seqadmin, 07-19-2024, 07:20 AM
|
0 responses
159 views
0 likes
|
Last Post
by seqadmin
07-19-2024, 07:20 AM
|
||
Started by seqadmin, 07-16-2024, 05:49 AM
|
0 responses
127 views
0 likes
|
Last Post
by seqadmin
07-16-2024, 05:49 AM
|
Leave a comment: