SEQanswers

Go Back   SEQanswers > Search Forums


Showing results 1 to 25 of 500
Search took 0.03 seconds.
Search: Posts Made By: mastal
Forum: Metagenomics 05-22-2018, 03:01 PM
Replies: 2
Views: 858
Posted By mastal
I think they look OK. The thing with 16S reads is...

I think they look OK. The thing with 16S reads is that they are amplicons, and that all the
sequences are very similar.

I had seen a document on the website of an Italian company named IGATech...
Forum: Bioinformatics 04-05-2018, 01:23 PM
Replies: 3
Views: 440
Posted By mastal
Are you sure that velveth produced the expected...

Are you sure that velveth produced the expected output?
What files did you get in the output directory after velveth ran?

If you think it didn't install properly the first time, you could...
Forum: Bioinformatics 02-19-2018, 02:25 PM
Replies: 2
Views: 473
Posted By mastal
A very low alignment rate could be a sign that...

A very low alignment rate could be a sign that the reads are not from the species you expected.
Forum: Bioinformatics 12-18-2017, 02:35 AM
Replies: 9
Views: 1,324
Posted By mastal
What does the ~ expand to on your computer? ...

What does the ~ expand to on your computer?

If it's something like /users/kay/ instead of /home/kay/
that might explain why you get the directory/file not found error.
Forum: Bioinformatics 12-17-2017, 01:32 PM
Replies: 9
Views: 1,324
Posted By mastal
The error message suggests that you have too many...

The error message suggests that you have too many parameters, and trimmomatic is expecting a trimming command at that position, it doesn't really care too much what you call your file. Do you have...
Forum: Bioinformatics 11-13-2017, 01:39 AM
Replies: 6
Views: 1,194
Posted By mastal
Have you tried running trimmomatic with only the...

Have you tried running trimmomatic with only the ILLUMINACLIP command and not the rest of the quality trimming, just to see how the Illuminaclip part works?

Are you sure you are using the right...
Forum: Bioinformatics 11-12-2017, 02:51 PM
Replies: 6
Views: 1,194
Posted By mastal
Your command looks OK, why do you think it's not...

Your command looks OK, why do you think it's not working?

Is trimmomatic running, and what output do you get when it finishes?
Forum: Bioinformatics 11-12-2017, 07:25 AM
Replies: 6
Views: 1,194
Posted By mastal
You have 'PE' in your command, but you appear to...

You have 'PE' in your command, but you appear to be trying to trim SE data, because you are providing the names of one input file and one output file only, so you should change that to 'SE'.
Forum: RNA Sequencing 09-10-2017, 07:59 AM
Replies: 5
Views: 1,180
Posted By mastal
I haven't used hisat2, but with tophat, not all...

I haven't used hisat2, but with tophat, not all of the stages of the analysis were multi-threaded. So unless you look at top while hisat2 is running one of the multi-threaded stages, you might only...
Forum: Bioinformatics 06-03-2017, 10:46 AM
Replies: 1
Views: 865
Posted By mastal
You don't need to run 'make' if you download the...

You don't need to run 'make' if you download the binary, only if you are trying to build from the source code.

What version of HISAT2 are you using?

I just downloaded the Linux binary for...
Forum: Bioinformatics 05-20-2017, 03:58 AM
Replies: 4
Views: 888
Posted By mastal
I wouldn't worry about the two peaks in the per...

I wouldn't worry about the two peaks in the per sequence quality scores although it looks a bit unusual.

Are the reads all the same length if it is Illumina data, or have they already been trimmed...
Forum: Bioinformatics 05-20-2017, 02:16 AM
Replies: 10
Views: 808
Posted By mastal
You could get a rough estimate, although it could...

You could get a rough estimate, although it could turn out wrong. Use the genome sizes of the most closely related bird species with known genomes as a guide.

You said you had made several...
Forum: Bioinformatics 05-20-2017, 01:56 AM
Replies: 6
Views: 758
Posted By mastal
Yes. that is correct.

Yes. that is correct.
Forum: Bioinformatics 05-19-2017, 09:28 AM
Replies: 6
Views: 758
Posted By mastal
I meant that you should check whether, for genes...

I meant that you should check whether, for genes on the - strand, the exon closest to the rightmost end of the gene is labelled as exon1 or not.

Usually the chromosomal coordinates for a gene are...
Forum: Illumina/Solexa 05-19-2017, 09:20 AM
Replies: 5
Views: 728
Posted By mastal
It's data from the SOLiD platform. They use...

It's data from the SOLiD platform.

They use colorspace rather than basespace,
so you get a sequence of 0,1,2,3 rather than A,T,G,C.

The . indicate missing base calls.
Forum: Bioinformatics 05-19-2017, 06:12 AM
Replies: 10
Views: 808
Posted By mastal
What assembler or assemblers have you used, and...

What assembler or assemblers have you used, and what coverage do you have? For some assemblers too high coverage also leads to problems.
Forum: Bioinformatics 05-19-2017, 05:54 AM
Replies: 6
Views: 758
Posted By mastal
In your gtf file, compare the exon annotations of...

In your gtf file, compare the exon annotations of genes that are on the + strand with those that are on the - strand. The gene you have shown above is on the + strand.
Forum: Bioinformatics 05-19-2017, 05:49 AM
Replies: 10
Views: 808
Posted By mastal
Can you try reference-guided assembly with a...

Can you try reference-guided assembly with a related bird genome?
Forum: Bioinformatics 05-19-2017, 05:19 AM
Replies: 10
Views: 808
Posted By mastal
That actually looks OK, as long as the %GC...

That actually looks OK, as long as the %GC content is what you expect for that genome.
Forum: Bioinformatics 05-19-2017, 04:47 AM
Replies: 10
Views: 808
Posted By mastal
What do the other FastQC plots look like, the per...

What do the other FastQC plots look like, the per base sequence quality, and the per base sequence content?

Have you done any adapter or quality trimming before doing the assembly?
Forum: Bioinformatics 04-29-2017, 02:15 PM
Replies: 189
Views: 58,287
Posted By mastal
Why do you think this indicates failure? You...

Why do you think this indicates failure? You asked trimmomatic to keep both reads of a pair when it finds adapter in palindrome mode, and it appears to have done that.
This does not mean that...
Forum: Bioinformatics 04-29-2017, 01:59 PM
Replies: 189
Views: 58,287
Posted By mastal
Unless you are using a very old version of...

Unless you are using a very old version of trimmomatic, you can use the -basein and -baseout options to specify the prefix of the input and output file names.
Forum: Bioinformatics 04-29-2017, 01:52 PM
Replies: 189
Views: 58,287
Posted By mastal
One point to note might be that the thresholds...

One point to note might be that the thresholds you are using in the ILLUMINACLIP command are 30 for the palindrome mode and 10 for the simple mode, and you are not allowing for any mismatches. So...
Forum: Bioinformatics 04-29-2017, 04:26 AM
Replies: 2
Views: 679
Posted By mastal
Did you check whether any part of the amino acid...

Did you check whether any part of the amino acid sequence corresponds to your nucleotide sequence?

Perhaps blastx translated the entire nucleotide sequence, 1-548, and not just the part annotated...
Forum: Bioinformatics 04-22-2017, 12:42 PM
Replies: 10
Views: 3,980
Posted By mastal
The error message suggests that there is still...

The error message suggests that there is still something not quite right with the syntax.

It looks like you have not added a name for the 2U output file.

I don't think there is anything wrong...
Showing results 1 to 25 of 500

 


All times are GMT -8. The time now is 08:09 AM.


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