Question: New To DiscoSnp-RAD
0
gravatar for Giovanni.madrigal12
9 months ago by
Giovanni.madrigal1220 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 • 329 views
ADD COMMENTlink modified 9 months ago by pierre.peterlongo830 • written 9 months ago by Giovanni.madrigal1220
2
gravatar for pierre.peterlongo
9 months ago by
France
pierre.peterlongo830 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 9 months ago by pierre.peterlongo830
1
gravatar for h.mon
9 months ago by
h.mon21k
Brazil
h.mon21k 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:

coli_muted_n_30_genome_0_reads.fasta
coli_muted_n_30_genome_1_reads.fasta
coli_muted_n_30_genome_2_reads.fasta
...

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 9 months ago by h.mon21k
Please log in to add an answer.

Help
Access

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