Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 1,815,227 | T→C | 9.3% | V2436A (GTT→GCT) | pksL → | polyketide synthase of type I (bacillaene synthesis) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 1,815,227 | 0 | T | C | 9.3% | 119.0 / 4.7 | 43 | V2436A (GTT→GCT) | pksL | polyketide synthase of type I (bacillaene synthesis) |
Reads supporting (aligned to +/- strand): ref base T (20/19); new base C (2/2); total (22/21) | |||||||||||
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.61e-01 |
AGACTTCAGAAACGCGCGGAGCGATTGCTGCATGCCATTCGTGAACAAACATATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTCTATCGCGGACAAGTAAAGAGCAACAAAGAAGCCATTGCTGCTTTTGCTGCTGAT > BS168_reference_genome/1815097‑1815354 | aGACTTCAGAAACGCGCGGAGCGATTGCTGCATGCCATTCGTGAACAAACATATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGaaa < 7:222978/139‑1 (MQ=255) gAAACGCGCGGAGCGATTGCTGCATGCCATTCGTGAACAAACATATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCaa < 8:22723/139‑1 (MQ=255) aaaGGCGCGGAGCGATTGCTGCATGCCATTGGTTAACAAACTTATGTTTAAGCATATTTACACCGAATTGCCTATACGGTGCATGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAag < 2:187324/139‑1 (MQ=255) cgcgGAGCGATTGCTGCATGCCATTCGTGAACAAACATATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGCAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGACACGATGCAAGAGCtt > 8:21403/1‑139 (MQ=255) cgcgGAGCGATTGCTGCATGCCATTCGTGAACAAACATATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAACCGATGCAAGAGCtt > 6:4157/1‑139 (MQ=255) cgcgGAGCGATTGCTGCATGCCATTCGTGAACAAACATATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCtt > 2:31016/1‑139 (MQ=255) cGATTGCTGCATGCCATTCGTGAACAAACATATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGaaa < 3:180007/139‑1 (MQ=255) gcATGCCATTCGTGAACAAACATATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTAt < 6:31402/139‑1 (MQ=255) gcATGCCATTCGTGAACAAACATATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTAt > 5:201634/1‑139 (MQ=255) cATGCCATTCGTGAACAAACATATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTAt < 3:87133/138‑1 (MQ=255) tGCCATTCGTGAACAAACATATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGaa > 8:36911/1‑139 (MQ=255) ttCGTGAACAAACATATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGa > 5:157249/1‑100 (MQ=255) ttCGTGAACAAACATATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGa < 6:157249/100‑1 (MQ=255) aacaaacaTATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCg > 5:171112/1‑138 (MQ=255) aacaaacaTATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCg > 5:51144/1‑138 (MQ=255) aaGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTg > 3:184702/1‑138 (MQ=255) aGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGaa < 6:171112/139‑1 (MQ=255) aGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGaa < 5:45836/139‑1 (MQ=255) gCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGaaa < 7:63450/79‑1 (MQ=255) gCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGaaa > 8:63450/1‑79 (MQ=255) cacCGAATTGCGTATACTCTGCAAGTAGGCCTCGAAGAGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTCTAt < 8:81574/138‑1 (MQ=255) cacCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTCTAt < 7:36911/138‑1 (MQ=255) aaTTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAAGGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGaa > 2:159001/1‑124 (MQ=255) aaTTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGaa < 1:159001/124‑1 (MQ=255) aaTTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTCTATCGCGGa > 1:7037/1‑139 (MQ=255) aaTTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTCTATCGCGGa > 6:129762/1‑139 (MQ=255) aaTTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGCTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCgg > 8:104635/1‑108 (MQ=255) aaTTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGCTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCgg < 7:101668/108‑1 (MQ=255) aaTTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGCTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCgg < 7:104635/108‑1 (MQ=255) aaTTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGCTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCgg > 8:101668/1‑108 (MQ=255) aTTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTCTATCGCGGAc < 7:139149/139‑1 (MQ=255) tGCAAGTCGGCCGGGAAGGGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTCTATCGCGGACAAGTAAAGAGcaa > 4:139872/1‑139 (MQ=255) gCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGaa < 8:118074/109‑1 (MQ=255) gCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGaa > 7:118074/1‑109 (MQ=255) gCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTCTATCGCGGACAAGTAAAGAGcaa > 7:180462/1‑138 (MQ=255) gCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTCTATCGCGGACAAGTAAAGAGcaa > 6:88852/1‑138 (MQ=255) cAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTCTATCGCGGACAAGTAAAGAGcaac > 4:188995/1‑138 (MQ=255) cAAGTCGGCCGGGAAGCGAAGAAAGAACGGCTGGCGTTCGTTGCTGAAACGATGCAAGAGCTTGAAGAAACGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTATATCGCGGACAAGTAAAGAGcaac > 5:223999/1‑138 (MQ=255) tCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTCTATCGCGGACAAGTAAAGAGCAACAaaga < 5:88852/139‑1 (MQ=255) gATGAAAGAACGGCTGGCGTTCGTTGAGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGaa < 5:128908/82‑1 (MQ=255) gATGAAAGAACGGCTGGCGTTCGTTGAGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGaa > 6:128908/1‑82 (MQ=255) tGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTCTATCGCGGACAAGTAAAGAGCAACAAAGAAGCCATTGCTGCttt < 6:51144/139‑1 (MQ=255) gCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTCTATCGCGGACAAGTAAAGAGCAACAAAGAAGCCATTGCTGCTTTTGCTGCtgat < 6:87825/139‑1 (MQ=255) | AGACTTCAGAAACGCGCGGAGCGATTGCTGCATGCCATTCGTGAACAAACATATGTTGAAGCAGATTTACACCGAATTGCCTATACGCTGCAAGTCGGCCGGGAAGCGATGAAAGAACGGCTGGCGTTCGTTGCGGAAACGATGCAAGAGCTTGAAGAAAAGCTGTATGAATGTATCTCCGGAACAGAAAATCGTGAATATGTCTATCGCGGACAAGTAAAGAGCAACAAAGAAGCCATTGCTGCTTTTGCTGCTGAT > BS168_reference_genome/1815097‑1815354 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |