SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Error in IlluminaBasecallsToSam.jar kga1978 Bioinformatics 0 01-30-2013 12:21 PM
picard error: Mismatch between read length and quals length writing read shawpa Bioinformatics 0 08-20-2012 05:52 AM
picard add read groups HGENETIC Bioinformatics 0 01-25-2012 03:19 AM
PicardTools IlluminaBaseCallsToSam Protaeus Bioinformatics 0 05-10-2011 10:06 AM
Read trimming and Picard JohnK Bioinformatics 4 11-02-2010 09:20 PM

Reply
 
Thread Tools
Old 02-26-2016, 10:05 AM   #1
user 31888
Member
 
Location: earth

Join Date: Nov 2015
Posts: 19
Default Picard IlluminaBasecallsToSam: Read Structure

Hi,

I am trying to convert MiSeq .bcl files into uSAM with Picard tools 'CheckIlluminaDirectory', 'ExtractIlluminaBarcodes' and 'IlluminaBasecallsToSam'.

These commands require to provide the READ_STRUCTURE option (described here)

Do you have any idea what the READ_STRUCTURE would be, knowing that I made a library with a TruSeq Exome kit (6 plex) and sequenced it on a MiSeq (2 x 75 bp - 76 cycles)?
Can we get this info from some of the sequencer run log files maybe?
user 31888 is offline   Reply With Quote
Old 02-26-2016, 10:11 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,059
Default

This information is in the RunInfo.xml file included in the flowcell_ID directory.

Looking at the help link for Picard above the code would be: 75T1S6B75T1S (Illumina's software generally drops the last base from READ1/2). If you need all 76 basepairs then 76T6B76T.
GenoMax is offline   Reply With Quote
Old 02-26-2016, 10:51 AM   #3
user 31888
Member
 
Location: earth

Join Date: Nov 2015
Posts: 19
Default

Thanks GenoMax !

Quote:
75T1S6B75T1S (Illumina's software generally drops the last base from READ1/2)
If Illumina's software drops the last base, it should not be sequenced, right?
So, why requesting Picard to skip 1 base then?
I read several times that people apply 1 extra cycle during the run.
Since I am interested to get 75bp long reads, does setting up MiSeq to 76 cycles really necessary?
user 31888 is offline   Reply With Quote
Old 02-26-2016, 11:00 AM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,059
Default

The last base (depending on the sequencing length) is always sequenced. As that last base does not have supporting phasing information it has been generally dropped by illumina's processing pipelines (CASAVA, bcl2fastq). There are options available to have the pipelines report that base, if you want all n+1 bp.

Some facilities supply n-1 bases, some will automatically setup n+1 runs (like in your case) or some may report all n bases.

Setting up an n+1 bp run ensures that resulting "n" bases have valid phasing information.
GenoMax is offline   Reply With Quote
Reply

Tags
illumina reads, picard

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 05:22 AM.


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