Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 3,947,046 | G→A | 10.4% | intergenic (+137/‑112) | tyrZR → / → tyrZ | transcriptional regulator controlling expression of tyrZ/tyrosyl‑tRNA synthetase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 3,947,046 | 0 | G | A | 10.4% | 110.6 / 3.6 | 39 | intergenic (+137/‑112) | tyrZR/tyrZ | transcriptional regulator controlling expression of tyrZ/tyrosyl‑tRNA synthetase |
Reads supporting (aligned to +/- strand): ref base G (20/14); new base A (2/2); total (22/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 = 6.06e-01 |
AAGCTTTGACGGGAATCCAGTAGCGGTTTGTCTCCCTGTTTCAGAGAGCCGGCGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACCGCGTGCATTGAGCCACGTCCCTTATTGGGATGGGCTCTTTTTTGTGTTTGTAAAAAAGGAGGAAAAATGATGAGAACATTTGAGCAGCTCA > BS168_reference_genome/3946908‑3947182 | aaGCTTTGACGGGAATCCAGTAGCGGTTTGTCTCCCTGTTTCAGAGAGCCGGCGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCg > 3:100201/1‑139 (MQ=255) gACGGGAATCCAGTAGCGGTTTGTCTCCCTGTTTCAGAGAGCCGGCGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGCACCTTAGTAGCAAAGAACGGGAAACCGTTATTc > 7:43058/1‑138 (MQ=255) gACGGGAATCCAGTAGCGGTTTGTCTCCCTGTTTCAGAGAGCCGGCGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTc > 8:20837/1‑138 (MQ=255) gACGGGAATCCAGTAGCGGTTTGTCTCCCTGTTTCAGAGAGCCGGCGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTc > 7:176512/1‑138 (MQ=255) gggAATCCAGTAGCGGTTTGTCTCCCTGTTTCAGAGAGCCGGCGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAAt < 8:176512/139‑1 (MQ=255) aaTCCAGTAGCGGTTTGTCTCCCTGTTTCAGAGAGCCGGAGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGACTTACCTCCCCAATCTTTCTTTGTGAACCGTAGTAGCAAAGAACGGGAAACCGTTATTCGatagag > 8:72365/1‑134 (MQ=255) gCGGTTTGTCTCCCTGTTTCAGAGAGCCGGCGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACa < 2:134780/139‑1 (MQ=255) gCGGTTTGTCTCCCTGTTTCAGAGAGCCGGCGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACa < 6:65156/139‑1 (MQ=255) cGGTTTGTCTCCCTGTTTCTGAGTGCCGGCGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACa < 4:30192/138‑1 (MQ=255) tgctcCCGGTGTCAGAGTGCCGGCGGGAGGTGCGATCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCt < 5:136186/136‑1 (MQ=255) agagCCGGCGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACaa < 4:192874/118‑1 (MQ=255) agagCCGGCGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTATCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACaa > 3:192874/1‑118 (MQ=255) agCCGGCGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACcg > 7:54111/1‑139 (MQ=255) agCCGGCGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACcg > 7:71223/1‑139 (MQ=255) ggcggTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACCGCGTg < 8:54111/139‑1 (MQ=255) aGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGATCGGAAGACAATGTCTTAAAACAGGGTGGTACCACGTGCAtt > 2:213423/1‑137 (MQ=255) gCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACCGCGTGCATTGAGCCACGTc < 7:72365/139‑1 (MQ=255) ccGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAg > 4:160656/1‑87 (MQ=255) ccGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAg < 3:160656/87‑1 (MQ=255) cgcACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACcgc > 8:221684/1‑117 (MQ=255) cgcACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACcgc < 7:221684/117‑1 (MQ=255) aaCGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTAc > 7:116864/1‑104 (MQ=255) aaCGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTAc < 8:116864/104‑1 (MQ=255) aTTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACCGCGTGCATTGAGCCACGTCCCTTATTGGGATGGGCTCtttt > 8:159077/1‑139 (MQ=255) ccTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACCGCGTGCATTGAGCCACGTCCCTTATTGGGATGGGCTCTTTTTtgt > 8:171417/1‑139 (MQ=255) ccccAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACcgc > 7:96969/1‑93 (MQ=255) ccccAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACcgc < 8:96969/93‑1 (MQ=255) ctttGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACCGCGTGCATTGAGCCACGTCCCTTATTGGGATGGGCTCTTTTTTGTGTTTGTAAAAAag > 3:193903/1‑138 (MQ=255) ctttGTGAACCTTAGTAGCAAAGAACGGGAAACCATTATTCGAATGAGCGGAAGACAATGTCTTCAAc < 8:94264/68‑1 (MQ=255) ctttGTGAACCTTAGTAGCAAAGAACGGGAAACCATTATTCGAATGAGCGGAAGACAATGTCTTCAAc < 8:97165/68‑1 (MQ=255) ctttGTGAACCTTAGTAGCAAAGAACGGGAAACCATTATTCGAATGAGCGGAAGACAATGTCTTCAAc > 7:97165/1‑68 (MQ=255) ctttGTGAACCTTAGTAGCAAAGAACGGGAAACCATTATTCGAATGAGCGGAAGACAATGTCTTCAAc > 7:94264/1‑68 (MQ=255) ttGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACCGCGTGCATTGAGCCACGTCCCTTATTGGGATGGGCTCTTTTTTGTGTTTGTAAAAAaggag > 7:84581/1‑139 (MQ=255) tttcaaCCTTAGTAGCAATCAACTGGAAACCTTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACCGCGTGCATTGAGCCACGTCCCTTATTGGGATGGGCTCTTTTTTGTGTTTGTAAAAAaggag < 1:213423/134‑1 (MQ=255) agtagCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACCGCGTGCATTGAGCCACGTCCCTTATTGGGATGGGCTCTTTTTTGTGTTTGTAAAAAAGGAGGAAAAatgat < 8:71223/138‑1 (MQ=255) aaGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACCGCGTGCATTGAGCCACGTCCCTTATTGGGATGGGCTCTTTTTTGTGTTTGTAAAAAAGGAGGAAAAATGATGAGAACAt > 6:191149/1‑139 (MQ=255) aaGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACCGCGTGCATTGAGCCACGTCCCTTATTGGGATGGGCTCTTTTTTGTGTTTGTAAAAAAGGAGGAAAAATGATGAGAACAt > 2:111309/1‑139 (MQ=255) cGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCCACCAGGGTGGTACCGCGTGCATTGAGCCACGTCACTTATTGGGATGGGCTCTTTTTTGTGTTTGTAAAAAAGGAGGAAAAATGATGAGAACATTTGag < 2:110572/139‑1 (MQ=255) cGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACCGCGTGCATTGAGCCACGTCCCTTATTGGGATGGGCTCTTTTTTGTGTTTGTAAAAAAGGAGGAAAAATGATGAGAACATTTGAGCAGCTca > 8:127751/1‑138 (MQ=255) cGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACCGCGTGCATTGAGCCACGTCCCTTATTGGGATGGGCTCTTTTTTGTGTTTGTAAAAAAGGAGGAAAAATGATGAGAACATTTGAGCAGCTca > 2:58544/1‑138 (MQ=255) | AAGCTTTGACGGGAATCCAGTAGCGGTTTGTCTCCCTGTTTCAGAGAGCCGGCGGTAGGTGCGAGCCGGCGCACAGGCAAACGTGAATTACCTCCCCAAGCTTTCTTTGTGAACCTTAGTAGCAAAGAACGGGAAACCGTTATTCGAATGAGCGGAAGACAATGTCTTCAACCAGGGTGGTACCGCGTGCATTGAGCCACGTCCCTTATTGGGATGGGCTCTTTTTTGTGTTTGTAAAAAAGGAGGAAAAATGATGAGAACATTTGAGCAGCTCA > BS168_reference_genome/3946908‑3947182 |
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: ‑ |