Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | BS168_reference_genome | 1,012,003 | G→T | 7.3% | A419E (GCA→GAA) | lytF ← | gamma‑D‑glutamate‑meso‑diaminopimelate muropeptidase (major autolysin, vegetative) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | BS168_reference_genome | 1,012,003 | 0 | G | T | 7.3% | 182.1 / 4.5 | 55 | A419E (GCA→GAA) | lytF | gamma‑D‑glutamate‑meso‑diaminopimelate muropeptidase (major autolysin, vegetative) |
Reads supporting (aligned to +/- strand): ref base G (22/29); new base T (2/2); total (24/31) | |||||||||||
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.72e-01 |
TTGGCATTTATGAATTCCCCGTTGCCTAGGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAAATGAATCCGCTGCAGTCAAACCCGGACGGTGTCGTGCCTCCCCAGCGATACGGAACCCCAAGCTGCGCTTTAGCGGCCGAAATCATTGTGT > BS168_reference_genome/1011871‑1012138 | ttGGCATTTATGAATTCCCCGTTGCCTAGGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGtt > 2:85640/1‑138 (MQ=255) tttATGTCTTCCCCGTTGCCTATGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAg < 6:186703/139‑1 (MQ=255) aTGAATTCCCCGTTGCCTAGGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGAt < 4:57585/138‑1 (MQ=255) tCCCCGTTGCCTAGGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGAtg > 6:190742/1‑139 (MQ=255) gTTGCCTAGGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACt < 1:126921/139‑1 (MQ=255) gTTGCCTAGGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACt < 5:259355/139‑1 (MQ=255) gTTGCCTAGGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACt < 2:19599/139‑1 (MQ=255) cTATGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGtt < 2:176367/139‑1 (MQ=255) cTAGGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGtt < 1:85640/139‑1 (MQ=255) cTAGGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGtt < 6:128462/139‑1 (MQ=255) cTAGGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGtt < 2:109270/139‑1 (MQ=255) aGGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGtt > 1:62188/1‑137 (MQ=255) aGGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGtt > 1:242220/1‑137 (MQ=255) tAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTa < 5:190742/139‑1 (MQ=255) tAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTa < 3:238713/139‑1 (MQ=255) aaatCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACa < 4:124231/136‑1 (MQ=255) aaTTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGAGAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACa < 6:290543/139‑1 (MQ=255) ttCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATa < 2:242220/139‑1 (MQ=255) ttCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATa < 1:188127/139‑1 (MQ=255) aCGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACt < 5:146378/120‑1 (MQ=255) aCGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACt > 6:146378/1‑120 (MQ=255) gCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATACCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAAATGAATCCGCTGCAGt < 2:62188/139‑1 (MQ=255) ttATATGTTGAAAAGAAAACGAAATCTCCGACAGCTGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTTTTAAGTACATAGTAAATGAATCCGCTCCAGTCa < 7:71341/138‑1 (MQ=255) atatGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTTCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTa < 3:276461/110‑1 (MQ=255) atatGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTTCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTa < 3:280960/110‑1 (MQ=255) atatGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTTCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTa > 4:280960/1‑110 (MQ=255) atatGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTTCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTa > 4:276461/1‑110 (MQ=255) ttGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTTATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAAATGAAGCCGCTGCAGTCAAACCCGGa > 7:249507/1‑139 (MQ=255) tGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGtt < 1:112223/73‑1 (MQ=255) tGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGtt > 2:112223/1‑73 (MQ=255) tGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCt > 2:130771/1‑78 (MQ=255) tGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCt < 1:130771/78‑1 (MQ=255) tGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACa < 6:297880/106‑1 (MQ=255) tGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACa > 5:297880/1‑106 (MQ=255) tGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAAATGAATCCGCTGCAGTCAACCCCGGAc > 2:311265/1‑139 (MQ=255) tGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAAATGAATCCGCTGCAGTCAAACCCGGAc > 3:282615/1‑139 (MQ=255) tGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAAATGAATCCGCTGCAGTCAAACCCGGAc > 1:142198/1‑139 (MQ=255) aTCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAAATGAATCCGCTGCAGTCAAACCCGGACGGTGTCGTGCCTccc < 2:142198/139‑1 (MQ=255) cagGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATa > 2:88798/1‑57 (MQ=255) cagGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATa < 1:120031/57‑1 (MQ=255) cagGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATa < 1:88798/57‑1 (MQ=255) cagGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATa > 2:120031/1‑57 (MQ=255) aaCTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAAATGAATCCGCTGCAGTCAAACCCGGACGGTGTCGTGCCTCCCCAGCGATACGGAACCCCAAGc > 3:122256/1‑139 (MQ=255) aTTGTGCTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAACTGAATCCGCTGCAGTCAAACCCGGACGGGGTCGTGCCTCCCCAGCGATACGGAACCCCAAGCTGCGCTTTAg < 5:230948/139‑1 (MQ=255) ttCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGtt > 5:13013/1‑47 (MQ=255) ttCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGtt < 6:13013/47‑1 (MQ=255) ttCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTa > 2:40830/1‑52 (MQ=255) ttCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTa < 1:40830/52‑1 (MQ=255) tatCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAAATGAATCCGCTGCAGTCAAACCCGGACGGTGTCGTGCCTCCCCAGCGATACGGAACCCCAAGCTGCGCTTTAGCGGCCGAAATCa < 5:77050/139‑1 (MQ=255) ccGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAAATGAATCCGCTGCAGTCAAACCCGGACGGTGTCGTGCCTCCCCAGCGATACGGAACCCCAAGCTGCGCTTTAGCGGCCGAAATCATtgt > 5:111326/1‑139 (MQ=255) cTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTa > 4:51033/1‑39 (MQ=255) cTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTa < 3:51033/39‑1 (MQ=255) cTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAAATGAATCCGCTGCAGTCAAACCCGGACGGTGTCGTGCCTCCCCAGCGATACGGAACCCCAAGCTGCGCTTTAGCGGCCGAAATCATTgtgt > 5:95937/1‑138 (MQ=255) cTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAAATGAATCCGCTGCAGTCAAACCCGGACGGTGTCGTGCCTCCCCAGCGATACGGAACCCCAAGCTGCGCTTTAGCGGCCGAAATCATTgtgt > 3:226128/1‑138 (MQ=255) cTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAAATGAATCCGCTGCAGTCAAACCCGGACGGTGTCGTGCCTCCCCAGCGATACGGAACCCCAAGCTGCGCTTTAGCGGCCGAAATCATTgtgt > 2:188521/1‑138 (MQ=255) | TTGGCATTTATGAATTCCCCGTTGCCTAGGTAAATTCCGACGTGGGAAGGGCCGGCTTTATATGTTGAAAAGAAAACGAAATCTCCGACAGCAGGCTGGCTAACTGATTTCATTGTGTTCCAATATCCCGCTGCCGTTAATCTAGATACGGATGTGACTTTGTTCAGTACATAGTAAATGAATCCGCTGCAGTCAAACCCGGACGGTGTCGTGCCTCCCCAGCGATACGGAACCCCAAGCTGCGCTTTAGCGGCCGAAATCATTGTGT > BS168_reference_genome/1011871‑1012138 |
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: ‑ |