User: abascalfederico
abascalfederico • 1.1k
- Reputation:
- 1,080
- Status:
- Trusted
- Location:
- Spain
- Last seen:
- 10 hours ago
- Joined:
- 4 years, 4 months ago
- Email:
- a**************@gmail.com
Posts by abascalfederico
<prev
• 178 results •
page 1 of 18 •
next >
0
votes
0
answers
48
views
0
answers
... Anyone knows of a VCF containing all human-chimp divergent sites mapped on human genome coordinates? It would be so nice if this was available :-)
Thanks!
Federico
...
written 4 weeks ago by
abascalfederico • 1.1k
0
votes
1
answer
549
views
1
answers
... For a given position, I want to read the pileup and when there is a match to the variant, I want to get the mapping coordinates of the read matching the variant. I do not have the exact code yet, but it would imply something like:
for pileupcolumn in samfile.pileup(chr,int(pos)-1,int(pos)):
...
written 20 months ago by
abascalfederico • 1.1k
0
votes
1
answer
549
views
1
answers
... Thanks, Devon. I don't want the programmatic burden ;-) I will instead make a call to python from R - not what I wanted though ...
written 20 months ago by
abascalfederico • 1.1k
0
votes
1
answer
549
views
1
answers
... Thanks Devon, I already had a look at that function, but it just provide a table of counts. I need to identify which reads support a given variant. This is relatively easy with pysam but I find no way to do this in Rsamtools. Maybe I missed something from the documentation? ...
written 20 months ago by
abascalfederico • 1.1k
2
votes
1
answer
549
views
1
answer
... For a given position, in which there is a variant called by some algorithm, I want to get the pileup and analyse the reads supporting that variant. Any idea on how to do this in R?
I know how to do this in Python using pysam but I need an R implementation.
Thanks!
...
written 20 months ago by
abascalfederico • 1.1k
• updated
20 months ago by
Devon Ryan ♦ 92k
0
votes
0
answers
585
views
0
answers
... Can anyone recommend a tool to segment the genome according to SNP density (using HMM-based segmentation, for example)?
Thanks!
...
written 2.0 years ago by
abascalfederico • 1.1k
• updated
2.0 years ago by
Biostar ♦♦ 20
0
votes
2
answers
1.5k
views
2
answers
... That's exactly what I needed, thanks! ...
written 2.3 years ago by
abascalfederico • 1.1k
3
votes
2
answers
1.5k
views
2
answers
... Does anybody know of a program to add the NM tag (number of mismatches) to a BAM file?
Usually aligned BAM files come with that tag, but I've found one without it and need that information.
Thanks! ...
written 2.3 years ago by
abascalfederico • 1.1k
• updated
10 months ago by
pammybiotech • 0
0
votes
0
answers
854
views
0
answers
... I'll send them an email, thanks for the suggestion.
Regarding the use of bwa aln, I found it in the supplementary material of:
Lek, M., Karczewski, K. J., Minikel, E. V., Samocha, K. E., Banks, E., Fennell, T., ... & Tukiainen, T. (2016). Analysis of protein-coding genetic variation in 60,706 hu ...
written 2.3 years ago by
abascalfederico • 1.1k
4
votes
0
answers
854
views
0
answers
... Anyone knows why in ExAC they used bwa `aln` instead of `mem`? I am not sure of the read lengths but somewhere said something about 150 bps x 2. Does it make sense to use bwa aln here? I though for long reads mem was the best choice
The actual parameters used in ExAC were:
bwa aln q 5 -l 32 -k ...
written 2.3 years ago by
abascalfederico • 1.1k
Latest awards to abascalfederico
Voter
4 weeks ago,
voted more than 100 times.
Good Answer
4 months ago,
created an answer that was upvoted at least 5 times.
For A: samtools view - F 4 removes sam header and truncates file??
Appreciated
5 months ago,
created a post with more than 5 votes.
For A: samtools view - F 4 removes sam header and truncates file??
Popular Question
5 months ago,
created a question with more than 1,000 views.
For Uniquely mapped PAIRS of reads
Popular Question
5 months ago,
created a question with more than 1,000 views.
For samtools mpileup: position of each base within their reads
Popular Question
5 months ago,
created a question with more than 1,000 views.
For Recode VCF file according to ancestral alleles (AA)
Popular Question
8 months ago,
created a question with more than 1,000 views.
For strange issues with absent EOF markers in BAM files
Popular Question
11 months ago,
created a question with more than 1,000 views.
For strange issues with absent EOF markers in BAM files
Popular Question
13 months ago,
created a question with more than 1,000 views.
For strange issues with absent EOF markers in BAM files
Popular Question
13 months ago,
created a question with more than 1,000 views.
For samtools mpileup: position of each base within their reads
Popular Question
14 months ago,
created a question with more than 1,000 views.
For strange issues with absent EOF markers in BAM files
Popular Question
19 months ago,
created a question with more than 1,000 views.
For strange issues with absent EOF markers in BAM files
Guru
20 months ago,
received more than 100 upvotes.
Popular Question
21 months ago,
created a question with more than 1,000 views.
For strange issues with absent EOF markers in BAM files
Popular Question
24 months ago,
created a question with more than 1,000 views.
For strange issues with absent EOF markers in BAM files
Popular Question
24 months ago,
created a question with more than 1,000 views.
For Duplicated read names in Fastq file
Popular Question
2.0 years ago,
created a question with more than 1,000 views.
For Duplicated read names in Fastq file
Popular Question
2.3 years ago,
created a question with more than 1,000 views.
For List of CNV (copy number variation) regions in the human genome
Student
2.3 years ago,
asked a question with at least 3 up-votes.
For Samtools' mpileup with pysam: enabling and disabling BAQ probabilistic realignment
Commentator
2.4 years ago,
created a comment with at least 3 up-votes.
For C: Is it good practice to convert lowewst quality genotyes to missing?
Popular Question
2.4 years ago,
created a question with more than 1,000 views.
For List of CNV (copy number variation) regions in the human genome
Popular Question
2.6 years ago,
created a question with more than 1,000 views.
For F3 test for admixture
Commentator
2.7 years ago,
created a comment with at least 3 up-votes.
For C: Is it good practice to convert lowewst quality genotyes to missing?
Scholar
2.7 years ago,
created an answer that has been accepted.
For A: How to properly merge multiple runs into one bam file ?
Teacher
2.7 years ago,
created an answer with at least 3 up-votes.
For A: distance between two genes
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 2013 users visited in the last hour