User: manuel.belmadani
manuel.belmadani • 30
- Reputation:
- 30
- Status:
- New User
- Location:
- Canada
- Twitter:
- pragmatwit
- Scholar ID:
- Google Scholar Page
- Last seen:
- 7 hours ago
- Joined:
- 2 years, 9 months ago
- Email:
- m***************@gmail.com
Posts by manuel.belmadani
<prev
• 6 results •
page 1 of 1 •
next >
0
votes
3
answers
94
views
3
answers
... It seems like the initial report of those flanking CTCF sites are from this [paper][1]:
"Recent studies have identified binding sites for CTCF as essential components of vertebrate insulator elements [21,22,23,24]. **CTCF-binding sites are about 50 bp and variable**, possibly because CTCF can use ...
written 5 days ago by
manuel.belmadani • 30
1
vote
1
answer
120
views
1
answers
... I got a reply from an SRA curator:
> The PATH2 is intended to be the full SRA Run accession, and is not
> restricted to a character limit.
So it appears that doing `PATH2=${SRR:0:10}` is not necessary or requested by any convention, and would probably indeed break if 11-character long identi ...
written 11 days ago by
manuel.belmadani • 30
0
votes
1
answer
120
views
1
answers
... Thanks for the reply!
I agree that what NCBI choses to do is speculative, however, wonderdump.sh *explicitly* truncates the SRR portion at 10 characters to create PATH2, when currently there's no difference between doing this and returning the full SRR. I'm assuming there's a good reason for this p ...
written 13 days ago by
manuel.belmadani • 30
2
votes
1
answer
120
views
1
answer
... I've been using [wonderdump.sh][1] from the Biostars handbook for some time. I'm now curious about the part that builds the ftp site url:
PATH1=${SRR:0:6}
PATH2=${SRR:0:10}
URL="ftp://ftp-trace.ncbi.nih.gov/sra/sra-instant/reads/ByRun/sra/SRR/${PATH1}/${PATH2}/${SRR}.sra"
I've seen SRR ...
written 13 days ago by
manuel.belmadani • 30
2
votes
3
answers
183
views
3
answers
... Keep them separated, like Nicolas said, most modern NGS software should handle paired-end reads. One they're aligned, you should have a single SAM/BAM file containing reads from both ends.
Depending on your purpose, you may need to choose different tools. For variant calling using whole-genome-sequ ...
written 8 weeks ago by
manuel.belmadani • 30
1
vote
1
answer
1.2k
views
1
answers
Answer:
A: Annovar using VCF file
... As @Sam pointed out, you're missing the -out before myanno. The -vcfinput argument (which you had in the second command but not your first) is also necessary, and omitting it wouldn't cause the "Syntax error" message, but instead all your variants would get discarded (and written to myanno.invalid_ ...
written 2.4 years ago by
manuel.belmadani • 30
Latest awards to manuel.belmadani
Scholar
11 days ago,
created an answer that has been accepted.
For A: wonderdump.sh and FTP site URL conventions for SRR identifiers
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 964 users visited in the last hour