SEQanswers

Go Back   SEQanswers > Search Forums


Showing results 1 to 25 of 500
Search took 0.47 seconds.
Search: Posts Made By: GenoMax
Forum: Illumina/Solexa Today, 03:16 AM
Replies: 1
Views: 27
Posted By GenoMax
Which ever side starts first likely grabs the...

Which ever side starts first likely grabs the next numeric number.
Forum: Bioinformatics Yesterday, 06:59 AM
Replies: 14
Views: 10,550
Posted By GenoMax
Use the index for blast. Nothing else needs to be...

Use the index for blast. Nothing else needs to be done beyond uncompressing files and making sure that all of them are in one folder. When you specify the index you only need to use the "basename" so...
Forum: Bioinformatics 03-23-2019, 04:51 AM
Replies: 14
Views: 10,550
Posted By GenoMax
Have you tried to use the index? It should be...

Have you tried to use the index? It should be fine to have one copy of nt.nal for sure. If there is anything wrong you will get an error.
Forum: Bioinformatics 03-21-2019, 11:27 AM
Replies: 12
Views: 485
Posted By GenoMax
I don't think you can split to a fasta format...

I don't think you can split to a fasta format file directly. Can you try following?

/home/hd55218/BBSplit/bbmap/bbsplit.sh -Xmx100g threads=2 in=/home/hd55218/BBSplit/QualTrimmed_bran11.fastq...
Forum: Bioinformatics 03-21-2019, 11:07 AM
Replies: 12
Views: 485
Posted By GenoMax
Now we have a different error. Is that...

Now we have a different error.



Is that file actually in fastq format and is it in that location?

Post the output from head -4 /home/hd55218/BBSplit/QualTrimmed_bran11.fastq
Forum: Illumina/Solexa 03-20-2019, 01:04 PM
Replies: 2
Views: 148
Posted By GenoMax
Not in immediate future AFAIK. Reagents should be...

Not in immediate future AFAIK. Reagents should be available into early 2020's. Someone may have saved that letter from Illumina and will provide specifics.
Forum: Bioinformatics 03-18-2019, 03:50 AM
Replies: 665
Views: 150,961
Posted By GenoMax
If your OS does not support shellscripts, replace...

If your OS does not support shellscripts, replace 'bbmap.sh' like this:
java -XmxNNg -cp /path/to/current align2.BBMap in=reads.fq out=mapped.sam

(NN will be a real number on your system).
Forum: Illumina/Solexa 03-13-2019, 06:49 AM
Replies: 2
Views: 241
Posted By GenoMax
It should basically be very similar to other...

It should basically be very similar to other Illumina sequencers. Quality scores on NovaSeq are binned by default so they would not be directly comparable. BBMap may be able to recalibrate the scores...
Forum: Bioinformatics 03-12-2019, 07:08 AM
Replies: 1
Views: 160
Posted By GenoMax
Cross-posted on Biostars, see my answer there:...

Cross-posted on Biostars, see my answer there: https://www.biostars.org/p/369000/#369008
Forum: General 03-11-2019, 06:37 AM
Replies: 5
Views: 224
Posted By GenoMax
If this is an assignment then use what you have...

If this is an assignment then use what you have to but comm should work (as long as you have enough RAM available). Since you are working with only read names (if you are not then you should).
Forum: General 03-11-2019, 06:06 AM
Replies: 5
Views: 224
Posted By GenoMax
You could simply get the names (field 1 as you...

You could simply get the names (field 1 as you already note, sort | uniq them in bash) and do a "comm" comparison of the three results. If your aim is just to find which reads are present in all...
Forum: Illumina/Solexa 03-07-2019, 03:59 AM
Replies: 4
Views: 464
Posted By GenoMax
If you use your MiSeq consistently (at least once...

If you use your MiSeq consistently (at least once a week) investing in the higher level contract with PM may be useful. Consider it cost of doing business and roll that into your operating costs.
Forum: RNA Sequencing 03-05-2019, 08:54 AM
Replies: 5
Views: 252
Posted By GenoMax
You can use DESeq2 without replicates...

You can use DESeq2 without replicates (https://support.bioconductor.org/p/101210/). It is not recommended. Your analysis would not have any statistical significance.

Hopefully you have some other...
Forum: RNA Sequencing 03-05-2019, 07:49 AM
Replies: 5
Views: 252
Posted By GenoMax
You will get raw counts from featureCounts. You...

You will get raw counts from featureCounts. You will need to use DESeq2/edgeR etc to actually do normalization and analysis.
Forum: RNA Sequencing 03-05-2019, 07:37 AM
Replies: 5
Views: 252
Posted By GenoMax
featureCounts part of Subread package...

featureCounts part of Subread package (http://subread.sourceforge.net/). You can feed all your BAM files at the same time to get a matrix of counts (rows as genes and columns as your samples).
Forum: RNA Sequencing 03-05-2019, 03:59 AM
Replies: 1
Views: 274
Posted By GenoMax
You can use Integrated Genome Viewer (Broad...

You can use Integrated Genome Viewer (Broad Institute) to do batch analyses/exports of images. Details are included here (https://software.broadinstitute.org/software/igv/batch).
...
Forum: Bioinformatics 03-04-2019, 05:53 AM
Replies: 9
Views: 439
Posted By GenoMax
Are you running this on a cluster/remote server?...

Are you running this on a cluster/remote server? Perhaps your process is encountering disk quota limit?
Forum: Bioinformatics 03-04-2019, 04:18 AM
Replies: 9
Views: 439
Posted By GenoMax
Are your data files very large? If the SAM file...

Are your data files very large? If the SAM file being made gets very large then you could possibly run into a space issue.
Forum: Bioinformatics 03-03-2019, 04:22 PM
Replies: 9
Views: 439
Posted By GenoMax
That is odd because the error says that hisat is...

That is odd because the error says that hisat is unable to read this file.


Warning: Could not open read file "/home/RNA-seq/Human/eynden/clb-bar-ctrl.fastq" for reading; skipping...
Error: No...
Forum: Bioinformatics 03-03-2019, 12:35 PM
Replies: 9
Views: 439
Posted By GenoMax
Are you trying to pass more than one set of files...

Are you trying to pass more than one set of files in your command line? Can you post your command line?
Forum: Illumina/Solexa 03-03-2019, 04:54 AM
Replies: 2
Views: 416
Posted By GenoMax
Have you consulted Illumina tech support? You are...

Have you consulted Illumina tech support? You are not doing something totally custom that is unsupported by Illumina?
Forum: Illumina/Solexa 02-28-2019, 02:25 PM
Replies: 3
Views: 483
Posted By GenoMax
Interesting. Let us know what happens.

Interesting. Let us know what happens.
Forum: General 02-28-2019, 02:25 PM
Replies: 2
Views: 364
Posted By GenoMax
If you know the barcode sequences you may be able...

If you know the barcode sequences you may be able to use bbduk.sh from BBMap suite to bin the reads you are interested in.
Forum: Illumina/Solexa 02-28-2019, 07:27 AM
Replies: 3
Views: 483
Posted By GenoMax
Sounds like bcl2fastq experienced a software...

Sounds like bcl2fastq experienced a software issue. I suggest that you re-run the demultiplexing. I have seen this posted rarely and if I recall had experienced it one time. bcl2fastq re-run fixed...
Forum: General 02-28-2019, 05:20 AM
Replies: 4
Views: 893
Posted By GenoMax
It depends. If you want a quick scan you could...

It depends. If you want a quick scan you could use something like fastq screen (https://www.bioinformatics.babraham.ac.uk/projects/fastq_screen/). You will need to select the database carefully.
...
Showing results 1 to 25 of 500

 


All times are GMT -8. The time now is 05:38 AM.


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