Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,928,928 | C→T | R238C (CGC→TGC) | rho → | transcription termination factor Rho |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,928,928 | 0 | C | T | 96.8% | 110.2 / NA | 31 | R238C (CGC→TGC) | rho | transcription termination factor Rho |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (16/14); total (16/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.84e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGAC > REL606/3928858‑3928993 | agaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCAGCATCTTGCCAc < 2:603120/76‑1 (MQ=255) agTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGACACCGCATCTTGCCACGtt > 2:815645/1‑76 (MQ=255) gTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTc > 1:175450/1‑76 (MQ=255) gTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTc > 2:55033/1‑76 (MQ=255) gTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTGTGACGAACCCGCATCTTGCCACGTTc < 1:963180/76‑1 (MQ=255) gagaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTg < 2:5203/76‑1 (MQ=255) aTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCgg > 1:82938/1‑76 (MQ=255) cAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGaaa < 2:680167/76‑1 (MQ=255) aGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAAt < 1:383106/76‑1 (MQ=255) cGTCTGGTAAAAGGTGAATTTCTTGCTTCTACCTTTGACGAACACGCATCTTGCCACGTTCAGGTTGCGGAAATgg > 2:976700/1‑76 (MQ=255) tCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTg < 2:345438/76‑1 (MQ=255) cTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAANGGTGa > 2:875727/1‑76 (MQ=255) cTGGCAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGa > 2:221183/1‑76 (MQ=255) tAAAAGGTGAAGTTCTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCga < 2:676256/76‑1 (MQ=255) aaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCgaga > 1:303696/1‑76 (MQ=255) aaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGCTTGCGGAAATGGTGATCgaga > 2:905379/1‑76 (MQ=255) aaaCGTGAAGGTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGNTCAGGTTACGGAAATGGTGATCgcga < 2:937429/76‑1 (MQ=255) aaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGaa > 2:237314/1‑76 (MQ=255) aGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAg < 1:197635/76‑1 (MQ=255) tggaGTTGTTGCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGcc > 2:631720/4‑76 (MQ=255) aGTTGTTGCTTCTACCCTTGACGAACCCGCATCGTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCaaa < 2:724902/76‑1 (MQ=255) ttgttgCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACg > 1:727513/1‑76 (MQ=255) ttgttgCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACg < 2:568218/76‑1 (MQ=255) ttgCTTCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCt < 1:1072461/76‑1 (MQ=255) tCTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTg > 2:390899/1‑76 (MQ=255) cTACCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAACAGCCTGGTTGa > 1:204982/1‑76 (MQ=255) aCCTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGc < 1:429818/76‑1 (MQ=255) ccTTTGACGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGca > 2:1058557/1‑76 (MQ=255) tttGACGAACCCACATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcaca > 1:619914/1‑76 (MQ=255) cGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaa > 1:552230/1‑76 (MQ=255) cGAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaa < 2:750836/76‑1 (MQ=255) gAACCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaag < 1:586922/76‑1 (MQ=255) aCCCGCATCTTGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGTCAAACGCCTGGTTGAGCACAAGAAAGAc > 2:492528/1‑76 (MQ=255) | AGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGAC > REL606/3928858‑3928993 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |