Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 649,601 | T→C | I401V (ATC→GTC) | 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,601 | 0 | T | C | 97.0% | 121.6 / ‑3.4 | 33 | I401V (ATC→GTC) | mrdA | transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2) |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (17/15); total (18/15) | |||||||||||
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.05e-01 |
CGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGA > REL606/649528‑649671 | cGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGAccc > 2:773069/1‑76 (MQ=255) cGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGAccc < 1:282747/76‑1 (MQ=255) gTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATGCGGAGAGGCGATCGAcccc < 2:1267927/76‑1 (MQ=255) cttcCGCAAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCat > 2:1220267/1‑76 (MQ=255) ccAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTa > 1:1103911/1‑76 (MQ=255) gTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGcc > 1:390697/1‑76 (MQ=255) tCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCa < 2:1009091/76‑1 (MQ=255) cGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCAc > 2:200443/1‑76 (MQ=255) aTACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCt < 1:1543324/76‑1 (MQ=255) aTACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCt > 2:1378117/1‑76 (MQ=255) tACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTg > 2:710347/1‑76 (MQ=255) tccGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGa < 2:1790981/75‑1 (MQ=255) ccGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGAt > 1:1211700/1‑76 (MQ=255) cGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATa > 2:349656/1‑76 (MQ=255) ggTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATag < 1:1726737/76‑1 (MQ=255) tAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATagaaga < 1:1430205/76‑1 (MQ=255) aTGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGaagaag > 2:164591/1‑76 (MQ=255) aTGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGaagaag < 2:1352868/76‑1 (MQ=255) aTGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGaagaag > 1:1694707/1‑76 (MQ=255) gACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGt < 1:629383/76‑1 (MQ=255) gACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGt < 1:1014200/76‑1 (MQ=255) cATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATc > 2:1413165/1‑76 (MQ=255) aTAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATcc > 2:1647262/1‑76 (MQ=255) aTAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATcc > 1:886528/1‑76 (MQ=255) cGAATTTACCCATCCATGCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGa < 1:224851/76‑1 (MQ=255) aaTTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGAtt > 2:909552/1‑76 (MQ=255) tttACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGAttct < 1:1585530/76‑1 (MQ=255) tACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGAttcttc > 1:281778/1‑76 (MQ=255) cccATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCa > 1:1542573/1‑76 (MQ=255) cccATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCa < 1:1196808/76‑1 (MQ=255) catccatTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGc < 2:1488558/76‑1 (MQ=255) ccatTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGCATCCGCAGATTCTTCCAGCGAt > 2:1208263/1‑76 (MQ=255) ccatGCGGGGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGAt < 1:303033/76‑1 (MQ=255) ttCGGAGATGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCtt < 2:472265/76‑1 (MQ=255) ggAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTtgtg < 1:1755127/76‑1 (MQ=255) gagaGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGGATCCGCAGATTCTTCCAGCGATCCTGTGa > 1:599238/1‑76 (MQ=255) agGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACAt < 1:948964/76‑1 (MQ=255) cgatcgacCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCa > 1:707651/1‑76 (MQ=255) gatcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAg > 1:1503738/1‑76 (MQ=255) atcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGa < 2:1101169/76‑1 (MQ=255) atcgaccCCCAGAGCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGa < 2:25484/76‑1 (MQ=255) | CGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGA > REL606/649528‑649671 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |