User: peri.tobias

gravatar for peri.tobias
peri.tobias10
Reputation:
10
Status:
New User
Location:
Last seen:
5 days, 7 hours ago
Joined:
3 years, 9 months ago
Email:
p**********@gmail.com

Posts by peri.tobias

<prev • 18 results • page 1 of 2 • next >
0
votes
0
answers
777
views
0
answers
Comment: C: RATT annotation transfer tool usage.
... Hello Mehmet, I am late to this thread, but would appreciate some guidance please if you are able to help. I have spent the whole day trying to get this software running after installing PAGIT locally and also using EMBLmyGFF3 to build .embl from augustus output. I try to run this script: `start. ...
written 7 weeks ago by peri.tobias10
0
votes
0
answers
237
views
0
answers
Why is the pacbio bam 10X the original bam file size?
... Running the genomicconsensus arrow on aligned raw reads against a de novo pacbio assembly (1 Gbase assembly). Used minimap2 for the aligned.bam and then pbbamify to format for use with arrow. Can't find information on what this reformatting does but the initial aligned.bam = 35G and the pb_align ...
arrow pbbamify genomicconsensus pacbio written 9 months ago by peri.tobias10
0
votes
3
answers
2.5k
views
3
answers
Comment: C: Alternative to BLASR ?
... Thanks for these tips @harish. I did not know that I could directly use subreads with Pbmm2. Actually this all looks more promising as a way forward. Pbbamify worked after 130 hours of walltime. It output a very large pb.bam (351G) from aligned.bam (36G). Arrow still needs a pbi indexed file accord ...
written 10 months ago by peri.tobias10
0
votes
3
answers
2.5k
views
3
answers
Comment: C: Alternative to BLASR ?
... Yes, you are correct. I should have been clearer. This is the script I used. minimap2 -ax map-pb genome.fasta *.fastq \ | samtools view -hF 256 - \ | samtools sort -@ 8 -o genome_aligned.bam I should also say that I am still having some trouble with pbbamify. So I haven't got this whole step wor ...
written 10 months ago by peri.tobias10
1
vote
3
answers
2.5k
views
3
answers
Comment: C: Alternative to BLASR ?
... Hello and thanks for all these suggestions. These further comments might help someone else. I have used minimap2 to map PacBio raw reads (converted to fastq as required for minimap2) against my de novo assembly. I recommend this approach over blasr. Use samtools to get the aligned.bam. Minimap2 doe ...
written 10 months ago by peri.tobias10
0
votes
1
answer
533
views
1
answer
inputs to genomicconsensus arrow algorithm
... Firstly apologies, this is a cross-post from here as I was not sure if I had sent to the correct forum. If I get a useful answer I will make sure it is across both platforms. https://github.com/PacificBiosciences/pbcore/issues/118 I have assembled a de novo genome (1.98 Gb) with canu v1.6 using pa ...
assembly polish arrow pacbio written 11 months ago by peri.tobias10 • updated 6 months ago by harish230
0
votes
2
answers
4.4k
views
2
answers
Comment: C: BWA-MEM using long PacBio reads
... thanks, appreciated! ...
written 20 months ago by peri.tobias10
0
votes
2
answers
4.4k
views
2
answers
Comment: C: BWA-MEM using long PacBio reads
... Thanks very much for your comment, genomax. I initially used dextract to output fasta, quiver and arrow files from bam and hdf5. https://dazzlerblog.wordpress.com/tag/arrow/ Quiver files are fastq-like but not exactly fastq format - though a simple conversion is probably possible. I intend to polis ...
written 20 months ago by peri.tobias10
0
votes
2
answers
4.4k
views
2
answers
Comment: C: BWA-MEM using long PacBio reads
... Just wanted to add that bbmap only seems able to take one in=file. My above trials at listing files does not work. Perhaps you could concatenate but my fasta files are huge so not going to risk. Will run them individually and output multiple unmapped and mapped .fasta files. ...
written 20 months ago by peri.tobias10
0
votes
2
answers
4.4k
views
2
answers
Comment: C: BWA-MEM using long PacBio reads
... another way to list all the fasta files is: mapPacBio.sh in=`echo ./foo/*.fasta` outu=unmapped.fasta outm=mapped.fasta maxlen=6000 ...
written 20 months ago by peri.tobias10

Latest awards to peri.tobias

Popular Question 5 months ago, created a question with more than 1,000 views. For Making bed files from fasta
Popular Question 15 months ago, created a question with more than 1,000 views. For Making bed files from fasta
Popular Question 2.5 years ago, created a question with more than 1,000 views. For Making bed files from fasta
Scholar 3.4 years ago, created an answer that has been accepted. For A: Making bed files from fasta

Help
Access

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