Python Forum

HELLO. I AM TRYING TO COUNT G'S AND C'S IN A LIST OF STRINGS AND HAVE THE LOOP PUT OUT THE GC COUNT AND THEN CALCULATE A GC PERCENTAGE. I HAVE TRIED VARIOUS CONFIGURATIONS OF A NESTED LOOP AS WELL AS TRYING TO FIGURE OUT HOW TO SOLVE IT USING A FUNCTION, AND AM JUST STUCK. ANY HELP WOULD BE GREATLY APPRECIATED. I AM STUCK ON PART 2/3.

print("\n########## PART 3 ##########\n")

#PART 3 -  Iterate through a list of strings, calculate value, write to file

#(1) Iterate through the list of sequences
#(2) Count G's and C's for each sequence number in a list.
#(3) Calculate the GC Percentage for the sequence nmer.
#(4) Write the sequence and its GC percentage (separated by a "Tab")
#    to a file named:
#       gc_output.txt

#For example if given the following sequence list
#seq_list=["AAUG","AGCG","GCUA"]

#The loop would write the following to the file gc_output.txt
#AAUG    25.0
#AGCG    75.0
#GCUA    50.0

#== PART 3 ==

#seq_list=["AAUG","AGCG","GCUA"]
#seq_list=["AAUGGAGACU","AGCGCCCC","GCUAAAAAAU"]

#(1)
seqs=["AAUG","AGCG","GCUA"]
for element in seqs:
    print(element)

#(2)
gc_count=0
total=len(seqs)
for seq in seqs:
    for x in seq:
        if x == "G" or x == "C":
            gc_count+=1
            print(gc_count)

#c_count=0
#seq=["AAUG"]
#for x in seq:
#    if x =="G" or x =="C":
#        gc_count+=1
#        print(gc_count)

#def gc_count(string)

#Doing a for loop for one of the three sequences
seq="AGGC"
gc_count=0
for b in seq:
    if b=="G" or b=="C":
        gc_count+=1
    else:
        gc_count+=0
print(gc_count)

gc_count=0
seq_list=["AAUG", "AGCG", "GCUA"]
def gc_calc(seq_list):
    for b in seq_list:
        if re.search(g,c):
            gc_count+=1
        else:
            pass
    return gc_count

            
#Lab 04 Part 3 Work
#(3)
#gc_count=0
#seqs=["AAUG","AGCG","GCUA"]
#total=len(seqs)
#for seq in seqs:
#    for x in seq:
 #       if x == "G" or x == "C":
#            gc_count+=1
#            gc_percentage = gc_count/float(total)
#            print("gc_percentage", gc_percentage)
    
#gc_percentage = gc_count/float(total)
#print("gc_percentage", gc_percentage)


#g=0
#c=0
#t=0
#a=0
#gc_count=0
#a_count=0
#t_count=0
#seq_list=["AAUG","AGCG","GCUA"]
#total=len(seq_list)
#print("total = ", total) #this refers to how many strings are in the list
#for seq in seq_list:
#    if i == "G" or i == "C":
#        gc_count+=1
#        print("gc=", gc_count)
#    elif i == "A":
#        a_count+=1
#        print(a_count)
#    elif i == "T":
#        t_count+=1
#        print(t_count)

#gc_count=0
#seq=["AAUG"]
#for char in seq:
 #   if i == "G" or i =="C":
 #       gc_count+=1
 #       print(gc_count)


#for x in seq_list:
#    if "C" in x:
#        c+=1    
#    elif "G" in x:
#        g+=1
#    elif "A" in x:
#        a+=1    
#    elif "T" in x:
#        t+=1

#print ("C=%total, G=%total, A=%total, T=%total" %(c,g,a,t))

#gc_content=(g+c)*100/(a+t+g+c)

#print ("gc_content= %f" %(gc_content))

Please, don't use ALL CAPS! It's yelling