Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,202,004 | A→G | Y186H (TAT→CAT) | iclR ← | DNA‑binding transcriptional repressor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,202,004 | 0 | A | G | 96.0% | 191.9 / ‑1.4 | 50 | Y186H (TAT→CAT) | iclR | DNA‑binding transcriptional repressor |
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (30/18); total (32/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.30e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01 |
GAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACC > REL606/4201931‑4202076 | gAATAACCCCGTTTGCGCGTTTGTGGCAGAATTTCTTTTAAAAGCTCAGGAGACACCAGCGTTGCGTGGGTATAgg < 2:588881/76‑1 (MQ=255) aaTAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGc > 2:1562353/1‑76 (MQ=255) aTAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCa > 1:1028562/1‑76 (MQ=255) aaCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATg > 2:985295/1‑76 (MQ=255) ccccGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTa > 2:1680163/1‑76 (MQ=255) cccGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTaa < 2:395897/76‑1 (MQ=255) cGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAAcc > 1:1553851/1‑76 (MQ=255) ttGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCtt < 1:102895/76‑1 (MQ=255) ttGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCtt > 1:765620/1‑76 (MQ=255) tGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCttt < 2:1442100/76‑1 (MQ=255) tGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGAGGGTATGGGCATGTAACCCttt > 2:642161/1‑76 (MQ=255) gcgcGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTg > 2:1424144/1‑76 (MQ=255) cgcgTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGc > 2:834084/1‑76 (MQ=255) cgcgTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGc > 1:1197316/1‑76 (MQ=255) cgTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCgg < 1:1611323/76‑1 (MQ=255) cgTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCgg < 1:285581/76‑1 (MQ=255) cgTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCgg < 1:563086/76‑1 (MQ=255) cgTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCgg < 1:648584/76‑1 (MQ=255) gTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGt > 1:44362/1‑76 (MQ=255) tttGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTg < 1:1055277/76‑1 (MQ=255) ttGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTgc < 2:1566443/76‑1 (MQ=255) tGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTgca > 2:506855/1‑76 (MQ=255) tGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTgca < 1:1384060/76‑1 (MQ=255) ggCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTgcagc < 2:139080/76‑1 (MQ=255) cGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGcagcag > 2:822137/1‑76 (MQ=255) agaTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCt < 1:850786/76‑1 (MQ=255) aTCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGCGCAGCAGCtcc > 2:1053436/1‑74 (MQ=255) tcttTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCa > 2:89388/1‑76 (MQ=255) tcttTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCa > 2:450546/1‑76 (MQ=255) tttAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCt > 1:1611088/1‑76 (MQ=255) tttAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCt > 1:591664/1‑76 (MQ=255) ttAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTg < 2:1404539/76‑1 (MQ=255) tAAATGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGt > 2:854520/1‑76 (MQ=255) aaTGCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGttc < 2:1508631/76‑1 (MQ=255) gCACAGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGttcttc > 2:331748/1‑76 (MQ=255) cacaGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCg > 2:890254/1‑76 (MQ=255) caGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCt < 1:812664/76‑1 (MQ=255) aGGAGACACCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTc > 2:285275/1‑76 (MQ=255) acacCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTg > 2:1370631/1‑76 (MQ=255) acacCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTg > 2:1232143/1‑76 (MQ=255) cacCAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTgg > 2:1108098/1‑76 (MQ=255) ccAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGc > 2:923978/1‑76 (MQ=255) cAGCGTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCt < 1:1079894/76‑1 (MQ=255) gTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTaaaa < 1:626457/76‑1 (MQ=255) gTTGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTaaaa < 1:700638/76‑1 (MQ=255) tGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAg > 1:339496/1‑76 (MQ=255) tGCGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAg > 1:176882/1‑76 (MQ=255) cGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGc > 2:1671140/1‑76 (MQ=255) cGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGc > 1:982517/1‑76 (MQ=255) cGTGGGTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGc > 2:1214125/1‑76 (MQ=255) gTATGGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTAc > 1:375944/1‑76 (MQ=255) tatgGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTAcc > 2:880802/1‑76 (MQ=255) | GAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTACC > REL606/4201931‑4202076 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |