| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 3,575,608 | 0 | T | G | 35.1% | 213.2 / 33.2 | 114 | T89P (ACC→CCC) | prlC | oligopeptidase A |
| Reads supporting (aligned to +/- strand): ref base T (16/58); new base G (40/0); total (56/58) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.42e-18 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
Only 100 of 135 total aligned reads displayed.
CAGGTCGCGATACGCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCC > REL606/3575534‑3575683 | cAGGTCGCGATACGCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGtt > 2:2466991/1‑76 (MQ=255) aGGTCGCGATACGCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGttt < 2:181353/76‑1 (MQ=255) aGGTCGCGATACGCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGttt < 2:2177156/76‑1 (MQ=255) cgATACGCTTTATACAGCCCTTCATGTTGTCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGt < 1:1480425/76‑1 (MQ=255) cgATACGCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGt < 2:703920/76‑1 (MQ=255) cgATACGCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGt > 2:310321/1‑76 (MQ=255) cgATACGCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGt > 1:2417235/1‑76 (MQ=255) gATACGCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTa < 2:2921213/76‑1 (MQ=255) tACGCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACg > 1:3709812/1‑76 (MQ=255) aCGCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGc > 1:2652141/1‑76 (MQ=255) gCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGcct > 1:3505604/1‑74 (MQ=255) tttATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCa < 2:477737/76‑1 (MQ=255) tttATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCa < 2:1849756/76‑1 (MQ=255) tttATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCa > 2:3233052/1‑76 (MQ=255) ttATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCAc < 2:2984261/76‑1 (MQ=255) ttATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCAc < 2:3793538/76‑1 (MQ=255) tataCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTccac > 2:2717151/1‑73 (MQ=255) tataCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACg > 1:3087740/1‑76 (MQ=255) ataCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGc < 2:1786912/76‑1 (MQ=255) ataCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGc < 2:771111/76‑1 (MQ=255) ataCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGc < 2:5020626/76‑1 (MQ=255) ataCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGc < 2:4237761/76‑1 (MQ=255) ataCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCCTCAccc > 1:3458463/1‑74 (MQ=255) ataCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTATTGGTACGCTTCCCGc > 1:2379218/1‑76 (MQ=255) taCAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCa < 2:2100676/76‑1 (MQ=255) ccAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCATGCAg > 2:4495073/2‑76 (MQ=255) cAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAAGCTTCACGCAGt > 2:2170336/1‑76 (MQ=255) aGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGtt < 1:1040910/76‑1 (MQ=255) aGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGtt < 1:4318211/76‑1 (MQ=255) gCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCAGTTc > 2:3323596/1‑76 (MQ=255) cccTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCCTCACGCAGtttc > 2:3967508/1‑74 (MQ=255) ccTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCAGTTCCg > 1:1736113/1‑76 (MQ=255) ccTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCAGTTCCg > 1:1967306/1‑76 (MQ=255) cTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCgg < 1:945432/76‑1 (MQ=255) cTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCgg < 1:489341/76‑1 (MQ=255) cATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCt < 1:1590274/76‑1 (MQ=255) cATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTTACGCAGTTCCggggt > 1:3803856/1‑74 (MQ=255) cATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCAGTTCCGGGCt > 2:121684/1‑76 (MQ=255) aTGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTa < 1:4292090/76‑1 (MQ=255) aTGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCAGTTCCGGGCTa > 2:1111600/1‑76 (MQ=255) gTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCCCGCAGTTCCGGGCTAtt > 1:4957957/1‑76 (MQ=255) tGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTAtttt < 2:4608604/76‑1 (MQ=255) tGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTTGTACGCTTCACGCCCTTTCCGGGTAtttt > 1:2566885/1‑76 (MQ=255) tGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCAGTTCCGGGCTAtttt > 1:4835657/1‑76 (MQ=255) cccTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTc < 1:3334238/76‑1 (MQ=255) cTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCAc < 1:1158686/76‑1 (MQ=255) cTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTTTTCGTACGCTTCACGCAGTTCCGGGGTATTTTTccc > 1:4120412/1‑74 (MQ=255) cTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCCCGCAGTTCCGGGGTATTTTTCAc > 1:1530356/1‑76 (MQ=255) cTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTTAc > 1:2981642/1‑76 (MQ=255) tACCGAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCAcc > 2:3938969/1‑76 (MQ=255) tACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTTTTCGTACGCTTCACGCAGTTCCGGGGCATTTTTCAcc > 1:3892757/1‑76 (MQ=255) tACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTAAGCTTCACGCAGTTCCGGGCTATTTTTTAcc > 2:3265160/1‑76 (MQ=255) aCCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCg < 1:448077/76‑1 (MQ=255) cccAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGa < 2:2859476/76‑1 (MQ=255) cAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGt < 1:4434413/76‑1 (MQ=255) cAGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCCTTACGCCGTTCCGGGGTAATTTTTACCCAGt > 2:2551806/1‑76 (MQ=255) aGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCCGTTCCGGGCTATTTTTCACCCAggt > 1:3793945/1‑74 (MQ=255) aGGTGCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCAGTTCCGGGGTATTTTTTACCGAGtt > 2:4525141/1‑76 (MQ=255) ggTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTc < 1:3898492/76‑1 (MQ=255) gTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCCGTTTCGGGCTATTTTTCACCCGGTTcc > 1:564836/1‑75 (MQ=255) gCTGTATTCCGACAGCAGCGGCAGGGGTTGTTCCGACGCTTCACGCAGTTCCGGGGTACTTTTCAGCCAGTTCacg > 2:2316383/1‑74 (MQ=255) cTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGt < 2:1830213/76‑1 (MQ=255) cTGTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGt > 1:3371467/1‑76 (MQ=255) gTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTgg < 1:4417427/76‑1 (MQ=255) gTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTgg < 2:2898713/76‑1 (MQ=255) gTATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCAGGTCCGGGCTATTTTTCACCGAGTTCAGGTgg > 2:2353758/1‑76 (MQ=255) tATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGc < 1:3089397/76‑1 (MQ=255) tATTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGc > 2:1417030/1‑76 (MQ=255) aTTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGGGGCt > 2:4192314/1‑76 (MQ=255) aTTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACCCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCt > 1:2376722/1‑76 (MQ=255) aTTCCGACAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCt > 1:387620/1‑76 (MQ=255) ttCCGACAGCAGCGGCAGTGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTg < 2:3709305/76‑1 (MQ=255) ttCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTg < 2:4907133/76‑1 (MQ=255) tCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGa < 1:1161872/76‑1 (MQ=255) tCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGa < 2:369370/76‑1 (MQ=255) ccGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGGTCAGGTGGCTGAc > 1:358603/1‑76 (MQ=255) aCAGCAGCGGCAGGGGTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTTACCGAGTTCAGGTGGCTGACCgg > 2:4159027/1‑76 (MQ=255) cagcagCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCggg < 2:1396299/76‑1 (MQ=255) cagcagCGGCAGGGGTTGTTCGTACGCTTCCCCCAGTTCCGGGCTATTTTTTACCCAGTTCAGGTGGCTGGCCcgg > 2:1931543/1‑76 (MQ=255) agcagcGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCgggg < 2:3466941/76‑1 (MQ=255) agcagcGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCgggg < 2:5015751/76‑1 (MQ=255) agcagcGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCgggg < 2:5015752/76‑1 (MQ=255) cagcGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGgag < 1:1392914/76‑1 (MQ=255) agcGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGgaga < 1:807590/76‑1 (MQ=255) agcGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGgaga < 1:1931543/76‑1 (MQ=255) agcGGCAGGGGTTGTTTGTACGCTTCACGCCGTTCCGGGCTATTTTTTACCCAGTTCCGGTGGCTGACCGGGggga > 2:3451904/1‑76 (MQ=255) gcGGCAGGGTTTGTTCGTACGCTTCGCGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGagaa < 1:1670087/76‑1 (MQ=255) ggCAGGGTTTGTTCGTACGCTTCACGCAGCTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGagaaga > 1:37961/1‑76 (MQ=255) gCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGAt < 1:2531694/76‑1 (MQ=255) gCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGAt < 1:3292721/76‑1 (MQ=255) cAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATa > 2:1169940/1‑76 (MQ=255) cAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATa < 2:2343871/76‑1 (MQ=255) aGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATAc > 2:4103103/1‑76 (MQ=255) aGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATAc > 2:3367561/1‑76 (MQ=255) aGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATAc > 2:4623322/1‑76 (MQ=255) aGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATAc > 1:13962/1‑76 (MQ=255) aGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATAc > 1:2709015/1‑76 (MQ=255) gggTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACg < 1:2712105/76‑1 (MQ=255) ggTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGc < 1:768564/76‑1 (MQ=255) ggTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGc < 1:3230569/76‑1 (MQ=255) | CAGGTCGCGATACGCTTTATACAGCCCTTCATGTTGCCCTACCCAGGTGCTGTATTCCGACAGCAGCGGCAGGGTTTGTTCGTACGCTTCACGCAGTTCCGGGCTATTTTTCACCGAGTTCAGGTGGCTGACCGGGGAGAAGATACGCCC > REL606/3575534‑3575683 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |