Answers.Body: <p>Try this pipeline: </p>
  
  <ol>
  <li>Clean raw fastq files (trimmomatic)</li>
  <li>Align clean fastq files to the reference genome (dna-to-dna aligner of your choice)</li>
  <li>Convert the alignment sam file to bam (samtools)</li>
  <li>Call SNPs from bam file (samtools)</li>
  </ol>
  
  <p>Here's a start tutorial on this: <a href="http://ged.msu.edu/angus/tutorials-2013/snp_tutorial.html" rel="nofollow noreferrer">http://ged.msu.edu/angus/tutorials-2013/snp_tutorial.html</a></p>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-08-19T13:57:56.257Z
Answers.Text: Thanks, that's a great tutorial. I'm able to call snps and generate vcf files, but I'm looking to quickly query the vcf files for a specific list of snps, preferably listed by reference genome coordinates. I'm not sure on the wording so I haven't been able to find a tool.
Answers.UserDisplayName: Dan K
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-08-19T06:34:42.700Z
Answers.LastActivityDate: 2017-08-19T06:34:42.700Z
Answers.UserDisplayName: aechchiki
Answers.Body: <p>You could use <a href="http://www.biorxiv.org/content/early/2017/04/07/118000" rel="nofollow noreferrer">ARIBA</a>, which will map you reads to the regions of your interest (more on that in a second) and then assemble/call variants.</p>
  
  <p><a href="https://i.sstatic.net/ZHRVk.jpg" rel="nofollow noreferrer"><img src="https://i.sstatic.net/ZHRVk.jpg" alt="enter image description here"></a></p>
  
  <p>The tool is mostly used to find known variants related to AMR (antimicrobial resistance) from a set of databases, but <a href="https://github.com/sanger-pathogens/ariba/wiki/Task%3A-prepareref#user-provided-data---fasta-and-metadata" rel="nofollow noreferrer">according to the project's wiki</a> you can build your own database with your region of interest.</p>
  
  <pre><code>Additionally, reference sequences can be either of:
  
      Presence/absence sequences. ARIBA will look for these sequences in the input reads and report any that it finds, and also any variants between these sequences and the input reads.
  
      Variants only sequences. These should have known variant details specified in the metadata file (see below). ARIBA reports only when it finds at least one of the given variants in each of these these sequences. If a sample has one of these sequences, but does not have one of the supplied variants, then it is not reported. If you supply a variants only sequence, but no variant, then the sequence will be removed during sanity checks.
  </code></pre>
  
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-08-22T08:22:53.450Z
Answers.Text: @DanK, please consider accepting one of the answers if they helped you solve your problem, thanks
Answers.UserDisplayName: mgalardini
Answers.ContentLicense: CC BY-SA 3.0
Answers.CreationDate: 2017-08-20T22:16:29.723Z
Answers.LastActivityDate: 2017-08-20T22:16:29.723Z
Answers.UserDisplayName: mgalardini
Body: <p>Is there a tool that can scan fastq files without assembling them for a custom list of user defined snps?</p>
  
Comments.ContentLicense: CC BY-SA 3.0
Comments.CreationDate: 2017-08-19T06:35:51.143Z
Comments.Text: Where do fastq come from? How does the list look like? Do you have a reference genome? Please clarify.
Comments.UserDisplayName: aechchiki
Comments.ContentLicense: CC BY-SA 3.0
Comments.CreationDate: 2017-08-19T13:45:37.980Z
Comments.Text: They're bacteria sequences from a MiSeq, paired end fastq. They are all the same species and I do have a good reference genome. I'd like to provide a list of genomic coordinates I. The reference genome and find of those snps are present in the new sequences. I'm looking for an easy way just to find the small list I'd like instead of all variants from the reference.
Comments.UserDisplayName: Dan K
Comments.ContentLicense: CC BY-SA 3.0
Comments.CreationDate: 2017-08-29T20:18:59.287Z
Comments.Text: could you edit the question and add a sample of your vcf file and the reference genome coordinates file? would be cool to play with, maybe some sed/awk one-liner will do the job
Comments.UserDisplayName: aechchiki
Tags: <fastq><vcf><snp>
Title: Is there a way to quickly verify the presence of some SNPs in Fastq files?
ThreadText: # Is there a way to quickly verify the presence of some SNPs in Fastq files?
  
  (Asked by: username_0 on 2017-08-19)
  
  Is there a tool that can scan fastq files without assembling them for a custom list of user defined snps?
  
  **Tags:** fastq, vcf, snp
  
  > **Comments:**
  > **username_1**: Where do fastq come from? How does the list look like? Do you have a reference genome? Please clarify.  
  > **username_0**: They're bacteria sequences from a MiSeq, paired end fastq. They are all the same species and I do have a good reference genome. I'd like to provide a list of genomic coordinates I. The reference genome and find of those snps are present in the new sequences. I'm looking for an easy way just to find the small list I'd like instead of all variants from the reference.  
  > **username_1**: could you edit the question and add a sample of your vcf file and the reference genome coordinates file? would be cool to play with, maybe some sed/awk one-liner will do the job  
  
  ---
  
  ## 2 answers
  
  ### Answer by username_1 on 2017-08-19. Score: 4
  
  Try this pipeline:
  
  1. Clean raw fastq files (trimmomatic)
  2. Align clean fastq files to the reference genome (dna-to-dna aligner of your choice)
  3. Convert the alignment sam file to bam (samtools)
  4. Call SNPs from bam file (samtools)
  
  Here's a start tutorial on this: <http://ged.msu.edu/angus/tutorials-2013/snp_tutorial.html>
  
  > **Comments:**
  > **username_0**: Thanks, that's a great tutorial. I'm able to call snps and generate vcf files, but I'm looking to quickly query the vcf files for a specific list of snps, preferably listed by reference genome coordinates. I'm not sure on the wording so I haven't been able to find a tool.  
  
  ---
  
  ### Answer by username_2 on 2017-08-20. Score: 1
  
  You could use [ARIBA](http://www.biorxiv.org/content/early/2017/04/07/118000), which will map you reads to the regions of your interest (more on that in a second) and then assemble/call variants.
  
  [![enter image description here](https://i.sstatic.net/ZHRVk.jpg)](https://i.sstatic.net/ZHRVk.jpg)
  
  The tool is mostly used to find known variants related to AMR (antimicrobial resistance) from a set of databases, but [according to the project's wiki](https://github.com/sanger-pathogens/ariba/wiki/Task%3A-prepareref#user-provided-data---fasta-and-metadata) you can build your own database with your region of interest.
  
  ```
  Additionally, reference sequences can be either of:
  
      Presence/absence sequences. ARIBA will look for these sequences in the input reads and report any that it finds, and also any variants between these sequences and the input reads.
  
      Variants only sequences. These should have known variant details specified in the metadata file (see below). ARIBA reports only when it finds at least one of the given variants in each of these these sequences. If a sample has one of these sequences, but does not have one of the supplied variants, then it is not reported. If you supply a variants only sequence, but no variant, then the sequence will be removed during sanity checks.
  
  ```
  
  > **Comments:**
  > **username_2**: @username_0, please consider accepting one of the answers if they helped you solve your problem, thanks  
  
  ---
  
  