User: abbysue

gravatar for abbysue
abbysue10
Reputation:
10
Status:
New User
Location:
Last seen:
5 days, 5 hours ago
Joined:
1 week, 4 days ago
Email:
a******@cableone.net

Posts by abbysue

<prev • 7 results • page 1 of 1 • next >
0
votes
1
answer
59
views
1
answers
Comment: C: BWA index getting killed
... I realized that recently, stupid mistake by me. Thanks!! ...
written 5 days ago by abbysue10
1
vote
1
answer
59
views
1
answer
BWA index getting killed
... I'm trying to index fastq files using this command bwa index -p testbwaidx -a bwtsw test.fastq However, the process keeps getting killed after ~10 minutes despite having 53G of memory open. Here's the error message [bwa_index] Pack FASTA... Killed ...
fastq bam bwa index written 5 days ago by abbysue10 • updated 5 days ago by swbarnes26.0k
0
votes
0
answers
105
views
0
answers
Comment: C: Pipe SamToFastq, BWA mem not working
... Nope (this 20 character min... damn) ...
written 5 days ago by abbysue10
1
vote
0
answers
105
views
0
answers
Comment: C: Pipe SamToFastq, BWA mem not working
... I ended up running gatk SamToFastq -I markedadapters.bam -F f797.fastq --CLIPPING_ATTRIBUTE XT --CLIPPING_ACTION 2 --INTERLEAVE true -NON_PF true -TMP_DIR= path/hello which of course worked. I tried your command and it worked!! So should I include the " > file.fastq" blurb to get my piped ...
written 5 days ago by abbysue10
8
votes
0
answers
105
views
5 follow
0
answers
Pipe SamToFastq, BWA mem not working
... I'm going through the GATK pipeline to realign bam files. I'm using this command line: gatk SamToFastq -I markedadapters.bam -F /dev/stdout --CLIPPING_ATTRIBUTE XT --CLIPPING_ACTION 2 --INTERLEAVE true -NON_PF true -TMP_DIR= path/hello | bwa mem -M -t 12 -p /mydata/Homo_sapiens_assembly38.fasta ...
gatk sam bwa pipe written 5 days ago by abbysue10
1
vote
1
answer
119
views
1
answer
Stats on WGS variant calls
... What's the best way to collect some statistics on 4 mil variants identified with GATK pipeline (and annotated with Funcotator)? I'm having trouble parsing info in the 'Funcotator' annotation since it contains several pieces of info that I want, separated by pipe characters. I want to know # of Intro ...
wgs variants written 10 days ago by abbysue10 • updated 10 days ago by reza.jabal330
0
votes
0
answers
111
views
0
answers
Parsing through VCF annotated by Funcotator
... I have a variant set annotated by GATK's Funcotator... creates an extra field in the VCF called "Funcotator" with several [annotations][1] separated by pipe characters. I want to get some statistics on % intronic vs exonic and missense/nonsense mutations. Is the easiest way to do this using grep? Or ...
vcf gatk annotated funcotator written 11 days ago by abbysue10

Latest awards to abbysue

No awards yet. Soon to come :-)

Help
Access

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