mothur

make.sra command only generate submission.xml and submit.ready files


#1

Hi,

I’ve been running problems using make.sra command. It used to work great.

Now I have issues like make.sra only generates submission.xml and submit.ready files, without sequences for each samples.

I used raw sequencing data, make.sra (file=xx.file, project=xx.project, mimark=xx.txt)

my xx.file contains information

name_1 DMP04993_L1_name_1_1.fq.gz DMP04993_L1_name_1_2.fq.gz
name_2 DMP04993_L1_name_2_1.fq.gz DMP04993_L1_name_2_2.fq.gz

If I run make.sra(fastq=xx.fastq, oligos=xx.oligos, project=xx.project, mimark=xx.txt)

the fastq file is made by merging all the raw data(forward and reverse fastq) files together.

All the seqs end up in the scarp file.

How can I correct it, please?

Thanks you so much for your help


#2

With the following command and the setup of the xx.file, mothur should not create new files for each sample because the samples are already separated by sample.

mothur > make.sra (file=xx.file, project=xx.project, mimark=xx.txt)

my xx.file contains information

name_1 DMP04993_L1_name_1_1.fq.gz DMP04993_L1_name_1_2.fq.gz - sample name_1 forward and reverse files for upload.
name_2 DMP04993_L1_name_2_1.fq.gz DMP04993_L1_name_2_2.fq.gz - sample name_2 forward and reverse files for upload.


#3

Makes sense! Thank you so much!


#4

Any luck with your submission last year…?

I believe I’m currently experiencing the same problem. I.e. most of my seqs end up in the scrap file after running the make.sra command, either straight from the fastq with oligos file or using “file=xx.file”.

Based on some of your other post (make.sra generated the length of the quality reads (in submitted file) is longer than the sequence reads) I suspect you were trying to submit some of your old “Mr. DNA” fastq files, which is the case for me at the moment.

I followed your advice from the previous post, generating a merge fastq from: make.contigs(ffastq=xx.R1.fastq, rfastq=xx.R2.fastq), trim.seqs(), and make.fastq(fasta=current, qfile=“current”.qual).

Then using the merged fastq and an oligos file to generate sample/group specific fastq files.

Anyway, any tips would be great!
Thanks.