SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Perl or Python ETHANol Bioinformatics 19 08-29-2019 05:30 AM
running GATK Queue script under PBS environment? orionzhou Bioinformatics 1 12-09-2012 05:24 AM
DEXSeq python script error message shilez Bioinformatics 4 09-08-2011 11:55 PM
python dror Bioinformatics 0 11-29-2010 04:22 AM
python in ABySS dror RNA Sequencing 0 11-28-2010 06:19 AM

Reply
 
Thread Tools
Old 09-07-2011, 03:57 AM   #1
m_elena_bioinfo
Member
 
Location: Ospedali Riuniti di Bergamo, ITALY

Join Date: Oct 2009
Posts: 99
Default gatk python script

Hi NGS users,
anyone could say me where i can find the python script of GATK (ConvertTableToAnnotatorRod.py, GenerateTranscriptToInfo.py...)?

Thanks a lot,
ME
m_elena_bioinfo is offline   Reply With Quote
Old 09-07-2011, 05:34 AM   #2
raonyguimaraes
Member
 
Location: Belo Horizonte - Brazil

Join Date: Jun 2010
Posts: 38
Default

Yes I know where you can get, the problem is that it's not public yet

https://svn.broadinstitute.org/Sting...criptToInfo.py
raonyguimaraes is offline   Reply With Quote
Old 09-07-2011, 05:38 AM   #3
m_elena_bioinfo
Member
 
Location: Ospedali Riuniti di Bergamo, ITALY

Join Date: Oct 2009
Posts: 99
Default

thanks a lot.
I have to generate rod file for depth coverage in all exon of all genes in a whole-exome bam file (according to http://www.broadinstitute.org/gsa/wi...It_Was_Created instruction).
Any idea about fix this problem?
m_elena_bioinfo is offline   Reply With Quote
Old 09-07-2011, 05:56 AM   #4
raonyguimaraes
Member
 
Location: Belo Horizonte - Brazil

Join Date: Jun 2010
Posts: 38
Default

So I guess we are trying to do the same analysis

I think, you don't need rod files anymore, you should use VCF files instead

To calculate the depth of coverage:

http://www.broadinstitute.org/gsa/wi..._Coverage_v3.0

You will need:

-R /path/to/your/reference.fasta
-I /path/to/your/bam_file.bam
-o /path/to/your/output_file

And the intervals ...

-L /path/to/your/intervals.interval_list
raonyguimaraes is offline   Reply With Quote
Old 09-07-2011, 05:59 AM   #5
m_elena_bioinfo
Member
 
Location: Ospedali Riuniti di Bergamo, ITALY

Join Date: Oct 2009
Posts: 99
Default

You are right, we are performing the same analysis.
Ok, I'm already using the depth coverage v3.
As interval list what do you use?
I have a whole exome. Can I use the bed file?
How can i select the specific gene for calculate the coverage in each single exon?

Thanks a lot again!!!!
:-)
m_elena_bioinfo is offline   Reply With Quote
Old 09-07-2011, 07:06 AM   #6
svl
Member
 
Location: Netherlands

Join Date: Sep 2009
Posts: 43
Default

I'm still looking at an efficient way to do this too. I use samtools view and perl to create a table for each sample (.bam) with for each region (so in this case exon coordinates):

- average coverage in region
- percentage bases covered by at least 1
- percentage bases covered by at least 5
- etc

But still searching for a way to store this efficiently while being able to compare it between any sample (whether whole-exome or any other capture design). If you do it for each exon, there's around 600.000 of them to calculate the stats for..
svl is offline   Reply With Quote
Old 09-16-2011, 06:40 AM   #7
Seq84
Member
 
Location: Italy

Join Date: Feb 2011
Posts: 19
Default

Quote:
Originally Posted by m_elena_bioinfo View Post
You are right, we are performing the same analysis.
Ok, I'm already using the depth coverage v3.
As interval list what do you use?
I have a whole exome. Can I use the bed file?
How can i select the specific gene for calculate the coverage in each single exon?

Thanks a lot again!!!!
:-)
Yes sure you can use the .bed file. You should have a .bed file with all targeted regions (exons regions). Set in depthofcoverage the option -L exons.bed.
If you want to have a specific coverage (of exons) for each gene you can provide an interval list with custom coordinates of the exons of your gene.

Instead if you want only the coverage of all genes use the option -genelist.
Seq84 is offline   Reply With Quote
Old 11-20-2011, 07:40 AM   #8
frymor
Senior Member
 
Location: Germany

Join Date: May 2010
Posts: 149
Question

Quote:
Originally Posted by Seq84 View Post
If you want to have a specific coverage (of exons) for each gene you can provide an interval list with custom coordinates of the exons of your gene.
Hi,

is it also possible to use the GRanges object from the R package GenomicRanges as an interval list of exons/gernes?

Thanks

A.
frymor 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 03:03 AM.


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