User: bari.ballew

gravatar for bari.ballew
bari.ballew190
Reputation:
190
Status:
Trusted
Location:
USA/NIH
Twitter:
@BariBallew
Last seen:
6 days, 3 hours ago
Joined:
1 year, 9 months ago
Email:
b**********@gmail.com

Posts by bari.ballew

<prev • 26 results • page 1 of 3 • next >
0
votes
1
answer
78
views
1
answers
Answer: A: Selecting set of SNPs using vcftools
... Try altering chr11_snp_pos.txt to have a single tab-delimited position per line, as described in the vcftools documentation. E.g.: 1195855385 1172496148 Are there two different alt alleles at position 11:72496148? That could result in more than one line for a single position in the VCF. ...
written 5 weeks ago by bari.ballew190
1
vote
1
answer
178
views
1
answers
Answer: A: Snakemake FastQC MissingOutputFiles ErrorM
... FastQC names output files by taking the original filename, removing the `.fastq` extension, and replacing the extension with `_fastqc.html` or `_fastqc.zip`. So, your output should be `reports/{sample}.trimmed_fastq.html` and `reports/{sample}.trimmed_fastq.zip` Also, you only need to include the ...
written 5 weeks ago by bari.ballew190
0
votes
1
answer
112
views
1
answers
Answer: A: Query regarding somatic variant calling
... What exactly is "Non-Tumor"? To call somatic variants in your tumor sample, you should use a somatic caller (e.g. MuTect2) with the blood sample as the normal. If Non-Tumor is a nearby normal, that is, a biopsy of normal tissue taken adjacent or very close to the site of the tumor, I would not rec ...
written 8 weeks ago by bari.ballew190
0
votes
2
answers
154
views
2
answers
Answer: A: Removing white space from the beginning of the second field (sequence) in a fast
... Take a look at regex anchors, which tie your pattern to the beginning or end of a line. You can use sed to remove whitespace at the beginning of a line only using "^", which anchors the pattern to the beginning of the line. sed 's/^[\t ]*//' file.fa > secondFile.fa ...
written 8 weeks ago by bari.ballew190
2
votes
1
answer
104
views
1
answers
Answer: A: Not call some variants in some populations?
... To expand on what [RamRS][1] said, if you've merged VCFs and not gVCFs, note that VCFs only report a genomic location if there is a variant in that individual, so you are susceptible to a missing data problem. When a variant is reported in A.vcf, but not in B.vcf, the merged file will record the var ...
written 9 weeks ago by bari.ballew190
0
votes
1
answer
130
views
1
answers
Answer: A: Merging SNV and DIV vcfs of the same Individual
... You want to concatenate (add rows), not merge (add individuals). See my answer [here][1]. [1]: https://www.biostars.org/p/390264/#390538 ...
written 9 weeks ago by bari.ballew190
3
votes
2
answers
194
views
2
answers
Answer: A: skipping column names using awk
... You can get awk to behave differently for different rows by using conditional statements based on either line number: (NR==1) {do something} (NR>1) {do something else} Or by a characteristic of the line: $0~/^#CHROM/ {do something} $0!~/^#CHROM/ {do something else} HOWEVER, I strongly ...
written 10 weeks ago by bari.ballew190
0
votes
2
answers
4.3k
views
2
answers
Answer: A: Best way to query VCF for specific variants
... This question is rather old and has an accepted answer, but I'll add another option since [Ikirola](https://www.biostars.org/u/57233/) has found the post and is having trouble. Try [bcftools isec](https://samtools.github.io/bcftools/bcftools.html#isec). Make a minimal sites-only vcf containing y ...
written 11 weeks ago by bari.ballew190
3
votes
2
answers
248
views
2
answers
Answer: A: VCF merge or concatenate?
... Take a look at `bcftools concat` (https://samtools.github.io/bcftools/bcftools.html#concat). In general, concatenate means adding on rows to your vcf (e.g. re-combining split chromosomes), while merge means creating a superset of variant calls across multiple individuals. Concatenating is relati ...
written 3 months ago by bari.ballew190
0
votes
1
answer
170
views
1
answers
Comment: A: MutSig2CV with self-generated "maf file:
... It sounds like MutSig2CV is looking for, but not finding, a gene name column called either "gene," "Hugo_Symbol," or "Gene_name." Can you show us the first few lines of your input file? ...
written 3 months ago by bari.ballew190

Latest awards to bari.ballew

Scholar 9 weeks ago, created an answer that has been accepted. For A: tool to detect large deletions
Scholar 10 weeks ago, created an answer that has been accepted. For A: tool to detect large deletions
Teacher 10 weeks ago, created an answer with at least 3 up-votes. For A: tool to detect large deletions
Teacher 10 weeks ago, created an answer with at least 3 up-votes. For A: tool to detect large deletions
Scholar 6 months ago, created an answer that has been accepted. For A: tool to detect large deletions
Teacher 6 months ago, created an answer with at least 3 up-votes. For A: tool to detect large deletions

Help
Access

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