SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Searching for and extracting sequences within a fastQ file rtleenay Bioinformatics 3 10-27-2015 06:45 AM
Searching for a sequence in WGS / RNA-seq data ErikFas Bioinformatics 17 10-01-2015 03:53 AM
NCBI WGS submission: Need to trim sequences of various length from scaffolds lac302 Bioinformatics 2 04-10-2014 06:43 AM
Question: Searching FASTA file for specific IDs aw90 Bioinformatics 1 07-19-2013 03:14 AM
Patchwork - allele-specific somatic CNAs in cancer (WGS) mackan Bioinformatics 2 04-19-2012 04:54 AM

Reply
 
Thread Tools
Old 10-21-2016, 04:21 PM   #1
Kogan
Junior Member
 
Location: Moscow

Join Date: Oct 2016
Posts: 1
Default Searching for specific sequences in WGS data

Hi everyone.

I have never worked with WGS data and have a question. I have several thousands of sequences that I can convert in any format (e.g. fasta file) and I want to know whether each of the sequences in presented in my genome. Just 'yes' or 'no' for each sequence and, if it is possible, number of repeats for each sequence if it is presented in genome (e.g. presented twice or three times). I have both fastq and bam files. It is human genome.

I went through several presentations, but would be grateful if someone who has an experience in WGS data could suggest me the procedure/pipeline.

Thanks!
Kogan is offline   Reply With Quote
Old 10-24-2016, 08:33 PM   #2
finswimmer
Member
 
Location: Europe

Join Date: Oct 2016
Posts: 60
Default

What do you mean by "I have both fastq and bam files"? Is this realated to your "several thousands of sequences" or your genome?

If it's the first than you already have a file were all your sequenced were tried to mappe to a reference (was this your genome?). You can use than use samtools to asked how many reads are mapped.

Code:
samtools view -F 0x4 -c mybamfile.bam
or

Code:
samtools flagstat mybyfile.bam
fin swimmer
finswimmer is offline   Reply With Quote
Reply

Tags
wgs, whole genome sequencing

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 12:41 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO