Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 1,294,024 | A→G | 23.7% | D225G (GAT→GGT) | truB → | tRNA pseudouridine synthase B |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 1,294,024 | 0 | A | G | 23.7% | 19.4 / 5.1 | 17 | D225G (GAT→GGT) | truB | tRNA pseudouridine synthase B |
Reads supporting (aligned to +/- strand): ref base A (10/3); new base G (2/2); total (12/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.38e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01 |
GCTAGCAACTGATATTGGTGTGAAATTAGGCTTTCCGGCACATATGTCGAAATT‑AACACGAATCGAGTCTGGTGGATTTGTGTTGAAAGATAGCCTTACATTAGAACAAATAAAAGAACTTCATGAGCAGGATTCATTGCAAAATAAATTGTTTCCTTTAGAATATGGATTAAAGGGTTTGCCAAGCATTAAAATTAAAGATTCGCACATAAAAAAACGTATTTTAAATGGGCAGAAATTTAATAAAAA > USA300TCH1516_ALE/1293893‑1294141 | gCTAGCAACTGATATTGGTGTGAAATTAGGCTTTCGGACAAAAATCGCGAATTTCAACAGGAACAGGGGATTCT‑GTTTTGTGTTGAAAGATAGCCTTACATTAGAACAAATAAAAGAACTTCATGAGCAGGATTCATTGCa > 1:246026/1‑141 (MQ=255) tAGCAACTGATATTGGTGTGAAATTAGGCTTTCCGGCACATATGTCGAAATT‑AACACGAATCGAGTCTGGTGGATTTGTGGTGAAAGATAGCCTTACATTAGAACAAATAAAAGAACTTCATGAGCAGGATTCATTGCaaa > 1:81522/1‑141 (MQ=255) cTGATATTGGTGTGAAATTAGGCTTTCCGGCACATATGTCGAAATT‑AACACGAATCGAGTCTCGTGGATTTGTGTTGAAAGATAGCCTTACATTAGAACAAATAAAAGAACTTCATGAGCAGCATTCAATACCagataaat > 2:109734/1‑141 (MQ=255) tatTGGTGTGAAATTAGGCTTTCCGGCACATATGTCGAAATT‑AACACGAATCGAGTCTGGTGGATTTGTGTTGAAAGATAGCCTTACATTAGAACAAATAAAAGAACTTCATGAGCAGGATTCATTGCAAAATAAATTGtt > 2:216899/1‑141 (MQ=255) ttGGTGTGAAATTAGGCTTTCCGGCACATATGTCGAAATT‑AACACGAATCGAGTCTGGTGGATTTGTGTTGAAAGATAGCCTTACATTAGAACAAATAAAAGAACTTCATGAGCAGGATTCATTGCAAAATAAATTGTTTc > 2:49567/1‑141 (MQ=255) tttATGCTTTCCTCCATATATGTCCAAATT‑AAGACGGATCGAGTCTGGTGGTTGTATGGTGAATGATAGCCTTACATTACAACACATAAAAGAACTTCATGAGAAGGATTCATTGCAAAATAAATTGTTTCCTTTAGAata < 2:126574/140‑1 (MQ=255) aGGCTTTCCGGCACATATGTCGAAATT‑AACACGAATCGAGTCTGGTGGATTTGTGTTGAAAGATAGCCTTACATTAGAACAAATAAAAGAACTTCATGAGCAGGATTCATTGCAAAGTAAATTGTTTCCTTTAGAATTTgg > 1:71316/1‑141 (MQ=255) caTATGTCGAAATT‑AACACGAATCGAGTCTGGTGGATTTGTGTTGAAAGATAGCCTTACATTAGAACAAATAAAAGAACTTCATGAGCAGGGTTCATTGCAAAATAAATTGTTTCCTTTAGAATATGGATTAAAGGGTTTg > 1:116956/1‑141 (MQ=255) caTATGTCGAAATT‑AACACGAATCGAGTCTGGTGGATTTGTGTTGAAAGATAGCCTTACATTAGAACAAATAAAAGAACTTCATGAGCAGGGTTCATTGCAAAATAAATTGTTTCCTTTAGAATATGGATTAAAGGGTTTg > 1:116954/1‑141 (MQ=255) tCGAAATT‑AACACGAATCGAGTCTGGTGGATTTGTGTTGAAAGATAGCCTTACATTAGAACAAATAAAAGAACTTCATGAGCAGGATTCATTGCAAAATAAATTGTTTCCTTTAGAATATGGATTAAAGGGTTTGCCAAGc > 1:205717/1‑141 (MQ=255) gTCTGGTGGATTTGTGTTGAAAGATAGCCTTACATTAGAACAAATAAAAGAACTTCATGAGCAGGATTCATTGCAAAATAAATTGTTTCCTTTAGAATATGGATTAAAGGGTTTGCCAAGCATTAAAATTAAAGATTCGca > 1:236701/1‑141 (MQ=255) aaaGATAGCCTTACATTAGAACAAATAAAAGAACTTCATGAGCAGGATTCATTGCAAAATAAATTGTTTCCTTTAGAATATGGATTAAAGGGTTTGCCAAGCATTAAAGCTACATGTGCGCACATAAAAAAACGTATTTTa > 2:184828/1‑141 (MQ=255) gAACAAATAAAAGAACTTCATGAGCAGGGTTCATTGCAAAATAAATTGTTTCCTTTAGAATATGGATTAAAGGGTTTGCCAAGCATTAAAATTAAAGATTCGCACATAAAAAAACGTATTTTAGATGGGCAGAAATTtaat < 2:116956/141‑1 (MQ=255) gAACAAATAAAAGAACTTCATGAGCAGGGTTCATTGCAAAATAAATTGTTTCCTTTAGAATATGGATTAAAGGGTTTGCCAAGCATTAAAATTAAAGATTCGCACATAAAAAAACGTATTTTAAATGGGCAGAAATTtaat < 2:116954/141‑1 (MQ=255) cAAATAAAAGAACTTCATGAGCAGGATTCATTGCAAAATAAATTGTTTCCTTTAGAATATGGATTAAAGGGTTTGCCAAGCATTAAAATTAAAGATTCGCACATAAAAAAACGTATTTTAAATGGGCAGAAATTTAATaaa < 1:237558/141‑1 (MQ=255) cAAATAAAAGAACTTCATGAGCAGGATTCATTGCAAAATAAATTGTTTCCTTTAGAATATGGATTAAAGGGTTTGCCAAGCATTAAAATTAAAGATTCGCACATAAAAAAACGTATTTTAAATGGGCAGAAATTTAATaaa < 1:37542/141‑1 (MQ=255) aaTAAAAGAACTTCATGAGCAGGATTCATTGCAAAATAAATTGTTTCCTTTAGAATATGGATTAAAGGGTTTGCCAAGCATTAAAATTAAAGATTCGCACATAAAAAAACGTATTTTAAATGGGCAGAAATTTAATaaaaa > 2:21494/1‑141 (MQ=255) | GCTAGCAACTGATATTGGTGTGAAATTAGGCTTTCCGGCACATATGTCGAAATT‑AACACGAATCGAGTCTGGTGGATTTGTGTTGAAAGATAGCCTTACATTAGAACAAATAAAAGAACTTCATGAGCAGGATTCATTGCAAAATAAATTGTTTCCTTTAGAATATGGATTAAAGGGTTTGCCAAGCATTAAAATTAAAGATTCGCACATAAAAAAACGTATTTTAAATGGGCAGAAATTTAATAAAAA > USA300TCH1516_ALE/1293893‑1294141 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |