Skip to content
Snippets Groups Projects
Commit 18259c64 authored by nina.marthe_ird.fr's avatar nina.marthe_ird.fr
Browse files

modified output form

parent af142979
No related branches found
No related tags found
No related merge requests found
...@@ -259,11 +259,12 @@ for feature in Features.values(): # add the sequence of all features ...@@ -259,11 +259,12 @@ for feature in Features.values(): # add the sequence of all features
# analysing the variations for all the mrna : # analysing the variations for all the mrna :
for mrna_id in mrna_var.keys(): for mrna_id in mrna_var.keys():
print("analysis of the variations in the mRNA",mrna_id,":\n")
cds_var=mrna_var[mrna_id] cds_var=mrna_var[mrna_id]
for cds_id in cds_var: for cds_id in cds_var:
cds=Features[cds_id] cds=Features[cds_id]
print("analysis of the variations in the CDS",cds_id,":\n") print("\nCDS",cds_id,":")
frameshift=0 frameshift=0
for index, var in enumerate(cds_var[cds_id]): # for each variation in the current cds : for index, var in enumerate(cds_var[cds_id]): # for each variation in the current cds :
type_var=var[8] type_var=var[8]
...@@ -280,7 +281,7 @@ for mrna_id in mrna_var.keys(): ...@@ -280,7 +281,7 @@ for mrna_id in mrna_var.keys():
else: else:
length_alt=len(var[10]) length_alt=len(var[10])
print("variation",index+1, ":") print("\nvariation",index+1, ":")
if posVar[0]<3: if posVar[0]<3:
print("variation of the start codon, mRNA most likely wont be translated") print("variation of the start codon, mRNA most likely wont be translated")
...@@ -354,5 +355,3 @@ for mrna_id in mrna_var.keys(): ...@@ -354,5 +355,3 @@ for mrna_id in mrna_var.keys():
# possible that it prints too many variations : for ex if we have a snp on the first and the last base of a codon, # possible that it prints too many variations : for ex if we have a snp on the first and the last base of a codon,
# while printing the effect of the first snp we aso use the second one. # while printing the effect of the first snp we aso use the second one.
print("\n")
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment