User: steven.davis
steven.davis • 10
- Reputation:
- 10
- Status:
- New User
- Location:
- Last seen:
- 1 year, 8 months ago
- Joined:
- 4 years, 2 months ago
- Email:
- s***********@fda.hhs.gov
Posts by steven.davis
<prev
• 7 results •
page 1 of 1 •
next >
0
votes
0
answers
869
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 22 months ago by
steven.davis • 10
• updated
7 months ago by
Biostar ♦♦ 20
0
votes
1
answer
2.2k
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 2.4 years ago by
steven.davis • 10
0
votes
1
answer
1.4k
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 3.8 years ago by
steven.davis • 10
0
votes
4
answers
14k
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 3.8 years ago by
steven.davis • 10
1
vote
1
answer
1.4k
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 3.8 years ago by
steven.davis • 10
• updated
3.8 years ago by
Pierre Lindenbaum ♦ 133k
0
votes
1
answer
2.5k
views
1
answers
... Solved it with help from a colleague. Problem solved with "--merge all" command line option. ...
written 4.2 years ago by
steven.davis • 10
1
vote
1
answer
2.5k
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 4.2 years ago by
steven.davis • 10
Latest awards to steven.davis
Popular Question
3.8 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: 1695 users visited in the last hour