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

Moderator: Jorge Amigo

gravatar for Jorge Amigo
Jorge Amigo ♦ 10k
Reputation:
10,310
Status:
Trusted
Location:
Santiago de Compostela, Spain
Website:
https://www.researchga...
Scholar ID:
Google Scholar Page
Last seen:
23 hours ago
Joined:
7 years, 11 months 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 • 749 results • page 1 of 75 • next >
0
votes
1
answer
170
views
1
answers
Answer: A: filtering a VCF file based on genotype
... from your question I understand that you have 1 single vcf file with 2 samples (germline and tumor) in it, and that you want to filter out all variants where 1 or the 2 samples have an empty genotype (./.) if this is the case, a simple `grep` should be enough: grep -vP '\t\./\.' file.vcf > ...
written 8 days ago by Jorge Amigo ♦ 10k
0
votes
2
answers
151
views
2
answers
Comment: C: Is indel realigning necessary for INDEL discovery?
... I must agree with you both. I've updated my answer to be more precise on what GATK states and how I personally have always considered it. thank you for the clarification. ...
written 5 weeks ago by Jorge Amigo ♦ 10k
2
votes
2
answers
151
views
2
answers
Answer: A: Is indel realigning necessary for INDEL discovery?
... GATK's HaplotypeCaller is both capable of detecting SNVs and InDels using a method that performs local *de novo* assembly (kind of a local realignment) to call variants, although it doesn't output any realigned bam. so, in summary, there's no need to use IndelRealigner if you are going to call varia ...
written 5 weeks ago by Jorge Amigo ♦ 10k
3
votes
3
answers
1.8k
views
3
answers
Answer: A: find positions of a short sequence in a genome
... here are my 2 cents for a one-liner perl solution to look for a particular motif in a genome (change the `$seqMotif` variable to look for your motif of interest; note that it accepts a `regex` pattern) that prints the positions where the motif happens to appear: time perl -ne 'BEGIN { $seqM ...
written 6 months ago by Jorge Amigo ♦ 10k
0
votes
6
answers
11k
views
6
answers
Comment: C: Splitting A Vcf File
... thanks a lot for this idea. it is indeed much faster (3x at least on a local test splitting 3 exomes), and it can be condensed as follows: bcftools query -l MyData.vcf.gz | parallel -a - \ bcftools view -c1 -s {} -Oz --threads 8 -o {}.vcf.gz MyData.vcf.gz ...
written 6 months ago by Jorge Amigo ♦ 10k
0
votes
5
answers
16k
views
5
answers
Comment: C: How To Split Multiple Samples In Vcf File Generated By Gatk?
... sure. since you will be generating a file per sample you could just simply use `-o $sample.vcf.gz`, but I personally prefer to keep the original file name to know where that data came from. for that reason I use `.o ${file/.vcf*/.$sample.vcf.gz}`, which uses a bash string manipulation function to su ...
written 7 months ago by Jorge Amigo ♦ 10k
1
vote
11
answers
119k
views
11
answers
Comment: C: Tools To Calculate Average Coverage For A Bam File?
... that's because `NR` is an `awk` variable that stores the number of rows, therefore the `c` variable on my code is not wrong, but it is superfluous. ...
written 7 months ago by Jorge Amigo ♦ 10k
0
votes
1
answer
3.9k
views
1
answers
Comment: C: Filter Bam File Based On Coverage
... first, you have to detect those low coverage regions. then, you have to filter your bam file with those regions. I would suggest you open a new question rather than commenting a previous related one. ...
written 8 months ago by Jorge Amigo ♦ 10k
0
votes
11
answers
119k
views
11
answers
Comment: C: Tools To Calculate Average Coverage For A Bam File?
... in addition, if you have already plotted depths it means that you have got a temporal R data.frame (or similar type) in there containing chromosome, position and depth information. splitting such data.frame by chromosome and calculating average values should be straight-forward. I would just suggest ...
written 10 months ago by Jorge Amigo ♦ 10k
0
votes
0
answers
565
views
0
answers
Comment: C: RNA-SeQC gatkFlags option reading first word only
... it seems to be a problem with my `java` calling script. I also use aliases for `java7` and `java8`, although I didn't include them in the original question. it works fine when calling java directly and not through the aliases, so I'll have to correct my aliases for the future. thanks anyway ;) by t ...
written 12 months ago by Jorge Amigo ♦ 10k

Latest awards to Jorge Amigo

Teacher 14 days ago, created an answer with at least 3 up-votes. For A: Which Version Of Gatk Do People Use
Good Answer 28 days 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 6 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 6 weeks ago, created a post with more than 5 votes. For A: Order Of Gatk Commands
Appreciated 8 weeks ago, created a post with more than 5 votes. For A: Order Of Gatk Commands
Teacher 8 weeks ago, created an answer with at least 3 up-votes. For A: Which Version Of Gatk Do People Use
Teacher 10 weeks ago, created an answer with at least 3 up-votes. For A: Which Version Of Gatk Do People Use
Teacher 11 weeks ago, created an answer with at least 3 up-votes. For A: Which Version Of Gatk Do People Use
Appreciated 12 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
Appreciated 3 months ago, created a post with more than 5 votes. For A: Order Of Gatk Commands
Popular Question 4 months ago, created a question with more than 1,000 views. For force reads to map to a particular genome region
Good Answer 4 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 4 months ago, created a post with more than 5 votes. For A: Order Of Gatk Commands
Good Answer 6 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 6 months ago, created an answer with at least 3 up-votes. For A: Which Version Of Gatk Do People Use
Appreciated 6 months ago, created a post with more than 5 votes. For A: Order Of Gatk Commands
Good Answer 7 months ago, created an answer that was upvoted at least 5 times. For A: Phased And Unphased Genotypes In Vcf Files: Does The Order Of Alleles Matter?
Appreciated 7 months ago, created a post with more than 5 votes. For A: Order Of Gatk Commands
Great Question 8 months ago, created a question with more than 5,000 views. For LinkedIn PubMed Importer
Scholar 9 months ago, created an answer that has been accepted. For A: Filtration of bam file but with header
Teacher 10 months ago, created an answer with at least 3 up-votes. For A: Which Version Of Gatk Do People Use
Teacher 10 months ago, created an answer with at least 3 up-votes. For A: Moore'S Law
Teacher 10 months ago, created an answer with at least 3 up-votes. For A: Which Version Of Gatk Do People Use
Teacher 10 months ago, created an answer with at least 3 up-votes. For A: Which Version Of Gatk Do People Use
Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1661 users visited in the last hour