Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In

Moderator: Rob

gravatar for Rob
Rob ♦ 3.4k
Reputation:
3,360
Status:
Trusted
Location:
United States
Twitter:
@nomad421
Last seen:
34 minutes ago
Joined:
8 years, 6 months ago
Email:
r********@gmail.com

Posts by Rob

<prev • 221 results • page 1 of 23 • next >
3
votes
1
answer
174
views
1
answers
Answer: A: Does salmon use phred scores while pseudoaligning (and/or validating mapping)?
... Hi Anthony, No; salmon does not use the PHRED score, either while mapping or aligning. The main reason for not using it during alignment or mapping is that it just tends not to be very useful for this purpose. If a nucleotide inside a reasonably long exact match is bad quality (and therefore lik ...
written 11 weeks ago by Rob ♦ 3.4k
1
vote
2
answers
259
views
2
answers
Answer: C: do paired end reads both mapped to the same gene appear consecutive to each othe
... You might consider using `--no-discordant` and `--no-mixed` to avoid the complicated patterns that can occur when one alignment record (e.g. for read 1) is paied with multiple alignments (for read 2, or vice versa). Basically, for concordant alignments, Bowtie2 will fillow your expected format. Wi ...
written 3 months ago by Rob ♦ 3.4k
5
votes
1
answer
337
views
1
answers
Answer: C: Salmon (quant.sf) to txtimport
... Hi Morris, In general, you should *not* modify the structure of the salmon directories as tximport relies on it. The easiest thing to do is to have a separate quant directory for each sample, and to provide those each to tximport using whatever naming convention you want for the quant *folders*. ...
written 4 months ago by Rob ♦ 3.4k
0
votes
1
answer
476
views
1
answers
Comment: C: Understanding salmon codes
... How large are they? You could use a file-sharing website, or if it is easier, you should be able to upload the data to [this gdrive link](https://drive.google.com/drive/folders/17N35X4_-6yuR4JMfLpbcreDZhQow1jpM?usp=sharing). ...
written 4 months ago by Rob ♦ 3.4k
0
votes
1
answer
476
views
1
answers
Comment: C: Understanding salmon codes
... That's a very surprising difference if the annotations are really the same. Would you be able to share the reads from one of your samples? We'd be interested to take a look. ...
written 4 months ago by Rob ♦ 3.4k
1
vote
1
answer
476
views
1
answers
Comment: C: Understanding salmon codes
... Hi Morris, Can I ask how you verified that these are "human"? I took the sequences you gave above, and searched against the gencode transcriptome using Bowtie2 as the aligner. Both of these reads came back as completely unmapped. I did get a hit with blast (https://www.ncbi.nlm.nih.gov/nucleot ...
written 4 months ago by Rob ♦ 3.4k
1
vote
1
answer
476
views
1
answers
Comment: C: Understanding salmon codes
... Hi Morris, It looks like you have a big space between the `--` and the flag name that should follow immediately after. That is, your command line should be: ``` salmon quant -i gencode_v29_idx -l A -1 Treated_1_m1.fastq -2 Treated_1_m2.fastq -p 4 --validateMappings --writeUnmappedNames -o salmo ...
written 4 months ago by Rob ♦ 3.4k
1
vote
1
answer
476
views
1
answers
Comment: C: Understanding salmon codes
... Actually, salmon can accept gzipped FASTA files for reference index creation as well as gzipped FASTA/Q files for quantification. I believe that Morris was able to build the index, but is running into problems while quantifying. ...
written 4 months ago by Rob ♦ 3.4k
2
votes
1
answer
476
views
1
answers
Answer: C: Understanding salmon codes
... Hi Morris, You should add the index *directory* as the parameter to `-i`, so, try: ``` (salmon) [@ws7910 RNAseq]$ salmon quant -i gencode_v29_idx -l A -1 Treated_1_m1.fastq -2 Treated_1_m2.fastq -p 4 --validateMappings -o salm_treated1 ``` I've also removed a spurious `/` before `gencode_v29_idx ...
written 4 months ago by Rob ♦ 3.4k
0
votes
1
answer
437
views
1
answers
Comment: C: Salmon transcriptome generation
... Hi Morris, The `-k` argument has a default value (31) that is used if `-k` is not provided. I you wish to use the default, you don't have to pass that option explicitly. If you want to use another value of k, then you can pass that to the index command. ...
written 4 months ago by Rob ♦ 3.4k

Latest awards to Rob

Teacher 3 days ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Epic Question 10 days ago, created a question with more than 10,000 views. For Fragment Size: TLEN vs. isize
Appreciated 6 weeks ago, created a post with more than 5 votes. For A: How to handle Reverse complement in De-Bruijn Graph
Teacher 11 weeks ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Appreciated 3 months ago, created a post with more than 5 votes. For A: How to handle Reverse complement in De-Bruijn Graph
Appreciated 3 months ago, created a post with more than 5 votes. For A: How to handle Reverse complement in De-Bruijn Graph
Appreciated 4 months ago, created a post with more than 5 votes. For A: How to handle Reverse complement in De-Bruijn Graph
Teacher 4 months ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Teacher 4 months ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Appreciated 5 months ago, created a post with more than 5 votes. For A: How to handle Reverse complement in De-Bruijn Graph
Good Question 5 months ago, asked a question that was upvoted at least 5 times. For Suffix Array Vs. Compressed Suffix Array Vs. Bwt Search
Appreciated 5 months ago, created a post with more than 5 votes. For A: How to handle Reverse complement in De-Bruijn Graph
Teacher 5 months ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Scholar 5 months ago, created an answer that has been accepted. For A: Can pairwise global alignment be done in linear memory with affine gap penalties
Scholar 5 months ago, created an answer that has been accepted. For A: Can pairwise global alignment be done in linear memory with affine gap penalties
Teacher 5 months ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Good Answer 5 months ago, created an answer that was upvoted at least 5 times. For A: Why samtools sort input file is 4.3 Gb but output become 3.5Gb....
Scholar 6 months ago, created an answer that has been accepted. For A: Can pairwise global alignment be done in linear memory with affine gap penalties
Commentator 6 months ago, created a comment with at least 3 up-votes. For C: best workflow for RNA data analysis without replications
Good Answer 6 months ago, created an answer that was upvoted at least 5 times. For A: Why samtools sort input file is 4.3 Gb but output become 3.5Gb....
Teacher 6 months ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Commentator 8 months ago, created a comment with at least 3 up-votes. For C: best workflow for RNA data analysis without replications
Appreciated 8 months ago, created a post with more than 5 votes. For A: How to handle Reverse complement in De-Bruijn Graph
Teacher 8 months ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Scholar 8 months ago, created an answer that has been accepted. For A: Can pairwise global alignment be done in linear memory with affine gap penalties
Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1041 users visited in the last hour