User: Pin.Bioinf
Pin.Bioinf • 260
- Reputation:
- 260
- Status:
- Trusted
- Location:
- Malaga
- Last seen:
- 7 months, 1 week ago
- Joined:
- 2 years, 6 months ago
- Email:
- i*************@gmail.com
Im a Health Engineer with a Bioinformatics specialization doing epigenetics research.
Posts by Pin.Bioinf
<prev
• 134 results •
page 1 of 14 •
next >
0
votes
0
answers
351
views
0
answers
...
I came up with this solution that a contributor in this page proposed to calculate the percentage of reference genome covered :
#Determine number of bases at 0 read depth
zero=$(bedtools genomecov -ibam BAM -g hg38.fasta -bga | awk '$4==0 {bpCountZero+=($3-$2)} {print bpCountZero}' | tail ...
written 8 months ago by
Pin.Bioinf • 260
0
votes
0
answers
351
views
0
answers
... As I asked in a similar question that was closed (thought it was different enough to be asked separately), could I determine the transcriptome size by multiplying uniquely mapped reads by read size?
Here is an example of the report
Number of input reads | 57456000
...
written 8 months ago by
Pin.Bioinf • 260
0
votes
0
answers
164
views
0
answers
... Hello,
I mapped my .fastq files against genome and I got the mapping quality reports. I have that 60M - 70M reads were uniquely mapped and that my reads were 100bp in size.
Can I say that **60M * 100bp (6 000 000 000 bp = 6000 Mb)** is the size of the transcriptome in that sample or am I complet ...
written 8 months ago by
Pin.Bioinf • 260
0
votes
0
answers
351
views
0
answers
... I got this:
1891663623, so these are bp? Do these numbers seem reasonable?
What I did:
bedtools genomecov -ibam input.bam -bg | bedtools merge -i - > covered_36.bed
cat covered_36.bed | awk -F'\t' 'BEGIN{SUM=0}{ SUM+=$3-$2 }END{print SUM}'
1891663623 ...
written 8 months ago by
Pin.Bioinf • 260
0
votes
0
answers
351
views
0
answers
... Thank you! could you help me with the command? I am new to this, how do I keep all regions with at least one read and sum up the size of the intervals with bedtools? Here is my idea:
bedtools genomecov -ibam input.bam -bg to get the regions with at least one read, and then how can I calculate tota ...
written 8 months ago by
Pin.Bioinf • 260
5
votes
0
answers
351
views
0
answers
... Hello,
I want to calculate TMB (mutational load), i know it will not be accurate because its RNA-Seq, but what is the target region generally for RNA-Seq? Similar to WES (30Mb)?
We sequenced with 100b length reads.
I need this to divide non synonymous mutations by this region size number.
Thank ...
written 8 months ago by
Pin.Bioinf • 260
0
votes
0
answers
300
views
0
answers
... Hello,
I want to calculate the Tumor mutational burden of my responders and non responders to treatment.
I have many question regarding this matter:
1. Should I calculate a TMB for responders and another TMB for non responders or only a global TMB for all samples?
2. **samtools mpileup , bcftoo ...
written 8 months ago by
Pin.Bioinf • 260
0
votes
0
answers
278
views
0
answers
... Hello, I have a normalized expression matrix that has many genes. I also have clinical data for the same samples of the matrix. Some data is numeric (age, levels of LDL, tumor size...), and some is categorical (sex, response to therapy, subtype of tumour...)
What kind of tests can I use to assess c ...
written 9 months ago by
Pin.Bioinf • 260
0
votes
0
answers
285
views
0
answers
... I would like to build a phylogenetic tree for around 40-50 patients with breast cancer based on SNP and CNAs I was asked for it by a fellow biologist, I am completely new to phylogenetics and I have no clue where to start. Do you have any recommendations on pipelines or brief steps I should follow ...
written 9 months ago by
Pin.Bioinf • 260
0
votes
1
answer
378
views
1
answers
... Yes, I agree. I was using default parameters, which were "complete" and "euclidean". I tried euclidean with ward.D2 and still, the differentiation between tumor and controls is not clear, it is a blur in the middle of both, and some tumors are the same color as controls as in my original heatmap. ...
written 10 months ago by
Pin.Bioinf • 260
Latest awards to Pin.Bioinf
Popular Question
9 months ago,
created a question with more than 1,000 views.
For Is it possible to find a human genome STAR index to not have to create it on my laptop?
Teacher
9 months ago,
created an answer with at least 3 up-votes.
For A: STAR --quantMode GeneCounts function
Centurion
12 months ago,
created 100 posts.
Supporter
14 months ago,
voted at least 25 times.
Scholar
14 months ago,
created an answer that has been accepted.
For A: STAR --quantMode GeneCounts function
Autobiographer
2.5 years ago,
has more than 80 characters in the information field of the user's profile.
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 1858 users visited in the last hour