Multiple sequence alignment - TraesCS4A01G488500

Loading Multiple Alignment...

Download Multiple Sequence Alignment: (Fasta) (Clustal)



BLAST Results


Shown here are the results of a BLAST search against the genome using the input sequence. The pipeline uses this to extract flanking regions for cloning. If the result of the multiple sequence alignment is not as expected, these results could be used for further investigation.

qseqid sseqid percentage.identical alignment.length no.mismatch no.gap.openings qstart qend sstart send evalue bitscore
0 TraesCS4A01G488500 chr4A 100.000 2852 0 0 1 2852 740064767 740067618 0.000000e+00 5267.0
1 TraesCS4A01G488500 chr4A 87.727 1100 102 12 759 1850 733931161 733932235 0.000000e+00 1253.0
2 TraesCS4A01G488500 chr4A 88.254 928 95 7 924 1850 739828888 739827974 0.000000e+00 1098.0
3 TraesCS4A01G488500 chr4A 84.567 1121 140 14 759 1865 742273777 742274878 0.000000e+00 1081.0
4 TraesCS4A01G488500 chr4A 90.039 773 36 18 1 764 739829844 739829104 0.000000e+00 963.0
5 TraesCS4A01G488500 chr4A 82.932 1037 140 22 985 2004 740872978 740873994 0.000000e+00 900.0
6 TraesCS4A01G488500 chr4A 83.867 781 69 19 1 761 742272991 742273734 0.000000e+00 691.0
7 TraesCS4A01G488500 chr4A 83.877 521 46 19 1 521 733930482 733930964 2.000000e-126 462.0
8 TraesCS4A01G488500 chr4A 89.773 88 6 3 582 667 733930963 733931049 3.000000e-20 110.0
9 TraesCS4A01G488500 chr7A 89.306 879 92 1 989 1865 2570634 2571512 0.000000e+00 1101.0
10 TraesCS4A01G488500 chr7A 84.375 1024 127 20 989 2004 2859953 2860951 0.000000e+00 974.0
11 TraesCS4A01G488500 chr7A 81.363 939 110 36 945 1847 3007832 3008741 0.000000e+00 704.0
12 TraesCS4A01G488500 chr7A 81.385 881 134 12 985 1856 3324663 3323804 0.000000e+00 691.0
13 TraesCS4A01G488500 chr7A 92.308 65 3 1 2 66 2569626 2569688 1.090000e-14 91.6
14 TraesCS4A01G488500 chr7A 89.062 64 4 2 4 67 2856708 2856768 3.050000e-10 76.8
15 TraesCS4A01G488500 chr7D 89.294 878 93 1 989 1865 3740582 3741459 0.000000e+00 1099.0
16 TraesCS4A01G488500 chr7D 80.470 1234 173 37 819 2032 3435229 3434044 0.000000e+00 881.0
17 TraesCS4A01G488500 chr7D 83.333 186 28 3 1584 1767 28884302 28884118 4.890000e-38 169.0
18 TraesCS4A01G488500 chr7D 89.394 66 4 2 2 67 3739662 3739724 2.360000e-11 80.5
19 TraesCS4A01G488500 chr5B 86.853 251 30 3 2337 2585 485664653 485664902 7.780000e-71 278.0
20 TraesCS4A01G488500 chr5B 89.552 67 7 0 2347 2413 485664897 485664831 5.060000e-13 86.1
21 TraesCS4A01G488500 chrUn 83.243 185 30 1 1584 1767 90469202 90469386 4.890000e-38 169.0
22 TraesCS4A01G488500 chr2D 82.353 136 23 1 2423 2557 615472362 615472497 1.800000e-22 117.0
23 TraesCS4A01G488500 chr6D 96.875 32 1 0 2368 2399 94779696 94779665 1.000000e-03 54.7


These results show the homologous regions that are presented in the multiple sequence alignment. HSPs from the BLAST search that are within 2000 bases of each other have been grouped together. Homologous regions that are purely upstream or downstream may be removed if they don't overlap both primers.

query scaffold start end length rev.comp avg.bitscore max.bitscore avg.percent.identical query.start query.end num_hsp groupid homo_length
0 TraesCS4A01G488500 chr4A 740064767 740067618 2851 False 5267.000000 5267 100.000000 1 2852 1 chr4A.!!$F1 2851
1 TraesCS4A01G488500 chr4A 739827974 739829844 1870 True 1030.500000 1098 89.146500 1 1850 2 chr4A.!!$R1 1849
2 TraesCS4A01G488500 chr4A 740872978 740873994 1016 False 900.000000 900 82.932000 985 2004 1 chr4A.!!$F2 1019
3 TraesCS4A01G488500 chr4A 742272991 742274878 1887 False 886.000000 1081 84.217000 1 1865 2 chr4A.!!$F4 1864
4 TraesCS4A01G488500 chr4A 733930482 733932235 1753 False 608.333333 1253 87.125667 1 1850 3 chr4A.!!$F3 1849
5 TraesCS4A01G488500 chr7A 3007832 3008741 909 False 704.000000 704 81.363000 945 1847 1 chr7A.!!$F1 902
6 TraesCS4A01G488500 chr7A 3323804 3324663 859 True 691.000000 691 81.385000 985 1856 1 chr7A.!!$R1 871
7 TraesCS4A01G488500 chr7A 2569626 2571512 1886 False 596.300000 1101 90.807000 2 1865 2 chr7A.!!$F2 1863
8 TraesCS4A01G488500 chr7A 2856708 2860951 4243 False 525.400000 974 86.718500 4 2004 2 chr7A.!!$F3 2000
9 TraesCS4A01G488500 chr7D 3434044 3435229 1185 True 881.000000 881 80.470000 819 2032 1 chr7D.!!$R1 1213
10 TraesCS4A01G488500 chr7D 3739662 3741459 1797 False 589.750000 1099 89.344000 2 1865 2 chr7D.!!$F1 1863


AutoCloner calculated primer pairs

These primer pairs were selected by Autocloner to minimize primer penalty (calculated by Primer3), whilst remaining within the specified product range where possible.

Download AutoCloner chosen primers as CSV: (Forward) (Reverse)
Position MSA Position Penalty Sequence TM GC Self any TH Self end TH Hairpin End Stability Orientation
563 3237 0.392193 AAGGCGAGCATCCTGTTCTG 60.392 55.0 0.0 0.0 33.43 3.02 F
Position MSA Position Penalty Sequence TM GC Self any TH Self end TH Hairpin End Stability Orientation
2561 5416 0.032267 GTATCTGGTGCTACCGAGGC 59.968 60.0 0.0 0.0 42.58 4.7 R


All possible primers

Listed here are all possible primers based on the multiple sequence alignment and SNPs generated previously. A lower penalty score indicates a better primer. If the AutoCloner primers did not work during PCR or sequencing, these could be used as alternatives.

Download all possible primers: (Forward) (Reverse)



Based at the University of Bristol with support from BBSRC. BBSRC icon Bristol icon

AutoCloner maintained by Alex Coulton.