Moderator: Jorge Amigo
Jorge Amigo ♦ 11k
- Reputation:
- 11,340
- Status:
- Trusted
- Location:
- Santiago de Compostela, Spain
- Website:
- https://www.researchga...
- Scholar ID:
- Google Scholar Page
- Last seen:
- 28 minutes ago
- Joined:
- 9 years, 10 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.
Posts by Jorge Amigo
<prev
• 771 results •
page 1 of 78 •
next >
1
vote
9
answers
40k
views
9
answers
Comment:
C: How to count fastq reads
... the only very slight problem I see is that very rare 1 base sequences (trimmed perhaps) with "+" quality would be counted twice, but I love the simplicity of this answer. ...
written 4 days ago by
Jorge Amigo ♦ 11k
0
votes
1
answer
160
views
1
answers
... although chr1 is one of the biggest ones, and considering that you've already done the largest effort which is calling variants on NA12878, I'd definitely go for the entire genome's numbers rather than being fine with a relatively small subset of statistics. ...
written 12 weeks ago by
Jorge Amigo ♦ 11k
0
votes
1
answer
160
views
1
answers
... you've chosen a reference sample to apply your variant calling pipeline to. well done, because you're almost there. just get a set of [high confidence NA12878 variants][1] such as the [GiB project][2] ones, use a comparison tool such as [RTG vcfeval][3] (the one recommended by the [GiB project][4]), ...
written 12 weeks ago by
Jorge Amigo ♦ 11k
2
votes
2
answers
137
views
2
answers
... it's difficult to say anything without knowing the content of your VCF file, but here are a couple of suggestions:
1. have you checked your VCF and your BED files refer to the same reference? one quick dirty check is making sure that you're either using or not using in both files the "chr" prefix, ...
written 12 weeks ago by
Jorge Amigo ♦ 11k
1
vote
1
answer
148
views
1
answers
... you should check if IGV is filtering any kind of reads. typically this is caused by IGV's default duplicate reads filtering.
just unselect everything you see under "View > Preferences > Alignments > Filter..." and you should see in IGV's coverage track the same profile than the Bedgraph ha ...
written 12 weeks ago by
Jorge Amigo ♦ 11k
0
votes
1
answer
165
views
1
answers
... if you already have exome data from each parent, if you have aligned it, and you have called variants, then you would have already ended up with 2 bam files and 2 vcf files.
if you need now to compare variants among them you don't have to map them again. you just have to compare the vcf files you a ...
written 4 months ago by
Jorge Amigo ♦ 11k
0
votes
2
answers
313
views
2
answers
... sorting such a big file only for being able to use the uniq function wouldn't be the more efficient way to do it, plus a simple `sort -u` wouldn't address the problem described in the question as the lines that need to be merged/skipped do vary. if sorting would be considered (and it would be really ...
written 4 months ago by
Jorge Amigo ♦ 11k
0
votes
2
answers
313
views
2
answers
... this is the same as an uniq function on entire lines, and it doesn't work with your example considering that positions qualities vary from line to line. you need to group (index) only the columns you expect to be repeated, as explained in the example I suggested previously. ...
written 4 months ago by
Jorge Amigo ♦ 11k
3
votes
2
answers
255
views
2
answers
... I don't know of a direct method to accomplish what you need, but have a look at this tentative proposal:
awk 'FS=OFS="\t"{print $1, 0, $2}' human_hg19.fa.fai \
| bedtools makewindows -b - -w 1000 \
| bedtools map -a - -b input.bedgraph -c 4 -o sum \
| grep -P "\d$"
what it does is: ...
written 4 months ago by
Jorge Amigo ♦ 11k
2
votes
2
answers
144
views
2
answers
Comment:
C: awk gff column change
... not caring about the number of total columns:
awk '{n=split($1,a,"_"); $4+=a[n]; $5+=a[n]; print}' input.txt
not caring about the number of total columns, plus considering that both input and output are tabulated:
awk 'FS=OFS="\t"{n=split($1,a,"_"); $4+=a[n]; $5+=a[n]; print}' input.txt ...
written 4 months ago by
Jorge Amigo ♦ 11k
Latest awards to Jorge Amigo
Commentator
7 weeks ago,
created a comment with at least 3 up-votes.
For C: How To Analyse Snp Data From Different Sources?
Good Answer
7 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?
Scholar
3 months ago,
created an answer that has been accepted.
For A: Filtration of bam file but with header
Good Answer
3 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
3 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
4 months ago,
created an answer that has been accepted.
For A: Filtration of bam file but with header
Teacher
4 months ago,
created an answer with at least 3 up-votes.
For A: Which Version Of Gatk Do People Use
Scholar
5 months ago,
created an answer that has been accepted.
For A: Filtration of bam file but with header
Scholar
5 months ago,
created an answer that has been accepted.
For A: Filtration of bam file but with header
Epic Question
8 months ago,
created a question with more than 10,000 views.
For LinkedIn PubMed Importer
Popular Question
8 months ago,
created a question with more than 1,000 views.
For Which Programs Are You Relying On For Solid Data Analysis?
Good Answer
10 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
12 months ago,
created an answer with at least 3 up-votes.
For A: Which Version Of Gatk Do People Use
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?
Scholar
13 months ago,
created an answer that has been accepted.
For A: Filtration of bam file but with header
Epic Question
15 months ago,
created a question with more than 10,000 views.
For How To Split A .Vcf.Gz File
Teacher
16 months ago,
created an answer with at least 3 up-votes.
For A: Which Version Of Gatk Do People Use
Great Question
16 months ago,
created a question with more than 5,000 views.
For Which Programs Are You Relying On For Solid Data Analysis?
Good Answer
16 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
16 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?
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 809 users visited in the last hour