Skip to content
Snippets Groups Projects
Commit 56050c50 authored by Rohmer Coralie's avatar Rohmer Coralie
Browse files

Fix tool_consensus.py

parent c9b7e8a9
No related branches found
No related tags found
No related merge requests found
......@@ -33,7 +33,7 @@ for key in nuc_to_iupac.keys():
def consensus_iupac(alignment,threshold):
summary_align = AlignInfo.SummaryInfo(alignment)
nb_reads = len(summary_align.get_column(0))
align_size=len(alignment[1].seq)
align_size=len(alignment[0].seq)
letters=summary_align._get_all_letters()
seq=""
#Run the alignment column by column
......@@ -103,7 +103,7 @@ def consensus_iupac(alignment,threshold):
def consensus(alignment,threshold):
summary_align = AlignInfo.SummaryInfo(alignment)
nb_reads = len(summary_align.get_column(0))
align_size = len(alignment[1].seq)
align_size = len(alignment[0].seq)
letters = summary_align._get_all_letters()
letters = letters.replace('-','')
seq=""
......@@ -149,7 +149,7 @@ def majority_consensus(alignment):
threshold=0.5
summary_align = AlignInfo.SummaryInfo(alignment)
nb_reads = len(summary_align.get_column(0))
align_size = len(alignment[1].seq)
align_size = len(alignment[0].seq)
seq=""
#Run the alignment column by column
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment