Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,928,304 | C→T | R30C (CGT→TGT) | rho → | transcription termination factor Rho |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,928,304 | 0 | C | T | 93.1% | 198.9 / 5.2 | 58 | R30C (CGT→TGT) | rho | transcription termination factor Rho |
Reads supporting (aligned to +/- strand): ref base C (1/3); new base T (21/33); total (22/36) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01 |
AATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCG > REL606/3928230‑3928376 | aaTTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCg < 1:523229/76‑1 (MQ=255) ttAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTa < 1:545649/76‑1 (MQ=255) tAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTaa < 2:943987/76‑1 (MQ=255) tAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTaa > 1:1503879/1‑76 (MQ=255) tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa > 1:869186/1‑76 (MQ=255) aaaGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAg > 2:222701/1‑76 (MQ=255) aaGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGc < 2:574665/76‑1 (MQ=255) aGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCa < 1:1010200/76‑1 (MQ=255) aaTACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAgg < 2:1548700/76‑1 (MQ=255) aTACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGa > 1:388377/1‑76 (MQ=255) tACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGAc > 2:1610773/1‑76 (MQ=255) tACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGAc > 2:1020336/1‑76 (MQ=255) gCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACatt > 1:140007/1‑76 (MQ=255) ggTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACattatt > 2:796694/1‑76 (MQ=255) gTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTAttt < 1:1133970/76‑1 (MQ=255) gTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTAttt < 2:297855/76‑1 (MQ=255) tttCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTAtttt > 1:427694/1‑76 (MQ=255) ttCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGGATGTGTAAGCAGGAAAtttttttt > 2:1360605/1‑76 (MQ=255) cTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGc < 2:167127/76‑1 (MQ=255) tGAGCTGATCACTCTCGGCGAAAATATGGGCGTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGcc < 2:699918/76‑1 (MQ=255) gAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCa > 1:684932/1‑76 (MQ=255) gAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCa < 1:673915/76‑1 (MQ=255) ngCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATATGTAAGCAGGACATTATTTTTGCCAt > 1:1357386/2‑76 (MQ=255) aGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCAt > 2:481700/1‑76 (MQ=255) tGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCt > 1:370277/1‑76 (MQ=255) gATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTg > 1:334560/1‑76 (MQ=255) cACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAag < 2:1632856/76‑1 (MQ=255) ctctCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAagca < 2:154709/76‑1 (MQ=255) tctcGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAagcag < 1:440788/76‑1 (MQ=255) ctcGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAagcagc < 1:1533942/76‑1 (MQ=255) tcGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAgcagca < 1:349020/76‑1 (MQ=255) cGAAAATATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCa > 2:83665/1‑76 (MQ=255) aaaTATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAag > 1:765822/1‑76 (MQ=255) aaaTATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAag < 2:427694/76‑1 (MQ=255) aaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAag < 1:103083/76‑1 (MQ=255) aaTATGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAaga < 1:1207234/76‑1 (MQ=255) atatGGGGCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAagag < 2:1530420/76‑1 (MQ=255) gggCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCg < 2:607600/76‑1 (MQ=255) ggCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGa < 1:66139/76‑1 (MQ=255) gCTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGaa < 2:775587/76‑1 (MQ=255) cTGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAg < 2:899681/76‑1 (MQ=255) tGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGa > 2:1514545/1‑76 (MQ=255) tGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGa < 2:1121563/76‑1 (MQ=255) tGGAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGa < 2:184987/76‑1 (MQ=255) ggAAAACCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGat < 1:1374321/76‑1 (MQ=255) aCCTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCtt < 1:1228618/76‑1 (MQ=255) ccTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCttt < 2:659026/76‑1 (MQ=255) cTGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTg > 1:1463519/1‑76 (MQ=255) tGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTgg < 1:1529788/76‑1 (MQ=255) tGGCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTgg > 1:1129976/1‑76 (MQ=255) ggCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGt < 1:261605/76‑1 (MQ=255) gCTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTg < 2:1478289/76‑1 (MQ=255) cTCGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGa < 2:4274/76‑1 (MQ=255) cGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATg > 1:442896/1‑76 (MQ=255) cGTATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATg > 2:1417956/1‑76 (MQ=255) gTATGTGTAAGCAGGACATTATTTTTGCCATTCTGAAGCAGCGCGAAAAGAGTGGCGAAGATATCTTTGGTGATgg < 2:1247504/76‑1 (MQ=255) gTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATgg < 2:970987/76‑1 (MQ=255) tATGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGc < 1:383118/76‑1 (MQ=255) aTGTGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCg < 1:1422264/76‑1 (MQ=255) | AATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCG > REL606/3928230‑3928376 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |