User: kamiljaron

gravatar for kamiljaron
kamiljaron70
Reputation:
70
Status:
Trusted
Location:
Last seen:
1 month, 2 weeks ago
Joined:
1 year, 9 months ago
Email:
k*********@gmail.com

Posts by kamiljaron

<prev • 15 results • page 1 of 2 • next >
0
votes
3
answers
418
views
3
answers
Comment: C: Kallisto pseudobam to IGV
... Hi, the author is here. I wrote this script for one use only (I was not trying to make a very general resource for everybody, but maybe I should). The problem is that you got differently formatted gtf file and I am using instead of some implemented parsers simply my own. These are lines parsing a g ...
written 11 weeks ago by kamiljaron70
0
votes
2
answers
23k
views
2
answers
Comment: C: Adding Read Groups To Bam Files
... I think this solution is invalid. 1. you need to escape @ character, so it would be at least : `perl -e 'print "\@RG\tID:ga\tSM:hs\tLB:ga\tPL:Illumina\n\@RG\tID:454\tSM:hs\tLB:454\tPL:454\n"' > rg.txt`. The second problem is that the line defined only the header and `-r` flag will assign readgrou ...
written 3 months ago by kamiljaron70
1
vote
4
answers
1.5k
views
4
answers
Comment: C: Snakemake vs. Nextflow: strengths and weaknesses
... And [BioMake][1] is off the game? It uses prolog (which is both the weakness and the strength...). [1]: https://academic.oup.com/bioinformatics/article-abstract/doi/10.1093/bioinformatics/btx306/3806980/BioMake-a-GNU-Make-compatible-utility-for?redirectedFrom=fulltext ...
written 4 months ago by kamiljaron70
0
votes
1
answer
543
views
1
answers
Comment: C: Big Genome Assembly of Medicinal Plant
... 1. How do you know your genome size? 2. Have you tried to estimate a genome size from reads? (for example using https://github.com/gmarcais/Jellyfish). ...
written 8 months ago by kamiljaron70
0
votes
1
answer
2.5k
views
1
answers
Comment: C: How to specify Illumina Single End data in the MaSuRCA Assembler config file
... This have not worked for me (version 11152016). Also, believe that two numbers after tag of the library are not avg/std_dev _read_length but average/std dev fragment length (i.e. the distance of the pair end reads). ...
written 11 months ago by kamiljaron70
0
votes
3
answers
674
views
3
answers
Comment: C: unpaired R1 and R2 sequences
... I am not sure If I understand correctly. Do you mean by "designations in the identifiers" just read IDs? If yes, then the fact that they are unpaired means that there will be always only read with unique ID in the concatenated file, therefore from the perspective of the software there can not be a d ...
written 13 months ago by kamiljaron70
1
vote
3
answers
674
views
3
answers
unpaired R1 and R2 sequences
... Hello, trimomatic reports unpaired R1 R2 after trimming pair-end reads (for those pairs where one of reads failed for some reason). Well, I kind of can not figure out, why it trimmomatic saves unpoired R1 and R2 separately. Is there any particular reason why not to concatenate those two files in on ...
trimmomatic trimming pair-end written 13 months ago by kamiljaron70 • updated 13 months ago by Brian Bushnell14k
0
votes
0
answers
762
views
0
answers
Comment: C: Population-wise decay of linkage disequilibrium (LD)
... It is not too helpful. Anyway, if you have some identification of populations in your file (lets say that CHR_A), than you can just use `subset` function `ld <- read.table("plink.ld",header=T)` `ld$a <- (ld$BP_B-ld$BP_A) #I recomed to use meaningful names, not "a"` `for(family in levels(ld$CH ...
written 19 months ago by kamiljaron70
0
votes
0
answers
762
views
0
answers
Comment: C: Population-wise decay of linkage disequilibrium (LD)
... May be you can solve it on the level of R... How does the plink.ld file looks like? ...
written 19 months ago by kamiljaron70
0
votes
0
answers
697
views
0
answers
Comment: C: pair-end, mate-pairs coding
... I found that this merged fastq files are called interleaved FASTQ (found at webpage of another trimmer [adapterremoval][1]). Using this, I ve googled a [bash script][2] which is probably solution for my problem, but not really answering a question - How it is coded. Is it just the order of sequences ...
written 19 months ago by kamiljaron70

Latest awards to kamiljaron

Appreciated 8 months ago, created a post with more than 5 votes. For C: Blast - Formatting Output
Commentator 19 months ago, created a comment with at least 3 up-votes. For C: Blast - Formatting Output

Help
Access

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