User: Illinu

gravatar for Illinu
Illinu80
Reputation:
80
Status:
Trusted
Location:
Belgium
Last seen:
2 weeks, 3 days ago
Joined:
2 years, 5 months ago
Email:
s****@hotmail.com

Posts by Illinu

<prev • 30 results • page 1 of 3 • next >
0
votes
1
answer
87
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 4 weeks ago by Illinu80
0
votes
1
answer
87
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 5 weeks ago by Illinu80
0
votes
4
answers
790
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 5 weeks ago by Illinu80
0
votes
1
answer
87
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 5 weeks ago by Illinu80
0
votes
2
answers
13k
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 6 weeks ago by Illinu80
0
votes
4
answers
790
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 7 weeks ago by Illinu80
0
votes
1
answer
370
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 3 months ago by Illinu80
1
vote
1
answer
281
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 9 months ago by Illinu80
4
votes
1
answer
281
views
1
answer
create a gene2accession file from a 2 column file
... I need to create a gene2accession file where one gene ID might have one or more Uniprot entries. My file looks like this: comp100002_c0 Q9FFI3 comp100004_c0 B9DHK3 comp100004_c0 F4J3J5 comp100005_c0 P54150 and I need it to look like this: comp100002_c0 Q9FFI3 c ...
cluego gene2accession go enrichment uniprot written 9 months ago by Illinu80 • updated 9 months ago by EagleEye3.4k
0
votes
3
answers
606
views
3
answers
Answer: A: grep unexpected behaviour
... Just to update, this form: for f in $(cat list.OE.txt); do grep -w $f Trinity_uniref_2015_02_filt_ann_out.txt; done > OE.annocript.txt works A LOT faster than this one: grep -wf list.OE.txt Trinity_uniref_2015_02_filt_ann_out.txt > OE.annocript.txt ...
written 2.0 years ago by Illinu80

Latest awards to Illinu

Epic Question 9 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 21 months 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 21 months ago, created a question with more than 1,000 views. For Python: count how many lines have a specific word
Popular Question 21 months 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: 1497 users visited in the last hour