Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Reproducing "Tutorial for setting up SnapperDB Instance" : Not getting the sames SNP Address #22

Closed
ScaonE opened this issue Nov 13, 2018 · 1 comment

Comments

@ScaonE
Copy link

ScaonE commented Nov 13, 2018

Tutorial cmd & output :
run_snapperdb.py get_strains -c ebg4_config.txt;

AM933172 1.1.1.1.1.1.1
SRR5194193_1 1.1.1.1.1.1.9
SRR5055288_1 1.1.2.2.2.2.2
SRR5815674_1 1.1.3.3.3.3.3
SRR5864444_1 1.1.3.5.5.5.5
SRR5850014_1 1.1.3.5.5.8.8
SRR6131972_1 1.1.4.4.4.4.4
ERR2200244_1 1.1.7.7.7.7.7
SRR5583186_1 2.2.6.6.6.6.6

My cmd & output :
$snapperdb get_strains -c custom_salmo.txt;

AM933172 1.1.1.1.1.1.1
SRR5194193_1 1.1.1.1.1.1.9
ERR2200244_1 1.1.2.2.2.2.2
SRR5055288_1 1.1.3.3.3.3.3
SRR5815674_1 1.1.4.4.4.4.4
SRR5864444_1 1.1.4.6.6.6.6
SRR5850014_1 1.1.4.6.6.8.8
SRR6131972_1 1.1.5.5.5.5.5
SRR5583186_1 2.2.7.7.7.7.7

ERR2200244_1 seems to have an unexpected SNP address

May I ask which version of each external tool you recommand ? (biopython, psycopg2, paramiko, hashids, joblib, Postgres, PHEnix, Samtools, Picard, GATK, BWA)

Another question : Does fastq_to_db accept single-end FASTQ as input ?

@timdallman
Copy link
Contributor

  1. SNP Address:
    SNP addresses are order dependent so the cluster number assigned at each threshold is dependent on the order the isolates were added to SnapperDB.
  2. Fastq_to_db does not except single-end FASTQ - If you produce a VCF independetly from SnapperDB based in single end reads you can ingest them into SnapperDB with vcf_to_db

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants