Mercurial > repos > charles_s_test > seqsero2

comparison libs/BWA_analysis_O_new_dependent.py @ 2:0d65b71ff8df draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 464b391afaa5819bc681452e85bea9d882730eb6
author charles_s_test
date Sun, 12 Nov 2017 02:27:32 -0500
parents 6895de35a263
children 38ad1130d077
comparison
equal deleted inserted replaced
1:8cd7fc65c3a7 2:0d65b71ff8df
5 from Bio import SeqIO 5 from Bio import SeqIO
6 import sys 6 import sys
7 from Initial_functions import Uniq 7 from Initial_functions import Uniq
8 from Bio.Blast import NCBIXML 8 from Bio.Blast import NCBIXML
9 9
10 BwaPath="/nfs/sw/apps/bwa/bwa-0.7.15/bwa" 10 BwaPath="bwa"
11 SamTlsPth="/nfs/sw/apps/samtools/samtools-1.3.1/bin/samtools" 11 SamTlsPth="samtools"
12 Makebltdb="/nfs/sw/apps/blast/ncbi-blast-2.6.0+/bin/makeblastdb" 12 Makebltdb="/nfs/sw/apps/blast/ncbi-blast-2.6.0+/bin/makeblastdb"
13 Blastnpth="/nfs/sw/apps/blast/ncbi-blast-2.6.0+/bin/blastn" 13 Blastnpth="/nfs/sw/apps/blast/ncbi-blast-2.6.0+/bin/blastn"
14 14
15 def BWA_O_analysis(sra_name,additional_file,database,mapping_mode,file_mode): 15 def BWA_O_analysis(sra_name,additional_file,database,mapping_mode,file_mode):
16 if file_mode=="1":#interleaved 16 if file_mode=="1":#interleaved