Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 2,753,394 | T→A | 5.9% | N110K (AAT→AAA) | yraF → | putative spore coat protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 2,753,394 | 0 | T | A | 5.9% | 222.3 / 4.0 | 68 | N110K (AAT→AAA) | yraF | putative spore coat protein |
Reads supporting (aligned to +/- strand): ref base T (37/27); new base A (2/2); total (39/29) | |||||||||||
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.97e-01 |
AGAGATTACGGAACTGATGATCAGCAAAAAGTGGTTCCATCCATATGAGCTGAGCGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAGGCATTAAAAATGTTGTTGTCAAAGATGTCCCCGATCCAAAGATTGAAAAATCC > BS168_reference_genome/2753259‑2753526 | agagATTACGGAACTGATGATCAGCAAAAAGTGGTTCCATCCATATGAGCTGAGCGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCg > 2:19581/1‑138 (MQ=255) aTTACGGAACTGATGATCAGCAAAAAGTGGTTCCATCCATATGAGCTGAGCGAACAGTATAAGCTTGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAg < 1:19581/139‑1 (MQ=255) atgatCAGCAAAAAGTGGTTCCATCCATATGAGCTGAGCGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGAcc > 4:189209/1‑139 (MQ=255) atgatCAGCAAAAAGTGGTTCCATCCATATGAGCTGAGCGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGAcc > 6:1240/1‑139 (MQ=255) gatCAGCATAAAGTGGTGCCATCCATATGAGCTGAGCGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCg < 4:45094/138‑1 (MQ=255) aGCAAAAAGTGGTTCCATCCATATGAGCTGAGCGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGacac < 5:27009/139‑1 (MQ=255) aaGTGGTTCCATCCATATGAGCTGAGCGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCtgatg < 5:1240/139‑1 (MQ=255) attGGTTCCATCCATATGAGCTGAGCTAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATCAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCtgatga < 1:219922/137‑1 (MQ=255) ttCCATCCATATGAGCTGAGCGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGaacac > 5:50797/1‑138 (MQ=255) tccatATGAGCTGAGCGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCa < 1:162292/101‑1 (MQ=255) tccatATGAGCTGAGCGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCa > 2:162292/1‑101 (MQ=255) atatGAGCTGAGCGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTgg < 3:189209/139‑1 (MQ=255) atatGAGCTGAGCGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATAATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTgg < 6:50797/139‑1 (MQ=255) agcGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCtgatg > 8:108584/1‑115 (MQ=255) agcGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCtgatg < 7:108584/115‑1 (MQ=255) agcGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTa > 1:157388/1‑137 (MQ=255) aCAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAAACGCa < 1:52117/83‑1 (MQ=255) aCAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAAACGCa < 1:48878/83‑1 (MQ=255) aCAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAAACGCa > 2:48878/1‑83 (MQ=255) aCAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAAACGCa > 2:52117/1‑83 (MQ=255) cAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGa < 1:204524/139‑1 (MQ=255) cAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGa < 7:189151/139‑1 (MQ=255) ataAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAgg > 2:160568/1‑138 (MQ=255) cTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTaa < 2:157388/139‑1 (MQ=255) ggATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGt < 7:7208/135‑1 (MQ=255) ggATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGt > 8:7208/1‑135 (MQ=255) gATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACTCTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACg > 5:213980/1‑138 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGTCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGaacact < 6:120717/102‑1 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCTCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTa > 4:141019/1‑139 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCa > 6:181014/1‑69 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCa < 5:181014/69‑1 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGa < 7:93105/82‑1 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGa > 8:93105/1‑82 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCtgatg < 1:126228/96‑1 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCtgatg > 2:126228/1‑96 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGaacact > 5:120717/1‑102 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAgg < 2:2779/129‑1 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAgg > 1:2779/1‑129 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGt > 8:118132/1‑133 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGt < 7:118132/133‑1 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTa > 3:157732/1‑139 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTa > 8:98548/1‑139 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTa > 1:202115/1‑139 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTa > 1:50254/1‑139 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTa > 3:148244/1‑139 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTa > 4:33415/1‑139 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTa > 4:64461/1‑139 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTa > 4:69822/1‑139 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTa > 2:2552/1‑139 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATAAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGaacact < 5:86955/102‑1 (MQ=255) aTCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATAAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGaacact > 6:86955/1‑102 (MQ=255) gCTCTCTGCCAAAAACACGATCATGGTCGTCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCaa < 4:137946/139‑1 (MQ=255) gCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTATCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCaa < 8:95558/139‑1 (MQ=255) ctctctGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAg > 1:139475/1‑139 (MQ=255) tGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAGGCAt > 7:166319/1‑138 (MQ=255) tGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAGGCAt > 7:74399/1‑138 (MQ=255) tGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAGGCAt > 1:88617/1‑138 (MQ=255) aaaaCACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGATACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAGGCATTaaaaa > 8:105896/1‑139 (MQ=255) cacGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAGGCATTAAAAAtgtt < 1:160568/139‑1 (MQ=255) acGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAGGCATTAAAAAtgttg < 2:50254/139‑1 (MQ=255) aTGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAGGCATTAAAAATGTTGTTGTCa < 7:98548/139‑1 (MQ=255) cAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAGGCATTAAAAATGTTGTTGTCAAAGATGTc > 6:24918/1‑139 (MQ=255) cAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAGGCATTAAAAATGTTGTTGTCAAAGATGTc > 8:78523/1‑139 (MQ=255) tctTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAGGCATTAAAAATGTTGTTGTCAAAGATGTCCCCGATCCa > 8:42681/1‑137 (MQ=255) atgaAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGt > 6:129611/1‑79 (MQ=255) atgaAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGt < 5:129611/79‑1 (MQ=255) atgaAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAGGCATTAAAAATGTTGTTGTCAAAGATGTCCCCGATCCATAGATTGaaa > 8:73224/1‑138 (MQ=255) aaCAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAGGCATTAAAAATGTTGTTGTCAAAGATGTCCCCGATCCAAAGATTGAAAAATcc < 8:166319/139‑1 (MQ=255) tCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTAc < 4:87845/55‑1 (MQ=255) tCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTAc > 3:87845/1‑55 (MQ=255) | AGAGATTACGGAACTGATGATCAGCAAAAAGTGGTTCCATCCATATGAGCTGAGCGAACAGTATAAGCTGGATCAGCTCTCTGCAAAAAACACGATCATGGTCGGCAATATGAACCTCTTTCCTGATGAAACAAATCGCAAAGGGATGTTTGACCGGACACCTGATGAACACTAACACTGGAGGTTTTACAGCATGAAGGCAGTAACGTATCAAGGCATTAAAAATGTTGTTGTCAAAGATGTCCCCGATCCAAAGATTGAAAAATCC > BS168_reference_genome/2753259‑2753526 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |