User: DNAngel

gravatar for DNAngel
DNAngel40
Reputation:
40
Status:
New User
Location:
Last seen:
3 weeks ago
Joined:
2 years, 7 months ago
Email:
k********************@hotmail.com

Posts by DNAngel

<prev • 81 results • page 1 of 9 • next >
0
votes
2
answers
1.4k
views
2
answers
Comment: C: Optimize bwa mem in parallel on a SLURM cluster
... How would you modify this if you have many reference files to process for many fastq files? If my reference.fasta files are saved in a large list in reflist.txt, and I have another directory with my fastq files to process how does one instead use parallel to process many reference.fasta files simul ...
written 24 days ago by DNAngel40
0
votes
1
answer
190
views
1
answers
Comment: C: How to parallelize bwa mem for single-end data within a larger bash script?
... I know but because I have multiple fastq ffiles to process for 1000s of ref genes (which the link you provided and others I've found here do not have this problem) I am unsure which part of my script to change. Do I use parallel to call the individual REF files of which I have thousands and in that ...
written 24 days ago by DNAngel40
1
vote
1
answer
190
views
1
answer
How to parallelize bwa mem for single-end data within a larger bash script?
... I've seen some posts here that have tried to parallelize bwa-mem for their sample data but all examples were paired-end data and the code was standalone (i.e. they only needed it for bwa mem and maybe converting SAM to BAM). But I also have to run my script for 1000s of different reference files and ...
bwa-mem gnu-parallel bash written 24 days ago by DNAngel40 • updated 24 days ago by Pierre Lindenbaum124k
0
votes
1
answer
150
views
1
answers
Comment: C: How do you extract longest transcript from fasta file using length and header ID
... It worked! Thank you :) ...
written 4 weeks ago by DNAngel40
2
votes
1
answer
150
views
1
answer
How do you extract longest transcript from fasta file using length and header IDs?
... Here is an example file that I made up to practise and it is called Medaka_TEST.fa where the header is GENE_ID | TRANSCRIPT_ID: $ cat Medaka_TEST.fa >ENSORLG000000000AA|ENSORLT000000000AD # This gene has two transcripts (next sequence) KEEPTHISONETHISISTHELONGEST >ENSORLG000 ...
fasta awk written 4 weeks ago by DNAngel40 • updated 4 weeks ago by JC9.1k
0
votes
3
answers
7.3k
views
3
answers
Comment: C: Generate Consensus sequence from BAM file
... Have you solved this? I am also trying to do the same with no luck :/ ...
written 5 weeks ago by DNAngel40
2
votes
0
answers
139
views
0
answers
How to loop over a directory to set each file as a variable?
... Hi all, I have a bash script that I can call and I supply command line arguments to feed specific files into it. It works fine, but I want to optimize it but looping over a new directory where each file is set as a new argument. For example my script below: #!/bin/bash REF=$1 NAME=$2 ...
bash written 5 weeks ago by DNAngel40
1
vote
1
answer
170
views
1
answers
Comment: C: Does Abyss work on single-end Illumina reads only?
... Ah okay I thought -se would be for single-end as well but it does say it requires a short and a long read assembly which I do not have. That's unfortuante! ...
written 5 weeks ago by DNAngel40
0
votes
1
answer
170
views
1
answers
Comment: C: Does Abyss work on single-end Illumina reads only?
... I am not working on microbial genomes, vertebrate exome data is what I have. I've asked a lot of questions and it seems that trying to assemble exome data is difficult becuase of missing intergenic regions. Unfortunately a lot of the papers I read in my field have used WGS but the data I HAVE to wor ...
written 5 weeks ago by DNAngel40
0
votes
1
answer
914
views
1
answers
Answer: A: How to get a VCF file for all variants in a specific gene and specific exon?
... Have you figured this out? I am also interested! ...
written 5 weeks ago by DNAngel40

Latest awards to DNAngel

Scholar 11 months ago, created an answer that has been accepted. For A: How to save trimmomatic terminal summary to an output file?

Help
Access

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