Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 649,489 | C→A | R438L (CGC→CTC) | 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,489 | 0 | C | A | 94.4% | 119.3 / 1.4 | 36 | R438L (CGC→CTC) | mrdA | transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2) |
Reads supporting (aligned to +/- strand): ref base C (1/1); new base A (13/21); total (14/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-01 |
GTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGA > REL606/649420‑649560 | gTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTcn < 1:1422556/76‑2 (MQ=255) gTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCt < 1:1404835/76‑1 (MQ=255) tCGCTGTCCAGTAACCCTGCCCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTg < 2:769569/76‑1 (MQ=255) cGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGt > 1:593141/1‑76 (MQ=255) tGTCCAGTAACCCTGACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTc < 2:532792/76‑1 (MQ=255) gTCCAGTAACCCTGACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTcc > 2:1660018/1‑76 (MQ=255) ccAGTAACCCTGACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCAt < 1:1680245/76‑1 (MQ=255) aaCCCTGACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTcgcg > 2:1307135/1‑76 (MQ=255) aaCCCTGACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTcgcg > 2:868168/1‑76 (MQ=255) cTGACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTa < 1:1061022/76‑1 (MQ=255) tGACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAg < 1:1006563/76‑1 (MQ=255) tGACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAg > 1:501656/1‑76 (MQ=255) gACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAgg > 2:705223/1‑76 (MQ=255) gACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAgg > 1:972144/1‑76 (MQ=255) ccGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCa < 1:777954/76‑1 (MQ=255) ccGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGAAGGCa < 2:462931/76‑1 (MQ=255) cGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCAt > 2:818420/1‑76 (MQ=255) cGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCAt < 2:1201181/76‑1 (MQ=255) cGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCAt < 1:1240354/76‑1 (MQ=255) tACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTa < 1:1374243/68‑1 (MQ=255) tACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTa > 2:1374243/1‑68 (MQ=255) tACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCATGtt < 2:998705/76‑1 (MQ=255) cGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGaa < 2:461015/76‑1 (MQ=255) tgtCACCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGttcttc > 2:1536924/1‑76 (MQ=255) aCCCTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGcc < 2:160335/76‑1 (MQ=255) cTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAgg < 1:1472171/76‑1 (MQ=255) cTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAgg < 2:142668/76‑1 (MQ=255) cTGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAgg < 1:1692043/76‑1 (MQ=255) tGATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGt > 1:1026306/1‑76 (MQ=255) gATACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTc < 2:1174387/76‑1 (MQ=255) aTACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCg < 1:373002/76‑1 (MQ=255) tACCACGGTTTTTTAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGa > 1:1088892/1‑76 (MQ=255) cACGGTTTTTTAAAGAGTTGCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATAc > 2:1264441/1‑76 (MQ=255) tAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATg > 2:618127/1‑76 (MQ=255) tAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATg < 1:582618/76‑1 (MQ=255) tAAAGAGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATg < 2:782069/76‑1 (MQ=255) aaaGAGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGa < 1:152585/76‑1 (MQ=255) | GTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGGTGTCACCCTGATACCACGGTTTTTTAAAGCGTTTCTGTTTCCATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGA > REL606/649420‑649560 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |