Moderator: Rob

gravatar for Rob
Rob2.8k
Reputation:
2,840
Status:
Trusted
Location:
United States
Twitter:
@nomad421
Last seen:
4 hours ago
Joined:
7 years, 10 months ago
Email:
r********@gmail.com

Posts by Rob

<prev • 205 results • page 1 of 21 • next >
3
votes
2
answers
158
views
2
answers
Answer: A: Average Insert Size When Using Transcript Quantification Tools Such as Salmon?
... Is the data paired-end? In this case, salmon estimated the fragment length distribution. In the quantification directory the file `libParams/flenDist.txt` contains the parameters of the normalized distribution of observed fragment lengths. Each of the 1001 numbers gives the probability of observi ...
written 26 days ago by Rob2.8k
1
vote
1
answer
235
views
1
answers
Comment: C: best pipeline to get gene and transcript counts
... You can reduce the memory required by star by building the index with the `--genomeSAsparseD` parameter set to a value greater than 1. This samples the suffix array. If the value is too large, STAR will become slower, but small values (e.g. 2/4) can provide a nice memory savings without slowing th ...
written 28 days ago by Rob2.8k
1
vote
1
answer
235
views
1
answers
Comment: C: best pipeline to get gene and transcript counts
... No. STAR has the ability to perform a single alignment against genome (and gtf), and to provide 2 separate outputs using the same alignments. One BAM file contains the spliced alignments to the genome. The other BAM file contains the unspliced alignments to the transcriptome (suitable to provide ...
written 28 days ago by Rob2.8k
1
vote
1
answer
235
views
1
answers
Comment: C: best pipeline to get gene and transcript counts
... In this case, as I suggested, the easiest path is to align with STAR and keep both the genome alignments and get the transcriptome alignments. As swbarnes2 mentions below, you can even use STARs `--quantMode` for the genome-level counts. Note, however, that you won't be able to get transcript leve ...
written 29 days ago by Rob2.8k
0
votes
1
answer
235
views
1
answers
Comment: C: best pipeline to get gene and transcript counts
... I see. In that case, perhaps consider aligning with STAR, which can give you both the genome and transcriptome alignments at the same time. The genome alignments can be counted using featureCounts, and the transcriptome alignments can be fed to salmon for transcript-level quantification. Also, it ...
written 29 days ago by Rob2.8k
0
votes
1
answer
235
views
1
answers
Comment: C: best pipeline to get gene and transcript counts
... Did you align against the genome or transcriptome? If you want to use a traditional aligner to provide input to salmon (rather than using its builtin mapping algorithm), then you should either (1) align to the genome with STAR and tell it to project the alignments to the transcriptome (it can do th ...
written 4 weeks ago by Rob2.8k
1
vote
2
answers
339
views
2
answers
Comment: C: Salmon (error: SAM file and fasta file have different transcript length)
... The effect of the `--geneMap` flag is to cause salmon to output, in addition to it's transcript level quantifications, aggregated gene-level quantifications. However, this flag exists only for backward compatibility and is deprecated in favor of tximport. Though it attempts to do the same thing, t ...
written 4 weeks ago by Rob2.8k
1
vote
2
answers
339
views
2
answers
Comment: C: Salmon (error: SAM file and fasta file have different transcript length)
... I'd recommend going with the quasi-mapping, and including the `--validateMappings` option among your command line flags (including others that seem relevant for your data given the documentation). ...
written 5 weeks ago by Rob2.8k
1
vote
2
answers
339
views
2
answers
Comment: C: Salmon (error: SAM file and fasta file have different transcript length)
... As you point out, whether or not you use traditional alignment pr the builtin mapping algorithm, you won't be finding novel isoforms. If you want to do that, you should perform assembly via a tool like Stringtie or Scallop prior to quantification. Regardless of if you do that, if you want to use ...
written 5 weeks ago by Rob2.8k
1
vote
2
answers
339
views
2
answers
Comment: C: Salmon (error: SAM file and fasta file have different transcript length)
... Wait; did you align against the transcriptome using TopHat2? I would highly recommend against doing that. In alignment-based mode, salmon requires alignments to the transcriptome. The best ways to obtain such alignments are either to (1) align to the genome using STAR and use it's built-in facilit ...
written 5 weeks ago by Rob2.8k

Latest awards to Rob

Appreciated 4 days ago, created a post with more than 5 votes. For A: How to handle Reverse complement in De-Bruijn Graph
Teacher 16 days ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Scholar 16 days ago, created an answer that has been accepted. For A: Can pairwise global alignment be done in linear memory with affine gap penalties
Good Answer 19 days 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 5 weeks 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 weeks ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Appreciated 5 weeks ago, created a post with more than 5 votes. For A: How to handle Reverse complement in De-Bruijn Graph
Appreciated 5 weeks ago, created a post with more than 5 votes. For A: How to handle Reverse complement in De-Bruijn Graph
Teacher 8 weeks ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Scholar 9 weeks ago, created an answer that has been accepted. For A: Estimating RPKM or TPM in RNA-Seq data
Teacher 9 weeks ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Teacher 10 weeks ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Teacher 3 months ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Teacher 3 months ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Scholar 3 months ago, created an answer that has been accepted. For A: Estimating RPKM or TPM in RNA-Seq data
Teacher 3 months ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Commentator 3 months ago, created a comment with at least 3 up-votes. For C: best workflow for RNA data analysis without replications
Popular Question 4 months ago, created a question with more than 1,000 views. For "Un-Projecting" genome alignments onto a transcriptome
Scholar 6 months ago, created an answer that has been accepted. For A: Estimating RPKM or TPM in RNA-Seq data
Teacher 6 months ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Scholar 6 months ago, created an answer that has been accepted. For A: Estimating RPKM or TPM in RNA-Seq data
Teacher 7 months ago, created an answer with at least 3 up-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
Teacher 8 months ago, created an answer with at least 3 up-votes. For A: How to handle Reverse complement in De-Bruijn Graph
Student 10 months ago, asked a question with at least 3 up-votes. For Fragment Size: TLEN vs. isize

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1726 users visited in the last hour