User: Seq225
Seq225 • 90
- Reputation:
- 90
- Status:
- Trusted
- Location:
- Last seen:
- 1 week, 4 days ago
- Joined:
- 2 years, 7 months ago
- Email:
- m***********@gmail.com
Posts by Seq225
<prev
• 66 results •
page 1 of 7 •
next >
0
votes
1
answer
187
views
1
answers
Comment:
C: How to explain the blastx output?
... Thank you so much!!!!!!!!!!
It's now working!! ...
written 11 days ago by
Seq225 • 90
0
votes
1
answer
187
views
1
answers
Comment:
C: How to explain the blastx output?
... Wow! Great. Thanks a ton. Almost there!!
I only know pieces about bash. So, confuse things.....
What I am using now is:
for i in seq_*
do
for j in *.dmnd
do
diamond blastp -p 30 --outfmt 6 --db $j --query $i -e 0.00001 -o ${i}_${j}.txt -k 5 --sensitiv ...
written 11 days ago by
Seq225 • 90
0
votes
1
answer
187
views
1
answers
Comment:
C: How to explain the blastx output?
... Thanks. I got the echo part.
I have been using diamond blastp (my db files are huge). However, I can now switch to regular blastp.
Thanks a ton for pointing out the e-value cut-off.
But, when I run the above loop, the output files are empty. ...
written 11 days ago by
Seq225 • 90
0
votes
1
answer
187
views
1
answers
Comment:
C: How to explain the blastx output?
... Thanks. I think I am putting the two variables "i", and "j" in a wrong way. I have added echo in my commands:
for i in seq_*
do
for j in *.dmnd
do
echo diamond blastp -p 30 -v --outfmt 6 --db echo "$j" --query echo "$i" -e 5 -o "$i".txt -k 5 --sen ...
written 12 days ago by
Seq225 • 90
0
votes
1
answer
187
views
1
answers
Comment:
C: How to explain the blastx output?
... Yes, it works!
However, my problem has an extra step. I have multiple queries (1067 to be exact). I want to run them in a loop. Like this:
**
for i in seq_*
do
for j in *.dmnd
do
diamond blastp -p 30 -v --outfmt 6 --db "$j" --query "$i" -e 5 -o "$i".txt ...
written 12 days ago by
Seq225 • 90
0
votes
1
answer
187
views
1
answers
Comment:
C: How to explain the blastx output?
... I have another question:
I want to run blastp against multiple db, and output 5 best hits (-max_target_seqs 5) from each of them (each db hits) . Suppose, I am running blastp against 5 db; I want all 25 hits into one output file.
I see that people use blastdb_aliastool. I just want your opinion. ...
written 12 days ago by
Seq225 • 90
1
vote
1
answer
187
views
1
answers
Comment:
C: How to explain the blastx output?
... Ya!!!!!
Thanks. Great that I can run command line to get all the sequences. I actually need to write a script with a loop to get sequences against my individual queries.
Thanks!! ...
written 15 days ago by
Seq225 • 90
0
votes
1
answer
187
views
1
answers
Comment:
C: How to explain the blastx output?
... Thank you very much!
I have used another program to get all the ORFs. Is have solved my concerns.
Now, I need something similar to qseq. Can I get the entire sequence of the subject? I want the fasta sequences of the top 5 hits for each of my queries. I do not see any options fro that on blast hel ...
written 15 days ago by
Seq225 • 90
0
votes
0
answers
156
views
0
answers
... Wow! Great!! Thanks a ton!!! ...
written 23 days ago by
Seq225 • 90
1
vote
0
answers
156
views
0
answers
... Hi,
I am doing horizontal gene transfer (HGT). My pipeline is based on blast bit-score. So far, I have used blastx (query is a set of my assembled transcripts). Now, I am thinking to use blastp (query is the set of all the translated-ORF from my assembled transcripts).
I am not sure if blastx and ...
written 23 days ago by
Seq225 • 90
Latest awards to Seq225
Popular Question
3 months ago,
created a question with more than 1,000 views.
For Bismark Methyl-Seq Analysis
Popular Question
16 months ago,
created a question with more than 1,000 views.
For HOMER annotatePeaks.pl problem
Popular Question
16 months ago,
created a question with more than 1,000 views.
For De Novo Assembly of small RNA Reads (Illumina)
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 734 users visited in the last hour