Multiple sequence alignment - TraesCS2D01G158400

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 TraesCS2D01G158400 chr2D 100.000 5363 0 0 1 5363 102086724 102081362 0.000000e+00 9904.0
1 TraesCS2D01G158400 chr2D 83.710 442 28 20 4907 5345 537490969 537491369 1.410000e-100 377.0
2 TraesCS2D01G158400 chr2B 90.717 2413 155 28 2536 4912 153225648 153223269 0.000000e+00 3151.0
3 TraesCS2D01G158400 chr2B 88.660 2037 130 50 564 2539 153227670 153225674 0.000000e+00 2388.0
4 TraesCS2D01G158400 chr2B 90.749 454 35 4 4912 5362 725732165 725732614 2.760000e-167 599.0
5 TraesCS2D01G158400 chr2B 90.030 331 16 8 90 410 153228088 153227765 3.870000e-111 412.0
6 TraesCS2D01G158400 chr2B 82.759 87 11 4 5171 5255 751333229 751333145 2.070000e-09 75.0
7 TraesCS2D01G158400 chr2A 90.746 2172 122 23 2539 4684 101096096 101093978 0.000000e+00 2824.0
8 TraesCS2D01G158400 chr2A 90.596 1276 61 21 1294 2539 101097371 101096125 0.000000e+00 1637.0
9 TraesCS2D01G158400 chr2A 84.350 754 40 29 485 1215 101098200 101097502 0.000000e+00 667.0
10 TraesCS2D01G158400 chr2A 91.414 198 7 2 4725 4912 101093975 101093778 4.120000e-66 263.0
11 TraesCS2D01G158400 chr2A 78.696 230 17 12 271 492 101098715 101098510 2.030000e-24 124.0
12 TraesCS2D01G158400 chr1B 91.247 457 33 5 4909 5362 416580296 416580748 2.740000e-172 616.0
13 TraesCS2D01G158400 chr7B 90.426 470 34 8 4900 5362 245853060 245853525 4.590000e-170 608.0
14 TraesCS2D01G158400 chr6B 90.949 453 36 3 4912 5362 490207611 490207162 5.940000e-169 604.0
15 TraesCS2D01G158400 chr6B 90.566 53 4 1 4897 4948 564989754 564989702 9.640000e-08 69.4
16 TraesCS2D01G158400 chrUn 88.136 472 30 11 4912 5361 9993411 9992944 6.110000e-149 538.0
17 TraesCS2D01G158400 chr7D 93.594 281 12 4 4912 5190 213792535 213792811 1.070000e-111 414.0
18 TraesCS2D01G158400 chr4A 90.036 281 22 4 4912 5190 231098371 231098647 5.110000e-95 359.0
19 TraesCS2D01G158400 chr5A 83.704 270 39 4 4945 5210 429122523 429122791 3.210000e-62 250.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 TraesCS2D01G158400 chr2D 102081362 102086724 5362 True 9904.000000 9904 100.000000 1 5363 1 chr2D.!!$R1 5362
1 TraesCS2D01G158400 chr2B 153223269 153228088 4819 True 1983.666667 3151 89.802333 90 4912 3 chr2B.!!$R2 4822
2 TraesCS2D01G158400 chr2A 101093778 101098715 4937 True 1103.000000 2824 87.160400 271 4912 5 chr2A.!!$R1 4641


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
516 882 0.036388 CCCCGAAAAGCTAGCACTGA 60.036 55.000 18.83 0.0 0.00 3.41 F
564 930 0.178992 CTGAGCTTTTCACCACCCCA 60.179 55.000 0.00 0.0 0.00 4.96 F
1536 2041 0.182299 GATCTAGCCTTGCCTTGCCT 59.818 55.000 0.00 0.0 0.00 4.75 F
1538 2043 0.322816 TCTAGCCTTGCCTTGCCTTG 60.323 55.000 0.00 0.0 0.00 3.61 F
1757 2263 0.445436 CGGTGCTGAACTGCTGATTC 59.555 55.000 0.00 0.0 0.00 2.52 F
2736 3278 1.202302 GCTTCAGGCCTTTATTTCCGC 60.202 52.381 0.00 0.0 34.27 5.54 F
3945 4504 0.250338 ACGCTGTGAAGGGTTTCTCC 60.250 55.000 0.00 0.0 34.31 3.71 F
Position MSA Position Penalty Sequence TM GC Self any TH Self end TH Hairpin End Stability Orientation
1814 2321 0.605319 TCCAAAGCAGTTGCCGGTAG 60.605 55.000 1.90 0.0 43.38 3.18 R
2109 2619 1.816224 GGCATTTACAATCCGTGGTGT 59.184 47.619 0.00 0.0 0.00 4.16 R
2772 3314 0.904394 AACGCATTCCTTGGGCCATT 60.904 50.000 7.26 0.0 41.83 3.16 R
2874 3417 1.562475 GGGCCTGGTGGTATCAAGTTA 59.438 52.381 0.84 0.0 35.27 2.24 R
3544 4099 2.970324 GGCAACACGTTCTGCGGA 60.970 61.111 11.43 0.0 46.52 5.54 R
3969 4528 1.069973 TGCGCGAAACAAGAGTTGATG 60.070 47.619 12.10 0.0 38.17 3.07 R
4912 5508 0.102481 ACGCACGGGCATTCTACTAG 59.898 55.000 11.77 0.0 41.24 2.57 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.