Question: can not get the alignment result from needle alignment
0
gravatar for Medhat
2.8 years ago by
Medhat7.6k
Poland
Medhat7.6k wrote:

I am using this script to compare to sequence and parse the result:

 

from Bio.Emboss.Applications import NeedleCommandline
from Bio import AlignIO

needle_cline = NeedleCommandline(asequence="a.fa",
                bsequence="b.fa",
                gapopen=10, gapextend=0.5, outfile="needle.txt")

filename= "needle.txt"
format = "emboss"


alignments = list(AlignIO.parse(filename, format))

 

when i try to print the results in the alignments list it gives me only part of the last line 

like that

 

SingleLetterAlphabet() alignment with 2 rows and 2258 columns
TTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAA...TCT ref
TGGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAA...T-- corr

 

any Idea?

alignment sequence • 713 views
ADD COMMENTlink modified 2.8 years ago by Kamil1.8k • written 2.8 years ago by Medhat7.6k
0
gravatar for Kamil
2.8 years ago by
Kamil1.8k
Boston
Kamil1.8k wrote:

Consider checking the examples in the documentation for the MultipleSeqAlignment class returned by AlignIO.

I guess this might do what you want:

print alignments[:,:]
ADD COMMENTlink written 2.8 years ago by Kamil1.8k
1

I figured it I need to add another loop to get the sequence

 

for alignment in AlignIO.parse(filename, format):
  for sequence in alignment:
   print(sequence.seq)

ADD REPLYlink modified 2.8 years ago • written 2.8 years ago by Medhat7.6k

I changed the method from parse to read like this 

alignments = AlignIO.read(filename, format)
print alignments[:,:]

 

but, the result is the same 

SingleLetterAlphabet() alignment with 2 rows and 2258 columns
TTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAA...TCT ref
TGGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAA...T-- corr

any Idea?

ADD REPLYlink written 2.8 years ago by Medhat7.6k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1736 users visited in the last hour