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

Moderator: Damian Kao

gravatar for Damian Kao
Damian Kao ♦ 15k
Reputation:
15,300
Status:
Trusted
Location:
USA
Website:
http://blog.nextgeneti...
Last seen:
1 day, 7 hours ago
Joined:
8 years, 2 months ago
Email:
d*********@gmail.com

Bioinformatician at Janelia Research Campus.

Posts by Damian Kao

<prev • 912 results • page 1 of 92 • next >
0
votes
0
answers
216
views
0
answers
Comment: C: Right MongoDB approach for a multigenomic database?
... My recommendation is to not over-optimize unless this is a personal project for learning purposes. There are three main considerations here all related to scale: 1. How many concurrent users do you expect to be using this database? 2. Is this pre-dominantly a read-only database where users query ...
written 6 months ago by Damian Kao ♦ 15k
0
votes
1
answer
312
views
1
answers
Comment: C: Is there a preferred RNA-seq aligner for (very) short reads and potential homolo
... Try setting --outFilterMatchNmin to 20 to see if you can get more mapping. However, that means it only requires 20 bases to map, which is pretty low. ...
written 8 months ago by Damian Kao ♦ 15k
2
votes
2
answers
406
views
2
answers
Answer: A: Does Gene length corrected TMM [GeTMM] violate any assumptions of TMM normalizat
... Technically, RPK values do not violate assumptions of TMM. TMM is just a technique that tries to find the non-DE portion of the expression distribution by very liberally trimming off outliers. It doesn't matter what kind of expression units you are using. However, RPK values do violate assumptio ...
written 8 months ago by Damian Kao ♦ 15k
3
votes
1
answer
401
views
1
answers
Answer: A: Trimmomatic output file Issue
... It looks like you are not specifying an output file for the -trimlog parameter. So it thinks your `output_forward_paired.fq.gz` is an input. I hope you still can still redownload `sample1_R1_001.fastq.gz`, because it might have been overwritten. ...
written 8 months ago by Damian Kao ♦ 15k
1
vote
8
answers
646
views
8
answers
Answer: A: Code golf: detecting homopolymers of length N in the (human) genome
... Takes in fasta file and a second parameter for homopolymer length. It streams through each line to find homopolymers. Outputs chromosome, start, end, homopolymer base, length of homopolymer. I tested it out on this fasta file: >A AGTCAAAA GGGGTTTTCCCC >B AGTCCCCCTTTTAAA ...
written 8 months ago by Damian Kao ♦ 15k
5
votes
4
answers
608
views
4
answers
Answer: C: Getting the number of SNPs in some ranges
... You have a .csv file with chromosome, start, end coordinates. You can change that into a .bed file pretty easily. Same with your SNP csv file. Convert those csv files to bed and then use bedtools intersect. ...
written 9 months ago by Damian Kao ♦ 15k
0
votes
0
answers
447
views
0
answers
Comment: C: Strange Depth of Coverage distribution
... Yeah that could be it. You can check this by looking at the insert size distribution of your PE reads. See if it is smaller than 2 * average read length. The 9th column of your .sam/.bam should be the insert size. ...
written 9 months ago by Damian Kao ♦ 15k
2
votes
0
answers
523
views
0
answers
Comment: C: Getting read depth for normal and tumour
... I am not sure what program you are using to generate the .vcfs. But you should look into the manual of the program and see if it outputs depths using the format fields in the vcf. For example, you have "DP" field in your vcf that shows the depth of the individual samples. Perhaps one of the other f ...
written 9 months ago by Damian Kao ♦ 15k
0
votes
0
answers
523
views
0
answers
Comment: C: Getting read depth for normal and tumour
... Can you post a couple of lines of your vcf file? ...
written 9 months ago by Damian Kao ♦ 15k
5
votes
1
answer
415
views
1
answers
Answer: A: Samtools segmentation fault when extracting exons from bam file.
... There is a -L flag that you can use with samtools to output alignments overlapping intervals defined by a .bed file. That's probably your best option. You'll have to convert your exons.txt to a bed file. ...
written 9 months ago by Damian Kao ♦ 15k

Latest awards to Damian Kao

Teacher 7 weeks ago, created an answer with at least 3 up-votes. For A: Removing Redundant Amino Acid Sequences From Fasta - *But Also Give The Groups O
Popular Question 12 weeks ago, created a question with more than 1,000 views. For Interpreting Trinity components
Popular Question 4 months ago, created a question with more than 1,000 views. For segment duplication vs heterozygosity
Teacher 5 months ago, created an answer with at least 3 up-votes. For A: Removing Redundant Amino Acid Sequences From Fasta - *But Also Give The Groups O
Teacher 6 months ago, created an answer with at least 3 up-votes. For A: Removing Redundant Amino Acid Sequences From Fasta - *But Also Give The Groups O
Appreciated 7 months ago, created a post with more than 5 votes. For A: A Farewell To Bioinformatics
Teacher 7 months ago, created an answer with at least 3 up-votes. For A: Removing Redundant Amino Acid Sequences From Fasta - *But Also Give The Groups O
Scholar 7 months ago, created an answer that has been accepted. For C: What's the correct way to withdraw a published database/server ?
Scholar 8 months ago, created an answer that has been accepted. For C: What's the correct way to withdraw a published database/server ?
Appreciated 8 months ago, created a post with more than 5 votes. For A: A Farewell To Bioinformatics
Teacher 9 months ago, created an answer with at least 3 up-votes. For A: Removing Redundant Amino Acid Sequences From Fasta - *But Also Give The Groups O
Scholar 9 months ago, created an answer that has been accepted. For C: What's the correct way to withdraw a published database/server ?
Appreciated 9 months ago, created a post with more than 5 votes. For A: A Farewell To Bioinformatics
Teacher 9 months ago, created an answer with at least 3 up-votes. For A: Removing Redundant Amino Acid Sequences From Fasta - *But Also Give The Groups O
Good Answer 11 months ago, created an answer that was upvoted at least 5 times. For A: A Farewell To Bioinformatics
Commentator 12 months ago, created a comment with at least 3 up-votes. For C: Shall We Go Back To Stackexchange?
Popular Question 12 months ago, created a question with more than 1,000 views. For deseq2 user defined size factors
Teacher 13 months ago, created an answer with at least 3 up-votes. For A: Removing Redundant Amino Acid Sequences From Fasta - *But Also Give The Groups O
Popular Question 14 months ago, created a question with more than 1,000 views. For deseq2 user defined size factors
Scholar 17 months ago, created an answer that has been accepted. For C: Does trinity discard reads shorter than 25-mer?
Scholar 17 months ago, created an answer that has been accepted. For C: What's the correct way to withdraw a published database/server ?
Good Answer 17 months ago, created an answer that was upvoted at least 5 times. For A: A Farewell To Bioinformatics
Appreciated 17 months ago, created a post with more than 5 votes. For A: A Farewell To Bioinformatics
Teacher 17 months ago, created an answer with at least 3 up-votes. For A: Removing Redundant Amino Acid Sequences From Fasta - *But Also Give The Groups O
Appreciated 17 months ago, created a post with more than 5 votes. For A: A Farewell To Bioinformatics
Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1991 users visited in the last hour