User: mforthman
mforthman • 30
- Reputation:
- 30
- Status:
- New User
- Location:
- Last seen:
- 4 months ago
- Joined:
- 3 years, 3 months ago
- Email:
- m********@ufl.edu
Posts by mforthman
<prev
• 48 results •
page 1 of 5 •
next >
0
votes
0
answers
586
views
0
answers
... Running fastq-generated SRA files from NCBI through prinseq-lite. The program generates an error:
Use of uninitialized value $qual in scalar chomp at /apps/prinseq/0.20.4/bin/prinseq-lite.pl line 2583, line 26875674.
It still continues on until it has finished, but clearly there is somethin ...
written 16 months ago by
mforthman • 30
1
vote
2
answers
480
views
2
answers
... I have been trying to download SRA data from NCBI and putting it in fastq format using fastq-dump. A colleague and I have been trying to figure out why the resulting fastq files are causing some errors when inputted into prinseq-lite.
My collaborator has been using this fastq-dump command:
fas ...
0
votes
4
answers
1.2k
views
4
answers
... I've still been tinkering with the script and feel as though I might be getting closer to the solution:
#!/usr/bin/env python
import sys
import re
original_fn = sys.argv[1]
company_fn = sys.argv[2]
pattern = '(uce.+$|ENSOFAS.+$|[AB]_[0-9]+$)'
map = {} ...
written 22 months ago by
mforthman • 30
1
vote
4
answers
1.2k
views
4
answers
... You are correct, that was my mistake. Thanks for catching that! And I appreciate you taking a further look at this later. ...
written 22 months ago by
mforthman • 30
0
votes
4
answers
1.2k
views
4
answers
... Correct, I kept the original question simple in original post. In reality, the company.fasta file had many differently formatted headers. I used the script you had provided to process just the `uce` headers, which worked. I already had the `ENSOFAS` headers. This leaves the others to be formatted. T ...
written 22 months ago by
mforthman • 30
0
votes
4
answers
1.2k
views
4
answers
... Thanks for the info. I've modified the regex pattern to also skip over other headers that shouldn't be modified, e.g., like the ones in the original posting (they are all in one file). If I tested correctly, the new pattern is `(uce | ENSOFAS | _[AB]_[0-9]+$)`. Regarding the last part of your reply, ...
written 22 months ago by
mforthman • 30
0
votes
4
answers
1.2k
views
4
answers
... How might this script be modified to deal with other headers that are formatted differently in the company.fasta file? Specifically, I'm targeting the following example headers now:
>Clavigralla_tomentosicollis_gi_512427358_gb_GAJX01007276.1_0_rc
TGAGAAGCTCTCGCAGTCACACTGGCCCACGACCACTCTCA ...
written 22 months ago by
mforthman • 30
0
votes
0
answers
619
views
0
answers
... I had a similar situation posted [here][1], which was resolved, but I'm now dealing with differently formatted headers. I have two DNA sequence files I'm working with. One file (e.g., 'company.fasta') has headers that I need to change based on corresponding headers from another file (e.g., 'original ...
0
votes
4
answers
1.2k
views
4
answers
... PERFECT! Thank you and everyone else for all of your help. ...
written 23 months ago by
mforthman • 30
0
votes
4
answers
1.2k
views
4
answers
... Example original headers from one dataset source:
>uce-2015_p10 |design:hemiptera-v1,designer:faircloth,probes-locus:uce-2015,probes-probe:10,probes-source:oncfas1,probes-global-chromo:Scaffold80,probes-global-start:539805,probes-global-end:539925,probes-local-start:40,probes-local-end:160
...
written 23 months ago by
mforthman • 30
Latest awards to mforthman
Popular Question
16 months ago,
created a question with more than 1,000 views.
For Several exonerate options; which to choose?
Popular Question
16 months ago,
created a question with more than 1,000 views.
For Select set of genes present in all species tested and unique within each genome/transcriptome
Popular Question
16 months ago,
created a question with more than 1,000 views.
For FrameDP installation issue
Popular Question
16 months ago,
created a question with more than 1,000 views.
For How to extract exon sequences from annotated genome
Popular Question
16 months ago,
created a question with more than 1,000 views.
For Modifying perl script to not print redundant sequences and to fix length of sequences outputted from blast
Popular Question
16 months ago,
created a question with more than 1,000 views.
For search text in one file and then replace with text from another file
Popular Question
3.1 years ago,
created a question with more than 1,000 views.
For How to extract exon sequences from annotated genome
Popular Question
3.1 years ago,
created a question with more than 1,000 views.
For Reciprocal Best Hits across more than 2 species
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 1151 users visited in the last hour