Return to list

View alignment of MAF (12-files) vs LASTZ top hits (20-flies)

View all possible LASTZ hits (20-flies)

dme-mir-4957

chrX:15213986-15214070 -

Known_mirbase

[View on miRBase] [View on UCSC Genome Browser {Cornell Mirror}]

intergenic


Legend:

   mature strand at position 15214004-15214024
   star strand at position 15214034-15214054
   mutations
   read mismatches

20-way Fast Statistical Aligner (FSA) re-alignment

Visually pleasing best hits from LASTZ and MAF


Species Coordinate Source Confidence Alignment
dm3 chrX:15213971-15214085 - 1 CTTTTTT-TTCT----CTG----CTCG--------------------------TTGTATTTTACGCACTCTGCTGGCAGCGCTGTCG------GCGTCGCTGCCAGCTGCCAGCAGC------------------------------------GAGCG--------------GT--ACACTTGCAACTA---------TCCTGTGCACCCT---ACAATT
droSim1 chrX:11736931-11737042 - LASTZ 4 CTTTTTT-TTTTTTCTCTG----CTCG--------------------------TTGCATTTTACGCACTCTGCTGGCAGCGACGCCG------ACAGCGCTGCCAGCA-------GC------------------------------------GAGCG--------------GT--ACACTTGCAACTA---------TCCTGCGCACCCT---ACAATT
droSec1 super_56:101687-101796 - LASTZ 4 CTTTTTT-TTTT--CTCTG----CTCG--------------------------TTGCATTTTACGCACTCTGCTGGCAGCGACGCCG------ACAGCGGTGCCAGCA-------GC------------------------------------GAGCG--------------GT--ACACTTGCAACTA---------TCCTGCACACCCT---ACAATT
droYak2 chrX:9414156-9414285 - LASTZ 5 TTTTTTT-------CTCTG----CTCGTTTT---TCTCTCTTTTTCTTCTCGCCTGCATTTTACGCACTCTGCTAGCAGCGACGCCG------ACTGCGCTGCCAACA-------GC------------------------------------AAGCG--------------GT--ACACTTGCAACTA---------TTCTATGCACCCT---ACAATT
droEre2 scaffold_4690:11807313-11807440 + LASTZ 5 TTTTTTC-T-------CCG----CTCGTTTT---TCTCTCCTTTTCTGCTCGCCTGCATTTTACGCACTCTGCTGGCAGCGACGCCG------ACTGCGCTGCAAACA-------GC------------------------------------GAGCG--------------GT--ACACTTGCAACTA---------TTCTGTGCACCCT---ACAATT
droEug1 scf7180000409533:350907-351023 - LASTZ 5 TTTTTTC-TCT-----CTG----CTAGTTTTTTTT----------------GCTTGCATTTTACGCACTCTGCTGGCAGCGACGTCG------ACTGCGCTGCCAAGT-------GC------------------------------------AAGCG--------------GT--ACACTTGCAACTA---------TTCGTTGTACTCC----GAATT
droBia1 scf7180000299226:181007-181106 + LASTZ 5 CTTCCCC-TCCC----CCG----CCCG------------------CCCCGCGGCAGCGTTTTACGCACTCTGCCGGCAGCGACGTCG------ACCGC-----------------------------------------------------------G--------------GT--ACACTTGCCACAA---------TTCCTTGTACCCC---CGAGAG
droTak1 scf7180000415381:822078-822183 - LASTZ 5 TTTTTTT-------CTCAG----CC---------------------------CCTGCATTTTACGCACTCCGCTGGCAGCGACGTCG------ACTGCGCTGCCGACG-------TC------------------------------------AGGC-----------GGCGGT--ACACTTGCCACTA---------TTCTTTGTACCC------CCGT
droEle1 scf7180000491023:1386855-1386967 + LASTZ 5 TTTTTTT-TTTT----TTGGGG-----TTTTTTTT----------------GCCTGCATTTTACGCACTCTGCTGGCAGCGACGTCG------ACTGCGACGCCGACT-------GC------------------------------------A-GCC--------------GT--ACACTTGCAACTA---------TTCGTTGTACCC------AATT
droRho1 scf7180000777097:282566-282669 + LASTZ 5 TTTTTTT-TGCT---------TTTTC--------------------------GCTGCATTTTACGCACTCTGCTGGCAGCGACGTCG------ACTGCGACGGCGACT-------GC------------------------------------AGGCG--------------GT--ACACTTGCAACTA---------TTCTGTGTACCC------AATT
droFic1 scf7180000454038:448768-448871 - LASTZ 5 TTTTTTG-TTGC----T----TTT----------T----------------GCCTGCATTTTACGCACACTGCTGGCAGCGACGTCG------ACTGCGACGCCGACC-------GC------------------------------------GAGT---------------GT--ACACTTTCAACTA---------TTTCATGTACCCC---ACA---
droKik1 scf7180000302816:65604-65724 - LASTZ 0 CTTTTGC-TTTG--CTATG----CCTGTGG-CTCTGC-----CCGCAATTTGTCCGTGCCGTACGCAC--TGCGGGCAGAGACGCCG------ACTGCGTAGCA------------------------------------------------------------------------CCACTTGCAACTGGGTTATTTCTT--TTGTATT-G---TCTTTT
droAna3 scaffold_13117:1427273-1427395 - LASTZ 0 ATTTTTC-TTCC----C-A---TG----------CAACCTGTGTCCAT----TGTACGTTTTACGCACTCTCC--CCAGCAACACCG------CCGCCGATAGCGACT-------GC------------------------------------GACGAGTACGTTGGCAGCGG----CACTTGACTCT--------------TTTT--GCG---GCAATT
droBip1 scf7180000396425:686944-687063 + LASTZ 0 GCATTCA-ACCT----GTG----TCCA--------------------T----TTTACGTTTTACGCACTCAGCC--CAGCAACAACGCCGTCGAGTGC------GACT-------GC------------------------------------GACGAGTACGCTGGCAGCGG----CACTTGACTCTT-------TTTGC--GGCAATTT---TCGACA
dp4 chrXL_group1e:12152468-12152588 - LASTZ 0 CTTTAGC-TTTT-----TGCGT-----TTGCGTTTGCCCCGATTTCAATGCA-ATGTGTTGTACGCACACGGGTT---TTTATCTCA------GCTGTGGGATGAACG-------GC--------------------------------------CCG--------------AG--ACACTTGC--CGA---------TGCGGCGACCCC------AGGC
droPer1 super_14:1456573-1456693 + LASTZ 0 CTTTAGC-TTTT-----TGCGT-----TTGCGTTTGCCCCGATTTCAATGCA-ATGTGTTGTACGCACACGGGTT---TTTATCTCA------GCTGTGGGATGAACG-------GC--------------------------------------CCG--------------AG--ACACTTGC--CGA---------TGCGGCGACCCC------AGGC
droWil1 Unknown 0 ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
droVir3 scaffold_12928:4095173-4095244 - LASTZ 0 CTCTTTA-TTGC----A----ATG----------C----------------GAC-------TACGCATG-------------------------------------AT-------TT------------------------------------A-GGCGCAGGTTGGT----CG--ACACTTGCATTGG---------TTGCAAGT----------ACTT
droMoj3 scaffold_6680:3750941-3751063 - LASTZ 0 TTTGTTT-TTGT----TTT----CTTT--T----------------TGTAGCCTAGGAATTTACC---TCTTGTGGCATATCAGCTG---TTGGCTTTGGTTGCAAAC-------AC------------------------------------ACGTT-------GACT---TTTCACACACAGAATTC---------TGT--TGAAATGTTTAACACTT
droGri2 scaffold_14853:566008-566117 + LASTZ 0 CTCTTTTATTGC----A----ATG----------A----------GACTGCG-ACGCAT----------------------------------------------GAT-------TTGGGCGCAGGTCTATCACCACTTGCACTGACGACTGCTTTTG--------------CT--TTTGCTACTACTT---------TTACTTGCCAGCG---GAACTT

Legend

0- The ortholog is absent. A better synthenic sequence candidate could not be found.
1- Reads w/ good pileup on mature and star strand & folds into hairpin structure.
2- Reads w/ good pileup on mature and star strand BUT does not fold into hairpin structure.
3- Sporadic reads & folds into a hairpin structure.
4- No reads but folds into a hairpin structure.
5- Ortholog is confident but miRNA appears absent (no reads and does not fold into a hairpin)

Re-alignment after removal of non-confident orthologs

Species Alignment
dm3 CTTTTTT----TTCTCTGCTCG----------------------------TTGTATTTTACGCACTCTGCTGGCAGCGCTGTCGGCGTCGCTGCCAGCTGCCAGCAGCGA---GCGGTACACTTGCAACTATCCTGTGCACCCTACAATT
droSim1 CTTTTTTTTTTTTCTCTGCTCG----------------------------TTGCATTTTACGCACTCTGCTGGCAGCGACGCCGACAGCGCTGCCAGCAG-------CGA---GCGGTACACTTGCAACTATCCTGCGCACCCTACAATT
droSec1 CTTTTTT--TTTTCTCTGCTCG----------------------------TTGCATTTTACGCACTCTGCTGGCAGCGACGCCGACAGCGGTGCCAGCAG-------CGA---GCGGTACACTTGCAACTATCCTGCACACCCTACAATT
droYak2 TTTTTT------TCTCTGCTCGTTTTTCTCTCTTTTTCTTCTC-----GCCTGCATTTTACGCACTCTGCTAGCAGCGACGCCGACTGCGCTGCCAACAG-------CAA---GCGGTACACTTGCAACTATTCTATGCACCCTACAATT
droEre2 TTTTT-------TCTCCGCTCGTTTTTCTCTCCTTTTCTGCTC-----GCCTGCATTTTACGCACTCTGCTGGCAGCGACGCCGACTGCGCTGCAAACAG-------CGA---GCGGTACACTTGCAACTATTCTGTGCACCCTACAATT
droEug1 TTTTTTC-----TCTCTGCTAGTTT------------------TTTTTGCTTGCATTTTACGCACTCTGCTGGCAGCGACGTCGACTGCGCTGCCAAGTG-------CAA---GCGGTACACTTGCAACTATTCGTTGTACTCC-GAATT
droBia1 CTTCCCC----TCCCCCGCCCG--C------------------CCCGCGGCAGCGTTTTACGCACTCTGCCGGCAGCGACGTCGACCGC--------------------------GGTACACTTGCCACAATTCCTTGTACCCCCGAGAG
droTak1 TTTTTT------TCTCAGCC-----------------------------CCTGCATTTTACGCACTCCGCTGGCAGCGACGTCGACTGCGCTGCCGACGT-------CAGGCGGCGGTACACTTGCCACTATTCTTTGTACCCC---CGT
droEle1 TTTTTTT----TTTTTTG-GGGTTT------------------TTTTTGCCTGCATTTTACGCACTCTGCTGGCAGCGACGTCGACTGCGACGCCGACTG-------C-A---GCCGTACACTTGCAACTATTCGTTGTACCC---AATT
droRho1 TTTTTTT----TGC-----------------------------TTTTTCGCTGCATTTTACGCACTCTGCTGGCAGCGACGTCGACTGCGACGGCGACTG-------CAG---GCGGTACACTTGCAACTATTCTGTGTACCC---AATT
droFic1 TTTTTTG----TTG----------C------------------TTTTTGCCTGCATTTTACGCACACTGCTGGCAGCGACGTCGACTGCGACGCCGACCG-------CGA---GT-GTACACTTTCAACTATTTCATGTACCCCACA---

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 None
UNAFold Image None

Small RNA-seq read pileup

Species Alignment
droSim1
CTTTTTT-TTTTTTCTCTG----CTCG--------------------------TTGCATTTTACGCACTCTGCTGGCAGCGACGCCG------ACAGCGCTGCCAGCA-------GC------------------------------------GAGCG--------------GT--ACACTTGCAACTA---------TCCTGCGCACCCT---ACAATTSizeHit CountTotal NormV044
Embryo
M023
Head
M024
Male body
M025
Embryo
..........................................................TTTTACGCACTCTGCTGGCAGCG...........................................................................................................................................2311.000.001.000.000.00
.................................................................................................CGCTGCCAGCA-------GC------------------------------------GAGCG--------------GT--AC..........................................2211.000.000.000.001.00
...........................................................................................................A-------GC------------------------------------GAGCG--------------GT--ACACGGGC....................................1811.000.001.000.000.00
...............TCTG----CGCG--------------------------TTGCATTTTACGCA.........................................................................................................................................................2211.001.000.000.000.00
...........................................................................................................................................................................................................................C1830.330.000.330.000.00
............................................................TTACGCACTCTGCTGGCA..............................................................................................................................................1811.000.001.000.000.00
.............................................................TACGCACTCTGCTGGCAGCG...........................................................................................................................................2011.000.001.000.000.00
...........................................................TTTAAGCACTCTGATGGC...............................................................................................................................................1811.000.000.001.000.00
.............................................................TACGCACTCTGCTGGCAGCGA..........................................................................................................................................2111.000.001.000.000.00
............................................................................................................................................................CG--------------GT--ACACTTGCAACTA---------TCCT..................2111.000.000.001.000.00
.TTTTTT-TTTTTTCTGTG----ATC..................................................................................................................................................................................................2011.000.000.000.001.00
............................................................TTACGCACTCTGCTGGCAGCG...........................................................................................................................................2111.000.000.001.000.00
..........................................................................GGCAGCGACGTCG------ACTGC..........................................................................................................................1811.000.000.000.001.00
.............................................................TACGCACTCTGCTGGCAGCGACG........................................................................................................................................2311.000.001.000.000.00
droSec1
CTTTTTT-TTTT--CTCTG----CTCG--------------------------TTGCATTTTACGCACTCTGCTGGCAGCGACGCCG------ACAGCGGTGCCAGCA-------GC------------------------------------GAGCG--------------GT--ACACTTGCAACTA---------TCCTGCACACCCT---ACAATTSizeHit CountTotal NormV113
Male body
V114
Embryo
V115
Head
...................................................................................................GTGCCAGCA-------GC------------------------------------xCGCG--------------GT--A...........................................1814.001.002.001.00
...........................................................................GCAGCGACGCCG------ACAGCGGTG......................................................................................................................2112.000.000.002.00
.............................................................TACGCACTCTGCTGGCAGCGAC.........................................................................................................................................2211.000.001.000.00
..........................................................................GGCAGCGACGCCG------ACAGCGGT.......................................................................................................................2111.000.000.001.00
......................................................................................................CCAGGA-------GC------------------------------------GAGGG--------------GT--ACA.........................................1811.000.001.000.00
...........................................................................GCAGCGACGTCG------ACAGCGGTGCCG...................................................................................................................2411.001.000.000.00
.........................................................................TGGCAGCGACGCCG------ACAGCGT........................................................................................................................2111.001.000.000.00
.............................................................TACGCACTCTGCTGGCAGAGACGCCG------ACA............................................................................................................................2911.001.000.000.00
.........................................................................TGGCAGCGACGCCG------ACAGCGG........................................................................................................................2111.000.000.001.00
.............................................................TACGCACTCTGCTGGCAGCGA..........................................................................................................................................2111.001.000.000.00
.............................................................................................ACAGCAGTGCGAGCA-------GC------------------------------------G..................................................................1830.330.000.000.33
droYak2
TTTTTTT-------CTCTG----CTCGTTTT---TCTCTCTTTTTCTTCTCGCCTGCATTTTACGCACTCTGCTAGCAGCGACGCCG------ACTGCGCTGCCAACA-------GC------------------------------------AAGCG--------------GT--ACACTTGCAACTA---------TTCTATGCACCCT---ACAATTSizeHit CountTotal NormV046
Embryo
V052
Head
V120
Male body
V058
Head
M026
Head
droEre2
TTTTTTC-T-------CCG----CTCGTTTT---TCTCTCCTTTTCTGCTCGCCTGCATTTTACGCACTCTGCTGGCAGCGACGCCG------ACTGCGCTGCAAACA-------GC------------------------------------GAGCG--------------GT--ACACTTGCAACTA---------TTCTGTGCACCCT---ACAATTSizeHit CountTotal NormV040
Embryo
V060
Head
V047
Embryo
V108
Head
.......................................................................................................AAACA-------GC------------------------------------GAGCG--------------GT--ACACTTG.....................................2111.000.001.000.000.00
........T-------CCG----TTTGTTTT---TCTCTCCTTT................................................................................................................................................................................2211.001.000.000.000.00
............................................................................CTGCGACGCCG------GCTGCGCT.......................................................................................................................1911.000.000.000.001.00
...........................................................................................................A-------GC------------------------------------GAGCG--------------GT--ACACGGGC....................................1811.000.000.000.001.00
droEug1
TTTTTTC-TCT-----CTG----CTAGTTTTTTTT----------------GCTTGCATTTTACGCACTCTGCTGGCAGCGACGTCG------ACTGCGCTGCCAAGT-------GC------------------------------------AAGCG--------------GT--ACACTTGCAACTA---------TTCGTTGTACTCC----GAATTSizeHit CountTotal Norm
droBia1
CTTCCCC-TCCC----CCG----CCCG------------------CCCCGCGGCAGCGTTTTACGCACTCTGCCGGCAGCGACGTCG------ACCGC-----------------------------------------------------------G--------------GT--ACACTTGCCACAA---------TTCCTTGTACCCC---CGAGAGSizeHit CountTotal Norm
droTak1
TTTTTTT-------CTCAG----CC---------------------------CCTGCATTTTACGCACTCCGCTGGCAGCGACGTCG------ACTGCGCTGCCGACG-------TC------------------------------------AGGC-----------GGCGGT--ACACTTGCCACTA---------TTCTTTGTACCC------CCGTSizeHit CountTotal Norm
droEle1
TTTTTTT-TTTT----TTGGGG-----TTTTTTTT----------------GCCTGCATTTTACGCACTCTGCTGGCAGCGACGTCG------ACTGCGACGCCGACT-------GC------------------------------------A-GCC--------------GT--ACACTTGCAACTA---------TTCGTTGTACCC------AATTSizeHit CountTotal Norm
droRho1
TTTTTTT-TGCT---------TTTTC--------------------------GCTGCATTTTACGCACTCTGCTGGCAGCGACGTCG------ACTGCGACGGCGACT-------GC------------------------------------AGGCG--------------GT--ACACTTGCAACTA---------TTCTGTGTACCC------AATTSizeHit CountTotal Norm
droFic1
TTTTTTG-TTGC----T----TTT----------T----------------GCCTGCATTTTACGCACACTGCTGGCAGCGACGTCG------ACTGCGACGCCGACC-------GC------------------------------------GAGT---------------GT--ACACTTTCAACTA---------TTTCATGTACCCC---ACA---SizeHit CountTotal Norm
droKik1
CTTTTGC-TTTG--CTATG----CCTGTGG-CTCTGC-----CCGCAATTTGTCCGTGCCGTACGCAC--TGCGGGCAGAGACGCCG------ACTGCGTAGCA------------------------------------------------------------------------CCACTTGCAACTGGGTTATTTCTT--TTGTATT-G---TCTTTTSizeHit CountTotal Norm
droAna3
ATTTTTC-TTCC----C-A---TG----------CAACCTGTGTCCAT----TGTACGTTTTACGCACTCTCC--CCAGCAACACCG------CCGCCGATAGCGACT-------GC------------------------------------GACGAGTACGTTGGCAGCGG----CACTTGACTCT--------------TTTT--GCG---GCAATTSizeHit CountTotal NormV105
Male body
V106
Head
V055
Unknown
V039
Embryo
......................................................................................................GCGGCG-------GC------------------------------------GACGAGTACG.........................................................1811.000.001.000.000.00
........................................................................................................GACT-------GC------------------------------------GACGCGTAGGTT.......................................................1811.000.000.001.000.00
droBip1
GCATTCA-ACCT----GTG----TCCA--------------------T----TTTACGTTTTACGCACTCAGCC--CAGCAACAACGCCGTCGAGTGC------GACT-------GC------------------------------------GACGAGTACGCTGGCAGCGG----CACTTGACTCTT-------TTTGC--GGCAATTT---TCGACASizeHit CountTotal Norm
dp4
CTTTAGC-TTTT-----TGCGT-----TTGCGTTTGCCCCGATTTCAATGCA-ATGTGTTGTACGCACACGGGTT---TTTATCTCA------GCTGTGGGATGAACG-------GC--------------------------------------CCG--------------AG--ACACTTGC--CGA---------TGCGGCGACCCC------AGGCSizeHit CountTotal NormV112
Male body
V043
Embryo
V051
Head
M022
Male body
...............................................................................TTATCTCA------GCTGTGGGATGA...................................................................................................................2011.000.001.000.000.00
......................................CCGATTTCAATGCA-ATGTGTT................................................................................................................................................................2111.001.000.000.000.00
............................TGCGTTTGCCCCGATTTCAAT...........................................................................................................................................................................2111.001.000.000.000.00
......C-TTTT-----TGCGT-----TTGCGTTTGCC......................................................................................................................................................................................2111.000.000.001.000.00
......................................................................................A------GCTGTGGGATGAACC-------GA--------------------------------------.................................................................1811.000.001.000.000.00
..................................................................................TGTCx------GCTGTGGGATGAAC.................................................................................................................1811.000.000.000.001.00
droPer1
CTTTAGC-TTTT-----TGCGT-----TTGCGTTTGCCCCGATTTCAATGCA-ATGTGTTGTACGCACACGGGTT---TTTATCTCA------GCTGTGGGATGAACG-------GC--------------------------------------CCG--------------AG--ACACTTGC--CGA---------TGCGGCGACCCC------AGGCSizeHit CountTotal NormV111
Male body
V042
Embryo
V050
Head
V057
Head
M021
Embryo
..................................................................................TGTCx------GCTGTGGGATGAAC.................................................................................................................1813.002.000.000.000.001.00
....................................................................................TCA------GCTGTGGGATGAACG-------G........................................................................................................1911.001.000.000.000.000.00
..................................................................................TCTCx------GCTGTGGGATGAACC-------.........................................................................................................1911.001.000.000.000.000.00
......................................................................................A------GCTGTGGGATGAACG-------GC--------------------------------------C................................................................1911.000.000.000.001.000.00
droWil1
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------SizeHit CountTotal NormV117
Male body
V118
Embryo
V119
Head
M020
Head
droVir3
CTCTTTA-TTGC----A----ATG----------C----------------GAC-------TACGCATG-------------------------------------AT-------TT------------------------------------A-GGCGCAGGTTGGT----CG--ACACTTGCATTGG---------TTGCAAGT----------ACTTSizeHit CountTotal NormV116
Male body
V053
Unknown
V047
Embryo
M027
Male body
M028
Head
...........................................................................................................T-------TT------------------------------------A-GGCGCAGGTTGGT----CG--AC..........................................2112.000.000.000.002.000.00
....................................................................................................................T------------------------------------G-GGCGCAGGTTGGT----CG--ACACT.......................................2211.001.000.000.000.000.00
droMoj3
TTTGTTT-TTGT----TTT----CTTT--T----------------TGTAGCCTAGGAATTTACC---TCTTGTGGCATATCAGCTG---TTGGCTTTGGTTGCAAAC-------AC------------------------------------ACGTT-------GACT---TTTCACACACAGAATTC---------TGT--TGAAATGTTTAACACTTSizeHit CountTotal NormV110
Male body
V056
Unknown
V049
Head
V041
Embryo
droGri2
CTCTTTTATTGC----A----ATG----------A----------GACTGCG-ACGCAT----------------------------------------------GAT-------TTGGGCGCAGGTCTATCACCACTTGCACTGACGACTGCTTTTG--------------CT--TTTGCTACTACTT---------TTACTTGCCAGCG---GAACTTSizeHit CountTotal NormV109
Male body

Generated: 04/20/2012 at 05:17:17 EDT