Editing
Daniel:Notebook/ComboLock/2016-9-27
(section)
Jump to navigation
Jump to search
Warning:
You are not logged in. Your IP address will be publicly visible if you make any edits. If you
log in
or
create an account
, your edits will be attributed to your username, along with other benefits.
Anti-spam check. Do
not
fill this in!
=Positive Control X3 (Started [[Daniel:Notebook/ComboLock/2016-8-29|08/29/2016]])= [[Daniel:Notebook/ComboLock|Back to Calendar]] Sequencing data in Genome miner at /media/NAS3_volume2/Dan/combolock/miseq_20160912 ==Data Analysis== As a reminder, here is a table of the three samples, as well as the size select gels for the library prep. {| class="wikitable" <hiddentext>generated with [[:de:Wikipedia:Helferlein/VBA-Macro for EXCEL tableconversion]] V1.8</hiddentext> |- style="background-color:#D8E4BC;font-size:12pt;font-weight:bold" align="center" | width="72" height="25" | Sample | width="65" | Index | width="102" | Approx bp Size |- style="font-size:12pt" | height="15" valign="bottom" | Sample0 | align="center" align="center" valign="bottom" | 34 | align="center" align="center" valign="bottom" | 130 |- style="background-color:#D9D9D9;font-size:12pt" | height="15" valign="bottom" | Sample1 | align="center" align="center" valign="bottom" | 35 | align="center" valign="bottom" | >500 |- style="font-size:12pt" | height="15" valign="bottom" | Sample2 | align="center" align="center" valign="bottom" | 36 | align="center" align="center" valign="bottom" | 100 |} <gallery perrow=2 heights=300px widths=300px> File:2016-09-05-X3Dilution-SizeSelect.png|Before Image File:2016-09-05-X3Dilution-SizeSelect-After.png|After Image </gallery> ===Methods-Manual Inspection=== I used several methods to analyze this data set. For the simplest, I just took the raw sequences and put them into excel charts. Example for one read below: {| class="wikitable" <hiddentext>generated with [[:de:Wikipedia:Helferlein/VBA-Macro for EXCEL tableconversion]] V1.8</hiddentext> |- style="background-color:#B7DEE8;font-size:12pt;font-weight:bold" align="center" | width="226" height="28" | Name | width="65" | Read 1/2 | width="283" | Sequence | width="65" | Offset | width="114" | Seq 1 | width="65" | Seq1 Match | width="82" | Index34? (R1) | width="108" | Seq2 | width="146" | Seq 2 Match | width="165" | Seq3 |- style="font-size:12pt" valign="bottom" | height="30" | @M00159:51:000000000-ATU5U:1:1101:15771:2370 1:N:0:AACTGA | R1 | GTGTTTGCGGCTAAGCGTTCCTAGCNANNGAATCTCGTATGCCGTCTTCTGCNNGCACAAAAAAACTCTTTTTNC | GT | GTTTGCGGCTAAGCGTTCCTAGC | AmpR6.3-RevComp | NANNGA | ATCTCGTATGCCGTCTTCTGCNNGC | AmpR6.3-RevComp bp: 33-53 | ACAAAAAAACTCTTTTTNC |- style="background-color:#D9D9D9;font-size:12pt" valign="bottom" | height="30" | @M00159:51:000000000-ATU5U:1:1101:15771:2370 2:N:0:AACTGA | R2 | ACTTCGGACCTCGATAACATTTGAGAGTGGTGTAGATCTTGATGTTCGCCGAATCNTTAAAAANAAATGTNTTCNT | AC | TTCGGACCTCGATAACATTTGAG | AmpF6.4-RevComp | AGTGGT | GTAGATCTTGATGTTCGCCGAATCN | First 11 match AmpF6.4Rev Comp 30-40 | TTAAAAANAAATGTNTTCN |} These allow me to manually assess the sequences, but are also not great for larger scale processing (I only did ~4 for any given index). For larger scale analyses, I used blast. ===Methods-Blast=== To do blast, I used a database made of the various primer, adapter, and index sequences. I set a threshold of 80% alignment length match to get rid of small alignments (that is, the read must contain >=80% of the bases matching to the primer/adapter/index). The table below summarizes the sequences targeted and how many bp a match requires. [[Media:Analysis.miseq20160909.txt|Shell script for blast analysis]] {| class="wikitable" <hiddentext>generated with [[:de:Wikipedia:Helferlein/VBA-Macro for EXCEL tableconversion]] V1.8</hiddentext> |- style="background-color:#B7DEE8;font-size:12pt;font-weight:bold" align="center" | width="86" height="45" | Name | width="65" | Length | width="101" | Number of bp Correct for 80% match |- style="font-size:12pt" | height="15" valign="bottom" | primer2 | align="center" align="center" valign="bottom" | 18 | align="center" align="center" valign="bottom" | 15 |- style="background-color:#D9D9D9;font-size:12pt" | height="15" valign="bottom" | primer4 | align="center" align="center" valign="bottom" | 18 | align="center" align="center" valign="bottom" | 15 |- style="font-size:12pt" | height="15" valign="bottom" | primer6 | align="center" align="center" valign="bottom" | 18 | align="center" align="center" valign="bottom" | 15 |- style="background-color:#D9D9D9;font-size:12pt" | height="15" valign="bottom" | primer12 | align="center" align="center" valign="bottom" | 18 | align="center" align="center" valign="bottom" | 15 |- style="font-size:12pt" | height="15" valign="bottom" | IlluminaP5 | align="center" align="center" valign="bottom" | 20 | align="center" align="center" valign="bottom" | 16 |- style="background-color:#D9D9D9;font-size:12pt" | height="15" valign="bottom" | IlluminaP7 | align="center" align="center" valign="bottom" | 24 | align="center" align="center" valign="bottom" | 20 |- style="font-size:12pt" | height="15" valign="bottom" | LinkerL | align="center" align="center" valign="bottom" | 23 | align="center" align="center" valign="bottom" | 19 |- style="background-color:#D9D9D9;font-size:12pt" | height="15" valign="bottom" | LinkerR | align="center" align="center" valign="bottom" | 25 | align="center" align="center" valign="bottom" | 20 |- style="font-size:12pt" | height="15" valign="bottom" | Latch1Barcode | align="center" align="center" valign="bottom" | 8 | align="center" align="center" valign="bottom" | 7 |- style="background-color:#D9D9D9;font-size:12pt" | height="15" valign="bottom" | Latch2Barcode | align="center" align="center" valign="bottom" | 8 | align="center" align="center" valign="bottom" | 7 |- style="font-size:12pt" | height="15" valign="bottom" | Latch3Barcode | align="center" align="center" valign="bottom" | 8 | align="center" align="center" valign="bottom" | 7 |- style="background-color:#D9D9D9;font-size:12pt" | height="15" valign="bottom" | Latch6Barcode | align="center" align="center" valign="bottom" | 8 | align="center" align="center" valign="bottom" | 7 |- style="font-size:12pt" | height="15" valign="bottom" | Latch8Barcode | align="center" align="center" valign="bottom" | 8 | align="center" align="center" valign="bottom" | 7 |- style="background-color:#D9D9D9;font-size:12pt" | height="15" valign="bottom" | Latch9Barcode | align="center" align="center" valign="bottom" | 8 | align="center" align="center" valign="bottom" | 7 |} ===Interior Blast=== I also performed a blast on the templates in combolock (padlocks, positive control templates, and full sequencing primers) against the individual sequences, using the templates as the db. In the table below are the matches that were not expected (off-target matches). For this one I did not do any thresholding since I want to see all possible off-target bindings, although blast seems to have a threshold of at least 7 bases. Highlighted in the table below are the most concerning sequence matches (it's all the same similarity). That is the LinkerR region matches to somewhere on Primer2. It is a 9bp match, and while it is not very likely, is certainly possible to cause problems. blastn -db combolockseq.fa -query primerseqs.fa -task blastn-short -outfmt 7 {| class="wikitable" <hiddentext>generated with [[:de:Wikipedia:Helferlein/VBA-Macro for EXCEL tableconversion]] V1.8</hiddentext> |- style="background-color:#B7DEE8;font-size:12pt;font-weight:bold" align="center" | width="127" height="30" | query id | width="98" | subject id | width="65" | % identity | width="64" | alignment length | width="70" | mismatches | width="65" | gap opens | width="65" | q. start | width="65" | q. end | width="65" | s. start | width="65" | s. end | width="65" | evalue | width="65" | bit score | width="89" | Notes |- style="background-color:#FFFF00;font-size:12pt" valign="bottom" | height="15" | primer2 | Padlock0001-68 | align="center" | 100 | align="center" | 9 | align="center" | 0 | align="center" | 0 | align="center" | 5 | align="center" | 13 | align="center" | 84 | align="center" | 76 | align="center" | 0.014 | align="center" | 18.3 | LinkerR* |- style="background-color:#FFFF00;font-size:12pt" valign="bottom" | height="15" | primer2 | AmpF6.4 | align="center" | 100 | align="center" | 9 | align="center" | 0 | align="center" | 0 | align="center" | 5 | align="center" | 13 | align="center" | 54 | align="center" | 46 | align="center" | 0.014 | align="center" | 18.3 | LinkerR* |- style="font-size:12pt" valign="bottom" | height="15" | primer6 | Amp3 | align="center" | 100 | align="center" | 8 | align="center" | 0 | align="center" | 0 | align="center" | 4 | align="center" | 11 | align="center" | 44 | align="center" | 51 | align="center" | 0.057 | align="center" | 16.4 | L3BC* |- style="background-color:#D9D9D9;font-size:12pt" valign="bottom" | height="15" | IlluminaP5 | Amp3 | align="center" | 100 | align="center" | 7 | align="center" | 0 | align="center" | 0 | align="center" | 4 | align="center" | 10 | align="center" | 61 | align="center" | 55 | align="center" | 0.26 | align="center" | 14.4 | Primer6* |- style="font-size:12pt" valign="bottom" | height="15" | IlluminaP5 | Amp1 | align="center" | 100 | align="center" | 7 | align="center" | 0 | align="center" | 0 | align="center" | 4 | align="center" | 10 | align="center" | 61 | align="center" | 55 | align="center" | 0.26 | align="center" | 14.4 | Primer6* |- style="background-color:#FFFF00;font-size:12pt" valign="bottom" | height="15" | LinkerR | Amp3 | align="center" | 100 | align="center" | 9 | align="center" | 0 | align="center" | 0 | align="center" | 15 | align="center" | 23 | align="center" | 13 | align="center" | 5 | align="center" | 0.023 | align="center" | 18.3 | Primer2* |- style="background-color:#FFFF00;font-size:12pt" valign="bottom" | height="15" | LinkerR | Amp1 | align="center" | 100 | align="center" | 9 | align="center" | 0 | align="center" | 0 | align="center" | 15 | align="center" | 23 | align="center" | 13 | align="center" | 5 | align="center" | 0.023 | align="center" | 18.3 | Primer2* |- style="font-size:12pt" valign="bottom" | height="15" | Latch3Barcode | Amp3 | align="center" | 100 | align="center" | 7 | align="center" | 0 | align="center" | 0 | align="center" | 2 | align="center" | 8 | align="center" | 71 | align="center" | 65 | align="center" | 0.057 | align="center" | 14.4 | Primer6* |- style="background-color:#D9D9D9;font-size:12pt" valign="bottom" | height="15" | Latch3Barcode | Amp1 | align="center" | 100 | align="center" | 7 | align="center" | 0 | align="center" | 0 | align="center" | 2 | align="center" | 8 | align="center" | 71 | align="center" | 65 | align="center" | 0.057 | align="center" | 14.4 | Primer6* |} ==Discussion== So index34 (~130bp product) and index 36 (~100bp) product have similar content. For both, the R1 reads contain (~80% and ~70% respectively) the Illumina P7 and LinkerL regions, which are both part of the AmpR6.3 sequencing primer. The R2 reads contain (~50%) the Illumina P5 and LinkerR, which are the components of AmpF6.4Sol. This means that a good number of the reads are probably AmpF6.4 and AmpR6.3 dimers. I think the difference might be between primer dimers and off-target binding to primer 2. [[Image:Miseq20160912-BlastCountsHistogram.png|800px]] [[Category:ComboLock]] [[Category:20160829]]
Summary:
Please note that all contributions to ZhangLabWiki may be edited, altered, or removed by other contributors. If you do not want your writing to be edited mercilessly, then do not submit it here.
You are also promising us that you wrote this yourself, or copied it from a public domain or similar free resource (see
ZhangLabWiki:Copyrights
for details).
Do not submit copyrighted work without permission!
Cancel
Editing help
(opens in new window)
Navigation menu
Personal tools
Not logged in
Talk
Contributions
Create account
Log in
Namespaces
Page
Discussion
English
Views
Read
Edit
View history
More
Search
Navigation
Main Page
Current events
Recent changes
Random page
Investigators
Matt Cai
Song Chen
Eric Chu
Dinh Diep
Elizabeth Duong
Shicheng Guo
Alan Fung
Daniel Jacobsen
Blue Lake
Huy Lam
Alice Li
Andrew Richards
Brandon Sos
Chris Wei
Yan Wu
Kun Zhang
Tools
What links here
Related changes
Special pages
Page information