sblock154733

chr3R:20518218-20518323 +

Candidate

[View on UCSC Genome Browser {Cornell Mirror}]

intron


Legend:

   mature strand at position 20518242-20518265
   star strand at position 20518278-20518305
   mutations
   read mismatches
   coordinate mismatch (possible incorrect ortholog)

12-way Fast Statistical Aligner re-alignment (From MAFs)

Species Coordinate Alignment
dm3 chr3R:20518203-20518338 + TACACCATGAACAACATAAGATT--------------TATTCGGTCT--ACAAACTGAAGCTTAATTGTGATTTATGGATGTGTATTTTTTAGCAGTTAAAAATCACAATTTAGCCACT-TTGTAG-GTCGATGA-------A--ACAGCT----------------CACAACTCTCAC
droSim1 chr3R_random:984803-984937 + TACACCATGAACAACATAAGATA--------------TATTCGGTCT--ACAAACTGAAGCTTAATTGTGATTTAAGGATGTGTATTTTTTAGCACTAAAAAATCACAATTTAGCCAAT-TTGTAG-GTCGAT-A-------A--AAACCA----------------CACAACTCTCAC
droSec1 super_0:20852363-20852498 + TACACCATGAACAACATAAGATA--------------TATTCGGTCT--ACAAACTGAAGCTTAATTGTGATTTAAGGATGTGTATTTTTTTGCACTAAAAAATCACAATTTAGCCAAT-TTGTAG-GTCGATGA-------A--AAACCA----------------CACAACTCTCAC
droYak2 chr3R:25978425-25978557 - TACACCA---ACAACATAAGATA--------------TATTCGGCAT--ACAAACTTGAGCTTAATTGTGATTTAAGGATGTGAATT-TCAAGCAACGAAAAATCACAATGAAGCCAATTTTGTAG-GCCGATGG-------A--ACACCT----------------CACAACTCTCAC
droEre2 scaffold_4820:7539900-7540034 - CACACCATGAACAACATAAGATA--------------TATTCGGCCT--ACAAACATAAGCATAATTGTGATTTAAGGATGTATATT-TCAAGTAGTAATACATCACAATTTAGCCAAT-TTGTAG-GCCGATGA-------A--AAACCT----------------CACAACTCTCAC
droAna3 scaffold_13340:4657405-4657499 - TACACCATGTCTGACTTGA--------------------------------------AAATTTAATTGTGATTTGGTTT------------------TTTGAATCACAATTAAGTTT---TTTTAA-GCCGGA-A-------AACACACAT----------------TTCAACTCTCAC
dp4 chr2:1714385-1714528 - TAGACCATGAACAACATAAGAAAAACTTACACACGTTTCTCTGGTTTATACAAAATCAAGTTTCATTGTGATTGGAGCCT---TATATTTTGGGATAATCCAATCACAATCAAGCTTC--CTGTACAACCCATGA-------G--AG---------------------ACAATTCACAC
droPer1 super_7:1868921-1869065 - TACACCATGAACAACATAAGAAAAACTTACACACGTTTCTCTGGTTTATACAAAATGAAGTTTCATTGTGATTGGAGCCTATAT--GTTCTGGGATAATCCAATCACAATCAAGCTTC--CTGTACAACCCATGA-------G--AG---------------------ACAATTCACAC
droWil1 scaffold_181089:9576497-9576520 - caaagac---------------a---------------------------------------------------------------------------------------------------------------------------c-act----------------cacaatccacac
droVir3 scaffold_13047:824149-824226 + TAC--CAACAACAA------ATA--------------Caaaaaaaaa--aaa-------------------------------------------aaaaaacaaa------------------------acaaaa-------a--aaaCAATACAATATGTAacactcgc-ataaacac
droMoj3 scaffold_6540:4819885-4819951 + TAC--------CAAAACAAAATA--------------CACTCAAACA--ACA-------------------------------------------ACAAAAAAACACAA------------------------GAATTATTCA--AAACAA-------------------TATAAacac
droGri2 scaffold_14906:5486805-5486840 - TACAAAAAAACCAACTA-----------------------------------------------------------------------------------------------------------------------------A--AAACAC----------------AAATACACATAC

20-way Fast Statistical Aligner re-alignment (Top LASTZ hits)

Species Coordinate Lastz Id Alignment
dm3 chr3R:20518203-20518338 + NA TACACCATGAA----CAACAT--AAGATTTA---------T-T-------------------------CGGTCT--ACAAACTGAAGCTTAATTGTGATTTATG-G------------------------------------ATGTGTATTTT--T--TA---------------------GC--A--GTTAAAAA-----------------TCACAATTTAGCCAC-------------TTTGTAG-GTCGATGAAA-----------------------------CAGC----------TCAC-AACTCTCAC
droSim1 chr3R_random:984803-984937 + 64 TACACCATGAA----CAACAT--AAGATATA---------T-T-------------------------CGGTCT--ACAAACTGAAGCTTAATTGTGATTTAAG-G------------------------------------ATGTGTATTTT--T--TA---------------------GC--A--CTAAAAAA-----------------TCACAATTTAGCCAA-------------TTTGTAG-GTCGAT-AAA-----------------------------AACC----------ACAC-AACTCTCAC
droSec1 super_0:20852363-20852498 + 83 TACACCATGAA----CAACAT--AAGATATA---------T-T-------------------------CGGTCT--ACAAACTGAAGCTTAATTGTGATTTAAG-G------------------------------------ATGTGTATTTT--T--TT---------------------GC--A--CTAAAAAA-----------------TCACAATTTAGCCAA-------------TTTGTAG-GTCGATGAAA-----------------------------AACC----------ACAC-AACTCTCAC
droYak2 chr3R:25978425-25978557 - 899 TACACCA---A----CAACAT--AAGATATA---------T-T-------------------------CGGCAT--ACAAACTTGAGCTTAATTGTGATTTAAG-G------------------------------------ATGTGAATTTC-----AA---------------------GC--A--ACGAAAAA-----------------TCACAATGAAGCCAA------------TTTTGTAG-GCCGATGGAA-----------------------------CACC----------TCAC-AACTCTCAC
droEre2 scaffold_4820:7539900-7540034 - 2876 CACACCATGAA----CAACAT--AAGATATA---------T-T-------------------------CGGCCT--ACAAACATAAGCATAATTGTGATTTAAG-G------------------------------------ATGTATATTTC-----AA---------------------GT--A--GTAATACA-----------------TCACAATTTAGCCAA-------------TTTGTAG-GCCGATGAAA-----------------------------AACC----------TCAC-AACTCTCAC
droEug1 scf7180000409881:243948-244075 + 465 TACACCATGAA----CAACAT--AAGCTATA---------T-C-----------------------TCCGGTCT--ACAAGTTGCAGATTTATTA--------------C------------------------TGATTTCT--------A--------T------------------ATAGA--A--A---CGAA-----------------TAATAATTGCGAAAA----------TCACTTGTAG-GCCGGAAAAA-----------------------------CCAC----------TCAC-AACTCTCAC
droBia1 scf7180000302402:4298359-4298484 + 251 TACACCATGAA----CAACAT--ATG-TATA---------T-T-------------------------CGGCCT--ACAGGCTGGAGCTCAATAGTGATTCGAG--------------------------------------ATTTGA-----------T---------------------TT--A--GTAGAGAA-----------------TCACTACTGGGTT--------------TCCTGTAG-GCCGGAGGAA-----------------------------CACC----------TCAC-AACTCTCAC
droTak1 scf7180000415310:621230-621363 - 131 TACACCATGAA----CAACAT--AAGATATG---------T-T-------------------------CGGCCT--ACAAGTTGAAGGTCAACAGTGATTTGAG-A------------------------------------ATAAATGTTTG-----AA---------------------AT--G--GTAGAAAA-----------------TCACTGTTGGGTT--------------TCCTGTAG-GCCGGAAAAA-----------------------------CACC----------TCAC-AACTCTCAC
droEle1 scf7180000491194:574957-575084 + 690 TACACCATGAA----CAACAT--AAGATATA---------T-T-------------------------CGGCCT--ACAAGAAGAGGCTCAATTGTGATTTAAAGG------------------------------------AAATGT-----------G---------------------GA--A--ATGGAAAG-----------------TCACAATTGGGTC--------------TCCTGT-G-GCCGGAAAAA-----------------------------CACC----------TCAC-AACTCTCAC
droRho1 scf7180000779488:616949-617082 + 160 TACACCATGAA----CAACAT--AAGATATA---------T-C-------------------------CGG-CT--ACAAGACGAGGCTCAATTGTGATTTGAA-A------------------------------------AGAA-TAGT-CAAA--TG---------------------GA-TA--ATGGTAAA-----------------TCACAATTGGGTC--------------TCCTGT-G-GCCGGAAAAA-----------------------------CACC----------TCAC-AACTCTCAC
droFic1 scf7180000454106:1794197-1794333 - 430 TACACCATGAA----CAACAT--AAGAAAAA---------T-C-------------------------TGATCT--ACAAACCGAAGCTAACTTGTGATTTAAAGGATAA---------------------------------TAAATAATT---T--TA---------------------AT-----ATGGAAAA-----------------TCACAATTTAGCT--------------TCCTGTAG-GTTGGAAAGA-----------------------------CACC----------CCAC-AACTCTCAC
droKik1 scf7180000302610:1723162-1723292 - 359 CAAG--------------------AGTTTTA---------T-T-------------------------TTATTT--ACAAGCTAAAATCTAATTTTAAGCTAAA-G------------------------------------ATTTTAGTTAT--T--TA---------------------AA--A--G-------------------------------------------------------------TTAAAAAAATATTTAAAAATAATTTGAAACTAAATAAAAACCGAATTTCAAGTCTA-AGTTTACAA
droAna3 scaffold_12929:3171299-3171434 + 53 AATACCAGACA----CAAAAT--ACACAAAA---------TACTAAAAATGGTATACCAACGTTATTTTTAGCTATCGATCATTATTTCTC----------------CACCAAATCCTCAGTATTTTGATTTGGT-T-----CTAT---------------------------------------------------------------------------------------------------TAT-TTTC---GAA-----------------------------AAAC----------ACAC-CATACACAC
droBip1 scf7180000396696:391796-391915 + 164 AACACAAA-------CAAAAT--GGA----A---------C-T-------------------------CAAGTC--AAAAATATAAATGAAACAGAGTATTAGA-A------------------------------------ATGGGCATTCG-----A----------------------GT-TGAACAAGAAAA---TCGAAAAATTGAAAAGACGATTAAA-------------------------------------------------------------------C----------TTTT-GAGTGAGGC
dp4 chr2:1714385-1714528 - 545 TAGACCATGAA----CAACAT--AAGAAAAACTTACACAC---------------------GTTTCTCTGGTTTATACAAAATCAAGTTTCATTGTGATTGGAG-C------------------------------------CTTA-T--ATT--T--TG---------------------GGATA--A--TCCAA-----------------TCACAATCAAGCT--------------TCCTGTACAACCCATGAGA--------------------------------------------GAC-AATTCACAC
droPer1 super_7:1868921-1869065 - 1214 TACACCATGAA----CAACAT--AAGAAAAACTTACACAC---------------------GTTTCTCTGGTTTATACAAAATGAAGTTTCATTGTGATTGGAG-C------------------------------------CTATATGTT----C--TG---------------------GGATA--A--TCCAA-----------------TCACAATCAAGCT--------------TCCTGTACAACCCATGAGA--------------------------------------------GAC-AATTCACAC
droWil1 scaffold_180955:2505632-2505764 + 1167 AAAAGAAAAAA----AAACAT--AATAATAT---------T-T-------------------------ATGTGT--ACACTGTGAGTTTTTATTATGTGTTTTG-T------------------------------------TTTTGCTTTTA--T--TT-------------------------A--T------T-----------------TATGCGTCTAGGCGCATTGCATTTC--TCC-------AAATCCAAA-----------------------------CAC-----------ACACACACACACAC
droVir3 scaffold_12928:1372331-1372464 - 7630 TACT--ATGAAGATACAACAT--ATA-TATA----TA---T-A--T----------------------ATATAG--ATATATAGAGATTTATGTGAAACCTATT-G------------------------------------TAGTGTATTCT--T--GA---------------------GC--T--GGTTAAAACAAAGTTTTAGCTGGCA-----------TT--------------TGGTATATGACCTATG----------------------------------------------------TAGTAGTT
droMoj3 scaffold_6496:17122051-17122183 + 3688 TATATTTTTA-----TAATAC--ATA-AAAA---------T-TCAT---------------------ATTACATATATTAATATTATTATTAATAT----TATT-ATTAT--------------------------------CTAATTTTT----T--TA---------------------AT-TTA-ATAAAATA-----------------TTACAAAT---TA--------------TAATTTATAATATATGTAA--------------------------------------------TAGATTTTTAAAT
droGri2 scaffold_15081:1647423-1647554 + 3061 AACAAAAAAAA----CAAGAAATAACAACTA---------T-T-------------------------CAGAAA--GGAGAAAACAAATTT---------TGTG-GTTGC------------------------TTAGTTCA--------A----TTTCAGTGCAAATGTTTTTATTT---AT--G--CTTG-------------------------------------------------GCAGTCA-GTCATACACA-----------------------------CAC-----------TCACACACTCACAC

[View All Possible Orthologs]

Original 12-way multiz alignment with reads super-imposed

Species Alignment
dm3
TACACCATGAACAACATAAGATT--------------TATTCGGTCT--ACAAACTGAAGCTTAATTGTGATTTATGGATGTGTATTTTTTAGCAGTTAAAAATCACAATTTAGCCACT-TTGTAG-GTCGATGA-------A--ACAGCT----------------CACAACTCTCACSizeHit CountTotal Norm
droSim1
TACACCATGAACAACATAAGATA--------------TATTCGGTCT--ACAAACTGAAGCTTAATTGTGATTTAAGGATGTGTATTTTTTAGCACTAAAAAATCACAATTTAGCCAAT-TTGTAG-GTCGAT-A-------A--AAACCA----------------CACAACTCTCACSizeHit CountTotal NormV044
Embryo
M023
Head
M024
Male body
M025
Embryo
....................ATA--------------TATTCGGTCT--ACAAAC............................................................................................................................19141.000.0028.0013.000.00
...................GATA--------------TATTCGGTCT--ACAAAC............................................................................................................................2014.000.003.001.000.00
.......................................................TGAAGCTTAATTGTGATTTAAGGAT...................................................................................................2513.000.001.001.001.00
..............................................................................................................TTAGCCAAT-TTGTAG-GTC.................................................1812.000.000.001.001.00
.....................TA--------------TATTCGGTCT--ACAAAC............................................................................................................................1812.000.002.000.000.00
.........................................................................TAACGATGTGTTTTTTTT........................................................................................1811.001.000.000.000.00
..............CATAAGATA--------------TATTCGGTCT--ACA...............................................................................................................................2211.000.000.000.001.00
....................ATA--------------TATTCGGTCT--ACAGAC............................................................................................................................1911.000.001.000.000.00
.......................................................TGAAGCTTAATTGTGATTTAAGGG....................................................................................................2411.000.000.001.000.00
.................AAGATA--------------TATTCGGTCT--ACAAAC............................................................................................................................2211.000.000.001.000.00
..........................................................................AAGGATGTGTATTTTTTAGCA....................................................................................2111.000.001.000.000.00
...............ATAAGATA--------------TATTCGGTCT--ACAAAC............................................................................................................................2411.000.001.000.000.00
droSec1
TACACCATGAACAACATAAGATA--------------TATTCGGTCT--ACAAACTGAAGCTTAATTGTGATTTAAGGATGTGTATTTTTTTGCACTAAAAAATCACAATTTAGCCAAT-TTGTAG-GTCGATGA-------A--AAACCA----------------CACAACTCTCACSizeHit CountTotal NormV113
Male body
V114
Embryo
V115
Head
........................................................................TTTAGGATTTGTATTTTT.........................................................................................1812.000.002.000.00
............................................................................................................ATTTAGCCAAT-TTGTAG-GTC.................................................2011.000.001.000.00
....................................................................................ATTTGTTCGCACTAAAAA.............................................................................1811.000.001.000.00
..........................................................................................................................GTGG-GTCGATGA-------A--GAACC.............................1811.000.000.001.00
droYak2
TACACCA---ACAACATAAGATA--------------TATTCGGCAT--ACAAACTTGAGCTTAATTGTGATTTAAGGATGTGAATT-TCAAGCAACGAAAAATCACAATGAAGCCAATTTTGTAG-GCCGATGG-------A--ACACCT----------------CACAACTCTCACSizeHit CountTotal NormV046
Embryo
V052
Head
V120
Male body
V058
Head
M026
Head
.......................................................TTGAGCTTAATTGTGATTTAAGG.....................................................................................................2311.000.000.001.000.000.00
.........................................................................TAAGGATGTGAGTT-TCAAGCAA...................................................................................2211.001.000.000.000.000.00
.......................................................TTGAGCTTAATTGTGATTTAAG......................................................................................................2211.000.000.001.000.000.00
droEre2
CACACCATGAACAACATAAGATA--------------TATTCGGCCT--ACAAACATAAGCATAATTGTGATTTAAGGATGTATATT-TCAAGTAGTAATACATCACAATTTAGCCAAT-TTGTAG-GCCGATGA-------A--AAACCT----------------CACAACTCTCACSizeHit CountTotal NormV040
Embryo
V060
Head
V047
Embryo
V108
Head
droAna3
TACACCATGTCTGACTTGA--------------------------------------AAATTTAATTGTGATTTGGTTT------------------TTTGAATCACAATTAAGTTT---TTTTAA-GCCGGA-A-------AACACACAT----------------TTCAACTCTCACSizeHit CountTotal NormV105
Male body
V106
Head
V055
Unknown
V039
Embryo
..............CTTGA--------------------------------------AAATTTxATGGTG.............................................................................................................1711.000.000.001.000.00
....................................................................................................................................A-A-------AACACACAT----------------TTxAAATCTCAC2211.000.000.000.001.00
dp4
TAGACCATGAACAACATAAGAAAAACTTACACACGTTTCTCTGGTTTATACAAAATCAAGTTTCATTGTGATTGGAGCCT---TATATTTTGGGATAATCCAATCACAATCAAGCTTC--CTGTACAACCCATGA-------G--AG---------------------ACAATTCACACSizeHit CountTotal NormV112
Male body
V043
Embryo
V051
Head
M022
Male body
.AGACCATGAACAACATAAGA..............................................................................................................................................................2011.000.001.000.000.00
...............................................................................................................AAGCTTC--CTGTACAACCCATGA-------G--A.................................2411.000.001.000.000.00
droPer1
TACACCATGAACAACATAAGAAAAACTTACACACGTTTCTCTGGTTTATACAAAATGAAGTTTCATTGTGATTGGAGCCTATAT--GTTCTGGGATAATCCAATCACAATCAAGCTTC--CTGTACAACCCATGA-------G--AG---------------------ACAATTCACACSizeHit CountTotal NormV111
Male body
V042
Embryo
V050
Head
V057
Head
M021
Embryo
............................................................ATTCATTGTGACTGGAGCC....................................................................................................1916.000.000.000.000.006.00
............................................................................................................ATCAAGCTTC--CTGTACAACCCATGA-------G--A.................................2712.000.000.000.000.002.00
......ATGAACAACGTAAAAAAAA..........................................................................................................................................................1912.000.000.000.000.002.00
droWil1
caaagac---------------a---------------------------------------------------------------------------------------------------------------------------c-act----------------cacaatccacacSizeHit CountTotal NormV117
Male body
V118
Embryo
V119
Head
M020
Head
droVir3
TAC--CAACAACAA------ATA--------------Caaaaaaaaa--aaa-------------------------------------------aaaaaacaaa------------------------acaaaa-------a--aaaCAATACAATATGTAacactcgc-ataaacacSizeHit CountTotal NormV116
Male body
V053
Unknown
V047
Embryo
M027
Male body
M028
Head
......................................AAAAAAAAA--AAA-------------------------------------------AAAAAAAAAA------------------------AAAAAA-------A--AAACAA............................3582.250.002.130.130.000.00
......................A--------------AAAAAAAAAA--AAA-------------------------------------------AAAAAAAAAA------------------------ACAAAA-------A--AAACA.............................3662.830.002.500.330.000.00
........................................AAAAAAA--AAA-------------------------------------------AAAAAAAAAA------------------------A.................................................21470.280.000.000.000.210.06
......................................AAAAAAAAA--AAA-------------------------------------------AAAAAAAAA...........................................................................21470.280.000.000.000.210.06
.......................................AAAAAAAA--AAA-------------------------------------------AAAAAAAAAA------------------------..................................................21470.280.000.000.000.210.06
......................................AAAAAAAAA--AAA-------------------------------------------AAAAAAAAAA------------------------A.................................................23360.330.000.000.000.080.25
......................................AAAAAAAAA--AAA-------------------------------------------AAAAAACAAA------------------------ACAAAA-------A--AAAAAAA...........................3619.000.001.008.000.000.00
......................................AAAAAAAAA--AAA-------------------------------------------AAAAAAAAAA------------------------ACAAAA-------A--AAACAAA...........................3633.000.002.330.670.000.00
.....................................AAAAAAAAAA--AAA-------------------------------------------AAAAAAAAAA------------------------ACAAAA-------A--AAACAA............................3651.600.001.000.600.000.00
....................AAA--------------AAAAAAAAAA--AAA-------------------------------------------AAAAAACAAA------------------------ACAAAA-------A--AAA...............................3624.000.002.501.500.000.00
......................................AAAAAAAAA--AAA-------------------------------------------AAAAAACAAA------------------------AAAAAA-------A--AAACAAA...........................3632.330.001.331.000.000.00
.......................................AAAAAAAA--AAA-------------------------------------------AAAAAAAAAA------------------------A.................................................22410.150.000.000.000.020.12
......................................AAAAAAAAA--AAA-------------------------------------------AAAAAAAAAA------------------------..................................................22410.150.000.000.000.020.12
.......................................AAAAAAAA--AAA-------------------------------------------AAAAAAAAAA------------------------ACAAAA-------A--AAACAAAA..........................3632.000.001.330.670.000.00
......................A--------------AAAAAAAAAA--AAA-------------------------------------------AAAAAACAAA------------------------AAAAAA-------A--AAACA.............................3641.250.001.250.000.000.00
.............A------AAA--------------AAAAAAAAAA--AAA-------------------------------------------AAAAAACAAA------------------------ACAAAA-------A--AA................................3622.500.002.000.500.000.00
......................................AAAAAAAAA--AAA-------------------------------------------AAAAAAAAAA------------------------ACAAAA-------A--AAAAAA............................3560.830.000.670.170.000.00
.......................................AAAAAAAA--AAA-------------------------------------------AAAAAACAAA------------------------AAAAAA-------A--AAACAAAA..........................3631.670.001.000.670.000.00
......................A--------------CAAAAAAAAA--AAA-------------------------------------------AAAAAAAAAA------------------------AAAAAA-------A--AAACA.............................3650.600.000.000.600.000.00
.......................................AAAAAAAA--AAA-------------------------------------------AAAAAAAAAA------------------------AAAAAA-------A--AAACAATA..........................3631.000.000.670.330.000.00
............AA------AAA--------------AAAAAAAAAA--AAA-------------------------------------------AAAAAACAAA------------------------ACAAAA-------A--A.................................3640.750.000.250.500.000.00
droMoj3
TAC--------CAAAACAAAATA--------------CACTCAAACA--ACA-------------------------------------------ACAAAAAAACACAA------------------------GAATTATTCA--AAACAA-------------------TATAAacacSizeHit CountTotal NormV110
Male body
V056
Unknown
V049
Head
V041
Embryo
droGri2
TACAAAAAAACCAACTA-----------------------------------------------------------------------------------------------------------------------------A--AAACAC----------------AAATACACATACSizeHit CountTotal NormV109
Male body

12-way Secondary Structure

dm3 droSim1 droSec1 droYak2 droEre2 droAna3 dp4 droPer1 droWil1 droVir3 droMoj3 droGri2
RNAfold Image
UNAFold Image None None None

20-way Secondary Structure

dm3 droSim1 droSec1 droYak2 droEre2 droEug1 droBia1 droTak1 droEle1 droRho1 droFic1 droKik1 droAna3 droBip1 dp4 droPer1 droWil1 droVir3 droMoj3 droGri2
RNAfold Image
UNAFold Image

Generated: 04/20/2012 at 05:19:57 EDT