Return to list

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

View all possible LASTZ hits (20-flies)

dme-mir-2490

chr2LHet:27336-27410 +

Known_mirbase

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

antisense_to_intron


Legend:

   mature strand at position 27344-27366
   star strand at position 27382-27402
   mutations
   read mismatches

Notes: The reads that map to sechelia and yakuba are not uniquely mapped reads.


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

Visually pleasing best hits from LASTZ and MAF


Species Coordinate Source Confidence Alignment
dm3 chr2LHet:27321-27425 + 1 ACAATCGTACGTTCTGAGCGAGTTGAAGCGATAAAGCAGGTTGCAAACTGCTCTCTCGC--TTTGCT----------------CTTACTTTATCGCTTAAACTCGCAACTCGCTT--------------------------------------------------------------------------CGCT-----CTCA
droSim1 chr4:281112-281214 + LASTZ 4 AAAATCGTACGCTCTGAGCGAGTTGAAGCAATAAAGCAGGTTGCAAACTGCTCTTTCGC--TTTGCT----------------CTCGCTATATCGCTTGAGCTCGCAACTCGCTT--------------------------------------------------------------------------CGCT-------CA
droSec1 super_30:85009-85111 + LASTZ 3 AAAATCGTACGCTCTGAGCGAGTTGAAGCAATAAAGCAGGTTGCAAACTGCTCTTTCGC--TTTGCT----------------CTCGCTATATCGCTTGAGCTCGCAACTCGCTT--------------------------------------------------------------------------CGCT-------CA
droYak2 chrU:9061077-9061181 + LASTZ 3 ACAATCGTACGCTCTGAGCGAGTTGAAGCGATAAAGCAGGTTGCAAACTGCTCTCTCGC--TTTGCT----------------CTCGCTTTATCGCTTGAGCTCGCTACTCGCCT--------------------------------------------------------------------------CGCT-----CTCA
droEre2 scaffold_4784:23192819-23192904 - LASTZ 4 ACTATCAGACGCTCTGAGCGAGTTGCAGCGATAATGCAGGTTGCAAACTGTTCTCT-------------------------------------CGCTTGAGCTCGCAACTCGTTT--------------------------------------------------------------------------CGCT-----CTCA
droEug1 Unknown 0 ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
droBia1 scf7180000301391:7073-7197 - LASTZ 0 ACTTTTGTAGCTTCTGCAGGCATTTGG----CTGTGTAGATTGTATGGTCTCCGAT-------------------------------------AAGTAGTAGTCGAAAGTCCCT--------------CACGATGACAACGTTTGTAACG--------TACCAGGGGAGCTCGAATGCTCGCCG-----A--G-------CT
droTak1 Unknown 0 ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
droEle1 scf7180000490727:38228-38321 - LASTZ 0 AGAA-CTTACGTTTTGATTAAGTTCAAGCGATAAAGCAAGTTCAAAACTGATCTCTTTC--TTTGCC----------------CTCAATTTTTCGCTTCGTT----------TTT--------------------------------------------------------------------------CGCA-----GCCA
droRho1 Unknown 0 ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
droFic1 scf7180000453955:804345-804445 - LASTZ 0 AAAATTTTTCGATCGG-----G--------------CTCCCCCCAAAATTTTTTCTGCG--TCCGC-----------------------------------------CCCTGACTG------------------------GCTTGTGACC--------TACCAGGGAAGCCCGGATGCCCGTCG-----A--GCAGCACGTG
droKik1 Unknown 0 ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
droAna3 Unknown 0 ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
droBip1 scf7180000396710:138523-138625 + LASTZ 0 ACAATCGAACGCTCTGAGCGAGTTGAAGCGATGAAGCTGGTTGCAAACTGCTCTCTC----TTTCCT----------------CCCGCTTTATCGCTTGAACTCGCAACTCGCTT--------------------------------------------------------------------------CGGT-----TTTA
dp4 chr4_group4:3101286-3101457 + LASTZ 0 AGCAT-GAGCG---TCCACGAGATCAAGCGATAAGACATTTGCAAAGCTACTCTCTTGC--TTCGCTTGCACTCGTTGCTGTCTTCGCTATAA-------------------CTTGTTTCCCGCTCTTCACTGTAAATCGTTTTGCACTCTTCGCTCAAACTAGCCGCGCTGCTCTCTTCGCTTATCTTTGCT-----CTCA
droPer1 Unknown 0 ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
droWil1 scaffold_181143:687820-687924 + LASTZ 0 ACACACGTATAACAGCAGCTGGATGAAGCGATAAAGCGA-TTCAGAATCGCTTTCTCGCAACTCGCT-----------------TCAACTTGTCGCTTTGGCTCTCTTCTCGCTT--------------------------------------------------------------------------CGCT-----TGCA
droVir3 Unknown 0 ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
droMoj3 scaffold_6498:432956-433055 + LASTZ 0 AACAC---ACTTTAT-CATGCGATTAAGCGAGAAAGCGA-TTCAGGAATGCTTTCTCGC--AACTCG----------------CTTAACTTCTCGCTTTGGGTCTCTTCTCGCTT--------------------------------------------------------------------------CGCT-----TGCA
droGri2 Unknown 0 ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

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 ACAATCGTACGTTCTGAGCGAGTTGAAGCGATAAAGCAGGTTGCAAACTGCTCTCTCGCTTTGCTCTTACTTTATCGCTTAAACTCGCAACTCGCTTCGCTCTCA
droSim1 AAAATCGTACGCTCTGAGCGAGTTGAAGCAATAAAGCAGGTTGCAAACTGCTCTTTCGCTTTGCTCTCGCTATATCGCTTGAGCTCGCAACTCGCTTCGCT--CA
droSec1 AAAATCGTACGCTCTGAGCGAGTTGAAGCAATAAAGCAGGTTGCAAACTGCTCTTTCGCTTTGCTCTCGCTATATCGCTTGAGCTCGCAACTCGCTTCGCT--CA
droYak2 ACAATCGTACGCTCTGAGCGAGTTGAAGCGATAAAGCAGGTTGCAAACTGCTCTCTCGCTTTGCTCTCGCTTTATCGCTTGAGCTCGCTACTCGCCTCGCTCTCA
droEre2 ACTATCAGACGCTCTGAGCGAGTTGCAGCGATAATGCAGGTTGCAAACTGTTCTCT-------------------CGCTTGAGCTCGCAACTCGTTTCGCTCTCA

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

Small RNA-seq read pileup

Species Alignment
droSim1
AAAATCGTACGCTCTGAGCGAGTTGAAGCAATAAAGCAGGTTGCAAACTGCTCTTTCGC--TTTGCT----------------CTCGCTATATCGCTTGAGCTCGCAACTCGCTT--------------------------------------------------------------------------CGCT-------CASizeHit CountTotal NormV044
Embryo
M023
Head
M024
Male body
M025
Embryo
..AATCGTACGCTCTGAGCGA.....................................................................................................................................................................................19130.080.000.000.080.00
..........................................................................................TATCGCTTGAGCTCGCAACTC...........................................................................................21110.090.000.090.000.00
..AATCGTACGCTCTGAGCGG.....................................................................................................................................................................................19130.080.000.000.080.00
..AATCGTACGCTCTGAGCGAGTTAAAGC.............................................................................................................................................................................27100.100.000.000.100.00
droSec1
AAAATCGTACGCTCTGAGCGAGTTGAAGCAATAAAGCAGGTTGCAAACTGCTCTTTCGC--TTTGCT----------------CTCGCTATATCGCTTGAGCTCGCAACTCGCTT--------------------------------------------------------------------------CGCT-------CASizeHit CountTotal NormV113
Male body
V114
Embryo
V115
Head
......................TTGAAGCAATAAAGCAGGTTG...............................................................................................................................................................2121.000.000.001.00
.......................TGAAGCAATAAAGCAGGTTGCA.............................................................................................................................................................2220.500.000.000.50
.........................AAGCGATAAAGCAGGTTGCAAACTG........................................................................................................................................................2550.200.000.000.20
droYak2
ACAATCGTACGCTCTGAGCGAGTTGAAGCGATAAAGCAGGTTGCAAACTGCTCTCTCGC--TTTGCT----------------CTCGCTTTATCGCTTGAGCTCGCTACTCGCCT--------------------------------------------------------------------------CGCT-----CTCASizeHit CountTotal NormV046
Embryo
V052
Head
V120
Male body
V058
Head
M026
Head
...............................................................TGCT----------------CTCGCTTAATCGCTTGAGA....................................................................................................2370.140.000.000.000.140.00
.....................GTTGAAGCGATAAAGCAGGTTGCA.............................................................................................................................................................24100.100.000.000.100.000.00
.......................TGAAGCGATAAAGCAGGTTGCA.............................................................................................................................................................22100.100.100.000.000.000.00
.....................GTTGAAGCGATAAAGCAGGTTGC..............................................................................................................................................................23100.100.000.000.100.000.00
.......................TGAAGCGATAAAGCAGGTTGCC.............................................................................................................................................................22110.090.090.000.000.000.00
droEre2
ACTATCAGACGCTCTGAGCGAGTTGCAGCGATAATGCAGGTTGCAAACTGTTCTCT-------------------------------------CGCTTGAGCTCGCAACTCGTTT--------------------------------------------------------------------------CGCT-----CTCASizeHit CountTotal NormV040
Embryo
V060
Head
V047
Embryo
V108
Head
droEug1
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------SizeHit CountTotal Norm
droBia1
ACTTTTGTAGCTTCTGCAGGCATTTGG----CTGTGTAGATTGTATGGTCTCCGAT-------------------------------------AAGTAGTAGTCGAAAGTCCCT--------------CACGATGACAACGTTTGTAACG--------TACCAGGGGAGCTCGAATGCTCGCCG-----A--G-------CTSizeHit CountTotal Norm
droTak1
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------SizeHit CountTotal Norm
droEle1
AGAA-CTTACGTTTTGATTAAGTTCAAGCGATAAAGCAAGTTCAAAACTGATCTCTTTC--TTTGCC----------------CTCAATTTTTCGCTTCGTT----------TTT--------------------------------------------------------------------------CGCA-----GCCASizeHit CountTotal Norm
droRho1
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------SizeHit CountTotal Norm
droFic1
AAAATTTTTCGATCGG-----G--------------CTCCCCCCAAAATTTTTTCTGCG--TCCGC-----------------------------------------CCCTGACTG------------------------GCTTGTGACC--------TACCAGGGAAGCCCGGATGCCCGTCG-----A--GCAGCACGTGSizeHit CountTotal Norm
droKik1
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------SizeHit CountTotal Norm
droAna3
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------SizeHit CountTotal NormV105
Male body
V106
Head
V055
Unknown
V039
Embryo
droBip1
ACAATCGAACGCTCTGAGCGAGTTGAAGCGATGAAGCTGGTTGCAAACTGCTCTCTC----TTTCCT----------------CCCGCTTTATCGCTTGAACTCGCAACTCGCTT--------------------------------------------------------------------------CGGT-----TTTASizeHit CountTotal Norm
dp4
AGCAT-GAGCG---TCCACGAGATCAAGCGATAAGACATTTGCAAAGCTACTCTCTTGC--TTCGCTTGCACTCGTTGCTGTCTTCGCTATAA-------------------CTTGTTTCCCGCTCTTCACTGTAAATCGTTTTGCACTCTTCGCTCAAACTAGCCGCGCTGCTCTCTTCGCTTATCTTTGCT-----CTCASizeHit CountTotal NormV112
Male body
V043
Embryo
V051
Head
M022
Male body
.................................................................CTTGCACTCGTTTCGGTC.......................................................................................................................1813.000.000.000.003.00
..................................................................TTGCACTCGTTTCGGTCT......................................................................................................................1811.000.000.000.001.00
....................................................................GCACTCGTTGCTGTTGTCGC..................................................................................................................2011.000.001.000.000.00
.................................................................CTTGCACTCGTTTCGGTCT......................................................................................................................1911.000.000.000.001.00
droPer1
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------SizeHit CountTotal NormV111
Male body
V042
Embryo
V050
Head
V057
Head
M021
Embryo
droWil1
ACACACGTATAACAGCAGCTGGATGAAGCGATAAAGCGA-TTCAGAATCGCTTTCTCGCAACTCGCT-----------------TCAACTTGTCGCTTTGGCTCTCTTCTCGCTT--------------------------------------------------------------------------CGCT-----TGCASizeHit CountTotal NormV117
Male body
V118
Embryo
V119
Head
M020
Head
........................................TTCAGAATCGCTTTCTCGCAACTCGCT-----------------......................................................................................................................27130.310.000.310.000.00
........................................TTCAGAATCGCTTTCTCGCAACTCGCT-----------------T.....................................................................................................................28130.150.000.150.000.00
..........................................................................................TGTCGCTTTGGCTCTCTTCTCGCTT--------------------------------------------------------------------------.............25160.130.000.130.000.00
..........................................................................................TGTCGCTTTGGCTCTCTTCTC...........................................................................................21160.060.000.060.000.00
..........AAACGCAGCTGGATGAAG..............................................................................................................................................................................1811.000.000.000.001.00
..........................................................................................TGTCGCTTTGGCTCTCTTCTCGCT........................................................................................24160.060.000.060.000.00
...............................................................CGCT-----------------TCAACTTGTCGCTTTGGCTCT.................................................................................................25150.070.000.070.000.00
...............................................TCGCTTTCTCGCAACTCGCT-----------------TCAAC.................................................................................................................25140.070.000.070.000.00
................................................CGCTTTCTCGCAACTCGCT-----------------TCAAC.................................................................................................................24140.070.000.070.000.00
................................................CGCTTTCTCGCAACTCACT-----------------TCAACT................................................................................................................25130.080.000.080.000.00
.........................................TCAGAATCGCTTTCTCGCAACTCACT-----------------......................................................................................................................26140.070.000.070.000.00
..........................................CAGAATCGCTTTCTCGCAACTCGCT-----------------T.....................................................................................................................26140.070.000.070.000.00
.........................................TCAGAATCGCTTTCTCGCAACTCGCT-----------------......................................................................................................................26130.080.000.080.000.00
...........................................................AACTCGCT-----------------TCAACTTGTCGCTTTTGCTC..................................................................................................28150.070.000.070.000.00
...............................................................CGCT-----------------TCAACTTGTCGCTTTGGCT...................................................................................................23150.070.000.070.000.00
.....................................................TCTCGCAACTCGCT-----------------TCAACTTGTCG...........................................................................................................25150.070.000.070.000.00
....................................................................................TCAACTTGTCGCTTTGGCTCTCTTCT............................................................................................26150.070.000.070.000.00
...........................GCGATAAAGCGA-TTCAGAATCG........................................................................................................................................................22130.080.000.080.000.00
.............................................................................................CGCTTTGGCTCTCTTCTCGCTT--------------------------------------------------------------------------CGC..........25160.060.000.060.000.00
.............................................................................................CGCTTTGTCTCTCTTCTCGCTT--------------------------------------------------------------------------CGCT-----....26160.060.000.060.000.00
droVir3
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------SizeHit CountTotal NormV116
Male body
V053
Unknown
V047
Embryo
M027
Male body
M028
Head
droMoj3
AACAC---ACTTTAT-CATGCGATTAAGCGAGAAAGCGA-TTCAGGAATGCTTTCTCGC--AACTCG----------------CTTAACTTCTCGCTTTGGGTCTCTTCTCGCTT--------------------------------------------------------------------------CGCT-----TGCASizeHit CountTotal NormV110
Male body
V056
Unknown
V049
Head
V041
Embryo
droGri2
----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------SizeHit CountTotal NormV109
Male body

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