SEQanswers

Go Back   SEQanswers > Search Forums


Showing results 1 to 25 of 48
Search took 0.00 seconds.
Search: Posts Made By: ssully
Forum: Bioinformatics 04-26-2019, 01:47 PM
Replies: 1
Views: 995
Posted By ssully
multiple rounds of Pilon polishing: need new bam file?

Pilon takes an assembly and one or more sorted .bam files of reads, as input.

When doing more than one round of Pilon polishing on an assembly, for the second round, do I first need to create a...
Forum: Bioinformatics 03-21-2019, 12:04 PM
Replies: 0
Views: 690
Posted By ssully
trimming fastq records based on a coordinates file

I have some fastq files that came with 'clip' files with instructions for quality trimming; a table of lines that are:

read ID -- left clip coordinate --right clip coordinate


I'd like to...
Forum: De novo discovery 03-20-2019, 10:25 AM
Replies: 1
Views: 3,141
Posted By ssully
genome assembly polisher that can use Sanger reads?

We have a PacBio final assembly that requires some polishing to fix e.g. erroneous stop codons. No Illumina reads at hand but we do have a set of high quality Sanger reads from an older version of...
Forum: Bioinformatics 03-17-2019, 06:41 PM
Replies: 678
Views: 199,861
Posted By ssully
How do I run the basic BBmap.sh script in...

How do I run the basic BBmap.sh script in Windows? What is the corresponding java command? I want to align four sets of paired-end Illumina RNA-Seq reads to a genome assembly. I am particularly...
Forum: RNA Sequencing 03-04-2019, 02:25 PM
Replies: 1
Views: 1,340
Posted By ssully
visualizing RNAseq mapping to reference

I'm using Geneious (for Windows), which has two RNA Seq mappers (Geneious RNA Mapper, and BBmap). I'm looking for introns in a set of individual gene sequences. Mapping results in Geneious are...
Forum: General 08-25-2016, 07:33 AM
Replies: 6
Views: 3,327
Posted By ssully
You are getting the 10 best target sequences for...

You are getting the 10 best target sequences for each query. For some of the queries, there is more than one hit (hsp) to the target. Those are returned too, so sometimes there are more than 10...
Forum: Illumina/Solexa 12-06-2014, 06:38 AM
Replies: 35
Views: 14,987
Posted By ssully
I worked out the correct orientation and order of...

I worked out the correct orientation and order of 454 paired reads input for SPAdes, and have corrected the reads with --iontorrent option (ionhammer). Btu now I have questions regarding ionhammer...
Forum: Illumina/Solexa 12-03-2014, 07:52 AM
Replies: 35
Views: 14,987
Posted By ssully
I don't know; the second variant seems to be...

I don't know; the second variant seems to be saying to me , 'the reads from the right side of the library read (post-linker, 454_2.fastq) belong at the right end of the genome fragment' -- which...
Forum: Illumina/Solexa 12-02-2014, 05:43 PM
Replies: 35
Views: 14,987
Posted By ssully
I have removed the linkers and split the 454 mate...

I have removed the linkers and split the 454 mate pair reads with sff_extract; I have them now as (after deinterlacing) a pair of fastq files (454_1.fastq and 454_2.fastq) containing reads _1 and...
Forum: Illumina/Solexa 12-01-2014, 02:38 PM
Replies: 35
Views: 14,987
Posted By ssully
But 454 paired end reads are two 'end' reads...

But 454 paired end reads are two 'end' reads connected by a linker sequence. Does the IonHammer corrector actually recognize those and split the reads before correcting? Or do the 454 PE reads...
Forum: Illumina/Solexa 12-01-2014, 09:54 AM
Replies: 35
Views: 14,987
Posted By ssully
re: attempting SPAdes with 454 reads -- if I...

re: attempting SPAdes with 454 reads -- if I understand correctly:

2) SPAdes does not natively support 454 reads

2) 454 reads resemble Ion Torrent reads (similar technology) but SPAdes will...
Forum: De novo discovery 11-21-2014, 09:20 AM
Replies: 8
Views: 2,924
Posted By ssully
That was the wrong way to do it -- instead I've...

That was the wrong way to do it -- instead I've joined the MIRA user forum, which is what the author recommends.
Forum: De novo discovery 11-19-2014, 02:31 PM
Replies: 8
Views: 2,924
Posted By ssully
The example in the mulitple platform manifest...

The example in the mulitple platform manifest *seems* to indicate the 454 PEs can be left as one fastq file. The insert size and SD need to be provided (which I can do) or 'autorefine' to let MIRA...
Forum: De novo discovery 11-19-2014, 01:38 PM
Replies: 8
Views: 2,924
Posted By ssully
That helps a lot (for MIRA) , thanks! So as...

That helps a lot (for MIRA) , thanks!

So as long as I convert the 454 Paired End sff to Fastq with sff_extract, MIRA will process and assemble them correctly *as paired ends*? It still...
Forum: De novo discovery 11-19-2014, 11:19 AM
Replies: 8
Views: 2,924
Posted By ssully
de novo assembly including Illumina and 454 paired-end reads

There's a bunch of hybrid de novo assembly threads but not sure I'm finding the answers to these questions, specifically about paired-end reads from two platforms (Illumina and 454).

454's 'paired...
Forum: Illumina/Solexa 11-18-2014, 12:47 PM
Replies: 145
Views: 318,987
Posted By ssully
Illumina helped sort me out. The seqs I posted...

Illumina helped sort me out. The seqs I posted last time are from an old , discontinued version of the kit (pp 20-21 of the July 2014 letter). The newer TruSeq adapters('universal' and 'indexed')...
Forum: Illumina/Solexa 11-17-2014, 06:18 PM
Replies: 145
Views: 318,987
Posted By ssully
I'm quite lost. I've got a set of paired end 100...

I'm quite lost. I've got a set of paired end 100 nt HiSeq reads (library insert size was ~300bp) from a colleague and I can't determine what adapter sequences to trim. The current Illumina letter ...
Forum: Bioinformatics 11-07-2014, 09:54 AM
Replies: 11
Views: 5,168
Posted By ssully
I'm trying (on Windows 7) to use Next...

I'm trying (on Windows 7) to use Next Workbench 3's SFF editor on some 454 sff files, but when I do end clipping and filtering, none of it is saved and the sequence itself never changes even...
Forum: Bioinformatics 10-06-2014, 01:26 PM
Replies: 14
Views: 4,488
Posted By ssully
mview (http://bio-mview.sourceforge.net/)

mview (http://bio-mview.sourceforge.net/)
Forum: Bioinformatics 10-02-2014, 01:18 PM
Replies: 14
Views: 4,021
Posted By ssully
When I wrote 'I don't have 140 Gb of RAM' I meant...

When I wrote 'I don't have 140 Gb of RAM' I meant on my workstation (it has 32Gb). But I do have access to a cluster where I can allocate >140 gb to a job . Should I try that or are you suggesting...
Forum: Bioinformatics 10-01-2014, 11:17 AM
Replies: 14
Views: 4,021
Posted By ssully
Hmm..the interleaved file is 140 million reads. ...

Hmm..the interleaved file is 140 million reads. I don't have 140 Gb of RAM. Guess I have to run this on a unix cluster. Is the corrected package (for running the dedup command) ready for download?
Forum: Bioinformatics 09-30-2014, 01:35 PM
Replies: 14
Views: 4,021
Posted By ssully
That works! Thanks.

That works! Thanks.
Forum: Bioinformatics 09-30-2014, 09:08 AM
Replies: 14
Views: 4,021
Posted By ssully
tried it but 'java' was not recognized as a...

tried it but 'java' was not recognized as a command. I added the path to 'my' java.exe
and reran, but then it threw this error

Error: Could not find or load main class jgi.ReformatReads
...
Forum: Bioinformatics 09-29-2014, 04:49 PM
Replies: 14
Views: 4,021
Posted By ssully
I installed BBmap on my Windows 7 machine (Java...

I installed BBmap on my Windows 7 machine (Java RE 1.7 is installed). But I can't seem to get dedup to read any of my input files. They're standard paired fastq files I named test1.f1 and...
Forum: Bioinformatics 09-29-2014, 09:05 AM
Replies: 2
Views: 2,119
Posted By ssully
the introduction of this pape...

the introduction of this pape (http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4027187/)r contains nice brief descriptions of many of the softwares used to find and analyse genomic repeats.
Showing results 1 to 25 of 48

 


All times are GMT -8. The time now is 03:19 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO