User: steven.davis
steven.davis • 10
- Reputation:
- 10
- Status:
- New User
- Location:
- Last seen:
- 2 weeks, 6 days ago
- Joined:
- 2 years, 6 months ago
- Email:
- s***********@fda.hhs.gov
Posts by steven.davis
<prev
• 7 results •
page 1 of 1 •
next >
0
votes
0
answers
126
views
0
answers
... What is the standard VCF INFO field, if any, to count the occurrences of variant genotypes having passed all filters?
Does AC count every occurrence of the alternate genotype, or only those occurrences where the call passes all filters?
Consider the haploid VCF snippet below. At position 300, the ...
written 9 weeks ago by
steven.davis • 10
0
votes
1
answer
1.1k
views
1
answers
... @jrj.healey This is 15 months after your original post, but may help others. There is a maf2fasta.pl conversion script included with mugsy. How did you solve the conversion problem? ...
written 9 months ago by
steven.davis • 10
0
votes
1
answer
883
views
1
answers
... Thank you Pierre. I'm writing a Python module to parse the sequence ids and extract various data elements. I will put it on GitHub sometime soon. ...
written 2.2 years ago by
steven.davis • 10
0
votes
4
answers
7.1k
views
4
answers
... What kind of instrument uses flowcells ending with BBXX? Here is the read sequence id line:
@ERR1417747.1 K00150:65:H77FVBBXX:7:1101:1570:1297 length=150 ...
written 2.2 years ago by
steven.davis • 10
1
vote
1
answer
883
views
1
answer
... Is there a way to extract instrument and flowcell from a read name like this:
>gnl|SRA|SRR1840614.1.1 FCC1KPRACXX:1:1101:1291:2172
FCC1KPRACXX : what is this?
1 = lane?
1101 = ?
1291 = x
2172 = y
I need to extract the flowcell if possible, so it can be assigned to reads ...
written 2.2 years ago by
steven.davis • 10
• updated
2.2 years ago by
Pierre Lindenbaum ♦ 120k
0
votes
1
answer
1.8k
views
1
answers
... Solved it with help from a colleague. Problem solved with "--merge all" command line option. ...
written 2.5 years ago by
steven.davis • 10
1
vote
1
answer
1.8k
views
1
answer
... I am merging multiple single-sample VCF files into one multi-sample VCF file with bcftools merge using this command:
bcftools merge --info-rules NS:sum -o outFilePath dir/*.gz
Previously, using bcftools 1.1, the output VCF file contained one row per position, like this:
#CHROM POS ID ...
written 2.5 years ago by
steven.davis • 10
Latest awards to steven.davis
Popular Question
2.2 years ago,
created a question with more than 1,000 views.
For VCF output when merging with bcftools merge
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 1799 users visited in the last hour