Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 2,401,489 | T→A | 11.5% | D157V (GAT→GTT) | ypdA ← | putative bacillithiol biosynthesis thiol disulfide oxidoreductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 2,401,489 | 0 | T | A | 11.5% | 85.0 / 5.5 | 35 | D157V (GAT→GTT) | ypdA | putative bacillithiol biosynthesis thiol disulfide oxidoreductase |
Reads supporting (aligned to +/- strand): ref base T (13/18); new base A (2/2); total (15/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.88e-01 |
CCGGCAGAATCCACGGCTTAATGCTCGGAGAATACTCATTGCCCCTATATAAAACTGTGACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGCCCATATAGTTTGGATGATCATAATAGCCTGTGGCGATGATACAATATGGCGTCGTGTAGGTCTCTTTTGACGT > BS168_reference_genome/2401352‑2401625 | ccGGCAGAATCCACGGCTTAAAGCTCGGAGAATACTCATTGCCCCTATATAAAACTGTGACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGACGATGACAACAACATCTTTAt < 8:128138/138‑1 (MQ=255) gCTTAATGCTCGGAGAATACTCATTGCCCCTATATAAAACTGTGACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGcc < 4:59216/139‑1 (MQ=255) cGGAGAATACTCATTGCCCCTATATAAAACTGTGACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTaaaa < 7:217323/139‑1 (MQ=255) ggAGAATACTCATTGCCCCTATATAAAACTGTGACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTaaaa > 8:217323/1‑138 (MQ=255) ggAGAATACTCATTGCCCCTATATAAAACTGTGACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTaaaa > 6:102338/1‑138 (MQ=255) ggAGAATACTCATTGCCCCTATATAAAACTGTGACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTaaaa < 5:102338/138‑1 (MQ=255) cATTGCCCCTATATAAAACTGTGACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATAc > 6:84765/1‑139 (MQ=255) tatataAAACTGTGACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCa > 8:47987/1‑138 (MQ=255) tgACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCCCGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCggg > 2:34985/1‑139 (MQ=255) tgACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCCCGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTAACAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCggg > 3:43006/1‑139 (MQ=255) tgACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCggg > 1:64544/1‑139 (MQ=255) tgACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCggg < 6:81535/139‑1 (MQ=255) tgACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCggg > 5:81535/1‑139 (MQ=255) tgACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCggg > 5:190928/1‑139 (MQ=255) tgACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCggg < 4:80413/139‑1 (MQ=255) tgACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCggg > 3:80413/1‑139 (MQ=255) tgACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTAACAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCggg > 3:46075/1‑139 (MQ=255) tgACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTACCGCCGATGACAACAACATCTTTATCAAAACACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCggg > 4:135935/1‑139 (MQ=255) gacagcgcCCCAGATTTCACGAGTTCAAGGGGGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGGGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTa < 7:67492/135‑1 (MQ=255) acacGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTa < 3:81981/139‑1 (MQ=255) cgcgCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGc > 3:18558/1‑139 (MQ=255) cgcgCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGAGGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGc < 2:78202/139‑1 (MQ=255) cccAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGCCCata < 4:18558/139‑1 (MQ=255) tCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATAATTATCAAAATACCGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCATCACCCGGTACGCCCATATAgctgt > 8:67492/1‑134 (MQ=255) cACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGCCCATATAGTTTgg < 7:47987/138‑1 (MQ=255) aCGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTGCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCGTCTGTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGGACGCCCATATAGTTTGgat > 8:91621/1‑139 (MQ=255) aCGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTAACAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGCCCATATAGTTTGgat < 4:46075/139‑1 (MQ=255) aCGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTAACAAAATACGGGGGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGCCCATATAGTTTGgat < 4:43006/139‑1 (MQ=255) aCGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTACCGCCGATGACAACAACATCTTTATCAAAACACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGCCCATATAGTTTGgat < 3:135935/139‑1 (MQ=255) gAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGCCCATATAGTTTGgatga < 1:8535/139‑1 (MQ=255) agatCAAGGGCGGCGTCCACGAGTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGCCCATATAGTTTGgatgat < 4:193399/136‑1 (MQ=255) gAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGCCCATATAGTTTGGATGATCATAATAGCCTGTGGCGATGATa < 5:84765/138‑1 (MQ=255) tttCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGCCCATATAGTTTGGATGATCATAATAGCCTGTGGCGATGATACAata < 6:190928/138‑1 (MQ=255) ccgccgATGACAACAACATCTTTATCAAACTACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGCCCATATAGTTTGGATGATCATAATAGCCTGTGGCGATGATACAATATGGcg < 4:18568/139‑1 (MQ=255) aTCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGCCCATATAGTTTGGATGATCATAATAGCCTGTGGCGATGATACAATATGGCGTCGTGTAGGTCTCTTTTGACGt < 2:144736/138‑1 (MQ=255) aTCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGCCCATATAGTTTGGATGATCATAATAGCCTGTGGCGATGATACAATATGGCGTCGTGTAGGTCTCTTTTGACGt < 1:34985/138‑1 (MQ=255) | CCGGCAGAATCCACGGCTTAATGCTCGGAGAATACTCATTGCCCCTATATAAAACTGTGACACGCGCCCCAGATTTCACGAGTTCAAGGGCGGCGTCCACGCTTGAATTTTTCCCGCCGATGACAACAACATCTTTATCAAAATACGGGTGGCCTTCTTTAAAATAATGAAATACTTTCGGCAAATCCTCACCGGGTACGCCCATATAGTTTGGATGATCATAATAGCCTGTGGCGATGATACAATATGGCGTCGTGTAGGTCTCTTTTGACGT > BS168_reference_genome/2401352‑2401625 |
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: ‑ |