Question: New To DiscoSnp-RAD
gravatar for Giovanni.madrigal12
3 months ago by
Giovanni.madrigal1210 wrote:

Dear all,

I am relatively new to bioinformatics and coding. Looking the DiscoSnp's supplementary data, the commands seem pretty straight forward; however, I just can't seem to wrap my head around the command where the fof.txt was created. For the experiment in the DIscoSnp++ paper, the command was

for((j=0;j<i;j++)); do echo coli_muted_n_30_genome_${j}_reads.fasta > fof.txt; done.

I am working with RAD-seq data and would like some advice on this. I appreciate any help and advice!

Sincerely, Giovanni Madrigal

snp discosnp • 228 views
ADD COMMENTlink modified 3 months ago by pierre.peterlongo820 • written 3 months ago by Giovanni.madrigal1210
gravatar for pierre.peterlongo
3 months ago by
pierre.peterlongo820 wrote:

Hi all,

This code comes from the supp mat of the original paper where we explained how experiments were conducted.

discoSnp (discoSnpRad or discoSnp++) takes as input read files organized in file of files. The idea (described in the documentation) is that a fof contains a list of read files.

The command above is a simple bash for loop, that creates automatically a fof file containing indeed all .fasta of n coli individuals.

More generally, the fof files enable to deal with pair end files, to virtually concatenate read sets, ... Have a look to the documentation.

Have a nice day, Pierre

ADD COMMENTlink written 3 months ago by pierre.peterlongo820
gravatar for h.mon
3 months ago by
h.mon15k wrote:

This is a bash for loop, contaminated with what I think is html code. This code just creates a file called fof.txt, with:


There will be $i names, last name is coli_muted_n_30_genome_(i-1)_reads.fasta.

Where did you get this code?

ADD COMMENTlink written 3 months ago by h.mon15k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 848 users visited in the last hour