User: Niek De Klein
Niek De Klein • 2.5k
- Reputation:
- 2,500
- Status:
- Trusted
- Location:
- Netherlands
- Last seen:
- 10 months, 4 weeks ago
- Joined:
- 10 years, 4 months ago
- Email:
- n**********@gmail.com
Posts by Niek De Klein
<prev
• 251 results •
page 1 of 26 •
next >
0
votes
1
answer
1.6k
views
1
answers
Comment:
C: Shapeit can't open existing file
... For me the file chr1.sample did not exist, so the error was pointing to the same file. If your problem is the same as mine you should check all input files if they exist, not only the input file that they show in the error. ...
written 11 months ago by
Niek De Klein • 2.5k
0
votes
2
answers
4.3k
views
2
answers
... Haven't tested for edge cases but can try something like this
def consensus(*sequences):
if(len(sequences) <= 1):
raise RuntimeError('Only '+str(len(sequences))+' sequences given, need at least 2')
consensus_sequence = ''
for index, amino_acid in ...
written 2.2 years ago by
Niek De Klein • 2.5k
0
votes
0
answers
1.6k
views
0
answers
... The problem was in the .dict file, which was apparently made using the gzipped reference genome fasta, while the CRAM was made using the unzipped reference genome fasta. I remade the .dict file with the unzipped genome fasta file and now it works.
Thanks for the help! ...
written 2.3 years ago by
Niek De Klein • 2.5k
0
votes
0
answers
1.6k
views
0
answers
... head -n30 /apps/data//ftp.ebi.ac.uk/pub/databases/gencode/Gencode_human/release_24/GRCh38.p5.genome.fa.fai
chr1 248956422 8 60 61
chr2 242193529 253105712 60 61
chr3 198295559 499335808 60 61
chr4 190214555 700936301 60 61
chr5 181538259 894321107 60 61
chr6 170805979 10788 ...
written 2.3 years ago by
Niek De Klein • 2.5k
0
votes
0
answers
1.6k
views
0
answers
... Yes the path exists and the double / is not a problem. If the reference file can not be found it gives this error: `Exception in thread "main" htsjdk.samtools.SAMException: Error opening file: /fake/path/genome.fa` ...
written 2.3 years ago by
Niek De Klein • 2.5k
0
votes
0
answers
1.6k
views
0
answers
... This is a cross-post from here: https://gatkforums.broadinstitute.org/gatk/discussion/13395/sequences-at-index-0-dont-match-but-using-the-same-reference-genome#latest
I am trying to run CollectMultipleMetrics from Picard tools on a CRAM file but I get an "Sequences at index 0 don't match, but usin ...
written 2.3 years ago by
Niek De Klein • 2.5k
0
votes
1
answer
4.3k
views
1
answers
... I made a stupid mistake, was counting lines of the gzipped fastq file and comparing to an unzipped fastq file. After fixing that there was still a small difference, as STAR by default doesn't put unmapped reads in the output file, so had to --outSAMunmapped Within option.
After that I got the corr ...
written 2.4 years ago by
Niek De Klein • 2.5k
0
votes
1
answer
4.3k
views
1
answers
... Thanks, I'm testing this now and seeing how it compares to removing secondary aligned reads ...
written 2.4 years ago by
Niek De Klein • 2.5k
7
votes
1
answer
4.3k
views
7 follow
1
answer
... I want to save my fastq files as CRAM files to save space. However, I do want to be able to get the reads back that I have in my fastq, in case I have to reallign in a different way in the future. Below I have the code that I use to do the conversions to BAM and back again. However, when I then chec ...
written 2.4 years ago by
Niek De Klein • 2.5k
• updated
2.4 years ago by
finswimmer ♦ 14k
2
votes
1
answer
1.1k
views
1
answer
... I have a file with RNA editing sites which I want to remove from my VCF. I can select all the RNA editing sites with:
java GenomeAnalysisTK.jar \
-T SelectVariants \
-R human_g1k_v37.fasta \
-V input.gg.vcf.gz \
-o output.gg.vcf.gz \
-L Human_AG_all_hg19_ ...
written 3.4 years ago by
Niek De Klein • 2.5k
• updated
3.4 years ago by
Pierre Lindenbaum ♦ 133k
Latest awards to Niek De Klein
Great Question
15 months ago,
created a question with more than 5,000 views.
For How to remove non-unique mapping reads with Hisat?
Popular Question
15 months ago,
created a question with more than 1,000 views.
For biomaRt LPR region to position giving different results than I expect (using R)
Popular Question
20 months ago,
created a question with more than 1,000 views.
For biomaRt LPR region to position giving different results than I expect (using R)
Popular Question
2.1 years ago,
created a question with more than 1,000 views.
For How to remove non-unique mapping reads with Hisat?
Great Question
2.1 years ago,
created a question with more than 5,000 views.
For How To Convert .Mgf (Mascot Generic File) To A Format (Mzxml?) So That They Can Be Used With Maxquant?
Student
2.1 years ago,
asked a question with at least 3 up-votes.
For How To Convert .Mgf (Mascot Generic File) To A Format (Mzxml?) So That They Can Be Used With Maxquant?
Popular Question
2.1 years ago,
created a question with more than 1,000 views.
For biomaRt LPR region to position giving different results than I expect (using R)
Popular Question
2.1 years ago,
created a question with more than 1,000 views.
For biomaRt LPR region to position giving different results than I expect (using R)
Great Question
2.2 years ago,
created a question with more than 5,000 views.
For What Is A Good Phylogenetic Tree Display Program For Large Data Sets?
Popular Question
2.3 years ago,
created a question with more than 1,000 views.
For How to remove non-unique mapping reads with Hisat?
Popular Question
2.5 years ago,
created a question with more than 1,000 views.
For How to remove non-unique mapping reads with Hisat?
Scholar
3.2 years ago,
created an answer that has been accepted.
For A: How To Install The Ucsc In The Linux Operating System
Great Question
3.2 years ago,
created a question with more than 5,000 views.
For What Is A Good Phylogenetic Tree Display Program For Large Data Sets?
Student
3.2 years ago,
asked a question with at least 3 up-votes.
For How To Parse Psiblast Results Using Biopython And Blast-2.2.24+?
Great Question
3.2 years ago,
created a question with more than 5,000 views.
For When To Choose Python Over Perl (And Vice Versa)?
Popular Question
3.4 years ago,
created a question with more than 1,000 views.
For How to remove non-unique mapping reads with Hisat?
Teacher
3.5 years ago,
created an answer with at least 3 up-votes.
For A: Is Conversion From Colorspace To Base Space Lossy?
Scholar
4.2 years ago,
created an answer that has been accepted.
For A: How To Install The Ucsc In The Linux Operating System
Popular Question
4.2 years ago,
created a question with more than 1,000 views.
For Biopython Ncbistandalone Blastall Gives Different Result Than Calling Blastall Directly From Cmd
Popular Question
4.2 years ago,
created a question with more than 1,000 views.
For What Is A Good Phylogenetic Tree Display Program For Large Data Sets?
Popular Question
4.2 years ago,
created a question with more than 1,000 views.
For How To Show More Blast Results Using Biopython?
Popular Question
4.2 years ago,
created a question with more than 1,000 views.
For What Is (Or How To Find Out) The Meaning Of Error Code -1073741510 From Local Blast 2.2.25+?
Popular Question
4.2 years ago,
created a question with more than 1,000 views.
For How to remove non-unique mapping reads with Hisat?
Good Question
4.2 years ago,
asked a question that was upvoted at least 5 times.
For What Is A Good Phylogenetic Tree Display Program For Large Data Sets?
Prophet
4.2 years ago,
created a post with more than 20 followers.
For What Is A Good Phylogenetic Tree Display Program For Large Data Sets?
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 1882 users visited in the last hour