Question: New To DiscoSnp-RAD
gravatar for Giovanni.madrigal12
13 days 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 • 122 views
ADD COMMENTlink modified 13 days ago by pierre.peterlongo710 • written 13 days ago by Giovanni.madrigal1210
gravatar for pierre.peterlongo
13 days ago by
pierre.peterlongo710 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 13 days ago by pierre.peterlongo710
gravatar for h.mon
13 days ago by
h.mon12k 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 13 days ago by h.mon12k
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: 1601 users visited in the last hour