SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
How to use Picard's MarkDuplicates cliff Bioinformatics 12 01-26-2015 10:56 PM
Picard, Markduplicates agraz Bioinformatics 0 11-18-2014 05:28 AM
picard MarkDuplicates READ_NAME_REGEX Lv Ray Bioinformatics 5 07-16-2014 07:30 PM
MarkDuplicates in picard bair Bioinformatics 3 12-23-2010 11:00 AM
Picard MarkDuplicates wangzkai Bioinformatics 2 05-18-2010 09:14 PM

Reply
 
Thread Tools
Old 04-02-2018, 11:02 PM   #1
sukesh14
Junior Member
 
Location: Chennai

Join Date: Sep 2016
Posts: 2
Default Reg Picard MarkDuplicates

Hi

I am beginner in python. I am trying to automate the process of variant calling for the whole genome sequencing data for plants in ubuntu. We wrote commands for each software (by giving path) and integrated all software's in python. I am having a problem when i execute the script, where the script file could not run picard tools Mark Duplicates. It did not give any error but process wont continue after that, whereas when i use same tool, picard markduplicates in terminal its working fine without any error. Can anyone help me in trouble shooting this.

Example script that we wrote

import subprocess,os,commands

m=raw_input("Enter the project destination")

os.environ["JAVA_HOME"]="/usr/java/jdk1.8.0_121/lib/visualvm/platform/modules"

os.chdir('/home/sukesh/picard-tools')

os.system('java -jar picard.jar MarkDuplicates I='+m+'/6SAMTOBAM/sortIRGSP.bam
O='+m+'/7MARKDUPS/dups.bam REMOVE_DUPLICATES=false METRICS_FILE=metricN.log ASSUME_SORTED=True VALIDATION_STRINGENCY=LENIENT')

os.chdir('/home/Downloads')

Or is there any other tool which can mark duplicates other than samtools which will be more helpful?
sukesh14 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 06:34 PM.


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