| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | BS168_reference_genome | 2,912,810 | T→G | 15.4% | T27P (ACT→CCT) | uvrC ← | excinuclease ABC (subunit C) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | BS168_reference_genome | 2,912,810 | 0 | T | G | 15.4% | 55.1 / 2.5 | 30 | T27P (ACT→CCT) | uvrC | excinuclease ABC (subunit C) |
| Reads supporting (aligned to +/- strand): ref base T (11/14); new base G (2/3); total (13/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGACGTCACAATATATTCAAAATCCTCGATTTCCGTCACAAGCCTTTGGGTTTTCGCGTCGTGAGAACCGGTGAAGTAGGAGCGCACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAGTTGTTTGTTCATATGTTACCCTTCCTTATGTTTAATCCTTTATTTCCAGTTTAACACGAACGTTTGTTTTT > BS168_reference_genome/2912682‑2912947 | ggaCGTCACAATATATTCAAAATCCTCGATTTCCGTCACAAGCCTTTGGGTTTTCGCGTCGTGAGAACCGGTGAAGTAGGAGCGCACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCgg < 2:219229/139‑1 (MQ=255) acttCACAAAATATTCAAAATCCTCGATTTCCGTCACAAGCCTGTGGGTTTTCGCGGCGTGAGAAGCGGTGAAGTAGGAGCGCACTCTGTTTTTCAGAACTTTTGCTTTCACTACGTAGATCACAGGCTGCTGCCGGTc < 8:48872/136‑1 (MQ=255) gTCACAATATATTCAAAATCCTCGATTTCCGTCACAAGCCTTTGGGTTTTCGCGTCGTGAGAACCGGTGAAGTAGGAGCGCACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGGct > 5:88789/1‑127 (MQ=255) gTCACAATATATTCAAAATCCTCGATTTCCGTCACAAGCCTTTGGGTTTTCGCGTCGTGAGAACCGGTGAAGTAGGAGCGCACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGGct < 6:88789/127‑1 (MQ=255) atTCAAAATCCTCGATTTCCGTCACAAGCCTTTGGGTTTTCGCGTCGTGAGAACCGGTGAAGTAGGAGCGCACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACATTCTGCTGCCCGTCCTGCATGAGAt > 2:71400/1‑138 (MQ=255) atTCAAAATCCTCGATTTCCGTCACAAGCCTTTGGGTTTTCGCGTCGTGAGAACCGGTGAAGTAGGAGCGCACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGAt > 6:230546/1‑138 (MQ=255) tCAAAATCCTCGATTTCCGTCACAAGCCTTTGGGTTTTCGCGTCGTGAGCACCGGTGAAGTAGGCGCGCACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGGCTGCTGCCGGGGCTTCATGAGATaac > 2:209160/1‑139 (MQ=255) tCAAAATCCTCGATTTCCGTCACAAGCCTTTGGGTTTTCGCGTCGTGAGAACCGGTGAAGTAGGAGCGCACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATaa > 4:216177/1‑138 (MQ=255) tCAAAATCCTCGATTTCCGTCACAAGCCTTTGGGTTTTCGCGTCGTGAGAACCGGTGAAGTAGGAGCGCACTCTGTTTTTCAGCACTTTTGCTTTCCCTAAGGAGATCACAGTCTGCTGCCGGTCCTTCAGGAGATaac > 6:16390/1‑139 (MQ=255) aaGCCTTTGGGTTTTCGCGTCGTGAGAACCGCTGAAGTAGGAGCGCACTCTGTTTTTCAGCACTTTTTCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCCGTTGATCAGGAaggagg > 6:54065/1‑139 (MQ=255) gtGAGAACCGGTGAAGTAGGAGCGCACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAgttgt < 4:19339/139‑1 (MQ=255) gtGAGAACCGGTGAAGTAGGAGCGCACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAgttgt < 1:124192/139‑1 (MQ=255) agaACCGGTGAAGTAGGAGCGCACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAGttgtttg < 4:114039/139‑1 (MQ=255) accgtggAAGTAGGAGCGCACTCTGTTTTTCAGCACTATTGCTTTCCCTACGTAGATCACAGGCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAGTTGTTTGTTc < 8:152439/133‑1 (MQ=255) tAGGAGCGCACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAGTTGTTTGTTCATATGTTAcc < 3:221455/139‑1 (MQ=255) gcgcACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTTCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGATGGCGAGTTTTTCTTTCAGTGGTTTGTTCATATGTTACccttcc > 6:12510/1‑139 (MQ=255) gcgcACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAGTTGTTTGTTCATATGTTAcc > 3:34389/1‑134 (MQ=255) gcgcACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAGTTGTTTGTTCATATGTTAcc < 4:34389/134‑1 (MQ=255) gcACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAGTTGTTTGTTCATATGTTACccttcctt < 5:230546/139‑1 (MQ=255) ctctGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGtt < 4:259990/78‑1 (MQ=255) ctctGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGtt > 3:259990/1‑78 (MQ=255) ctctGTTTTTCAGCACGTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACACCCCGGTTGAGCAGGAAGGAGGGCGAGTTTTTCTTTCTGTTGTTTGTTCATATGTTACcgttccgtttg > 1:79945/1‑134 (MQ=255) cTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACACCCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAgttgt > 1:5510/1‑98 (MQ=255) cTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAgttgt < 2:5510/98‑1 (MQ=255) ttGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAGTTGTTTGTTCATATGTTACCCTTCCTTATGTTTAATCCTTTATTTc < 5:16390/137‑1 (MQ=255) cGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAGTTGTTTGTTCATATGTTACCCTTCCTTATGTTTAATCCTTTATTTCCAGTTTAACACg < 1:71400/137‑1 (MQ=255) cGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAGTTGTTTGTTCATATGTTACCCTTCCTTATGTTTAATCCTTTATTTCCAGTTTAACACGaa > 8:272014/1‑139 (MQ=255) acaGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAGTTGTTTGTTCATATGTTACCCTTCCTTATGTTTAATCCTTTATTTCCAGTTTAACACGAACgtttgtt < 7:272014/139‑1 (MQ=255) acaGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAGTTGTTTGTTCATATGTTACCCTTCCTTATGTTTAATCCTTTATTTCCAGTTTAACACGAACgtttgtt < 1:209160/139‑1 (MQ=255) gTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAGTTGTTTGTTCATATGTTACCCTTCCTTATGTTTAATCCTTTATTTCCAGTTTAACACGAACGTTTGttttt < 5:54065/139‑1 (MQ=255) | GGACGTCACAATATATTCAAAATCCTCGATTTCCGTCACAAGCCTTTGGGTTTTCGCGTCGTGAGAACCGGTGAAGTAGGAGCGCACTCTGTTTTTCAGCACTTTTGCTTTCCCTACGTAGATCACAGTCTGCTGCCGGTCCTTCATGAGATAACAACCCGGTTGATCAGGAAGGAGGGCGAGTTTTTCTTTCAGTTGTTTGTTCATATGTTACCCTTCCTTATGTTTAATCCTTTATTTCCAGTTTAACACGAACGTTTGTTTTT > BS168_reference_genome/2912682‑2912947 |
| 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: ‑ |