Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,795,650 | 0 | A | C | 31.6% | 118.1 / 13.0 | 57 | V39G (GTG→GGG) | celF | cryptic phospho‑beta‑glucosidase, NAD(P)‑binding |
Reads supporting (aligned to +/- strand): ref base A (34/5); new base C (0/18); total (34/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCTGC > REL606/1795576‑1795725 | cGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATccac < 2:1917887/76‑1 (MQ=255) cGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATccac < 2:1928859/76‑1 (MQ=255) gCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATccacc > 1:1075207/1‑76 (MQ=255) gCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATccacc < 1:1561303/76‑1 (MQ=255) ccAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATccacca < 1:559480/76‑1 (MQ=255) cAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAg < 1:1852096/76‑1 (MQ=255) cAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTGTCGCTTTACCGCCTTCGACATCCACCAg < 1:1552507/76‑1 (MQ=255) gCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCAGCAGcc < 2:385878/76‑1 (MQ=255) cGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGccn < 1:1684651/76‑2 (MQ=255) tATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATaa > 2:549654/1‑76 (MQ=255) tATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATaa > 2:1911409/1‑76 (MQ=255) atcaatcaTCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAAt > 1:1516412/1‑76 (MQ=255) aatcaTCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCg > 2:1379573/1‑76 (MQ=255) aatcaTCCGTTGGCAGAGATCGAAAAAAATATCCATTTTCCCTTTTCCCCCTTCGACACCCCCCAGCCATAATTCg < 2:505546/76‑1 (MQ=255) atcCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGa > 1:862466/1‑76 (MQ=255) ccGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGAcc > 2:1719378/1‑76 (MQ=255) ttGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGc > 2:1671688/1‑76 (MQ=255) ttGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGc > 1:92568/1‑76 (MQ=255) tGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCa < 1:716516/76‑1 (MQ=255) tGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCa > 2:980900/1‑76 (MQ=255) gagaTCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCg > 2:1545991/1‑76 (MQ=255) gagaTCAAAAATAATATCCAGTTTCGCTTTCCCCCCTTCGCCACCCCCCAGCCATAATTCGCTGACCGGCAATTCg < 1:868360/76‑1 (MQ=255) gaTCGAAAATAATATCCAGTTTCGCTTTCCCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTg < 1:1545991/76‑1 (MQ=255) aTCGAAAATAATATCCAGTTTCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGa < 1:363312/76‑1 (MQ=255) tCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGAt > 1:625598/1‑76 (MQ=255) caaaaaTAATATCCAGTTTCGCTTTACCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATa < 1:666543/74‑1 (MQ=255) cGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa > 2:1370396/1‑76 (MQ=255) cGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATa > 2:1370385/1‑76 (MQ=255) aaaaataaaatCCAGTTCCCCTTTCCCCCCTCCGCCACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATaa < 2:1322405/75‑1 (MQ=255) aaaTAATATCCATTTTCGCTTTCCCCCCTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACg < 1:1943730/76‑1 (MQ=255) aaaTAATATCCAGTTTCGCTTTACCCCCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACg < 1:1866257/76‑1 (MQ=255) nataataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGc > 1:1722308/2‑76 (MQ=255) aataatATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGc > 2:351651/1‑76 (MQ=255) ataataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCt > 2:1831693/1‑76 (MQ=255) ataataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCt > 2:1796665/1‑76 (MQ=255) taataTCCAGTTTCGCTTTCCCACCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCtt < 1:348044/76‑1 (MQ=255) aataTCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTta > 2:432055/1‑76 (MQ=255) atatCCAGTTTCGCTTTACCCCCTCCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaa < 1:1538336/76‑1 (MQ=255) atCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaata > 2:1989992/1‑76 (MQ=255) tCCAGTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaataa < 1:1442824/76‑1 (MQ=255) cccAGTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTtaataa < 1:1907820/75‑1 (MQ=255) ccAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATaaa > 2:1855830/1‑76 (MQ=255) gTTTCGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATcc < 2:1149021/76‑1 (MQ=255) ttCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCtt > 1:1343286/1‑76 (MQ=255) ttCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCtt > 2:970091/1‑76 (MQ=255) tcngCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCtt > 2:1403494/4‑76 (MQ=255) ccGCTTTACCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTc < 2:813887/75‑1 (MQ=255) cGCTTTCCCCCCTTCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc < 1:524559/76‑1 (MQ=255) cGCTTTCCCCCCTTCCACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc < 2:1267889/76‑1 (MQ=255) cGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc > 1:357612/1‑76 (MQ=255) cGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc > 2:692217/1‑76 (MQ=255) cGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc > 2:662093/1‑76 (MQ=255) cGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTcc > 2:1294147/1‑76 (MQ=255) gCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCa > 2:130252/1‑76 (MQ=255) cTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAg > 2:464323/1‑76 (MQ=255) tttACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGt > 1:1946914/1‑76 (MQ=255) tttACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGt > 1:1354595/1‑76 (MQ=255) tcacccccTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTa < 1:155183/74‑1 (MQ=255) acccccTTCGACATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAAc < 1:63323/76‑1 (MQ=255) cccTTCGCCATCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTcc < 1:70457/75‑1 (MQ=255) cccTCCGACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTcc < 1:607955/75‑1 (MQ=255) tCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCgggg > 1:71480/1‑76 (MQ=255) gACACCCCCCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGta < 2:1766710/76‑1 (MQ=255) cATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGtata > 1:206014/1‑76 (MQ=255) cATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGtata > 2:1412545/1‑76 (MQ=255) caccaGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAgctg > 1:1518535/1‑76 (MQ=255) accaGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAgctgc > 2:685943/1‑76 (MQ=255) | CGCCAGCGTTATCAATCATCCGTTGGCAGAGATCGAAAATAATATCCAGTTTCGCTTTACCACCTTCGACATCCACCAGCCATAATTCGCTGACCGGCAATTCGTGATAACGCTTAATAAATCCTTCCAGTAACTCCGGGGTATAGCTGC > REL606/1795576‑1795725 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |