Return to list

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

View all possible LASTZ hits (20-flies)

dme-mir-276b-as

chr3L:10312885-10312983 -

Candidate

[View on UCSC Genome Browser {Cornell Mirror}]

intergenic


Legend:

   mature strand at position 10312906-10312927
   mutations
   read mismatches

Notes: TODO: This is bizarre. Why does the 12-way have perfect conservation vs. me not finding any with lastz


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

Visually pleasing best hits from LASTZ and MAF


Species Coordinate Source Confidence Alignment
dm3 chr3L:10312870-10312998 - 1 TT--------TTAA--CTAT--------A--TGCTGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATATAGGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTTTTT---------------------CG--TTC-GATGGTC
droSim1 chr3L:9702897-9703022 - LASTZ 5 TT--------TCAA--CTAT--------A--TGATGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATA--GGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTTTTC----------------------G--GTC-GATGGTC
droSec1 super_0:2526747-2526871 - LASTZ 5 TT--------TCAA--CTAT--------A--TGATGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATA--GGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTTTTC----------------------G--TTC--GTGGTC
droYak2 chr3L:10302128-10302256 - LASTZ 4 TT--------TTAA--CTTT--------A--TGATGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATA--GGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTTTTG-------------------TTCT--TTC-TGTTGTC
droEre2 scaffold_4784:10306167-10306292 - LASTZ 0 TT--------TCAA--CTTT--------A--TGATGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATA--GGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTTTTT----------------------A--TTC-GGTTGTC
droEug1 scf7180000409711:2174289-2174430 - LASTZ 0 TT--------TAAA--CTTT--------A--CGATGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----AAATA--GAAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGGCTTCGGTATTTTC------TTTGTTTTTCTTGGCCC--TTC-GATGGAC
droBia1 scf7180000302428:7436352-7436478 + LASTZ 0 TT--------TCAA--GTTT--------T--GGATGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATA--GGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTTTTG---------------------TA--TTC-GACGGTC
droTak1 scf7180000415209:361628-361756 - LASTZ 0 TT--------TCAA--CTTT--------A--TGATGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACG----GAATA--GGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGGCTTCGGTTATTTC------T-------------TCG--TTC-GATGGTC
droEle1 scf7180000490564:1145980-1146108 - LASTZ 0 TG--------TTAATTCTTT--------A--TGTTGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGATT----AAATG--GGAACGTAGGAACTCTATACCTCGCTGATGGTATAAAGACTTCGGTTTTCTG---------------------CT--TTG-AATGGTC
droRho1 scf7180000780105:56692-56817 + LASTZ 0 TG--------TAAA--CTTT--------A--TGTTGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATA--GGAACGTAGGAACTCTATACCTCGCTGATGGTAAA-AGACTTCGGATTTTTG---------------------GT--TTC-AATGGTC
droFic1 scf7180000453841:740991-741120 + LASTZ 0 TT--------TAGT--TTTT--------ATTTGATGGACGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATT--TGAACGTAGGAACTCTATACCTCGCTTATGGTAAAAAGGCTTCGGATTTGGG---------------------GT--TTCTGATGGTC
droKik1 scf7180000302334:94650-94789 + LASTZ 0 CTCTATAATAAACA--CTTT--------A--TGTTGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACTATTTG---A--GGAACGTAGGAACTCTATACCTCGCTGATGGCAAAAAGGCTTCGTTTGTTTT-----------------TGGCTT--TTC-AGTCTTC
droAna3 scaffold_13337:15233368-15233487 - LASTZ 0 CT--------T-------TC--------A--TGATGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACTAGTTTAAGT--GGCACGTAGGAACTCTATACCTCGCTGATGGCAAATAGGCTTCGGTTTCTTC----------------------A--------GAAAT
droBip1 scf7180000395933:88091-88219 - LASTZ 0 GA--------AAGA--CTTT--------C-TCGATGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACTAGTTAAAGT--GGCACGTAGGAACTCTATACCTCGCTGATGGTAAATAGGCTTCGGTTTCTTA--------------------------TA-GAATCTT
dp4 chrXR_group6:8706030-8706135 - LASTZ 0 --------------------------------GTTGATGGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACT-TT-GGATA-CGGCACGTAGGAACTCTATACCTCGCTGATGGTAAAGAGGCTTCGCTTT------------------------------------AGAT
droPer1 super_27:1146105-1146210 - LASTZ 0 --------------------------------GTTGATGGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACT-TT-GGATA-CGGCACGTAGGAACTCTATACCTCGCTGATGGTAAAGAGGCTTCGCTTT------------------------------------AGAT
droWil1 scaffold_181024:290397-290516 + LASTZ 0 TT-----------------------------TGTTGTTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGATTATTGCAATT--TGCACGTAGGAACTCTATACCTCGCTTATGGTAAATAGGCTTCGTTTATTTC----------------------A--A---ATTTGTT
droVir3 scaffold_13049:16093198-16093326 - LASTZ 0 TC--------TTTT--C-----------T--TGCTGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGATTATTGCAATT--TGCACGTAGGAACTCTATACCTCGCTGATGGTATATAGACTTCGGTATCTTT---------------------AA--ATCTTAAGATA
droMoj3 scaffold_6680:20678203-20678343 + LASTZ 0 TT--------TATT--CGTT--------T--TGTTGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACTATTGCAATT--TGCACGTAGGAACTCTATACCTCGCTGATGGTATATAGACTTCGTTATCTTTACAATCTCTG-----------TT--TCA-CTTGGTC
droGri2 scaffold_15110:12572308-12572448 + LASTZ 0 AC--------TTAA--TTGTTTTGATCAA--TGCTGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGATTATTGCAATT--TGCACGTAGGAACTCTATACCTCGCTAATGGTAAATAGACTTCGTTATCTTT---------------------CGGTTTT-GAGATTT

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 TTTTAACTATATGCTGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGACTGAATATAGGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTTT--TTCGTTCGATGGTC
droSim1 TTTCAACTATATGATGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGACTGAA--TAGGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTT---TTCGGTCGATGGTC
droSec1 TTTCAACTATATGATGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGACTGAA--TAGGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTT---TTCGTTCG-TGGTC
droYak2 TTTTAACTTTATGATGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGACTGAA--TAGGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTTTTGTTCTTTCTGTTGTC

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

Small RNA-seq read pileup

Species Alignment
droSim1
TT--------TCAA--CTAT--------A--TGATGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATA--GGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTTTTC----------------------G--GTC-GATGGTCSizeHit CountTotal NormV044
Embryo
M023
Head
M024
Male body
M025
Embryo
...............................................................................................TAGGAACTTCATACCTCG................................................................18120.170.000.170.000.00
......................................................................................................................................GTTTTTC----------------------G--GTC-GATCGTC1820.500.000.500.000.00
...............................................................................................TAAGAACTCTATACCTAG................................................................1880.130.000.130.000.00
droSec1
TT--------TCAA--CTAT--------A--TGATGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATA--GGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTTTTC----------------------G--TTC--GTGGTCSizeHit CountTotal NormV113
Male body
V114
Embryo
V115
Head
.....................................................AGCGCGGTATAAAGTTCCTACG......................................................................................................22160.310.190.000.13
.....................................................................................................CTCTATACCTCGCTGATGGCA.......................................................2180.130.000.000.13
.....................................................AGCGCGGTATTGAGTTCCTACG......................................................................................................22120.080.000.000.08
.....................................................AGCGCGGTATCAAGTTCCTACG......................................................................................................2280.130.000.000.13
.....................................................AGCGAGGTATTAAGTTCCT.........................................................................................................19120.080.000.000.08
.....................................................AGCAAGGTATAAAGTTCCTACG......................................................................................................22160.060.000.000.06
.....................................................AGCACGGTATAGAGTTCCTACG......................................................................................................22160.060.000.000.06
droYak2
TT--------TTAA--CTTT--------A--TGATGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATA--GGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTTTTG-------------------TTCT--TTC-TGTTGTCSizeHit CountTotal NormV046
Embryo
V052
Head
V120
Male body
V058
Head
M026
Head
.....................................................AGCGCGGTATAAAGTTCCTACG......................................................................................................22160.250.190.000.000.060.00
...............................................................................................TAGGAACTTCATACCTCGCT..............................................................20120.080.000.000.000.080.00
...............................................................................................TAGGAACTTCATACCTCG................................................................18120.080.000.000.000.080.00
...............................................................................................TAGGAACTTTATACCGCGCT..............................................................20160.060.000.000.000.060.00
.....................................................AGCGCGGTATTGAGTTCCTACG......................................................................................................22120.080.080.000.000.000.00
.....................................................AGCGCGGTATTAAGTTCCTACG......................................................................................................2280.130.000.000.000.130.00
.....................................................AGCAAGGTATAAAGTTCCTACG......................................................................................................22160.060.000.000.000.060.00
.....................................................AGCACGGTATAGAGTTCCTACG......................................................................................................22160.060.000.000.000.060.00
...............................................................................................TAGGAACTGTATACCGCGCT..............................................................2080.130.000.130.000.000.00
droEre2
TT--------TCAA--CTTT--------A--TGATGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATA--GGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTTTTT----------------------A--TTC-GGTTGTCSizeHit CountTotal NormV040
Embryo
V060
Head
V047
Embryo
V108
Head
...............................................................................................TAGGAACTTTATACCTCGCTCA............................................................2280.130.000.130.000.00
.....................................................AGCACGGTATAAAGTTCCTCCG......................................................................................................2280.130.130.000.000.00
.....................................................AGCGCGGTATCAAGTTCCTACG......................................................................................................2280.130.130.000.000.00
...............................................................................................TAGGAACTTCATACCTCGCTG.............................................................2180.130.000.130.000.00
droEug1
TT--------TAAA--CTTT--------A--CGATGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----AAATA--GAAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGGCTTCGGTATTTTC------TTTGTTTTTCTTGGCCC--TTC-GATGGACSizeHit CountTotal Norm
droBia1
TT--------TCAA--GTTT--------T--GGATGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATA--GGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGACTTCGGTTTTTTG---------------------TA--TTC-GACGGTCSizeHit CountTotal Norm
droTak1
TT--------TCAA--CTTT--------A--TGATGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACG----GAATA--GGAACGTAGGAACTCTATACCTCGCTGATGGTAAAAAGGCTTCGGTTATTTC------T-------------TCG--TTC-GATGGTCSizeHit CountTotal Norm
droEle1
TG--------TTAATTCTTT--------A--TGTTGGTCGACAGTCCTCCAAGAGCACGGTATTAAGTTCCTACGATT----AAATG--GGAACGTAGGAACTCTATACCTCGCTGATGGTATAAAGACTTCGGTTTTCTG---------------------CT--TTG-AATGGTCSizeHit CountTotal Norm
droRho1
TG--------TAAA--CTTT--------A--TGTTGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATA--GGAACGTAGGAACTCTATACCTCGCTGATGGTAAA-AGACTTCGGATTTTTG---------------------GT--TTC-AATGGTCSizeHit CountTotal Norm
droFic1
TT--------TAGT--TTTT--------ATTTGATGGACGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACT----GAATT--TGAACGTAGGAACTCTATACCTCGCTTATGGTAAAAAGGCTTCGGATTTGGG---------------------GT--TTCTGATGGTCSizeHit CountTotal Norm
droKik1
CTCTATAATAAACA--CTTT--------A--TGTTGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACTATTTG---A--GGAACGTAGGAACTCTATACCTCGCTGATGGCAAAAAGGCTTCGTTTGTTTT-----------------TGGCTT--TTC-AGTCTTCSizeHit CountTotal Norm
droAna3
CT--------T-------TC--------A--TGATGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACTAGTTTAAGT--GGCACGTAGGAACTCTATACCTCGCTGATGGCAAATAGGCTTCGGTTTCTTC----------------------A--------GAAATSizeHit CountTotal NormV105
Male body
V106
Head
V055
Unknown
V039
Embryo
.....................................................AGCGCGGTATAAAGTTCCTACG......................................................................................................22160.250.000.000.250.00
..........................................................................................................................ATTAGGCTTCGGTTTGTT.....................................1821.000.000.001.000.00
.....................................................AGCACGGTATAGAGTTCCTACG......................................................................................................22160.130.000.060.060.00
.....................................................AGCGCGGTATTGAGTTCCTACG......................................................................................................22120.170.080.000.080.00
...............................................................................................TAGGAACTCCATACCTCGCTC.............................................................2180.130.000.000.130.00
..........................................................................................................................ATTAGACTTCGGTTTCTT.....................................1820.500.000.000.000.50
.....................................................AGCACGGTATACAGTTCCTACG......................................................................................................2280.130.000.000.130.00
.....................................................AGCATGGTATAAAGTTCCTACG......................................................................................................2280.130.000.000.130.00
.....................................................AGCACGGTATAAAGTTCCTACG......................................................................................................2280.130.000.000.130.00
...............................................................................................TAGGAACTTAATACCTCGCT..............................................................20120.080.000.000.080.00
...............................................................................................TAGGAACTTTATACCGCGCT..............................................................20160.060.000.000.060.00
droBip1
GA--------AAGA--CTTT--------C-TCGATGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACTAGTTAAAGT--GGCACGTAGGAACTCTATACCTCGCTGATGGTAAATAGGCTTCGGTTTCTTA--------------------------TA-GAATCTTSizeHit CountTotal Norm
dp4
--------------------------------GTTGATGGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACT-TT-GGATA-CGGCACGTAGGAACTCTATACCTCGCTGATGGTAAAGAGGCTTCGCTTT------------------------------------AGATSizeHit CountTotal NormV112
Male body
V043
Embryo
V051
Head
M022
Male body
...............................................................................................TTGAAACTCTATACCTCGCTGA............................................................2280.130.000.130.000.00
droPer1
--------------------------------GTTGATGGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACT-TT-GGATA-CGGCACGTAGGAACTCTATACCTCGCTGATGGTAAAGAGGCTTCGCTTT------------------------------------AGATSizeHit CountTotal NormV111
Male body
V042
Embryo
V050
Head
V057
Head
M021
Embryo
.....................................................AGCACGGTATAGAGTTCCTACG......................................................................................................22160.190.000.000.000.190.00
.....................................................AGCGCGGTATTGAGTTCCTACG......................................................................................................22120.170.000.000.000.170.00
.....................................................AGCGCGGTATAAAGTTCCTACG......................................................................................................22160.130.000.000.000.130.00
.....................................................AGCACGGTATAGAGTTCCT.........................................................................................................19160.060.000.000.000.060.00
...............................................................................................TAGGAACTATATACCGCGCTG.............................................................2180.130.000.000.000.130.00
...............................................................................................TAGGAACTCCATACCTCGCTC.............................................................2180.130.000.000.000.130.00
...............................................................................................AAGGAGCTCTATACCTCGCTG.............................................................2180.130.000.000.000.130.00
...............................................................................................TAGGAACTTCATACCTCGCT..............................................................20120.080.000.000.000.080.00
.....................................................AGCACGGTATAAAGTTCCT.........................................................................................................1980.130.000.000.000.130.00
.....................................................AGCACGGTxTTAACTTCCTA........................................................................................................1940.250.000.000.000.250.00
...............................................................................................TAGGAACTTTATACCTCGCTC.............................................................2180.130.000.000.000.130.00
.....................................................AGCACGGTATCGAGTTCCTACG......................................................................................................2280.130.000.000.000.130.00
...............................................................................................TAGGAACTTTAGACCTCGCTGA............................................................2280.130.000.000.130.000.00
droWil1
TT-----------------------------TGTTGTTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGATTATTGCAATT--TGCACGTAGGAACTCTATACCTCGCTTATGGTAAATAGGCTTCGTTTATTTC----------------------A--A---ATTTGTTSizeHit CountTotal NormV117
Male body
V118
Embryo
V119
Head
M020
Head
.....................................................AGCGCGGTATAAAGTTCCTACG......................................................................................................22120.830.000.080.750.00
.....................................................AGCACGGTATAGAGTTCCTACG......................................................................................................22120.420.000.000.420.00
.....................................................AGCAAGGTATAAAGTTCCTACG......................................................................................................22120.170.000.000.170.00
.....................................................AGCGCGGTATTGAGTTCCTACG......................................................................................................22100.100.000.000.100.00
...............................................................................................TAGGAACTTAATACCTCG................................................................18100.100.000.000.000.10
droVir3
TC--------TTTT--C-----------T--TGCTGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGATTATTGCAATT--TGCACGTAGGAACTCTATACCTCGCTGATGGTATATAGACTTCGGTATCTTT---------------------AA--ATCTTAAGATASizeHit CountTotal NormV116
Male body
V053
Unknown
V047
Embryo
M027
Male body
M028
Head
.....................................................AGCGCGGTATTGAGTTCCTACG......................................................................................................22120.170.000.170.000.000.00
.....................................................AGCACGGTATAGAGTTCCTACG......................................................................................................22160.060.000.060.000.000.00
.....................................................AGCGCGGTATAAAGTTCCTACG......................................................................................................22160.060.000.060.000.000.00
...............................................................................................TAGGAACTTTATACCTxGCTGAT...........................................................2280.130.000.130.000.000.00
droMoj3
TT--------TATT--CGTT--------T--TGTTGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGACTATTGCAATT--TGCACGTAGGAACTCTATACCTCGCTGATGGTATATAGACTTCGTTATCTTTACAATCTCTG-----------TT--TCA-CTTGGTCSizeHit CountTotal NormV110
Male body
V056
Unknown
V049
Head
V041
Embryo
.....................................................AGCGCGGTATAAAGTTCCTACG......................................................................................................22160.310.000.190.130.00
...............................................................................................TAGGAACTTCATACCTCGCT..............................................................20120.170.000.170.000.00
...............................................................................................TAGGAACTCCATACCGCGCTG.............................................................2180.130.000.000.130.00
.....................................................AGCAAGGTATAAAGTTCCT.........................................................................................................19160.060.000.060.000.00
..................TT--------T--TGTTGGACGCCAGCC...................................................................................................................................1840.250.250.000.000.00
.....................................................AGCACGGTTTAAAGTTCCT.........................................................................................................1980.130.000.000.000.13
droGri2
AC--------TTAA--TTGTTTTGATCAA--TGCTGGTCGACAGCCCTCCAAGAGCACGGTATTAAGTTCCTACGATTATTGCAATT--TGCACGTAGGAACTCTATACCTCGCTAATGGTAAATAGACTTCGTTATCTTT---------------------CGGTTTT-GAGATTTSizeHit CountTotal NormV109
Male body

Generated: 04/20/2012 at 05:14:15 EDT