Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 2,854,034 | T→C | 5.9% | M176V (ATG→GTG) | spo0B ← | sporulation initiation phosphotransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 2,854,034 | 0 | T | C | 5.9% | 208.1 / 4.1 | 68 | M176V (ATG→GTG) | spo0B | sporulation initiation phosphotransferase |
Reads supporting (aligned to +/- strand): ref base T (34/30); new base C (2/2); total (36/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01 |
GTCGCCGCCTTTTACATATACTTTGACCTGATCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGAATCAGCTGTCTGTCAGGATGATCCGTTTGAAGCGAAACCGTTAAAT > BS168_reference_genome/2853906‑2854165 | gTCGCCGCCTTTTACATATACTTTGACCTGATCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACg < 6:369837/139‑1 (MQ=255) gTCGCCGCCTTTTACATATACTTTGACCTGATCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACg < 8:285400/139‑1 (MQ=255) gTCGCCGCCTTTTACATATACTTTGACCTGATCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACg < 2:352733/139‑1 (MQ=255) gTCGCCGCCTTTTACATATACTTTGACCTGATCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACg < 2:14842/139‑1 (MQ=255) taCTTTGACCTGATCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTc > 6:373381/1‑123 (MQ=255) taCTTTGACCTGATCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTc < 5:373381/123‑1 (MQ=255) taCTTTGACCTGATCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGc > 5:204101/1‑139 (MQ=255) cTTTGACCTGATCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGcc < 8:353355/138‑1 (MQ=255) tttGACCTGATCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCg > 7:173726/1‑138 (MQ=255) tGACCTGATCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGaa > 8:26224/1‑138 (MQ=255) tGACCTGATCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTACGCAAGTCCAACCCAATTTAAATCACACATTCGTGACTCGTGATTTAAAAACTCATGATATAAACATCCTCAAATACATTATGCAGaa > 2:49046/1‑138 (MQ=255) aCCTGATCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAata > 8:273659/1‑138 (MQ=255) aTCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAata < 4:388387/133‑1 (MQ=255) aTCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAata < 4:68549/133‑1 (MQ=255) aTCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAata > 3:388387/1‑133 (MQ=255) aTCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAata > 3:68549/1‑133 (MQ=255) aTCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATatcatc < 4:34871/138‑1 (MQ=255) aTCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATatcatc > 5:257972/1‑138 (MQ=255) aTCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATatcatc > 4:328075/1‑138 (MQ=255) aTCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATatcatc > 5:316243/1‑138 (MQ=255) aTAATGTCCTCCGTTCTAAACCGTTTAAAACTCCGCTATTCCAACCCAAATTCAATGAGAAATTCGTGACTCGTGATTTCATAACGCATGATATCAACGTCCTCTTATCCATTCAGCAGAATATCATCAAAAGCatatt > 3:403538/1‑138 (MQ=255) ctccGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCt > 2:185458/1‑95 (MQ=255) ctccGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCt > 2:179621/1‑95 (MQ=255) ctccGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCt < 1:179621/95‑1 (MQ=255) ctccGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCt < 1:185458/95‑1 (MQ=255) ctccGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGc > 7:229640/1‑138 (MQ=255) ctccGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGc > 7:57934/1‑138 (MQ=255) tccGTTCTAAACCGTTAAAAACTCAGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCa > 1:77225/1‑138 (MQ=255) ccGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCGACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCa > 7:224621/1‑137 (MQ=255) ccGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAata > 8:193083/1‑113 (MQ=255) ccGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAata < 7:193083/113‑1 (MQ=255) ccGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCACGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCaaa > 8:358543/1‑139 (MQ=255) ccGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCACGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCa > 2:361255/1‑137 (MQ=255) tAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGcc > 4:278724/1‑139 (MQ=255) aCCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTg < 7:89364/139‑1 (MQ=255) cGTTAAAAACTCCGCTAGTCCACCCCAATTTCAATCAGACGTTCGTGACTCGTGATTTAAAAACGCATGCTATCAACGTCCTCATATCCATTATGCCGAATATCTACAAAAGCAGAAGGAACTGCAAAGGCGCCGTGa > 6:206313/1‑138 (MQ=255) cGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCACGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGaa < 7:358543/139‑1 (MQ=255) cGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCACGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGaa < 1:361255/139‑1 (MQ=255) gTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATatcat < 8:138381/105‑1 (MQ=255) gTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATatcat > 7:138381/1‑105 (MQ=255) ttAAAACCTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGaaaa > 1:19323/1‑139 (MQ=255) ttAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACg < 7:249298/76‑1 (MQ=255) ttAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACg > 8:249298/1‑76 (MQ=255) ttAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATgg < 7:9858/117‑1 (MQ=255) ttAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATgg > 8:9858/1‑117 (MQ=255) ttAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGt < 6:189478/134‑1 (MQ=255) ttAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGt > 5:189478/1‑134 (MQ=255) ttAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGaaaa > 7:63737/1‑139 (MQ=255) ttAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGaaaa > 3:5922/1‑139 (MQ=255) ttAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGGGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGaaaa > 2:236158/1‑139 (MQ=255) aaaaCTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCa < 3:278724/139‑1 (MQ=255) aaaaCTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCa < 4:34696/139‑1 (MQ=255) cTCCGCTAGTCAAACCCAATTTCTATCATTCATTCGTGACTCTTGATTTAAAAACTAATGATATCAACGTCCTCATATCCATTCTGACGAATATCATCATATTCAGATGGATCGGCATTGGCGCAGTGAAAATCAAGGt < 1:49046/139‑1 (MQ=255) gTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGaa < 6:316243/137‑1 (MQ=255) cAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGAATCAg < 7:26224/138‑1 (MQ=255) aaCCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGAATCAGct < 3:328075/139‑1 (MQ=255) cccAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGAATCAGctgt > 8:53735/1‑139 (MQ=255) aaattCAATCAGACAATCGTGACTCGTGATTTCAAAACGCATGATAACAACGTCCTCATATCCCTTCTGCCGAATATCATCAAAAGCAGATGCATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGAATCAGctgtctg < 4:415797/136‑1 (MQ=255) tttCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCTGATGGATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGAATCAGCtgtctgtc < 6:204101/139‑1 (MQ=255) aTCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGAATCAGCTGTCTGTCAGgat > 7:237076/1‑139 (MQ=255) ttcggcaCTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGTATCGGCAAAGGCGCCGTGTAAATCAAGGTACAGAATCAGCTGTCTGTCAGGATGATCCGtt < 4:403538/133‑1 (MQ=255) ttCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGAATCAGCTGTCTGTCAGGATGATCCGtt < 8:21190/139‑1 (MQ=255) agTGACTCTTGATTTCAAAACGCATGATATCAACGTCCTCATATCGATTCTGCCTAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGAATCAGCTGTCTGTCAGGATGATCCGTTTg < 3:101794/138‑1 (MQ=255) cTCGTGATTTCAAAACGCATGATATCAACTTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGAGCCGTGAAAATCAAGGTACAGAATCAGCTGTCTGTCAGCATCAACCGTTTGAAGc > 2:209339/1‑138 (MQ=255) cGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGAATCAGCTGGCTGTCAGGATGCTCCGTTTGAAGCGa > 5:268750/1‑138 (MQ=255) ttCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGAATCAGCTGTCAGGCAGGATGAACCGTTAGAAGCGAAAcagct > 4:101794/1‑135 (MQ=255) aaaaCGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGAATCAGCTGTCTGTCAGGATGATCCGTTTGAAGCGAAACCGTTaa < 6:102964/138‑1 (MQ=255) aaaCGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAACAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGAATCAGCTGTCTGTCAGGATGATCCGTTTGAAGCGAAACCGTTAAAt < 4:423201/139‑1 (MQ=255) | GTCGCCGCCTTTTACATATACTTTGACCTGATCTACAAACATAATGTCCTCCGTTCTAAACCGTTAAAAACTCCGCTAGTCCAACCCAATTTCAATCAGACATTCGTGACTCGTGATTTCAAAACGCATGATATCAACGTCCTCATATCCATTCTGCCGAATATCATCAAAAGCAGATGGATCGGCAAAGGCGCCGTGAAAATCAAGGTACAGAATCAGCTGTCTGTCAGGATGATCCGTTTGAAGCGAAACCGTTAAAT > BS168_reference_genome/2853906‑2854165 |
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: ‑ |