User: murathong

gravatar for murathong
murathong90
Reputation:
90
Status:
Trusted
Location:
China
Last seen:
2 weeks, 1 day ago
Joined:
1 year, 9 months ago
Email:
m********@gmail.com

Posts by murathong

<prev • 16 results • page 1 of 2 • next >
3
votes
2
answers
471
views
2
answers
Answer: A: Difference between GTF file with CHR and without CHR. ENSEMBL
... The one without 'chr' contains annotations for genes on unplaced or unlocalized contigs, while the one with 'chr' only contains annotation for assembled chromosomes, both of them have no prefix 'chr' in chromosome name, see this example: zcat Danio_rerio.GRCz10.87.gtf.gz | cut -f1 | awk '{dict[ ...
written 6 weeks ago by murathong90
0
votes
1
answer
261
views
1
answers
Comment: C: why allelic count on forward and reverse strand do not add up to depth in the re
... the quality of base calling for those reads is ~30 at that position. I guess this is not the main problem. ...
written 9 weeks ago by murathong90
0
votes
1
answer
261
views
1
answer
why allelic count on forward and reverse strand do not add up to depth in the results of samtools mpileup?
... dmel_mitochondrion_genome 2878 . A <*> 0 . DP=15;I16=0,9,0,0,347,13419,0,0,180,3600,0,0,140,2570,0,0;QS=1,0;MQ0F=0 PL:ADF:ADR 0,27,86:0,0:9,0 See the above line for example. Here the depth is DP=15, however allelic count on forward strand and reve ...
vcf samtools mpileup written 3 months ago by murathong90 • updated 3 months ago by Tom_L150
0
votes
1
answer
228
views
1
answers
Comment: C: how to adjust size of intervals in bed12 format?
... I didn't intend to slop each block. For example, if first block is 8nt in length, I want to discard this block and then remove 2 nt from next block so that I removed 8 + 2 nt = 10nt in total. Anyway, I have written a script for this: [bed12_slop.py][1]. I haven't test it on all possible extreme cond ...
written 3 months ago by murathong90
0
votes
1
answer
228
views
1
answers
Comment: C: how to adjust size of intervals in bed12 format?
... Hi. You should consider the size of each block (column 11). if the first block is shorter than 10, then the start should adjusted based on the original start and the start position of the next block (start posn of each block is in column 12). ...
written 3 months ago by murathong90
1
vote
1
answer
228
views
1
answer
how to adjust size of intervals in bed12 format?
... I have a file containing genomic intervals in bed12 format. I want to adjust the size of each interval (specifically, remove N base from both end of each interval). `bedtools slop` seems only works with bed6. Is there a convenient way to achive this for bed12 intervals? ...
bedtools bed alignment written 3 months ago by murathong90 • updated 3 months ago by Carlo Yague2.8k
0
votes
1
answer
851
views
1
answers
Answer: A: Implementing McDonald-Kreitman test
... I met the same problem recently and I assign assign equal weight to all possible mutation path. For example, for TTC to CTA, TTC->TTA->CTA and TTC->CTC->CTA have the same probability, as I am dealing triplet change in noncoding region. The following article talks more sophiscated method ...
written 4 months ago by murathong90
3
votes
1
answer
895
views
1
answers
Answer: A: get case submitter ID for GDC miRNA quantification files
... following the instruction on https://gdc-docs.nci.nih.gov/API/Users_Guide/Search_and_Retrieval I figured out how to download the uuid and barcode associated with each miRNA quantification file. first determine how many files there are ```curl 'https://gdc-api.nci.nih.gov/files/ids?query=mirnas.qua ...
written 8 months ago by murathong90
6
votes
1
answer
895
views
1
answer
get case submitter ID for GDC miRNA quantification files
... I downloaded hundreds of miRNA quantifications files from GDC and wanted to analyze its association with clinical information. The problem is how to get corresponding sumitter ID (previous TCGA barcode) of each file. This is an example: https://gdc-portal.nci.nih.gov/files/b2804bb2-70f4-471a-b6db- ...
uuid tcga gdc mirna written 8 months ago by murathong90
0
votes
1
answer
473
views
1
answers
Comment: C: where can I get the parameters used by UCSC genomes to run RepeatMasker for dm6?
... hi, following the archives by genome@soe.ucsc.edu, I find the answer after looking up similar posts. [how ucsc do RepeatMasker][1] my $defaultSpecies = 'scientificName from dbDb'; my $species = $opt_species ? $opt_species : &HgAutomate::getSpecies($dbHost, $db); Thanks for your help! ...
written 13 months ago by murathong90

Latest awards to murathong

Student 6 weeks ago, asked a question with at least 3 up-votes. For get case submitter ID for GDC miRNA quantification files
Teacher 6 weeks ago, created an answer with at least 3 up-votes. For A: get case submitter ID for GDC miRNA quantification files
Teacher 6 months ago, created an answer with at least 3 up-votes. For A: get case submitter ID for GDC miRNA quantification files

Help
Access

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