This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Need Help to Identify my mistake in my code

Hey I have to find the FASTA file containing the sequence 1kb upstream of each gene on the X chromosome. Given: a GFF and a genome sequence file in fasta format, The code I used to solve is

from Bio import SeqIO
from Bio.SeqRecord import SeqRecord

fast_lgx_records = []
lgs_record_ids = []
records = list(SeqIO.parse(open("Midterm.fna"), 'fasta'))
for record in records:
    if "LGX" in record.description:
        fast_lgx_records.append(record)
        lgs_record_ids.appendrecord.id)

            positions_to_read = []
              with open("Midterm.gff") as f:
    for line in f:
        if not line.startswith("#"):
            split_line = line.split('\t')
            seq_id = split_line[0]
            feature_type = split_line[2]
            start = split_line[3]
            end = split_line[4]
            sign = split_line[6]
            if seq_id in lgs_record_ids and feature_type == "gene":
                if sign == "+":
                    start_index = int(start)
                    positions_to_read.append((start_index - 1001, start_index - 1))
                else:
                    end_index = int(end)
                    positions_to_read.append((end_index + 1, end_index + 1001))

                   # write the final sequence to a new fasta file
                final_data = []

My professor mentioned that : All the headers are identical so I can't match sequences to their corresponding genes or positions in the genome. There are incorrect sequences, but I can't be sure exactly the problem because I can't tell which gene they are supposed to belong to.

Please help.

gene python

What is in "Midterm.fna"? What is LGX? Do you have to use Python? If this is an assignment, shouldn't you figure this out on your own rather than having someone else work for your degree?

  1. Misterm.fna is a file which has all the sequences of the gene in a fasta form.
  2. LGX represents linkage group X
  3. yes I need to use python
  4. I already have been graded for this , not able to identify my mistake tats why I am seeking help.

0 answers

No answers yet.

Log in to answer this question.