I have a program that runs through a list of Amino Acid sequences for a protein in Influenza A that I'm analyzing and finds amino acid at position 627 in each amino acid sequence in my FASTA file for each protein.
My code works like this
with open(file, "r" ) as source:
for heading_and_lines in group_by_heading( source ):
heading= heading_and_lines[0]
lines= heading_and_lines[1:]
lines = ''.join(lines)
if lines[627-1] == 'K':
print "---------------MUTATION BELOW--------------"
print heading
print lines[627-1]
#print "-------------------------------------------"
print "end of file"
But my code does not work like this
with open(file, "r" ) as source:
for heading_and_lines in group_by_heading( source ):
heading= heading_and_lines[0]
lines= heading_and_lines[1:]
lines = ''.join(lines)
if lines[627-1] == 'K':
print "---------------MUTATION BELOW--------------"
print heading
print lines[627-1]
print "-------------------------------------------"
print "end of file"
For some reason printing the line below the mutation that I found returns an error, and unexpected indent. See I want the line beneath each mutation, just like the line on top, so I can clearly organize the mutations. Any reason why you think the printing statement below doesn't work?
Sorry if this sounded confusion, and as always thanks for your time.