Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,795,650 | 0 | A | C | 34.7% | 92.2 / 13.9 | 49 | V39G (GTG→GGG) | celF | cryptic phospho‑beta‑glucosidase, NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base A (28/4); new base C (0/17); total (28/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.21e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCT > REL606/1795575‑1795723 | aCGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATcca > 2:500826/1‑76 (MQ=255) cGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATccac < 2:130462/76‑1 (MQ=255) gCCAGCGTTATCAAACATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATccacc > 1:874987/1‑76 (MQ=255) cAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAg < 2:1352953/76‑1 (MQ=255) aGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGc > 2:878990/1‑76 (MQ=255) gTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCAt > 1:397389/1‑76 (MQ=255) atcaatcaCCCGTGGGCAGAGACCGAAAATAATATCCATTTCCGCTTTCCCTCCTTCGACACCCCCCAGCCATAAt < 2:864334/76‑1 (MQ=255) caatcaTCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTc < 1:60107/76‑1 (MQ=255) atcCGTTGGCGGAGATCGAAAATAATATCCAGTTTCGCTTTACCCCCTTCGACACCCCCCAGCCATAATTCGCTGa < 2:808816/76‑1 (MQ=255) atcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGa > 2:458475/1‑76 (MQ=255) atcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGa > 1:362718/1‑76 (MQ=255) atcCGGTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTTCCCCCTTCGCCAACCACCAGCCATAATTCGCTGa < 2:289745/76‑1 (MQ=255) ccGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGAcc > 1:1358034/1‑76 (MQ=255) cGTTGGCAGAGATCNAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCg > 1:886643/1‑76 (MQ=255) gTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCgg > 1:834549/1‑76 (MQ=255) gTTGGCAGAGACCAAAAATAATATCCAGTTTCGCTTTCCCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCgg < 1:233345/76‑1 (MQ=255) ttGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGGCCGGc > 2:1629680/1‑76 (MQ=255) tGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa > 2:540644/1‑76 (MQ=255) tGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa > 1:97894/1‑76 (MQ=255) gCAGAGAACAAAAAAAATATCCTTTTTCGCTTTCCCGCCTTGGCCATCCACCAGCCATAATTCGCTGACCGGCAAt < 2:1348002/76‑1 (MQ=255) agagATCAAAAATAATATCCAGTTTCGCTTTCCCACCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTc < 1:983053/76‑1 (MQ=255) aTCGAAAATAATATCCAGTTTCGCTTTACCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGa < 1:1167836/76‑1 (MQ=255) caaaaaaaACACCCGGTCTCGCTTTCCAACCTGGCACTCCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa < 2:642242/74‑1 (MQ=255) gAAAATAATATCCATTTCCGCTTTCCCCCCTCCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATaa < 1:1182349/76‑1 (MQ=255) aaaaaaaATATCCAGTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATaa < 1:727979/75‑1 (MQ=255) aaaaTAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAAc > 2:101557/1‑76 (MQ=255) aaaaTAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAAc > 1:1130655/1‑76 (MQ=255) aaaTAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACg > 2:951958/1‑76 (MQ=255) aataatATCCATTTCCGCTTTCCCCCCTTCGCCACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGc < 2:618711/76‑1 (MQ=255) taataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCtt > 1:479392/1‑76 (MQ=255) taataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCtt > 2:926493/1‑76 (MQ=255) aataTCCAGTTTCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTta < 2:1619181/76‑1 (MQ=255) aataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTta > 2:1148448/1‑76 (MQ=255) aataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTta > 2:678178/1‑76 (MQ=255) atgtCCCGTTTCCCTTCACCACCTGCGACACCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaa < 2:1374236/73‑1 (MQ=255) tatCCATTTCCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaat < 1:254039/76‑1 (MQ=255) tatCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaat > 2:1151911/1‑76 (MQ=255) atCCAGTTTCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaata < 1:619902/76‑1 (MQ=255) cccAGTTTCCCTTTCCCCCCTCCGCCACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaataa < 1:1143452/75‑1 (MQ=255) cattttCGCTTCACCCCCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAAt < 2:659916/73‑1 (MQ=255) ttCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCtt < 2:1090171/76‑1 (MQ=255) tCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTc > 2:963495/1‑76 (MQ=255) cGCTTTCCCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc < 2:1358034/76‑1 (MQ=255) cTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAg > 2:950360/1‑76 (MQ=255) tttACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGt > 1:924894/1‑76 (MQ=255) ttACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTa > 2:543652/1‑76 (MQ=255) ttACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTa > 1:558066/1‑76 (MQ=255) tACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTaa > 2:433673/1‑76 (MQ=255) cccTCCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTcc < 1:951958/75‑1 (MQ=255) accTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTcc > 1:1192249/1‑76 (MQ=255) ccTCCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCg < 2:631550/76‑1 (MQ=255) ttCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCggg < 2:96780/76‑1 (MQ=255) ttCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCggg < 2:1154347/76‑1 (MQ=255) cacccaccaGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGtata < 2:1567626/73‑1 (MQ=255) aTCCACCAGCCATAATTCGCTGACCGTCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAg > 1:224376/1‑76 (MQ=255) tCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAgc > 1:431850/1‑76 (MQ=255) tCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAgc > 2:557774/1‑76 (MQ=255) ccaccaGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAgct > 2:1376664/1‑76 (MQ=255) | ACGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCT > REL606/1795575‑1795723 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |