User: oconnwald
oconnwald • 20
- Reputation:
- 20
- Status:
- New User
- Location:
- Boulder, Colorado
- Last seen:
- 6 days, 13 hours ago
- Joined:
- 9 months, 2 weeks ago
- Email:
- a*************@yahoo.com
Posts by oconnwald
<prev
• 9 results •
page 1 of 1 •
next >
0
votes
2
answers
6.3k
views
2
answers
Comment:
C: bwa mem and multiply mapped reads
... Does anyone know where:
> When a read matches in its entirety, with an equal score in multiple
> locations, one of the locations is picked at random, is labeled as
> primary, will be given a mapping quality of zero and will have an XA
> tag that contains the alternative locations (this ...
written 6 weeks ago by
oconnwald • 20
0
votes
11
answers
54k
views
11
answers
... What about:
rand_records = random.sample(range(0, records), N)
rand_records.sort()
random.sample performs random selection without replacement, so the check you've added above would just be built-in to the call. ...
written 3 months ago by
oconnwald • 20
0
votes
0
answers
260
views
0
answers
... As an update, I've run the random-cased fastqs alignment with our reference modified to not include masking (i.e. I capitalized all lower-case bases), and I'm getting equivalent results as I did when I ran alignment of the fastqs with the reference with masking. So, I think the casing is here makin ...
written 5 months ago by
oconnwald • 20
0
votes
0
answers
260
views
0
answers
... That's interesting -- but the ordering of the reads isn't changing in each run (only casing), so I don't know if the first link applies. We are using a reference that has masked bases though. ...
written 5 months ago by
oconnwald • 20
0
votes
0
answers
260
views
0
answers
... Total reads includes supplemental reads -- you can see the offset on both is +12. If there were variation in the number of raw reads in the same, you'd see that reflected in the read1 and read2 counts ...
written 5 months ago by
oconnwald • 20
0
votes
0
answers
260
views
0
answers
... Just a simple Python script (I'm happy to post it if you'd like). I've confirmed that the random casing reads are identical in content (other than casing) by converting the random-cased bases back to the original upper casing and running a diff. ...
written 5 months ago by
oconnwald • 20
2
votes
0
answers
260
views
0
answers
... I was wondering how the casing of base calls affected BWA MEM performance.
I ran three sets of data to see if there was any difference at all:
- All base calls capitalized
- All base calls lower-cased
- Random casing for all base calls
I found that the first two sets (when casing was consistent) ...
written 5 months ago by
oconnwald • 20
0
votes
0
answers
212
views
0
answers
... When I use Samtools (v1.7 or v.1.10) mpileup with -l option to specify the positions I want, many REF bases are incorrect. When I use the -r command (in addition to -l), however, I get the correct bases.
For instance, I use these two commands to test consistency between using the -l vs -r options ...
written 7 months ago by
oconnwald • 20
1
vote
11
answers
167k
views
11
answers
... You can install the executable from their github release page: https://github.com/brentp/mosdepth/releases
Latest version at time of commenting (0.2.9):
(Download executable)
`wget https://github.com/brentp/mosdepth/releases/download/v0.2.9/mosdepth`
(Update permissions to run it)
`chmod a+x m ...
written 9 months ago by
oconnwald • 20
Latest awards to oconnwald
No awards yet. Soon to come :-)
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 1799 users visited in the last hour