Moderator: Jorge Amigo

gravatar for Jorge Amigo
Jorge Amigo11k
Reputation:
10,960
Status:
Trusted
Location:
Santiago de Compostela, Spain
Website:
https://www.researchga...
Scholar ID:
Google Scholar Page
Last seen:
1 week, 6 days ago
Joined:
9 years, 1 month ago
Email:
a****@yahoo.com

Scrutinizing genomic human variation by dealing with high throughput genotyping and next generation sequencing results, among many other things.

Bioinformatician @ Genomic Medicine Group

Hospital Clínico Universitario, Santiago de Compostela, Spain

Posts by Jorge Amigo

<prev • 758 results • page 1 of 76 • next >
1
vote
10
answers
44k
views
10
answers
Comment: C: Multiline Fasta To Single Line Fasta
... thanks for pointing it out. I've corrected my previous answer and tested thoroughly the new one. ...
written 4 months ago by Jorge Amigo11k
0
votes
3
answers
10k
views
3
answers
Comment: C: how to remove multiallelic from VCF
... note that the merging process when normalizing can be quite intense on such large files as gnomad's, so if you are interested in particular regions do not forget to indicate this in your normalization command. I've edited my answer to reflect this. ...
written 5 months ago by Jorge Amigo11k
1
vote
3
answers
10k
views
3
answers
Comment: C: how to remove multiallelic from VCF
... I have always used vcf.gz, although I've seen vcf.bgz probably to indicate that the compression has been performed with the `bgzip` algorithm rather than simply `gzip`. but this is simply informative for the human user, not for the software. ...
written 5 months ago by Jorge Amigo11k
1
vote
3
answers
10k
views
3
answers
Comment: C: how to remove multiallelic from VCF
... [bcftools' documentation][1] is very clear about this > Use -m2 -M2 -v snps to only view biallelic SNPs. so if it isn't working it must be either a problem with `bcftools` (that'd be odd) or either a problem with gnomad data. if you query gnomad data for that particular position you'll see th ...
written 5 months ago by Jorge Amigo11k
0
votes
0
answers
1.1k
views
0
answers
Comment: C: Merge hundreds bedgraph files with unionbedg in one bash line
... that's precisely what `-header` option does. the above command `bedtools unionbedg -i *.bedgraph -filler NA -header` does exactly what you need. ...
written 6 months ago by Jorge Amigo11k
0
votes
1
answer
5.9k
views
1
answers
Comment: C: bwa mem on a bam file
... from [bwa manual page][1]: > If -p is used, the command assumes the 2i-th and the (2i+1)-th read in reads.fq constitute a read pair (such input file is said to be interleaved) but unfortunately `samtools fastq` does not output interleaved reads. you can either try sorting reads by name BEFORE ...
written 7 months ago by Jorge Amigo11k
1
vote
3
answers
1.6k
views
3
answers
Comment: C: Remove a list of positions form a VCF file
... my answer was very simple. this one adds more *grep* functionality: `-F` option looks for fixed strings rather than regular expressions, and `-w` option looks for whole words rather than just matching patterns. I don't know how `-F` works in conjunction with `-w`, but it looks like an overall faster ...
written 7 months ago by Jorge Amigo11k
0
votes
10
answers
44k
views
10
answers
Comment: C: Multiline Fasta To Single Line Fasta
... just create a simple script.pl file containing this (do not use!) while (<>) { chomp unless /^>/; print } EDIT: the previous code is wrong. use this one into script.pl instead: while (<>) { $. > 1 and /^>/ ? print "\n" : chomp; print } and run it perl script.pl ...
written 8 months ago by Jorge Amigo11k
1
vote
2
answers
3.2k
views
2
answers
Answer: A: Identifying private and shared SNPs using VCFtools
... update: *bcftools isec* does exactly these kind of comparisons, and faster. this quick example bcftools isec -p isec_results file1 file2 generates in a *isec_results* folder vcf files with private variants from *file1*, private variants from *file2*, *file1* variants shared with *file2*, and * ...
written 11 months ago by Jorge Amigo11k
5
votes
3
answers
652
views
3
answers
Answer: A: h38 total genes ?
... locate refGene file, uncompress it, look for the 13th column containing gene names, make sure there are no repetitions, and count them: curl http://hgdownload.cse.ucsc.edu/goldenPath/hg38/database/refGene.txt.gz \ | zcat | cut -f13 | sort -u | wc -l 28054 ...
written 13 months ago by Jorge Amigo11k

Latest awards to Jorge Amigo

Appreciated 6 weeks ago, created a post with more than 5 votes. For A: Order Of Gatk Commands
Good Answer 8 weeks ago, created an answer that was upvoted at least 5 times. For A: Is It Ok To Use One End Of A Set Of Paired-End Reads As A Set Of Single Reads?
Appreciated 10 weeks ago, created a post with more than 5 votes. For A: Order Of Gatk Commands
Appreciated 3 months ago, created a post with more than 5 votes. For A: Order Of Gatk Commands
Teacher 4 months ago, created an answer with at least 3 up-votes. For A: Which Version Of Gatk Do People Use
Good Answer 5 months ago, created an answer that was upvoted at least 5 times. For A: Is It Ok To Use One End Of A Set Of Paired-End Reads As A Set Of Single Reads?
Scholar 5 months ago, created an answer that has been accepted. For A: Filtration of bam file but with header
Appreciated 5 months ago, created a post with more than 5 votes. For A: Order Of Gatk Commands
Epic Question 7 months ago, created a question with more than 10,000 views. For How To Split A .Vcf.Gz File
Teacher 7 months ago, created an answer with at least 3 up-votes. For A: Which Version Of Gatk Do People Use
Great Question 8 months ago, created a question with more than 5,000 views. For Which Programs Are You Relying On For Solid Data Analysis?
Good Answer 8 months ago, created an answer that was upvoted at least 5 times. For A: Is It Ok To Use One End Of A Set Of Paired-End Reads As A Set Of Single Reads?
Good Answer 8 months ago, created an answer that was upvoted at least 5 times. For A: Is It Ok To Use One End Of A Set Of Paired-End Reads As A Set Of Single Reads?
Teacher 8 months ago, created an answer with at least 3 up-votes. For A: Which Version Of Gatk Do People Use
Good Answer 9 months ago, created an answer that was upvoted at least 5 times. For A: Is It Ok To Use One End Of A Set Of Paired-End Reads As A Set Of Single Reads?
Appreciated 9 months ago, created a post with more than 5 votes. For A: Order Of Gatk Commands
Commentator 10 months ago, created a comment with at least 3 up-votes. For C: How To Analyse Snp Data From Different Sources?
Scholar 11 months ago, created an answer that has been accepted. For A: Filtration of bam file but with header
Great Question 11 months ago, created a question with more than 5,000 views. For LinkedIn PubMed Importer
Appreciated 11 months ago, created a post with more than 5 votes. For A: Order Of Gatk Commands
Teacher 11 months ago, created an answer with at least 3 up-votes. For A: Which Version Of Gatk Do People Use
Teacher 12 months ago, created an answer with at least 3 up-votes. For A: Which Version Of Gatk Do People Use
Teacher 13 months ago, created an answer with at least 3 up-votes. For A: Which Version Of Gatk Do People Use
Scholar 13 months ago, created an answer that has been accepted. For A: Filtration of bam file but with header
Good Answer 13 months ago, created an answer that was upvoted at least 5 times. For A: Is It Ok To Use One End Of A Set Of Paired-End Reads As A Set Of Single Reads?

Help
Access

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