Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,928,907 | A→G | T231A (ACC→GCC) | rho → | transcription termination factor Rho |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,928,907 | 0 | A | G | 98.2% | 218.1 / ‑4.2 | 55 | T231A (ACC→GCC) | rho | transcription termination factor Rho |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (30/24); total (31/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01 |
GCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTT > REL606/3928837‑3928978 | gctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCttt > 2:266813/1‑76 (MQ=255) ctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTg < 2:1468305/76‑1 (MQ=255) tgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGa < 2:50936/76‑1 (MQ=255) gATCGACGAACGTCCGGAAGAAGTAACCGGGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGAc > 2:73482/1‑76 (MQ=255) gATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGAc > 2:792287/1‑76 (MQ=255) tCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGa > 1:1517535/1‑76 (MQ=255) tCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGa > 2:60088/1‑76 (MQ=255) acgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAAcc < 1:884773/76‑1 (MQ=255) acgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAAcc > 2:1265061/1‑76 (MQ=255) gTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCAtc > 1:847721/1‑76 (MQ=255) tCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCAtct > 2:1306279/1‑76 (MQ=255) ggAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGc > 2:1054094/1‑76 (MQ=255) gaagaaGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGcc > 2:1448900/1‑76 (MQ=255) agaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCAc < 1:1181251/76‑1 (MQ=255) agaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCAc > 1:481521/1‑76 (MQ=255) aagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGt < 2:995889/76‑1 (MQ=255) aagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGt > 2:635083/1‑76 (MQ=255) gTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTc < 1:60088/76‑1 (MQ=255) gTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTc < 1:1216624/76‑1 (MQ=255) aaCCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAg > 1:260431/1‑76 (MQ=255) gagaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTg < 1:885220/76‑1 (MQ=255) gagaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTg < 1:674920/76‑1 (MQ=255) aTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCgg < 2:1579989/76‑1 (MQ=255) tGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGa < 1:108658/76‑1 (MQ=255) cAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGaaa > 1:871175/1‑76 (MQ=255) cAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGaaa < 1:363454/76‑1 (MQ=255) cAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGaaa > 1:55674/1‑76 (MQ=255) cGTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATgg > 2:940697/1‑76 (MQ=255) gTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGt < 1:253010/76‑1 (MQ=255) gTCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGt < 1:1447594/76‑1 (MQ=255) tCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTg > 2:1017938/1‑76 (MQ=255) tCTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTg < 2:1060425/76‑1 (MQ=255) cTGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTCGCGGAAATGGTGa > 1:1549486/1‑76 (MQ=255) tGGTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGAt > 2:1184400/1‑76 (MQ=255) ggTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATc < 1:1403885/76‑1 (MQ=255) gTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCg < 2:105632/76‑1 (MQ=255) gTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCg > 1:346293/1‑76 (MQ=255) tAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCga > 1:1259553/1‑76 (MQ=255) tAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCga > 1:638394/1‑76 (MQ=255) tAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCga < 1:68960/76‑1 (MQ=255) aaaaGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCgag < 2:878719/76‑1 (MQ=255) aaaaGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCgag < 1:767709/76‑1 (MQ=255) aaaaGGTGAAGTTGTTGCTTCCGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCgag > 1:735415/1‑76 (MQ=255) aaGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGaa > 1:425120/1‑76 (MQ=255) aaGGGGATGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGGTCAGGTTGCGGGAATGGTGATCGATaa > 2:789370/1‑76 (MQ=255) ggTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAgg < 1:418475/76‑1 (MQ=255) tGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGcc > 1:1000448/1‑76 (MQ=255) gttgttGCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAAc < 2:1320115/76‑1 (MQ=255) ttgttgCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCaaccc > 1:796375/1‑73 (MQ=255) ttgttgCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACg > 1:222119/1‑76 (MQ=255) tgttgCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGc > 2:1341756/1‑76 (MQ=255) ttgCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCt < 2:186524/76‑1 (MQ=255) ttgCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCt < 2:1250272/76‑1 (MQ=255) ttgCTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCt > 1:765734/1‑76 (MQ=255) cTTCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGt > 1:6068/1‑76 (MQ=255) ttCTGCCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGtt < 2:432085/76‑1 (MQ=255) | GCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTT > REL606/3928837‑3928978 |
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: ‑ |