Multiple sequence alignment - TraesCS4A01G016400

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 TraesCS4A01G016400 chr4A 100.000 6870 0 0 1 6870 11182282 11189151 0.000000e+00 12687.0
1 TraesCS4A01G016400 chr4D 92.193 2946 142 43 2969 5870 457504132 457501231 0.000000e+00 4085.0
2 TraesCS4A01G016400 chr4D 88.053 3005 198 69 1 2927 457507138 457504217 0.000000e+00 3411.0
3 TraesCS4A01G016400 chr4D 91.954 522 23 6 5922 6434 457501218 457500707 0.000000e+00 713.0
4 TraesCS4A01G016400 chr4D 91.160 362 22 6 6073 6434 457498148 457497797 3.730000e-132 483.0
5 TraesCS4A01G016400 chr4B 89.870 2458 157 39 1 2399 571474787 571472363 0.000000e+00 3075.0
6 TraesCS4A01G016400 chr4B 90.808 2252 98 41 4491 6701 571470360 571468177 0.000000e+00 2911.0
7 TraesCS4A01G016400 chr4B 91.731 1040 49 12 3423 4441 571471386 571470363 0.000000e+00 1410.0
8 TraesCS4A01G016400 chr4B 91.224 433 29 3 2969 3393 571471815 571471384 1.280000e-161 580.0
9 TraesCS4A01G016400 chr4B 94.624 279 9 3 2502 2780 571472254 571471982 1.770000e-115 427.0
10 TraesCS4A01G016400 chr4B 85.350 157 22 1 6702 6857 53190626 53190782 1.980000e-35 161.0
11 TraesCS4A01G016400 chr7A 84.069 408 39 11 2511 2899 170288198 170288598 3.030000e-98 370.0
12 TraesCS4A01G016400 chr7A 87.135 171 17 5 6702 6870 35593699 35593532 9.100000e-44 189.0
13 TraesCS4A01G016400 chr7B 88.235 170 16 4 6702 6870 49608898 49609064 4.200000e-47 200.0
14 TraesCS4A01G016400 chr1D 87.059 170 19 3 6702 6870 64115332 64115165 9.100000e-44 189.0
15 TraesCS4A01G016400 chr1D 85.294 170 23 2 6702 6870 86114839 86114671 2.550000e-39 174.0
16 TraesCS4A01G016400 chr1D 93.478 46 3 0 4444 4489 333492168 333492213 1.240000e-07 69.4
17 TraesCS4A01G016400 chr2B 86.047 172 18 6 6701 6870 45036407 45036240 5.470000e-41 180.0
18 TraesCS4A01G016400 chr3D 87.261 157 18 2 6716 6870 41428750 41428594 1.970000e-40 178.0
19 TraesCS4A01G016400 chr3D 84.211 171 24 3 6702 6870 47601005 47601174 5.510000e-36 163.0
20 TraesCS4A01G016400 chr7D 84.024 169 24 3 6702 6870 440005490 440005325 7.130000e-35 159.0
21 TraesCS4A01G016400 chr1B 84.177 158 11 1 3828 3971 6714117 6713960 2.580000e-29 141.0
22 TraesCS4A01G016400 chr1B 84.177 158 11 1 3828 3971 6718503 6718346 2.580000e-29 141.0
23 TraesCS4A01G016400 chr1B 83.544 158 12 1 3828 3971 6705284 6705127 1.200000e-27 135.0
24 TraesCS4A01G016400 chr1B 82.911 158 13 1 3828 3971 6709701 6709544 5.590000e-26 130.0
25 TraesCS4A01G016400 chr1B 82.278 158 14 1 3828 3971 6701107 6700950 2.600000e-24 124.0
26 TraesCS4A01G016400 chr1B 94.444 72 4 0 3828 3899 6721662 6721591 2.030000e-20 111.0
27 TraesCS4A01G016400 chr1B 94.444 72 4 0 3828 3899 6725258 6725187 2.030000e-20 111.0


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 TraesCS4A01G016400 chr4A 11182282 11189151 6869 False 12687.0 12687 100.0000 1 6870 1 chr4A.!!$F1 6869
1 TraesCS4A01G016400 chr4D 457497797 457507138 9341 True 2173.0 4085 90.8400 1 6434 4 chr4D.!!$R1 6433
2 TraesCS4A01G016400 chr4B 571468177 571474787 6610 True 1680.6 3075 91.6514 1 6701 5 chr4B.!!$R1 6700


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
612 624 0.037697 TACTACACAAGTGCGCCAGG 60.038 55.0 4.18 0.00 39.39 4.45 F
647 665 0.384309 GTGATGGATAGCACGGACGA 59.616 55.0 0.00 0.00 33.04 4.20 F
1823 1895 0.246635 TTCTCGTGTCAAGGCTAGCC 59.753 55.0 27.19 27.19 0.00 3.93 F
2633 2756 0.178981 ATGGGTAGGTGGATGCATGC 60.179 55.0 11.82 11.82 0.00 4.06 F
2707 2830 0.106569 TTGCCTGGCCATATGGAGTG 60.107 55.0 26.47 12.81 37.39 3.51 F
4453 4661 0.033306 ACTACTCCTTCCCTCCGTCC 60.033 60.0 0.00 0.00 0.00 4.79 F
4454 4662 0.033405 CTACTCCTTCCCTCCGTCCA 60.033 60.0 0.00 0.00 0.00 4.02 F
Position MSA Position Penalty Sequence TM GC Self any TH Self end TH Hairpin End Stability Orientation
1776 1848 0.243907 GAGAAGCTGCACGTCCTACA 59.756 55.0 1.02 0.00 0.00 2.74 R
2482 2605 0.454285 GTGTGTGTGTGTGTGTGTGC 60.454 55.0 0.00 0.00 0.00 4.57 R
2670 2793 0.106519 AAGTCCACTTGCCAAGCTGT 60.107 50.0 3.88 0.00 34.38 4.40 R
3491 3678 0.250513 GTCTGCGTCTTCCCCATCTT 59.749 55.0 0.00 0.00 0.00 2.40 R
4643 4862 0.173481 TCTCCTTGGCGAAGATCACG 59.827 55.0 10.40 7.74 0.00 4.35 R
5393 5637 0.035056 ATCAAGTGGCGGAAGGGAAG 60.035 55.0 0.00 0.00 0.00 3.46 R
6268 6560 0.110238 CATCCAAAACAGTGCGTCGG 60.110 55.0 0.00 0.00 0.00 4.79 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.