This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Exclude non coding sequences from multifasta

Hello there

I have a multifasta that have sequences with multiple stop codons and/or starts without AUG.

How do I exclude those sequences?

sequence gene
• 3,241 views • link
updated by Ram Ram • written by radha.jg • radha.jg

I've started a python program like this:

use Bio::SeqIO
sequenceFile= open(standarinput,'ro') #con bioseq
outFile =open(standarinput,'rw')
for seq in sequenceFile:
  if startAUG(seq):
    if stopAtTheEnd(seq):
      if NOTstopMiddle(seq):
        outfile.print(seq)

def startAUG(seq):

  #SOMETHING INSIDE HERE


def stopAtTheEnd(seq):

# SOMETHING INSIDE HERE

def if NOTstopMiddle(seq):

# SOMETHING INSIDE HERE

So, any ideas how to do this on python or perl?

• 0 views • link
updated by Ram Ram • written by radha.jg • radha.jg

the sequences are CDSs from tritrypdb and they usually have AUG as start codons. I'd like to keep it simple and can live without rare start codons because it's for use with codonW without sequences that trigger warnings

• 0 views • link
updated by Ram Ram • written by radha.jg • radha.jg

1 answer

Things done so far

from Bio import SeqIO

infile= input('input fasta file')
out = input('output fasta file')
outFile =open(out,'rw')

#Opens output file
outFile =open(standarinput,'rw')

########################   Main Program   ####################################
#opens and parse the fasta file and each sequense is seq_record

for seq_record in SeqIO.parse(infile, "fasta"):

  if startAUG(seq_record.seq):
    if stopAtTheEnd(seq_record.seq):
      if NOTstopMiddle(seq_record.seq):
        outfile.print(">",seq_record.id"\n")
        outfile.print(seq_record.seq,"\n")

#close exit file
close outFile

######################## End Main program  ####################################

######################     Functions       #######################################

########this one checks if there's AUG at the begining
def startAUG(seq):
  if seq[0]=='A' and seq[1]=='U' seq[2]=='G'
    return true
  else
    return false
#############################################

######## this one checks if there's stop codon at the end
def stopAtTheEnd(seq):
  if (seq[len(seq)-1]=='U' and seq[len(seq)-2]=='A' and seq[len(seq)-3]=='A') or (seq[len(seq)-1]=='U' and seq[len(seq)-2]=='A' and seq[len(seq)-3]=='G') or (seq[len(seq)-1]=='U' and seq[len(seq)-2]=='G' and seq[len(seq)-3]=='A')
    return true
  else
    return false
#############################################

######## stop codons in the middle
def if NOTstopMiddle(seq):

##############################################

Please I need some help here

• 0 views • link
updated by Ram Ram • written by radha.jg • radha.jg

Log in to answer this question.

More posts like this