User: rizoic

gravatar for rizoic
rizoic190
Reputation:
190
Status:
Trusted
Location:
Last seen:
20 minutes ago
Joined:
8 months, 1 week ago
Email:
r********@gmail.com

Posts by rizoic

<prev • 12 results • page 1 of 2 • next >
1
vote
1
answer
468
views
1
answers
Comment: C: How to remove singletons in a BAM file using samtools
... In your e.g. R1 is the singleton read since its mate is unmapped and thus R1 read will be removed when you filter for the flag `0x08`. ...
written 4 months ago by rizoic190
3
votes
1
answer
468
views
1
answers
Answer: C: How to remove singletons in a BAM file using samtools
... As per the [SAM flag guide][1] the flag that you are filtering for i.e. `0x04` is for unmapped reads. This explains the removal of unmapped reads. In case you want to filter for singletons i.e. reads for which the `mate is unmapped` you an select that as the condition. This gives a SAM flag value ...
written 4 months ago by rizoic190
3
votes
2
answers
281
views
2
answers
Answer: C: Filter DNA sequence in FASTA file by minimum length and redirect sequences that
... In case it is a requirement to do the solution without any additional packages then you can use the solution by Heather. However when parsing standard file formats it is always recommended to use specialized packages for handling them. You can use the `Biopython` package in python for this. An examp ...
written 4 months ago by rizoic190
0
votes
1
answer
294
views
1
answers
Comment: C: How to separate a strain of MTB drug wise?
... For Pre-XDR as per the condition you mention i.e. `Isoniazid, Rifampicin and (CIPROFLOXACIN|MOXIFLOXACIN|OFLOXACIN|LEVOFLOXACIN) any of these or one or two depends.` the command should be the same as in the answer right? i.e. tail -n +2 data.txt |\ awk '{a[$1] = a[$1]","$2};END{for(i in a)p ...
written 4 months ago by rizoic190
2
votes
2
answers
327
views
2
answers
Comment: C: What is the best way to compare transcriptome between different condition?
... If you are interested in just knowing which of the knockdowns i.e. B or C is close to lets say A, you can do hierarchical clustering on the counts post applying a transform like `vst() or rld()` in DEseq2. You can find an example [here][1]. [1]: http://people.duke.edu/~ccc14/duke-hts-2017/Statis ...
written 5 months ago by rizoic190
1
vote
1
answer
294
views
1
answers
Comment: C: How to separate a strain of MTB drug wise?
... A simplistic tweak would be to change the condition on these lines tail -n +2 data.txt |\ awk '{a[$1] = a[$1]","$2};END{for(i in a)print i"\t"a[i]}'|\ grep -P "\t,ISONIAZID,RIFAMPICIN$" Here the -P makes grep evaluate the expression as a perl compatible regex.The `\t and $`(indicating ...
written 6 months ago by rizoic190
0
votes
1
answer
358
views
1
answers
Comment: C: Error: Snakemake skiped a rule
... Changing rule bowtie2Aln: input: rules.trimmen.output To rule bowtie2Aln: input: rules.trimmen.output, rules.bowtie2Build.output Should solve this error for you. It would be better to state an explicit dependency than specifying a priority in this case. You can use the `--dry ...
written 6 months ago by rizoic190
0
votes
1
answer
358
views
1
answers
Answer: A: Error: Snakemake skiped a rule
... One of the problems I can see is that you have to add the output of rule `bowtie2Build` as input for rule `bowtie2Aln`. This way snakemake understands that it should perform the alignment only after the generation of the bowtie index. What it might do right now is execute the alignment step even b ...
written 6 months ago by rizoic190
4
votes
1
answer
294
views
1
answers
Answer: A: How to separate a strain of MTB drug wise?
... Incase you need a grep/awk based solution based on your tags. You can run something on these lines. If you paste the data into a file lets say `data.txt`. Example command for ` Pre-XDR` tail -n +2 data.txt |\ awk '{a[$1] = a[$1]","$2};END{for(i in a)print i"\t"a[i]}'|\ grep ISONIAZID| ...
written 6 months ago by rizoic190
3
votes
1
answer
353
views
1
answers
Answer: C: unable to load quant.sf files
... Can you add the output when you enter `files` on the terminal. You should not have both double/single quotes when you initialize files. Try removing the single quotes and changing the statement where you initialize files to `files <- c("/home/dell/Desktop/quant/quant3.sf", {same for other files}) ...
written 6 months ago by rizoic190

Latest awards to rizoic

Teacher 4 months ago, created an answer with at least 3 up-votes. For C: unable to load quant.sf files
Teacher 4 months ago, created an answer with at least 3 up-votes. For C: unable to load quant.sf files
Supporter 6 months ago, voted at least 25 times.
Scholar 6 months ago, created an answer that has been accepted. For C: unable to load quant.sf files
Teacher 6 months ago, created an answer with at least 3 up-votes. For C: unable to load quant.sf files
Scholar 6 months ago, created an answer that has been accepted. For C: unable to load quant.sf files
Teacher 6 months ago, created an answer with at least 3 up-votes. For C: unable to load quant.sf files

Help
Access

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