SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
after piccard MarkDuplicates drop in mapping percentage trotos Bioinformatics 1 05-19-2014 07:55 AM
Low mapping percentage of reads on assembled contigs morning latte Bioinformatics 20 03-23-2014 09:01 AM
Typical alignment mapping percentage with genome? obscurite RNA Sequencing 8 05-02-2013 09:39 AM
Low mapping percentage with TopHat2 DerSeb RNA Sequencing 3 06-05-2012 05:35 AM
Calculating Percentage of n's for Each Transcript in FASTA File KDS Bioinformatics 7 06-07-2011 11:06 PM

Reply
 
Thread Tools
Old 05-24-2018, 12:52 PM   #1
SuzuBell
Member
 
Location: U.S.A

Join Date: Nov 2013
Posts: 18
Default Determine mapping percentage from .fastq file

Is it possible to obtain the mapping percent from .fastq file?

If I do the following:

grep "@HISEQ" -c file.fastq

I believe I obtain the number of sequences. Can this be used to determine the mapping rate without needing a .bam or .sam file?
SuzuBell is offline   Reply With Quote
Old 05-24-2018, 06:00 PM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,747
Default

No. Unless the data is mapped to a reference you would not be able to determine the mapping percentage. The % will change depending on the reference you use and/or options you choose for the aligner.

If for some reason you do not wish to create an alignment file you could use bbmap.sh from BBMap suite to get alignment statistics by not providing an "out=" file that would normally store the alignments. You can also sample a set fraction of reads (using "reads=1000000") option to get an idea of what the alignment percentage may be like.
GenoMax is offline   Reply With Quote
Reply

Tags
fastq, mapping coverage

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 08:14 AM.


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