SEQanswers

Go Back   SEQanswers > Search Forums


Showing results 1 to 25 of 500
Search took 0.06 seconds.
Search: Posts Made By: GenoMax
Forum: Bioinformatics Today, 04:51 AM
Replies: 11
Views: 10,446
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: 476
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: 476
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: 138
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,840
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: 227
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: 150
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: 204
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: 204
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: 452
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: 241
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: 241
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: 241
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: 263
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: 432
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: 432
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: 432
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: 432
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: 407
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: 471
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: 351
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: 471
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: 882
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.
...
Forum: Bioinformatics 02-26-2019, 04:42 AM
Replies: 327
Views: 98,858
Posted By GenoMax
@aushev: Unfortunately Brian no longer has time...

@aushev: Unfortunately Brian no longer has time to participate on this forum. He would really be the only person who can authoritatively answer your questions. You could try to create a ticket on SF...
Forum: Bioinformatics 02-25-2019, 04:11 PM
Replies: 327
Views: 98,858
Posted By GenoMax
Any reads that align to the ribosomal repeat will...

Any reads that align to the ribosomal repeat will be identified and separated in a file. Isn't that what you are looking to do?
Showing results 1 to 25 of 500

 


All times are GMT -8. The time now is 06:32 AM.


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