Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,928,259 | A→T | I15F (ATC→TTC) | rho → | transcription termination factor Rho |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,928,259 | 0 | A | T | 95.1% | 231.5 / 3.2 | 61 | I15F (ATC→TTC) | rho | transcription termination factor Rho |
Reads supporting (aligned to +/- strand): ref base A (1/2); new base T (23/35); total (24/37) | |||||||||||
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.05e-01 |
TTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC > REL606/3928188‑3928327 | ttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCAc > 1:1054711/1‑76 (MQ=255) tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCAct > 2:308390/1‑76 (MQ=255) tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCAct < 1:512815/76‑1 (MQ=255) cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACtctc < 2:1330166/76‑1 (MQ=255) cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACtctc < 2:1615375/76‑1 (MQ=255) cAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCgg > 1:804928/1‑76 (MQ=255) aaGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGc < 2:493606/76‑1 (MQ=255) aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCg < 2:1239620/76‑1 (MQ=255) tttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGaa < 1:354244/76‑1 (MQ=255) tAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGaaaa > 1:803855/1‑76 (MQ=255) aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAat < 1:227678/76‑1 (MQ=255) aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAata > 2:1365004/1‑76 (MQ=255) aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAata < 1:974870/76‑1 (MQ=255) taaCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAatat > 2:190630/2‑76 (MQ=255) aaCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATg > 1:386779/1‑76 (MQ=255) aaCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATg > 1:386798/1‑76 (MQ=255) tCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATgggg < 1:950473/76‑1 (MQ=255) cacaCCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGc > 2:876990/1‑76 (MQ=255) cacaCCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGc > 2:200546/1‑76 (MQ=255) cacaCCACTATGAATCTTACCGAATTAAAGAATACGCCGATTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGc > 2:527725/1‑76 (MQ=255) acacCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCt < 1:668242/76‑1 (MQ=255) caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTg < 1:1239038/76‑1 (MQ=255) ccacTATGAATCTTACCAAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGa < 1:360440/76‑1 (MQ=255) cTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGaaaa < 1:97172/76‑1 (MQ=255) tCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCt < 2:1647518/76‑1 (MQ=255) tACCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt < 1:1008864/76‑1 (MQ=255) aCCGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTa < 2:139196/76‑1 (MQ=255) ccGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt < 2:1100724/76‑1 (MQ=255) cGAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATg < 2:804927/76‑1 (MQ=255) gAATTAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGc < 1:1705251/76‑1 (MQ=255) aaTTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCg < 1:785439/76‑1 (MQ=255) ttAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa < 2:758289/76‑1 (MQ=255) ttAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa < 1:18902/76‑1 (MQ=255) ttAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa < 1:308391/76‑1 (MQ=255) tAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGtga < 2:1575077/76‑3 (MQ=255) tAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa < 1:785653/76‑1 (MQ=255) tAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa < 2:1614569/76‑1 (MQ=255) tAAAGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa > 1:807673/1‑76 (MQ=255) aaaGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg > 1:1263558/1‑76 (MQ=255) aaGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc < 1:335157/76‑1 (MQ=255) aaGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc < 2:511921/76‑1 (MQ=255) aGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa < 2:368875/76‑1 (MQ=255) aGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa > 2:1271410/1‑76 (MQ=255) aGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa < 2:368861/76‑1 (MQ=255) aGAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa > 1:1032411/1‑76 (MQ=255) gAATACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg > 1:475415/1‑76 (MQ=255) aaTACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt < 1:319212/63‑1 (MQ=255) aaTACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt > 2:319211/1‑63 (MQ=255) aaTACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAgg < 2:930296/76‑1 (MQ=255) tACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAc < 1:567306/76‑1 (MQ=255) tACGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAc > 2:681456/1‑76 (MQ=255) cGCCGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACat < 1:1374018/76‑1 (MQ=255) cGGTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat < 1:1088032/76‑1 (MQ=255) ggTTTCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattatt > 1:1012908/1‑76 (MQ=255) tttCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt < 2:1576687/76‑1 (MQ=255) tttCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt < 2:1576682/76‑1 (MQ=255) tttCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt < 2:1038447/76‑1 (MQ=255) tttCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt > 1:329900/1‑76 (MQ=255) nttCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt > 1:173989/2‑76 (MQ=255) ttCTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttttt > 1:66523/1‑76 (MQ=255) cTGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGc > 2:957769/1‑76 (MQ=255) tGAGCTGTTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGcc > 1:266267/1‑76 (MQ=255) | TTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCC > REL606/3928188‑3928327 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |