Return to list

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

View all possible LASTZ hits (20-flies)

dme-mir-305-as

chr2L:7425972-7426044 -

Novel

[View on UCSC Genome Browser {Cornell Mirror}]

intergenic


Legend:

   mature strand at position 7425979-7426001
   star strand at position 7426016-7426036
   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 chr2L:7425957-7426059 - 1 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGA-GACACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droSim1 chr2L:7217527-7217627 - LASTZ 1 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAGCATGTGCCGGGTTACGA-GACACACCAGAGCA-CTGATGAAGTACAATAGACATGGGAGACAGT-G--------ATAC------------------------------------------
droSec1 super_3:2941418-2941520 - LASTZ 1 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGA-GACACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droYak2 chr2L:5473614-5473716 + LASTZ 1 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGG-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droEre2 scaffold_4929:16344683-16344785 - LASTZ 1 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGG-ACCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droEug1 scf7180000409463:614863-614965 - LASTZ 4 GT----------------ATGCAAATCACCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACAG-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droBia1 scf7180000302422:9125681-9125783 - LASTZ 4 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGG-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTCG--------ATAC------------------------------------------
droTak1 scf7180000415604:142269-142370 + LASTZ 4 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACG--GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droEle1 scf7180000491028:1624701-1624808 + LASTZ 4 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACTT-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATACACAAC-------------------------------------
droRho1 scf7180000766421:7519-7621 + LASTZ 4 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACTT-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droFic1 scf7180000454115:155833-155935 - LASTZ 4 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTGCGG-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droKik1 scf7180000302384:54674-54769 + LASTZ 4 AT----------------TTGGAAATCACCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGA-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAG---------------ATAC------------------------------------------
droAna3 scaffold_12916:15334141-15334251 - LASTZ 1 GT----------------TTGACATTCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGG-TGGACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACATATTTCTCAAAAATAC------------------------------------------
droBip1 scf7180000396572:1689256-1689352 + LASTZ 4 GA----------------TTGACATTCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGA-TAGACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACATT--------------------------------------------------------
dp4 chr4_group1:4175881-4175974 - LASTZ 4 GC----------------ACGAACTTCGCCACATATTGAGTGTACTTCAACATGTGCCGGGTTACAG-GACACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGAC-----------------------------------------------------------
droPer1 super_5:1086315-1086408 + LASTZ 1 GC----------------ACGAACTTCGCCACATATTGAGTGTACTTCAACATGTGCCGGGTTACAG-GACACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGAC-----------------------------------------------------------
droWil1 scaffold_180772:5257378-5257513 + LASTZ 1 ------------------ATGAAGATCGCCACATATTGAGTGTACTTCAACATGTGCCGGGTTACAAACGGACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAA--------------C--AACAATAAATTAAAAATGCCAATAAAAAAAAAAGAGAATC
droVir3 scaffold_12963:15419827-15419927 + LASTZ 1 AA----------------ATGCAAATCAACGCATATTGAGTGTACTTCAACATGTGCCGGGTTACAA-GTCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACA--TG--------AAAC------------------------------------------
droMoj3 scaffold_6500:18399247-18399346 + LASTZ 1 GC--------------AACTGCAAATCACCGCATATTGAGTGTACTTCAACATGTGCCGGGTTTCCA-CTGACACCAGAGCACCTGATGAAGTACAATAGACATGGGGGA-----G--------ACAC------------------------------------------
droGri2 scaffold_15126:5008541-5008648 - LASTZ 1 TTAATTAAAACGCAGCAAATGCCAATCACCGCATATTGAGTGTACTTCAACATGTGCCGGGTTACGA-GTGACACCAGAGCACCTGATGAAGTACAATAGACATGGGAG-------------------------------------------------------------

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 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGA-GACACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droSim1 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAGCATGTGCCGGGTTACGA-GACACACCAGAGCA-CTGATGAAGTACAATAGACATGGGAGACAGT-G--------ATAC------------------------------------------
droSec1 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGA-GACACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droYak2 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGG-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droEre2 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGG-ACCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droEug1 GT----------------ATGCAAATCACCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACAG-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droBia1 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGG-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTCG--------ATAC------------------------------------------
droTak1 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACG--GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droEle1 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACTT-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATACACAAC-------------------------------------
droRho1 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACTT-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droFic1 GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTGCGG-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------
droKik1 AT----------------TTGGAAATCACCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGA-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAG---------------ATAC------------------------------------------
droAna3 GT----------------TTGACATTCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGG-TGGACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACATATTTCTCAAAAATAC------------------------------------------
droBip1 GA----------------TTGACATTCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGA-TAGACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACATT--------------------------------------------------------
dp4 GC----------------ACGAACTTCGCCACATATTGAGTGTACTTCAACATGTGCCGGGTTACAG-GACACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGAC-----------------------------------------------------------
droPer1 GC----------------ACGAACTTCGCCACATATTGAGTGTACTTCAACATGTGCCGGGTTACAG-GACACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGAC-----------------------------------------------------------
droWil1 ------------------ATGAAGATCGCCACATATTGAGTGTACTTCAACATGTGCCGGGTTACAAACGGACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAA--------------C--AACAATAAATTAAAAATGCCAATAAAAAAAAAAGAGAATC
droVir3 AA----------------ATGCAAATCAACGCATATTGAGTGTACTTCAACATGTGCCGGGTTACAA-GTCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACA--TG--------AAAC------------------------------------------
droMoj3 GC--------------AACTGCAAATCACCGCATATTGAGTGTACTTCAACATGTGCCGGGTTTCCA-CTGACACCAGAGCACCTGATGAAGTACAATAGACATGGGGGA-----G--------ACAC------------------------------------------
droGri2 TTAATTAAAACGCAGCAAATGCCAATCACCGCATATTGAGTGTACTTCAACATGTGCCGGGTTACGA-GTGACACCAGAGCACCTGATGAAGTACAATAGACATGGGAG-------------------------------------------------------------

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
GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAGCATGTGCCGGGTTACGA-GACACACCAGAGCA-CTGATGAAGTACAATAGACATGGGAGACAGT-G--------ATAC------------------------------------------SizeHit CountTotal NormV044
Embryo
M023
Head
M024
Male body
M025
Embryo
......................................ATTGTACTTCATCATGTGC.................................................................................................................1921.500.001.500.000.00
......................................................................................ATGAAGxACAATAGGCATG.................................................................1820.500.000.500.000.00
......................................ATTGTACTTCATCATGTGCC................................................................................................................2020.500.000.500.000.00
droSec1
GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGA-GACACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------SizeHit CountTotal NormV113
Male body
V114
Embryo
V115
Head
...........................................................................CAGAGCACCTGATGAAGTACAAT........................................................................23210.006.000.004.00
...........................................................................CAGAGCACCTGATGAAGTACA..........................................................................2123.002.000.001.00
...........................................................................CAGAGCACCTGATGAAGTACAA.........................................................................2221.501.000.500.00
............................................................................AGAGCACCTGATGAAGTACAAT........................................................................2221.000.000.500.50
............................................................................AGAGCACCTGTTGAAGTACAAT........................................................................2220.500.000.000.50
...........................................................................CAGAGCACCTGATGAAGTTCAAT........................................................................2320.500.500.000.00
...........................................................................CAGAGCACCTGATGAAGTACAAA........................................................................2320.500.000.000.50
...........................................................................CAGAGCACCTGATGAAGTA............................................................................1920.500.000.000.50
droYak2
GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGG-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------SizeHit CountTotal NormV046
Embryo
V052
Head
V120
Male body
V058
Head
M026
Head
...........................................................................CAGAGCACCTGATGAAGTACA..........................................................................2123.000.000.001.501.000.50
...........................................................................CAGAGCACCTGATGAAGTACAAT........................................................................2322.000.000.001.500.500.00
...........................................................................CAGAGCACCTGATGAAGTACAA.........................................................................2221.000.000.001.000.000.00
...........................................................................CAGAGCACCTGATGAAGTAC...........................................................................2020.500.000.000.000.500.00
droEre2
GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGG-ACCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------SizeHit CountTotal NormV040
Embryo
V060
Head
V047
Embryo
V108
Head
...........................................................................CAGAGCACCTGATGAAGTACAAT........................................................................2322.000.000.000.002.00
...........................................................................CAGAGCACCTGATGAAGTACA..........................................................................2121.000.000.000.001.00
......................................ATTGTACTTCATCATGTGC.................................................................................................................1940.250.000.000.000.25
...........................................................................CATAGCACCTGATGAAGTAC...........................................................................2020.500.000.000.000.50
...........................................................................CAGAGCACCTGATGAAGTAC...........................................................................2020.500.000.000.000.50
...........................................................................CAGAGCACCTGATGAGGTA............................................................................1920.500.000.000.000.50
......................................ATTGTACTTCATCATGTG..................................................................................................................1840.250.000.000.000.25
droEug1
GT----------------ATGCAAATCACCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACAG-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------SizeHit CountTotal Norm
droBia1
GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGG-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTCG--------ATAC------------------------------------------SizeHit CountTotal Norm
droTak1
GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACG--GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------SizeHit CountTotal Norm
droEle1
GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACTT-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATACACAAC-------------------------------------SizeHit CountTotal Norm
droRho1
GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACTT-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------SizeHit CountTotal Norm
droFic1
GT----------------ATGCAAATCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTGCGG-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAGTTG--------ATAC------------------------------------------SizeHit CountTotal Norm
droKik1
AT----------------TTGGAAATCACCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGA-GCCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAG---------------ATAC------------------------------------------SizeHit CountTotal Norm
droAna3
GT----------------TTGACATTCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGG-TGGACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACATATTTCTCAAAAATAC------------------------------------------SizeHit CountTotal NormV105
Male body
V106
Head
V055
Unknown
V039
Embryo
...........................................................................CAGAGCACCTGATGAAGTACA..........................................................................2120.500.000.500.000.00
droBip1
GA----------------TTGACATTCGCCTCATATTGAGTGTACTTCAACATGTGCCGGGTTACGA-TAGACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACATT--------------------------------------------------------SizeHit CountTotal Norm
dp4
GC----------------ACGAACTTCGCCACATATTGAGTGTACTTCAACATGTGCCGGGTTACAG-GACACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGAC-----------------------------------------------------------SizeHit CountTotal NormV112
Male body
V043
Embryo
V051
Head
M022
Male body
.......................................GTGTACTTCAACCTGGGC.................................................................................................................1820.500.000.000.000.50
droPer1
GC----------------ACGAACTTCGCCACATATTGAGTGTACTTCAACATGTGCCGGGTTACAG-GACACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGAC-----------------------------------------------------------SizeHit CountTotal NormV111
Male body
V042
Embryo
V050
Head
V057
Head
M021
Embryo
...........................................................................CAGAGCACCTGATGAAGTACA..........................................................................2120.500.500.000.000.000.00
...........................................................................CAGAGCACCTGATGAAGTACAAT........................................................................2320.500.500.000.000.000.00
droWil1
------------------ATGAAGATCGCCACATATTGAGTGTACTTCAACATGTGCCGGGTTACAAACGGACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACAA--------------C--AACAATAAATTAAAAATGCCAATAAAAAAAAAAGAGAATCSizeHit CountTotal NormV117
Male body
V118
Embryo
V119
Head
M020
Head
...........................................................................CAGAGCACCTGATGAAGTACAAT........................................................................2321.000.000.001.000.00
....................................................................................................................................................TCAATAAAAGAAAAAGAGA...1940.250.000.250.000.00
....................GAAGATCxCGACATATTGAGT.................................................................................................................................2020.500.000.500.000.00
......................................................................................................................................................AAAAAAAAAAAAAGAGAA..18320.030.030.000.000.00
................................................................................................................................................AATGTTAATAAAAAAAAA........18150.070.000.070.000.00
..............................................................................................................................................AAAATGACATTAAAAAAAA.........1930.330.000.330.000.00
droVir3
AA----------------ATGCAAATCAACGCATATTGAGTGTACTTCAACATGTGCCGGGTTACAA-GTCACACCAGAGCACCTGATGAAGTACAATAGACATGGGAGACA--TG--------AAAC------------------------------------------SizeHit CountTotal NormV116
Male body
V053
Unknown
V047
Embryo
M027
Male body
M028
Head
........................................TGTACTTCAACATGTGCCGGGT............................................................................................................2220.500.000.000.000.500.00
.................................................................AA-GTCACACCAGAGCACCTGA...................................................................................2120.500.500.000.000.000.00
......................................AGTGTACTTCATCAGGTG..................................................................................................................1840.250.000.250.000.000.00
droMoj3
GC--------------AACTGCAAATCACCGCATATTGAGTGTACTTCAACATGTGCCGGGTTTCCA-CTGACACCAGAGCACCTGATGAAGTACAATAGACATGGGGGA-----G--------ACAC------------------------------------------SizeHit CountTotal NormV110
Male body
V056
Unknown
V049
Head
V041
Embryo
...........................................................................CAGAGCACCTGATGAAGTACA..........................................................................2120.500.500.000.000.00
droGri2
TTAATTAAAACGCAGCAAATGCCAATCACCGCATATTGAGTGTACTTCAACATGTGCCGGGTTACGA-GTGACACCAGAGCACCTGATGAAGTACAATAGACATGGGAG-------------------------------------------------------------SizeHit CountTotal NormV109
Male body
............................................................................AGAGCACCTGATGAAGTACAAT........................................................................2220.500.50

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