Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 380,714 | G→A | 7.3% | P1249P (CCG→CCA) | srfAA → | surfactin synthetase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 380,714 | 0 | G | A | 7.3% | 171.7 / 2.6 | 55 | P1249P (CCG→CCA) | srfAA | surfactin synthetase |
Reads supporting (aligned to +/- strand): ref base G (26/24); new base A (2/2); total (28/27) | |||||||||||
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.35e-01 |
TTTATAAAGGCGAACAGCTTCCTGATCTTCGTCTCCAGTATAAGGACTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGGCGCCACTTTATATA > BS168_reference_genome/380581‑380850 | tttATAAAGGCGAACAGCTTCCTGCTCTTCGTCTCCAGTATAAGGCCTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAATGAAGTCTTTTCAGGCGAGAGACCGGTGCt > 4:67413/1‑139 (MQ=255) aGGCGAACAGCTTCCTGATCTTCGTCTCCAGTATAAGGACTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCtt < 3:235290/139‑1 (MQ=255) gAACAGCTTCCTGATCTTCGTCTCCAGTATAAGGACTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGt > 1:186618/1‑139 (MQ=255) aGCTTCCTGATCTTCGTCTCCAGTATAAGGACTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGa < 3:237606/139‑1 (MQ=255) tataAGGACTACGCTGTATTGCAAAGCAGACAGGCTGCTGAAGGGTACAAGACGGACCAGGCTTATTTGAAGGAAGTCTTTGCAGGCGAGCTCCCTGTGCTTCAGCTTCTGTCCTATTACCCAAGACCACCTGTTCat < 5:224922/138‑2 (MQ=255) ggACTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAACGCttt > 2:126098/1‑139 (MQ=255) agtACGCTGTATGGCAAAGCAGTCAGGCTGCTGTAGGGTACATGTAGGAACTTGCTTATTGGAAGGCTGTCTTTGCAGGCGAGCTGCCGGTTCGTCAGCGGCTGTCCGATTACCGAATACCAGCTGTTCAAAGCTTTGa < 8:78890/137‑1 (MQ=255) aCTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCGTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCCGCGTCTGGCCGAGAACCCAAGCCCCCCTGGTCAAAGCTTTg > 4:96864/1‑138 (MQ=255) aCGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGGAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTc < 4:163542/93‑1 (MQ=255) aCGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGGAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTc > 3:163542/1‑93 (MQ=255) aCGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAgg > 6:224922/1‑139 (MQ=255) aCGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCAGTGCtt > 3:238280/1‑92 (MQ=255) aCGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCAGTGCtt > 3:234142/1‑92 (MQ=255) aCGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCAGTGCtt < 4:234142/92‑1 (MQ=255) aCGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCAGTGCtt < 4:238280/92‑1 (MQ=255) tATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGCCCACCTGTTCAAAGCTTTGAAGGGGATCg > 2:94252/1‑139 (MQ=255) tATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGCAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTAAGCTTCTGTCACATTCAGCAACACCACCTGTTCAAAGATTTGAAGGGGATcc > 7:78890/1‑138 (MQ=255) cAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGgtgt < 3:34964/139‑1 (MQ=255) aGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGtcc > 7:159049/1‑137 (MQ=255) aGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGtcaa > 2:221525/1‑139 (MQ=255) cAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGtcaat > 2:28257/1‑138 (MQ=255) ctgctgAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGa > 3:5472/1‑139 (MQ=255) gctgAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATg > 1:96116/1‑139 (MQ=255) ctgAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGc < 1:132907/139‑1 (MQ=255) aaGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCg < 1:221525/137‑1 (MQ=255) ggTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCggggg > 3:42891/1‑138 (MQ=255) ggTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACAACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGTATTCggggg > 5:225939/1‑138 (MQ=255) ggTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACAACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCggggg > 7:282869/1‑138 (MQ=255) tACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCg < 4:10071/132‑1 (MQ=255) tACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCg > 3:10071/1‑132 (MQ=255) gaagGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGAt < 2:96116/139‑1 (MQ=255) ggACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAGAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGc > 2:15448/1‑139 (MQ=255) aCCAGGCTTATTGGAATGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGttaa < 6:237502/82‑3 (MQ=255) aCCAGGCTTATTGGAATGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCa > 5:237502/1‑82 (MQ=255) aCCAGGCTTATTGGAAGGAAATCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCt > 7:79254/1‑139 (MQ=255) aCCAGGCTTATAGGAAGGAAGTCTTTGCAGCCGACCTCCCGGTGCTTCAGCTACAGTCCCAATACCCAAGACCACCTGTTAAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATCCGGGAGTAAAGTATCGcat > 8:35459/1‑137 (MQ=255) aGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCaa > 8:259238/1‑139 (MQ=255) aggaagTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGaccac < 2:134328/60‑1 (MQ=255) aggaagTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGaccac > 1:134328/1‑60 (MQ=255) aggaagTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGa > 3:205195/1‑76 (MQ=255) aggaagTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGa < 4:205195/76‑1 (MQ=255) aggaagTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTgg > 5:243219/1‑135 (MQ=255) aggaagTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTgg < 6:243219/135‑1 (MQ=255) gCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGgcg < 1:126098/139‑1 (MQ=255) ccggCGAGCGCCCGGGGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGgcgc < 1:245641/137‑1 (MQ=255) cAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTGCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGgcg < 6:121027/138‑1 (MQ=255) aGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTc < 7:236606/112‑1 (MQ=255) aGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTc > 8:236606/1‑112 (MQ=255) cgcgCTCCCGGGGCTTCAGCTGCTGTCCGATAACCCAAGACCACCTGTTCAACGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGGGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGGCGCCACt < 6:206360/136‑1 (MQ=255) cGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGGGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGGCGCCACt < 3:80143/139‑1 (MQ=255) cGAGCTCCCGGTGCTTCAGCTGCGGTCCGATTACCCAAGACCACCTGGTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGAGCGCCTCAATCGTTTGGCTGAACAAAACGGCGCCACt < 6:225939/139‑1 (MQ=255) cccccGGGGCGTCAGCGTCTGCCCGATACCCCACCACCACCTGTTCAAAGCTTTGACGGGGATCGGGTGTCAATCACGCGGGAGGCGGGGGTAACGGCTCGCCGCAATCGTTTGGCTCAAAAAAACGGCGCCACTTtat < 4:42891/137‑1 (MQ=255) ccGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAACGCTTTGAAGGGGATCGGGTGTCAATCAAGCGGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGGCGCCACTTtata < 4:84581/137‑1 (MQ=255) ccGGTGCTGCAGCTTCTGACCGATTACCCAACACCACCTGTTCAAAGCTTTGAAGGGGAGCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGGCGCCACTTtata < 4:5472/137‑1 (MQ=255) ccGGTGCTCCAGCTTCTGTCCGATTACCCACGACCACCTGTGCAAAGCTTTGAAGGGCATCGGGGGCCAAAAAAGCTGGATGGGGGGGTAACGGAGGGCCACAATCGTTTGCCCCAACAAAACGGCGCCACTTtatata < 8:159049/139‑1 (MQ=255) | TTTATAAAGGCGAACAGCTTCCTGATCTTCGTCTCCAGTATAAGGACTACGCTGTATGGCAAAGCAGACAGGCTGCTGAAGGGTACAAGAAGGACCAGGCTTATTGGAAGGAAGTCTTTGCAGGCGAGCTCCCGGTGCTTCAGCTTCTGTCCGATTACCCAAGACCACCTGTTCAAAGCTTTGAAGGGGATCGGGTGTCAATCAAGCTGGATGCGGGGGTAAAGGATCGCCTCAATCGTTTGGCTGAACAAAACGGCGCCACTTTATATA > BS168_reference_genome/380581‑380850 |
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: ‑ |