SeqMonk
Perhaps Seqmonk would be of use. SeqMonk is a program to enable the visualisation and analysis of mapped sequence data. The progam allows the user to divide the genome into sliding windows and then quantify each region.
Seqmonk homepage
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Originally posted by mathew View PostI am looking for a tool which can allow me to divide a bacterial genome into different sliding windows and count the average of read count. Am I have to go only with custom script or is there any software out there which can allow me to do that.
Thanks
This snippet will divide your genome file (format: chromosome-name<tab>chr size) in windows of size 1000bp by sliding them by 100bp
Code:bedtools makewindows -g mygenome.bed -w 1000 -s 100 -i srcwinnum > mygenome.windows.bed
Code:bedtools coverage -abam myreads.bam -b mygenome.windows.bed > mycov.bed
Good luck!
Dario
Leave a comment:
-
Yikes! Quite an unhelpful response there -- for what is really a pretty trivial program.
Once you have all the reads mapped in a BAM file, then this is easy to do with Perl and the Bio:B::Sam library for Perl (or Python and same or so forth). Probably about a page of Perl (or similar language);
Not speedy, but you can also do it with a series of samtools commands:
samtools view myaln.sam seq1:1000-2000 | wc
samtools view myaln.sam seq1:1500-2500 | wc
The first number of each line will be the number of reads found in that window. Generating this series of commands can even be run entirely on the command line.
Yet another approach would be to generate the pileup & parse that; the read ends are marked in the output and so could be used to do this counting (to depth at middle of window, add all the read start markers found within the window but to the right of the middle plus all the read end markers found within the window but to the left of the middle).
Leave a comment:
-
software for sliding window
I think all scientific fields go hand in hand. I am asking for a pointer if some one came across such tool in any software.
Thanks
Leave a comment:
-
1, If we developers do make such programs, anyone is willing to pay for it?
2, If not, will any journal publish a program like this?
3, If not, a free coffee? Too far away.
Leave a comment:
-
Software for sliding window
I am looking for a tool which can allow me to divide a bacterial genome into different sliding windows and count the average of read count. Am I have to go only with custom script or is there any software out there which can allow me to do that.
Thanks
Tags: None
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, 07-25-2024, 06:46 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
07-25-2024, 06:46 AM
|
||
Started by seqadmin, 07-24-2024, 11:09 AM
|
0 responses
26 views
0 likes
|
Last Post
by seqadmin
07-24-2024, 11:09 AM
|
||
Started by seqadmin, 07-19-2024, 07:20 AM
|
0 responses
160 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: