User: Illinu

gravatar for Illinu
Illinu90
Reputation:
90
Status:
Trusted
Location:
Belgium
Last seen:
5 months, 2 weeks ago
Joined:
2 years, 12 months ago
Email:
s****@hotmail.com

Posts by Illinu

<prev • 32 results • page 1 of 4 • next >
0
votes
1
answer
260
views
1
answers
Comment: C: too many adj pvalues = 1
... Hi Devon, It just seems weird that the test finds insignificant Gene A between T and C for genotype B, when it goes from 2000 to 8000 counts with small error bars. To me there is a clear induction of this gene by the treatment compared to genotype 1. But the test considers this difference as not si ...
written 5 months ago by Illinu90
1
vote
1
answer
260
views
1
answer
too many adj pvalues = 1
... Hi, While analyzing a set of DEGs resulting from DESeq2, I noticed for some genes that were I 'see' differences they are not significant, but then another gene with similar behaviour, the difference is significant. They both come from different comparisons, the former with 6 replicates and the late ...
p-value deseq2 rna-seq written 5 months ago by Illinu90 • updated 5 months ago by Devon Ryan71k
0
votes
1
answer
215
views
1
answers
Answer: A: emtpy bcf file
... This was corrected at the mpileup step with -A option samtools mpileup -ug -A -f thal.fasta $bname-thal.nodup.sorted.bam > $bname-thal.bcf -A, --count-orphans Do not skip anomalous read pairs in variant calling. Maybe this is because the sequences I used as reference are short and the p ...
written 7 months ago by Illinu90
0
votes
1
answer
215
views
1
answers
Comment: C: emtpy bcf file
... Thanks @Pierre, I normally run the pipeline for some commands, some of those intermediate files I need for other tasks. But in the script above I am testing this way to see at what step it fails. So far up to the *.nodup.sorted.bam it works well, that's how I know that it is the bcf file that is not ...
written 8 months ago by Illinu90
0
votes
4
answers
951
views
4
answers
Comment: C: comparing three genomes?
... and what is the difference for you? It still consists of fasta files with sequences right? The question is how do you do it for three sets of 'genes' (if you want) instead of two. The problem is that what you propose doesn't work when the gene names are not the same, and when there is gene expansion ...
written 8 months ago by Illinu90
0
votes
1
answer
215
views
1
answer
emtpy bcf file
... I am aligning genomic reads to a set of 6 genes all saved in one file in fasta format. The sam and bam files look ok but the mpileup step returns empty bcf files, with only the header until: #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT sample1-thal.nodup.sorted.bam Th ...
alignment samtools mpileup written 8 months ago by Illinu90
0
votes
2
answers
14k
views
2
answers
Comment: C: Venndiagram Using R
... Hi Ly, but with this method you cannot see up and down regulated genes right? You just plot significantly differentially expressed ...
written 8 months ago by Illinu90
0
votes
4
answers
951
views
4
answers
Comment: C: comparing three genomes?
... I don't see this for a set of three transcriptomes, it only presents option for comparing two nucleotide databases, can you explain how you do this if you have three databases? Thank you ...
written 8 months ago by Illinu90
0
votes
1
answer
592
views
1
answers
Comment: C: vcf-consensus Error: Fasta sequence does not match REF
... Hi microfuge, I am getting the same error and I mapping reads from several genotypes to one heterozygous reference and was hoping to get the consensus files for each genotype. Does this error mean it is not possible to get the consensus for each genotype individually? ...
written 10 months ago by Illinu90
1
vote
1
answer
400
views
1
answers
Comment: C: create a gene2accession file from a 2 column file
... Thanks a million @EagleEye, I asked in stackoverflow and my question was blocked because it was '*unclear*' what I was asking for. Note that to get rid of the last occurrence of '|' I used sed 's/BACK_SLASH(.*\)|/\1/'. Thanks again ...
written 16 months ago by Illinu90

Latest awards to Illinu

Epic Question 16 months ago, created a question with more than 10,000 views. For filter out multiple rows based on values in columns - perl and R attempt
Great Question 2.3 years ago, created a question with more than 5,000 views. For filter out multiple rows based on values in columns - perl and R attempt
Popular Question 2.3 years ago, created a question with more than 1,000 views. For Python: count how many lines have a specific word
Popular Question 2.3 years ago, created a question with more than 1,000 views. For filter out multiple rows based on values in columns - perl and R attempt

Help
Access

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