Hi all,
I know that maybe there's many threads that talked about a pieces of what I'm going to ask, but I'm really in a none issue situation and running out of time with this piece of project while I'm not used to use python.
the situation is the following: I have a set of exons sequences coming from several patients (forward and revers sequences for each exon) in the .abi files.
What I'm asked to do is to write a script to:
1. convert .abi files to fasta+qual or to fastq files.
2. merge this files into a multi fasta or multi fastq file (one file for each exon).
3. rebasecall with lucy or tracetuner or something equivalent in biopython.
4. align everything against the refseq of each exon.
Still have a question: should I put the forward and reverse sequences in the same place ?
Any help or advise is welcomed.
Thanks.
I know that maybe there's many threads that talked about a pieces of what I'm going to ask, but I'm really in a none issue situation and running out of time with this piece of project while I'm not used to use python.
the situation is the following: I have a set of exons sequences coming from several patients (forward and revers sequences for each exon) in the .abi files.
What I'm asked to do is to write a script to:
1. convert .abi files to fasta+qual or to fastq files.
2. merge this files into a multi fasta or multi fastq file (one file for each exon).
3. rebasecall with lucy or tracetuner or something equivalent in biopython.
4. align everything against the refseq of each exon.
Still have a question: should I put the forward and reverse sequences in the same place ?
Any help or advise is welcomed.
Thanks.