User: Niek De Klein

gravatar for Niek De Klein
Niek De Klein2.4k
Reputation:
2,420
Status:
Trusted
Location:
Netherlands
Last seen:
2 days, 13 hours ago
Joined:
8 years, 2 months ago
Email:
n**********@gmail.com

Posts by Niek De Klein

<prev • 252 results • page 1 of 26 • next >
0
votes
0
answers
130
views
0
answers
Comment: C: How can I fix picard error: Sequences at index 0 don't match
... 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 22 days ago by Niek De Klein2.4k
0
votes
0
answers
130
views
0
answers
Comment: C: How can I fix picard error: Sequences at index 0 don't match
... 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 23 days ago by Niek De Klein2.4k
0
votes
0
answers
130
views
0
answers
Comment: C: How can I fix picard error: Sequences at index 0 don't match
... 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 23 days ago by Niek De Klein2.4k
0
votes
0
answers
130
views
0
answers
How can I fix picard error: Sequences at index 0 don't match
... 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 ...
gatk picard written 24 days ago by Niek De Klein2.4k
0
votes
1
answer
405
views
1
answers
Comment: C: How to convert fastq to cram, and from cram back to fastq
... 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 9 weeks ago by Niek De Klein2.4k
0
votes
1
answer
405
views
1
answers
Comment: C: How to convert fastq to cram, and from cram back to fastq
... Thanks, I'm testing this now and seeing how it compares to removing secondary aligned reads ...
written 11 weeks ago by Niek De Klein2.4k
2
votes
1
answer
405
views
1
answer
How to convert fastq to cram, and from cram back to fastq
... 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 ...
fastq cram written 11 weeks ago by Niek De Klein2.4k • updated 11 weeks ago by finswimmer6.8k
2
votes
1
answer
533
views
1
answer
GaTK reverse select interval list
... 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_ ...
gatk intervallist written 14 months ago by Niek De Klein2.4k • updated 14 months ago by Pierre Lindenbaum114k
0
votes
1
answer
646
views
1
answers
Comment: C: Bedtools merge only part of feature that overlaps, not whole feature
... Thanks for the detailed explanation! For others that find this question, this is how I went from ensemble GTF to BED with CHR, start, stop, ensemble gene ID: awk '{if ($3 == "gene") print $0;}' Homo_sapiens.GRCh37.75.gtf \ | bedtools sort -i stdin \ | awk '{print $1 "\t" $4 "\t ...
written 15 months ago by Niek De Klein2.4k
1
vote
2
answers
502
views
2
answers
Answer: A: Get sequences from oligonucleotide probe ids
... The second column is an NCBI nucleotide ID, you can search the nucleotide database for the sequence. For example: https://www.ncbi.nlm.nih.gov/gquery/?term=XM_001122142.1 ...
written 15 months ago by Niek De Klein2.4k

Latest awards to Niek De Klein

Great Question 9 days 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 26 days ago, created a question with more than 1,000 views. For How to remove non-unique mapping reads with Hisat?
Popular Question 3 months ago, created a question with more than 1,000 views. For How to remove non-unique mapping reads with Hisat?
Student 11 months ago, asked a question with at least 3 up-votes. For How To Parse Psiblast Results Using Biopython And Blast-2.2.24+?
Scholar 11 months ago, created an answer that has been accepted. For A: How To Install The Ucsc In The Linux Operating System
Great Question 11 months ago, created a question with more than 5,000 views. For What Is A Good Phylogenetic Tree Display Program For Large Data Sets?
Great Question 11 months ago, created a question with more than 5,000 views. For When To Choose Python Over Perl (And Vice Versa)?
Popular Question 14 months ago, created a question with more than 1,000 views. For How to remove non-unique mapping reads with Hisat?
Teacher 16 months ago, created an answer with at least 3 up-votes. For A: Is Conversion From Colorspace To Base Space Lossy?
Popular Question 23 months ago, created a question with more than 1,000 views. For How To Show More Blast Results Using Biopython?
Popular Question 23 months ago, created a question with more than 1,000 views. For How to remove non-unique mapping reads with Hisat?
Good Question 23 months ago, asked a question that was upvoted at least 5 times. For What Is A Good Phylogenetic Tree Display Program For Large Data Sets?
Popular Question 23 months ago, created a question with more than 1,000 views. For What Is A Good Phylogenetic Tree Display Program For Large Data Sets?
Scholar 23 months ago, created an answer that has been accepted. For A: How To Install The Ucsc In The Linux Operating System
Prophet 23 months ago, created a post with more than 20 followers. For What Is A Good Phylogenetic Tree Display Program For Large Data Sets?
Popular Question 2.0 years ago, created a question with more than 1,000 views. For What Is A Good Phylogenetic Tree Display Program For Large Data Sets?
Teacher 2.3 years ago, created an answer with at least 3 up-votes. For A: Is Conversion From Colorspace To Base Space Lossy?
Popular Question 2.3 years ago, created a question with more than 1,000 views. For How To Show More Blast Results Using Biopython?
Appreciated 2.3 years ago, created a post with more than 5 votes. For A: Why Gene Regulatory Network Is Sparse?
Popular Question 2.5 years ago, created a question with more than 1,000 views. For Windows Cmd Not Responding With Python Blastall 2.2.17
Popular Question 2.5 years ago, created a question with more than 1,000 views. For When To Choose Python Over Perl (And Vice Versa)?

Help
Access

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