Hey Gireesh, The sequences are not just reverse but reverse complement, AND they are written ONLY in the 5' to 3' direction. Taking the example of the first 16 bases ( in red) are the reverse complement of the last 16 ( in pink). 5'-3':TTCATGTTTTAACACT 5'-3':AGTGTTAAAACATGAA Now flip the lower sequence around ( write is 3'-5' direction) and the beauty of the DNA strand becomes obvious with A's pairing with T and G with C 5'-3':TTCATGTTTTAACACT 3'-5':AAGTACAAAATTGTGA
On Fri, Sep 17, 2010 at 7:08 PM, John McPherson <John.McPherson@oicr.on.ca> wrote: I think that you have extracted the “+” strand for your comparison. Below is the “-” strand from UCSC hg18
hg18_dna range=chr10:123229360-123229525 5'pad=0 3'pad=0 strand=+ repeatMasking=lower TTCATGTTTTAACACTGCCGTTTATGTGTGGATACTGAGGAAGGCATGGT TCGTAAGGCATGGGGTCTGGAGAAAAAACAGAATCATCTCCTGAAGAACA AGAACTTCTTGTGTCAGGGTAACTAGGTGAATACTGTTCGAGAGGTTGGC TGAGGTCCAAGTATTC
hg18_dna range=chr10:123229360-123229525 5'pad=0 3'pad=0 strand=- repeatMasking=lower GAATACTTGGACCTCAGCCAACCTCTCGAACAGTATTCACCTAGTTACCC TGACACAAGAAGTTCTTGTTCTTCAGGAGATGATTCTGTTTTTTCTCCAG ACCCCATGCCTTACGAACCATGCCTTCCTCAGTATCCACACATAAACGGC AGTGTTAAAACATGAA On Sep 17, 2010, at 5:52 PM, gireesh bogu wrote:
No I extrcated -. But your sense antisense sequences should be the same but in reverse right? But your two sequences are completely different ? Why is that?
On Fri, Sep 17, 2010 at 7:08 PM, John McPherson <John.McPherson@oicr.on.ca> wrote: I think that you have extracted the “+” strand for your comparison. Below is the “-” strand from UCSC hg18
hg18_dna range=chr10:123229360-123229525 5'pad=0 3'pad=0 strand=+ repeatMasking=lower TTCATGTTTTAACACTGCCGTTTATGTGTGGATACTGAGGAAGGCATGGT TCGTAAGGCATGGGGTCTGGAGAAAAAACAGAATCATCTCCTGAAGAACA AGAACTTCTTGTGTCAGGGTAACTAGGTGAATACTGTTCGAGAGGTTGGC TGAGGTCCAAGTATTC
hg18_dna range=chr10:123229360-123229525 5'pad=0 3'pad=0 strand=- repeatMasking=lower GAATACTTGGACCTCAGCCAACCTCTCGAACAGTATTCACCTAGTTACCC TGACACAAGAAGTTCTTGTTCTTCAGGAGATGATTCTGTTTTTTCTCCAG ACCCCATGCCTTACGAACCATGCCTTCCTCAGTATCCACACATAAACGGC AGTGTTAAAACATGAA
On 10-09-16 11:28 PM, "gireesh bogu" <gireeshkbogu@gmail.com> wrote:
Dear Galaxy
I found a bug in extracting FATSA sequence tool in galaxy.
I extracted sequence of this location chr10 123229360 123229525 fgfr2 0 -
hg18_chr10_123229360_123229525_- GAATACTTGGACCTCAGCCAACCTCTCGAACAGTATTCACCTAGTTACCC TGACACAAGAAGTTCTTGTTCTTCAGGAGATGATTCTGTTTTTTCTCCAG ACCCCATGCCTTACGAACCATGCCTTCCTCAGTATCCACACATAAACGGC AGTGTTAAAACATGA
And I checked back in UCSC browser by pasting the same location. The sequences are completely different (hg18).
Could you please help me in this issue
Thanx Gireesh
_______________________________________________ galaxy-user mailing list galaxy-user@lists.bx.psu.edu http://lists.bx.psu.edu/listinfo/galaxy-user