Skip to content

Commit

Permalink
stripping of variant file reinstated and typos fixed
Browse files Browse the repository at this point in the history
  • Loading branch information
rbundschuh authored Jan 30, 2020
1 parent 7afe32b commit 1e2e503
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions SMaSH.py
Original file line number Diff line number Diff line change
Expand Up @@ -28,7 +28,7 @@
#
# SMaSH - Sample Matching using SNPs in Humans
#
# Maximilian Westpha, David Frankhouser, Carmine Sonzone, Peter G. Shields,
# Maximilian Westphal, David Frankhouser, Carmine Sonzone, Peter G. Shields,
# Pearlly Yan, and Ralf Bundschuh
#
# Usability and Python3 update by Alan Hoyle, UNC Lineberger Bioinformatics
Expand Down Expand Up @@ -170,7 +170,7 @@ def write_list(inList, outFile): outFile.write("\t".join( map( str, inList ) ) +
elif bam_comma_list == "ALL" or bam_comma_list == "*" or bam_comma_list == "*.bam" :
print (strftime("[%Y-%m-%d %H:%M:%S]"), "Adding *.bam from current directory...")
bams += glob.glob("*.bam")
elif len(bam_comma_list)>4: # has to be greater than for for the file to be x.bam
elif len(bam_comma_list)>4: # has to be greater than four for the file to be x.bam
bams += bam_comma_list.split(',')

elif len(bam_comma_list)>0:
Expand Down Expand Up @@ -247,7 +247,7 @@ def dd():

for lines in variants:
if lines[0] == "#": continue
cols = lines.split('\t')
cols = lines.strip().split('\t')
chrom = cols[chr_index]
# checks if the bam file's chromosome annotation matches the vcf file annotation
# and edits the chromosome accordingly to be able to properly retrieve reads from the bam file
Expand Down

0 comments on commit 1e2e503

Please sign in to comment.