User: DNAngel
DNAngel • 20
- Reputation:
- 20
- Status:
- New User
- Location:
- Last seen:
- 4 days, 15 hours ago
- Joined:
- 1 year, 9 months ago
- Email:
- k********************@hotmail.com
Posts by DNAngel
<prev
• 32 results •
page 1 of 4 •
next >
2
votes
1
answer
78
views
1
answer
... I'm creating consensus sequences using a custom script but I am unsure if lowercase bases (the bases that are on the reverse strand) are already reverse complemented, or they need to be reverse complemented?
In IGV, when I see reads aligned to a position and I have a mix of A (forward reads), and ' ...
written 9 days ago by
DNAngel • 20
• updated
9 days ago by
Devon Ryan ♦ 87k
0
votes
2
answers
904
views
2
answers
... Hi, this sounds like something I need as well. Does this count the actual number of reads that contain a particular base or just how many reads aligned at that position? For example at position 1 if I know I have 100 reads aligned there, would it tell me I have 50 A's, 25 a, 25 C? Something like tha ...
written 10 days ago by
DNAngel • 20
0
votes
1
answer
809
views
1
answers
... Did you find out how to create the consensus sequences without IUPAC codes? I'm having the same problem using samtools and bcftools. ...
written 14 days ago by
DNAngel • 20
0
votes
0
answers
83
views
0
answers
... I can create consensus sequences now, I just cannot get them to not include IUPAC codes so I do think they are different questions. ...
written 14 days ago by
DNAngel • 20
0
votes
0
answers
83
views
0
answers
... Tried it but it didn't work no output consensus file would even be created for me to check. It definitely runs (no weird errors pop up) and something happens but it does not give me the output. I didn't understand how to use mosdepth, it seemed like I would have to relearn quite a bit when I was at ...
written 14 days ago by
DNAngel • 20
1
vote
0
answers
83
views
0
answers
... The only code that appears to work for me when I want to generate consensus sequences have been:
samtools mpileup -Q 20 -q 20 -d 5000 -uf ref.fa sorted.bam | bcftools call -c | vcfutils.pl vcf2fq -d 10 -Q 20 -l 1 > consensus_output.fq
However, I don't need ambiguity and almost every spot in ...
written 14 days ago by
DNAngel • 20
0
votes
0
answers
269
views
0
answers
... Yes that is me! Thank you I will try this as well. Also, I think I may have gotten it with bcftools using call as you suggested but using -d 10 and -Q 20 settings to ensure high quality base calls are kept. And that (to me right now at least) makes the consensus only contain the high quality base ca ...
written 27 days ago by
DNAngel • 20
0
votes
3
answers
2.7k
views
3
answers
Comment:
C: Ngs Question ~ Consensus
... Coming in 4 year later, I've tried these lines of code for my single-end BAM file but the consensus sequence fills in gaps with bases from the reference.fasta file. Is there a way to just get a consensus of the BAM file or even the mpileup file to also keep gaps where there was no depth of coverage ...
written 28 days ago by
DNAngel • 20
0
votes
0
answers
269
views
0
answers
... Hi, I know you mentioned to me earlier in another site that bcftools consensus will fill in gaps with bases from the reference...but I know from other scripts where colleagues have worked with paired-end data that this is not true. They were able to get the consensus sequence of their species1.bam f ...
written 28 days ago by
DNAngel • 20
0
votes
0
answers
269
views
0
answers
... I've ran into a new issue, now I cannot create a consensus sequence at all and I don't understand why.
I have my indexed (from bwa) reference files and my sorted alignment bam files for 30 species in the same directory. I am trying to make consensus sequences of them all. I am even trying to do thi ...
written 4 weeks ago by
DNAngel • 20
Latest awards to DNAngel
Scholar
5 weeks ago,
created an answer that has been accepted.
For A: How to save trimmomatic terminal summary to an output file?
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 1438 users visited in the last hour