Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,928,260 | T→A | I15N (ATC→AAC) | rho → | transcription termination factor Rho |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,928,260 | 0 | T | A | 94.5% | 206.4 / 5.2 | 55 | I15N (ATC→AAC) | rho | transcription termination factor Rho |
Reads supporting (aligned to +/- strand): ref base T (1/2); new base A (27/25); total (28/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.11e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-02 |
TTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT > REL606/3928187‑3928332 | tttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACa < 2:1501134/76‑1 (MQ=255) ttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAacac > 1:3503/1‑76 (MQ=255) aCCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActc > 1:1443817/1‑76 (MQ=255) ccccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActct < 2:1310859/76‑1 (MQ=255) ccccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActct < 2:685654/76‑1 (MQ=255) cccAAGTTTAAGAACTCGCACCACTATGAATCTTGCCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc < 1:1731333/76‑1 (MQ=255) ccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCg > 2:466815/1‑76 (MQ=255) ccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCg < 2:473285/76‑1 (MQ=255) cAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCgg < 1:438056/76‑1 (MQ=255) aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCg > 1:1770794/1‑76 (MQ=255) tttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaa < 2:1772321/76‑1 (MQ=255) tttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaa > 1:1549524/1‑76 (MQ=255) ttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaa > 1:60219/1‑76 (MQ=255) tAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaaa > 2:1874697/1‑76 (MQ=255) gAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAatat < 2:723258/76‑1 (MQ=255) aaCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATg > 1:809123/1‑76 (MQ=255) aaCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATg < 1:630870/76‑1 (MQ=255) cTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCG‑AAATATgggg < 2:872592/76‑1 (MQ=255) tCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATgggg > 2:1955134/1‑76 (MQ=255) acacCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCt < 2:1807902/76‑1 (MQ=255) caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTg < 1:375946/76‑1 (MQ=255) accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTgg < 2:1344098/76‑1 (MQ=255) cacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaa < 1:205372/76‑1 (MQ=255) cacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaa < 1:77245/76‑1 (MQ=255) acTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaa < 1:1933328/76‑1 (MQ=255) cTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaaa > 2:1306411/1‑76 (MQ=255) tATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAc < 2:1193492/76‑1 (MQ=255) aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAcc > 2:1368485/1‑76 (MQ=255) aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAcc > 1:16274/1‑76 (MQ=255) aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAcc > 2:83070/1‑76 (MQ=255) aTCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGc > 1:681383/1‑76 (MQ=255) cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc < 2:1907849/76‑1 (MQ=255) cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc > 1:486740/1‑76 (MQ=255) ttACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCg < 2:1487282/76‑1 (MQ=255) ttACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCg > 2:1987332/1‑76 (MQ=255) tACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt < 1:1488833/76‑1 (MQ=255) ccGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt < 2:45024/76‑1 (MQ=255) aaTTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCg > 1:328154/1‑76 (MQ=255) ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa > 1:973874/1‑76 (MQ=255) tAAAGAATACGCCGGTTTCTGAGCTGAACACTCCCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa < 1:456057/76‑1 (MQ=255) aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg > 2:1665566/1‑76 (MQ=255) aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg > 2:1944470/1‑76 (MQ=255) aaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc < 1:261027/76‑1 (MQ=255) gAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg > 1:991884/1‑76 (MQ=255) gAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg < 1:328693/76‑1 (MQ=255) aaTACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAgg < 1:304032/76‑1 (MQ=255) cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTGGTATGCGTAAGCAGGACat < 1:1245537/76‑1 (MQ=255) cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACat < 1:1951183/76‑1 (MQ=255) ccGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACatta < 2:812945/76‑1 (MQ=255) cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTCAGCAGGCCattat > 2:719595/1‑76 (MQ=255) cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat > 2:372847/1‑76 (MQ=255) cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat > 2:719592/1‑76 (MQ=255) ngTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattatt > 1:345390/2‑76 (MQ=255) gTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttt < 1:1222418/76‑1 (MQ=255) gAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCa > 1:558913/1‑76 (MQ=255) aGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCAt > 1:7230/1‑76 (MQ=255) tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt > 1:970654/1‑76 (MQ=255) tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt > 2:1125389/1‑76 (MQ=255) | TTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT > REL606/3928187‑3928332 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |