Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 649,829 | C→A | V325F (GTT→TTT) | mrdA ← | transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 649,829 | 0 | C | A | 97.5% | 119.6 / ‑5.5 | 40 | V325F (GTT→TTT) | mrdA | transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2) |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base A (15/24); total (15/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01 |
TGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGAATAATCTTTGCTGGAGAT > REL606/649757‑649901 | tGCGCGTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAaccc < 1:1314567/76‑1 (MQ=255) cgcgTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCt < 1:1357987/76‑1 (MQ=255) cgcgTGATCCCCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCt < 1:185487/76‑1 (MQ=255) gTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCtgtg > 1:148294/1‑76 (MQ=255) tCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAAACCCTGTGTGgc > 2:270304/1‑76 (MQ=255) cACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGCTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGgcg > 2:695798/1‑76 (MQ=255) gcTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCAcc < 1:259393/76‑1 (MQ=255) gcTAAAGGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCAcc < 1:844560/76‑1 (MQ=255) tCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAg < 2:1235349/76‑1 (MQ=255) aaggCCAAACCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCg < 1:1012164/73‑1 (MQ=255) aaTGCCGAACCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCg < 2:26812/76‑1 (MQ=255) tGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGgt > 2:29727/1‑76 (MQ=255) ccGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGgtgt < 2:169554/76‑1 (MQ=255) ccGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGgtgt > 1:8515/1‑76 (MQ=255) ccGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGgtgt > 1:723993/1‑76 (MQ=255) ccGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGgtgt > 1:406567/1‑76 (MQ=255) cagAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGgtgt < 2:314738/74‑1 (MQ=255) cGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTa > 1:176871/1‑76 (MQ=255) cGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTa > 1:1501138/1‑76 (MQ=255) aaCCGTCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTc < 1:1347560/76‑1 (MQ=255) aaCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTc > 2:785903/1‑76 (MQ=255) aaCCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTCCACCAGCGGTGTATTc > 2:231728/1‑76 (MQ=255) cacaTAGTGTTTAACTGTAGGCGCGGGAGGATAAATCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCg < 2:362282/76‑1 (MQ=255) cacaTAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCg > 2:554766/1‑76 (MQ=255) caTAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGtt < 2:329237/76‑1 (MQ=255) aTAGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTc > 1:891868/1‑76 (MQ=255) aGGGTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTcaa > 2:842134/1‑76 (MQ=255) ggTTTAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTcaaca < 2:1014974/76‑1 (MQ=255) tAACTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGc > 1:940529/1‑76 (MQ=255) aaCTGGAGACGCGGGAGGATAAAGCCGCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCg < 2:1437601/76‑1 (MQ=255) aCTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCgg < 1:978167/76‑1 (MQ=255) aCTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCgg < 1:696856/76‑1 (MQ=255) cTGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGa < 1:819675/76‑1 (MQ=255) tGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaa < 1:966035/76‑1 (MQ=255) tGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaa < 1:774059/76‑1 (MQ=255) tGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaa < 1:376733/76‑1 (MQ=255) tGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaa < 1:248901/76‑1 (MQ=255) tGTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaa < 1:1293955/76‑1 (MQ=255) gTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaat < 1:774330/76‑1 (MQ=255) gTAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaat < 1:1326375/76‑1 (MQ=255) tAGACGCGGGAGGATAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGaata < 1:1473472/76‑1 (MQ=255) gcgGGAGGATAACCCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGAATAATCtt < 2:1424112/76‑1 (MQ=255) ggaTAAAACCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGAATAATCTTTGCTgg < 1:1365706/76‑1 (MQ=255) aTAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGAATAATCTTTGCTGgag < 1:583364/76‑1 (MQ=255) aaaaCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGAATAATCTTTGCTGGAGAt < 1:689084/76‑1 (MQ=255) | TGCGCGTGATCACCCCGGCGCTCAATGCCGAAACCGCCACATAGGGTTTAACTGTAGACGCGGGAGGATAAACCCCCTGTGTGGCGCGGTTCACCAGCGGTGTATTCGGGTCGTTCAACAAGGCGGAATAATCTTTGCTGGAGAT > REL606/649757‑649901 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |