Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,928,938 | A→T | Q241L (CAG→CTG) | rho → | transcription termination factor Rho |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,928,938 | 0 | A | T | 96.2% | 203.4 / ‑1.6 | 53 | Q241L (CAG→CTG) | rho | transcription termination factor Rho |
Reads supporting (aligned to +/- strand): ref base A (0/2); new base T (25/26); total (25/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.92e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01 |
ACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTG > REL606/3928865‑3929005 | aCCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTgg < 1:1172626/76‑1 (MQ=255) aCCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTgg > 1:904905/1‑76 (MQ=255) cGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGtt > 2:995895/1‑76 (MQ=255) cGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGtt < 1:785845/76‑1 (MQ=255) gagaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTg < 2:772980/76‑1 (MQ=255) gCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATg < 1:1046498/76‑1 (MQ=255) ggTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATc > 1:845916/1‑76 (MQ=255) tAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCga > 1:1021684/1‑76 (MQ=255) aaaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCgag < 1:733845/76‑1 (MQ=255) aaaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCgag < 1:138421/76‑1 (MQ=255) aaaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCgag > 2:609259/1‑76 (MQ=255) aaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCgaga > 1:12550/1‑76 (MQ=255) aaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCgaga > 1:1233760/1‑76 (MQ=255) aaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCgaga < 2:1259812/76‑1 (MQ=255) aaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCgaga > 1:568286/1‑76 (MQ=255) aGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAg < 2:988568/76‑1 (MQ=255) aGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAg < 1:67151/76‑1 (MQ=255) tGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGcc > 2:840179/1‑76 (MQ=255) tGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGcc < 1:769626/76‑1 (MQ=255) tGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGcc > 1:749660/1‑76 (MQ=255) aaGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCaa > 2:438312/1‑76 (MQ=255) tgttgCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGc > 1:762903/1‑76 (MQ=255) tttgCTTCTCCCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGcc < 2:384986/75‑1 (MQ=255) gttgCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGcc > 1:1482122/1‑76 (MQ=255) ttgCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCt > 1:883393/1‑76 (MQ=255) tgCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTg < 2:1188103/76‑1 (MQ=255) tCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTg < 1:732608/76‑1 (MQ=255) tCTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTg < 1:413006/76‑1 (MQ=255) cTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGa < 2:873557/76‑1 (MQ=255) cTACCTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGa > 1:280386/1‑76 (MQ=255) tACCTTTGACGAACCCGCATCTCGCCCCGTTCTGGTTGCGGAAATGGCGATCGAGAAGGCCAAACGCCTGGGTGAg > 1:1342377/1‑76 (MQ=255) ccTTTGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGca > 1:193306/1‑76 (MQ=255) tttGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcaca > 1:934611/1‑76 (MQ=255) tGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaag < 1:1317240/76‑1 (MQ=255) tGACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaag < 1:1317247/76‑1 (MQ=255) ngACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaag > 1:134292/2‑76 (MQ=255) gACGAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaaga < 1:1120070/76‑1 (MQ=255) gACGAACCCGCAGCTCGCCGCGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaaga < 1:585704/76‑1 (MQ=255) gAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaag < 1:204503/76‑1 (MQ=255) gAACCCGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaag < 2:810908/76‑1 (MQ=255) ccGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGt > 2:687138/1‑76 (MQ=255) ccGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGt < 1:412176/76‑1 (MQ=255) cGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGtt < 2:495474/76‑1 (MQ=255) cGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGtt > 2:1502409/1‑76 (MQ=255) cGCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGtt < 1:65394/76‑1 (MQ=255) cGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGtt < 2:111609/76‑1 (MQ=255) gCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTa < 2:762902/76‑1 (MQ=255) gCATCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTa < 1:1236839/76‑1 (MQ=255) cATCTCGCCACGTTCTGGTTGCGGAAATGGAGATCGAGAAGGCCAAACGCCCGGTTGAGCACAAGAAAGACGTTAt > 1:1343547/1‑76 (MQ=255) aTCTCGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATc > 2:1071938/1‑76 (MQ=255) cGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTc > 2:454388/1‑76 (MQ=255) cGCCACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTc > 1:1332144/1‑76 (MQ=255) ccACGTTCTGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTg < 2:710108/76‑1 (MQ=255) | ACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTG > REL606/3928865‑3929005 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |