Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 2,794,537 | A→T | 13.8% | Y180* (TAT→TAA) | yrrN ← | putative hydrolase small subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 2,794,537 | 0 | A | T | 13.8% | 67.4 / 5.6 | 29 | Y180* (TAT→TAA) | yrrN | putative hydrolase small subunit |
Reads supporting (aligned to +/- strand): ref base A (15/10); new base T (2/2); total (17/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01 |
ATCGTTCGGGCTCATGATGTGAGTGCCGTTTTCATCTTCGAAAATCGGATATTTATTATCGCGTTCTTTATCGTGTAAAAACATGCCTGATTCTTTTTTCTTTCTTTCAATGTCCATGACTTTGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAATCTCAACAATGCTGTCCATGTTTAATTCTCTGGCAAGCACAGAACG > BS168_reference_genome/2794399‑2794671 | aTCGTTCGGGCTCATGATGTGAGTGCCGTTTTCATCTTCGAAAATCGGATATTTATTATCGCGTTCTTTATCGTGTAAAAACATGCCTGATTCTTTTTTCTTTCTTTCAATGTCCATGACTTTGCCCTGATATTCaaaa > 2:9944/1‑139 (MQ=255) gTTCGGGCTCATGATGTGAGTGCCGTTTTCATCTTCGAAAATCGGATATTTATTATCGCGTTCTTTATCGTGTAAAAACATGCCTGATTCTTTTTTCTTTCTTTCAATGTCCATGACTTTGCCCTGATATTCAAAATAg > 8:77495/1‑139 (MQ=255) tgatgTGAGTGCCGTTTTAATCTTAGAAAATAGGAAATTTATAATAGCGTTCTTTATCGTGTAAAAACATGCCTGATTCTTTTTTCTTTATTTCAATTTACATGACTTTGCCCTGATAATGAAAATATTTTAAAATCa > 6:217782/1‑138 (MQ=255) tgtgAGTGCCGTTTTCATCTTCGAAAATCGGATATTTATTATCGCGTTCTTTATCGTGTAAAAACATGCCTGATTCTTTTTTCTTTCTTTCAATGTCCATGACTTTGCCCTGATATTCAAAATAGTTTCCAATCAAAGa < 1:9944/139‑1 (MQ=255) gtgAGTGCCGTTTTCATCTTCGAAAATCGGATATTTATTATCGGGTTCTTTATCGTGTAAAAACATGCCTGATTCTTTTTTCTTTCTCTCAATGTCCATGACTTTGCCCTGATATTCAAAATAGTTTCCAATCAAAGa < 5:100699/138‑1 (MQ=255) tgacGTTTTCACCTTCGAAAAGCGGATATTTATTATCGCGTTCTTTATCGTGTAAAAACAGGCCTCACTCTTTTTTCTTACTTTCAATGTCCATGACTGTCCCCTGATATACAAAATAGTTTCCAATCAAACACCGCtt < 6:141812/136‑1 (MQ=255) tGCCGTTTTCATCTTCGAAAATCGGATATTTATTATCGCGTTCTTTATCGTGTAAAAACATGCCTGATTCTTTTTTCTTTCTTTCAATGTCCATGACTTTGCCCTGATATTCAAATTAGTTTCCAATCAAAGAACGCtt > 1:8596/1‑139 (MQ=255) tGCCGTTTTCATCTTCGAAAATCGGATATTTATTATCGCGTTCTTTATCGTGTAAAAACATGCCTGATTCTTTTTTCTTTCTTTCAATGTCCATGACTTTGCCCTGATATTCAAATTAGTTTCCAATCAAAGAACGCtt > 1:12041/1‑139 (MQ=255) gTTTTCATCTTCGAAAATCGGATATTTATTATCTCGTTCTTTTTCGTGTAAAAACATGCCTGATTCTTTTTTCTTTCTTTCAATGTCCATGACTTTGCCCTGATATTCAAAATAGTTTCCAATCAACGAACGCTTCGAt < 5:58728/139‑1 (MQ=255) gAAAATCGGATATTTATTATCGCGTTCTTTATCGTGTAAAAACATGCCTGATTCTTTTTTCTTTCTTTCAATGTCCATGACTTTGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCa < 4:197817/138‑1 (MQ=255) tATCGTGTAAAAACATGCCTGATTCTTTTTTCTTTCTTTCAATGTCCATGACTTTGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAAtt < 7:77495/137‑1 (MQ=255) tgtAAAAACATGCCTGATTCTTTTTTCTTTCTTTCAATGTCCATGACTTTGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCTCGTCATTCCATGCACCTGGATTTCAATTTCGACtt > 6:56155/1‑139 (MQ=255) tAAAAACATGCCTGATTCTTTTTTCTTTCTTTCAATGTCCATGACTTTGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACtt > 8:143102/1‑137 (MQ=255) tacttGCCTGATTCTTTTTTATTTCTTTCATTGTCCTTGACTTTGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACTTAATTCCATGCACCTGGATTTCAATTTCGACTTCAGCAt < 4:213606/135‑1 (MQ=255) cTGATTCTTTTTTCTTTCTTTCAATGTCCATGACTTTGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCttt < 6:81432/139‑1 (MQ=255) cATGACTTTGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAATCTCAACAATGCTGTCCATGTTTAAt > 8:31560/1‑139 (MQ=255) cATGACTTTGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAATCTCAACAATGCTGTCCATGTTTAAt > 1:151620/1‑139 (MQ=255) tGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAAt > 5:75048/1‑106 (MQ=255) tGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAAt < 6:75048/106‑1 (MQ=255) tGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAATCTCAACAATg > 6:139129/1‑116 (MQ=255) tGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAATCTCAACAATg < 5:139129/116‑1 (MQ=255) tGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAATCTCAACAATGCTGTCCATGTTTAATTCTCTgg > 7:12246/1‑138 (MQ=255) tGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAATCTCAACAATGCTGTCCATGTTTAATTCTCTgg > 7:26235/1‑138 (MQ=255) tGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAATCTCAACAATGCTGTCCATGTTTAATTCTCTgg > 7:91286/1‑138 (MQ=255) tGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAATCTCAACAATGCTGTCCATGTTTAATTCTCTgg > 8:92671/1‑138 (MQ=255) tGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAATCTCAACAATGCTGGCCATGTTTAATTCTCTgg > 6:49348/1‑138 (MQ=255) gCCCTGATATTCAAATTAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAATCTCAACAATGCTGTCCATGTTTAATTCTCTGGc < 2:8596/138‑1 (MQ=255) gCCCTGATATTCAAATTAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAATCTCAACAATGCTGTCCATGTTTAATTCTCTGGc < 2:12041/138‑1 (MQ=255) atatTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAATCTCAACAATGCTGTCCATGTTTAATTCTCTGGAAAGca > 6:121492/1‑137 (MQ=255) cAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCAGTTTCTTTAATCTCAACAATGCTGTCCATGTTTACTTCTCTGGCAAGCACAGAAcg > 2:215100/1‑139 (MQ=255) | ATCGTTCGGGCTCATGATGTGAGTGCCGTTTTCATCTTCGAAAATCGGATATTTATTATCGCGTTCTTTATCGTGTAAAAACATGCCTGATTCTTTTTTCTTTCTTTCAATGTCCATGACTTTGCCCTGATATTCAAAATAGTTTCCAATCAAAGAACGCTTCGATTGGAACATGCACGTCATTCCATGCACCTGGATTTCAATTTCGACTTCAGCATTTTCTTTAATCTCAACAATGCTGTCCATGTTTAATTCTCTGGCAAGCACAGAACG > BS168_reference_genome/2794399‑2794671 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |