Commit d85901f6 authored by Tizian Schulz's avatar Tizian Schulz
Browse files

Changed source paths in test files and added graph building to workflow

parent 8cf71cc3
......@@ -10,3 +10,4 @@ bugRep*
*.txt
long*
small*
test/makefile
......@@ -7,7 +7,18 @@ def getGffFiles(sp):
return ["%s/annotations/%s.gff" %(sp, basename(g).split('.')[0]) for g in glob("%s/genomeSequences/*.fasta" %sp)]
def getAssemblies(wcs):
return ["%s/genomeSequences/%s" %(wcs.spec, basename(g)) for g in glob("%s/genomeSequences/*.fasta" %wcs.spec)]
gList = ["%s/genomeSequences/%s" %(wcs.spec, basename(g)) for g in glob("%s/genomeSequences/*.fas*" %wcs.spec)]
if wcs.gtype.find("wotSf") > -1:
gList.remove("arabidopsis/genomeSequences/sf_2.v7.PR_in_lowercase.fas")
return gList
def listGenomes(sp):
return ["%s/genomeSequences/%s.fasta" %(sp, a.strip()) for a in open("%s/genomeSequences/accessions.txt" %sp, 'r')]
def getReadSets(acc):
return ["arabidopsis/reads/" + basename(f) for l in open("arabidopsis/reads/filereport_read_run_PRJEB2457_tsv.txt", 'r') if l.find(acc.capitalize()) > -1 for f in l.split('\t')[6].split(';')]
rule all:
input:
......@@ -22,8 +33,8 @@ rule all:
180, 200, 220, 240]),
expand("arabidopsis/graphs/assemblyGraph_k21_Core_m18_d60.{gFSuf}", gFSuf=config['gFileSuffixes']),
"arabidopsis/benchmarks/benchmark_assemblyGraph_k21_Core_m18_d60.txt",
expand("arabidopsis/graphs/assemblyGraph_k21_wotSf_Core_m17_d60.{gFSuf}", gFSuf=config['gFileSuffixes']),
"arabidopsis/benchmarks/benchmark_assemblyGraph_k21_wotSf_Core_m17_d60.txt",
expand("arabidopsis/graphs/assemblyGraph_wotSf_k21_Core_m17_d60.{gFSuf}", gFSuf=config['gFileSuffixes']),
"arabidopsis/benchmarks/benchmark_assemblyGraph_wotSf_k21_Core_m17_d60.txt",
expand("arabidopsis/graphs/readGraph_k21_Core_m17_d60.{gFSuf}", gFSuf=config['gFileSuffixes']),
"arabidopsis/benchmarks/benchmark_readGraph_k21_Core_m17_d60.txt",
expand("{s}/panarooRes/gene_presence_absence.csv", s=config['species']),
......@@ -46,6 +57,84 @@ rule calculateCore:
" /usr/bin/time -v %s/Corer -q {params.qrm} -d {params.dlt} -i $(echo {input} | cut -d'.' -f1) " %config['corer_bin_dir'] +
"-o $(echo {output.grph} | cut -d'.' -f1) 2> {output.bnchmrk}"
rule buildPangenomeGraph:
input:
"{sp}/{gtype}List.txt"
threads:
8
output:
expand("{sp}/{g}_k{k}_s11.{gFSuf}", sp='{sp}', g='{gtype}', k='{k}', gFSuf=G_SUFS)
shell:
"mkdir -p {wildcards.sp}/graphs; Bifrost build -r {input.lst} -o $(echo {output} | cut -d'.' -f1) -t {threads} -v -c -k {wildcards.k}"
rule createReadGraphList:
input:
["arabidopsis/reads/%s.fasta" %basename(l).split('_')[0] for l in open("arabidopsis/genomeSequences/links.txt", 'r') if l.find("sf") < 0]
output:
"arabidopsis/readGraphList.txt"
shell:
"ls -l {input} | tr -s ' ' | cut -d' ' -f9 > {output}"
rule extractGraphSeqs:
input:
"arabidopsis/graphs/readGraph_{acc}_k21.gfa"
output:
"arabidopsis/reads/{acc}.fasta"
shell:
"python3 scripts/extractGFAseqs.py {input} > {output}"
rule buildGenomeGraph:
input:
"arabidopsis/{acc}List.txt"
threads:
8
output:
"arabidopsis/graphs/readGraph_{acc}_k21.gfa"
shell:
"mkdir -p arabidopsis/graphs; Bifrost build -s {input} -o $(echo {output} | cut -d'.' -f1) -t {threads} -v -k 21"
rule createReadList:
input:
getReadSets
output:
"arabidopsis/{acc}List.txt"
shell:
"ls -l {input} | tr -s ' ' | cut -d' ' -f9 > {output}"
rule downloadReadSets:
input:
"arabidopsis/reads/filereport_read_run_PRJEB2457_tsv.txt"
output:
["arabidopsis/reads/" + basename(f) for l in open("arabidopsis/reads/filereport_read_run_PRJEB2457_tsv.txt", 'r') for f in l.split('\t')[6].split(';')]
shell:
"wget {input}; mv *.fastq.gz arabidopsis/reads"
rule createAssemblyList:
input:
getAssemblies
output:
temp("{sp}/{gtype}List.txt")
wildcard_constraints:
gtype="[a]+"
shell:
"ls -l {input} | tr -s ' ' | cut -d' ' -f9 > {output}"
rule getNCBI_Genomes:
input:
"{sp}/genomeSequences/accessions.txt"
output:
listNCBI_Genomes
shell:
"python3 script/get_fastas_for_acc_list.py {input}"
rule: getArabidopsisGenomes:
input:
"arabidopsis/genomeSequences/links.txt"
output:
expand("arabidopsis/genomeSequences/{g}", g=[u.split('/')[-1].strip() for u in open("arabidopsis/genomeSequences/links.txt", 'r')])
shell:
"for u in $(cat {input}); do wget $u; done; mv *.fas arabidopsis/genomeSequences"
rule runPanaroo:
input:
getGffFiles
......
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/bur_0.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/can_0.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/ct_1.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/edi_0.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/hi_0.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/kn_0.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/ler_0.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/mt_0.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/no_0.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/oy_0.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/po_0.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/rsch_4.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/sf_2.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/tsu_0.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/wil_2.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/ws_0.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/wu_0.v7.PR_in_lowercase.fas
http://mtweb.cs.ucl.ac.uk/mus/www/19genomes/fasta/MASKED/zu_0.v7.PR_in_lowercase.fas
study_accession sample_accession experiment_accession run_accession tax_id scientific_name fastq_ftp submitted_ftp sra_ftp
PRJEB2457 SAMEA1027642 ERX011733 ERR031531 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031531/ERR031531_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031531/ERR031531_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031531/Bur_0_bur_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031531
PRJEB2457 SAMEA1027636 ERX011735 ERR031532 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031532/ERR031532_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031532/ERR031532_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031532/Can_0_can_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031532
PRJEB2457 SAMEA1027636 ERX011734 ERR031533 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031533/ERR031533.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031533/ERR031533_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031533/ERR031533_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031533/Can_0_can_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031533
PRJEB2457 SAMEA1027631 ERX011737 ERR031534 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031534/ERR031534_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031534/ERR031534_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031534/Ct_1_ct_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031534
PRJEB2457 SAMEA1027631 ERX011736 ERR031535 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031535/ERR031535.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031535/ERR031535_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031535/ERR031535_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031535/Ct_1_ct_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031535
PRJEB2457 SAMEA1027637 ERX011739 ERR031536 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031536/ERR031536_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031536/ERR031536_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031536/Edi_0_edi_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031536
PRJEB2457 SAMEA1027637 ERX011738 ERR031537 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031537/ERR031537.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031537/ERR031537_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031537/ERR031537_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031537/Edi_0_edi_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031537
PRJEB2457 SAMEA1027635 ERX011740 ERR031538 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031538/ERR031538_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031538/ERR031538_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031538/Hi_0_hi_amplified.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031538
PRJEB2457 SAMEA1027635 ERX011741 ERR031539 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031539/ERR031539_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031539/ERR031539_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031539/Hi_0_hi_nonamp.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031539
PRJEB2457 SAMEA1027635 ERX011742 ERR031540 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031540/ERR031540.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031540/ERR031540_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031540/ERR031540_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031540/Hi_0_hi_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031540
PRJEB2457 SAMEA1027629 ERX011744 ERR031541 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031541/ERR031541.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031541/ERR031541_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031541/ERR031541_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031541/Kn_0_kn_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031541
PRJEB2457 SAMEA1027629 ERX011743 ERR031542 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031542/ERR031542_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031542/ERR031542_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031542/Kn_0_kn_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031542
PRJEB2457 SAMEA1027626 ERX011746 ERR031543 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031543/ERR031543_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031543/ERR031543_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031543/Ler_0_ler_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031543
PRJEB2457 SAMEA1027626 ERX011745 ERR031544 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031544/ERR031544_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031544/ERR031544_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031544/Ler_0_ler_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031544
PRJEB2457 SAMEA1027638 ERX011748 ERR031545 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031545/ERR031545_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031545/ERR031545_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031545/Mt_0_mt_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031545
PRJEB2457 SAMEA1027638 ERX011747 ERR031546 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031546/ERR031546_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031546/ERR031546_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031546/Mt_0_mt_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031546
PRJEB2457 SAMEA1027641 ERX011750 ERR031547 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031547/ERR031547.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031547/ERR031547_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031547/ERR031547_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031547/No_0_no_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031547
PRJEB2457 SAMEA1027641 ERX011749 ERR031548 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031548/ERR031548_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031548/ERR031548_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031548/No_0_no_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031548
PRJEB2457 SAMEA1027627 ERX011752 ERR031549 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031549/ERR031549_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031549/ERR031549_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031549/Oy_0_oy_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031549
PRJEB2457 SAMEA1027627 ERX011751 ERR031550 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031550/ERR031550.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031550/ERR031550_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031550/ERR031550_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031550/Oy_0_oy_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031550
PRJEB2457 SAMEA1027640 ERX011754 ERR031551 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031551/ERR031551_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031551/ERR031551_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031551/Po_0_po_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031551
PRJEB2457 SAMEA1027640 ERX011753 ERR031552 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031552/ERR031552.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031552/ERR031552_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031552/ERR031552_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031552/Po_0_po_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031552
PRJEB2457 SAMEA1027634 ERX011756 ERR031553 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031553/ERR031553_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031553/ERR031553_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031553/Rsch_4_rsch_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031553
PRJEB2457 SAMEA1027634 ERX011755 ERR031554 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031554/ERR031554.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031554/ERR031554_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031554/ERR031554_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031554/Rsch_4_rsch_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031554
PRJEB2457 SAMEA1027630 ERX011758 ERR031555 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031555/ERR031555_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031555/ERR031555_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031555/Tsu_0_tsu_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031555
PRJEB2457 SAMEA1027630 ERX011757 ERR031556 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031556/ERR031556.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031556/ERR031556_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031556/ERR031556_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031556/Tsu_0_tsu_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031556
PRJEB2457 SAMEA1027628 ERX011760 ERR031557 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031557/ERR031557_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031557/ERR031557_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031557/Wil_2_wil2_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031557
PRJEB2457 SAMEA1027628 ERX011759 ERR031558 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031558/ERR031558.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031558/ERR031558_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031558/ERR031558_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031558/Wil_2_wil2_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031558
PRJEB2457 SAMEA1027632 ERX011762 ERR031559 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031559/ERR031559_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031559/ERR031559_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031559/Ws_0_ws_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031559
PRJEB2457 SAMEA1027632 ERX011761 ERR031560 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031560/ERR031560_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031560/ERR031560_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031560/Ws_0_ws_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031560
PRJEB2457 SAMEA1027632 ERX011763 ERR031561 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031561/ERR031561_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031561/ERR031561_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031561/Ws_0_ws_tsl.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031561
PRJEB2457 SAMEA1027633 ERX011765 ERR031562 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031562/ERR031562.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031562/ERR031562_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031562/ERR031562_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031562/Wu_0_wu_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031562
PRJEB2457 SAMEA1027633 ERX011764 ERR031563 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031563/ERR031563_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031563/ERR031563_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031563/Wu_0_wu_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031563
PRJEB2457 SAMEA1027639 ERX011767 ERR031564 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031564/ERR031564_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031564/ERR031564_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031564/Zu_0_zu_phaseI.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031564
PRJEB2457 SAMEA1027639 ERX011766 ERR031565 3702 Arabidopsis thaliana ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031565/ERR031565_1.fastq.gz;ftp.sra.ebi.ac.uk/vol1/fastq/ERR031/ERR031565/ERR031565_2.fastq.gz ftp.sra.ebi.ac.uk/vol1/run/ERR031/ERR031565/Zu_0_zu_PII.bam ftp.sra.ebi.ac.uk/vol1/err/ERR031/ERR031565
#!/usr/bin/env python3
from sys import argv
counter = 1
K = 21
for l in open(argv[1], 'r'):
if l.startswith('S'):
print('>%d' %counter)
print(l.split('\t')[2])
counter += 1
\ No newline at end of file
#include <gtest/gtest.h>
#include "BridgingTest.h"
#include "../Bridging.cpp"
#include "../src/Bridging.cpp"
//Tests for function const bool lCrTooFar(const size_t&, const list<pair<uint32_t, uint32_t>>&, const uint32_t&)//
// 1. The core list is empty and the unitig is (not) too long DONE
......
#ifndef CORE_TEST_HPP
#define CORE_TEST_HPP
#include "../Core.h"
#include "../src/Core.h"
#define DEFAULT_TEST_K 9
#define DEFAULT_TEST_G 4
......
......@@ -2,7 +2,7 @@
#include <bifrost/CompactedDBG.hpp>
#include <bifrost/ColoredCDBG.hpp>
#include "../Core.cpp"
#include "../src/Core.cpp"
#include "CoreTest.h"
//Tests for function const bool chkQrm(UnitigColorMap<CoreInfo>&, const uint32_t&)//
......
#include <gtest/gtest.h>
#include "../IO.cpp"
#include "../src/IO.cpp"
#include "IOtest.h"
#include "../Bridging.h"
#include "../src/Bridging.h"
//Tests for function const bool prsArgs(int&, char**, string&, uint32_t&, uint32_t&)//
// 1. Input graph prefix is (not) given DONE
......
#ifndef IO_TEST_HPP
#define IO_TEST_HPP
#include "../Traversal.h"
#include "../src/Traversal.h"
#include "CoreTest.h"
class PrsArgsTest : public ::testing::Test {
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment