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 | 87.1% | 85.4 / 5.5 | 31 | I401V (ATC→GTC) | mrdA | transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2) |
Reads supporting (aligned to +/- strand): ref base T (0/4); new base C (14/13); total (14/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.07e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01 |
GTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACG > REL606/649529‑649673 | gTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATccc < 1:937933/76‑1 (MQ=255) ttcttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACcccc > 2:460604/1‑76 (MQ=255) tcttcCGCCAGGTCNATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCa > 1:872933/1‑76 (MQ=255) ccGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATc > 2:1433560/1‑76 (MQ=255) gCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGt < 2:351735/76‑1 (MQ=255) ccAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTa < 2:1226554/76‑1 (MQ=255) cAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAg < 2:400377/76‑1 (MQ=255) aGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAgg < 1:1252366/76‑1 (MQ=255) tACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTg < 2:363594/76‑1 (MQ=255) gtgtAATGACCATAGCCGAATTGTCCCATCCATACGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATaga < 2:163009/76‑1 (MQ=255) gtgtAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATaga > 2:111514/1‑76 (MQ=255) tgtAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATagaa < 2:1315568/76‑1 (MQ=255) tAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATagaaga > 1:1362864/1‑76 (MQ=255) gACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGt < 1:482111/76‑1 (MQ=255) aTAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATcc > 2:379960/1‑76 (MQ=255) aTAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATcc > 1:984116/1‑76 (MQ=255) tAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCg < 2:93983/76‑1 (MQ=255) tAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGCATCCg > 1:171465/1‑76 (MQ=255) tAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGAATCCg > 2:41952/1‑76 (MQ=255) aaTTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATACGCAGAtt > 1:194859/1‑76 (MQ=255) tACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGAttcttc < 1:1261684/76‑1 (MQ=255) aCCCATCCATTCGGAGAGGCGATCGCCCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTcc < 2:555867/76‑1 (MQ=255) cccATCCATTCGGAGAGGTGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCa < 1:1144226/76‑1 (MQ=255) cccATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCa < 1:1404782/76‑1 (MQ=255) ccatTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGAt > 1:549213/1‑76 (MQ=255) cGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAATGTATCCGCAGATTCTTCCATCGATCTtgt < 2:194859/76‑1 (MQ=255) cGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTtgt > 1:938982/1‑76 (MQ=255) cGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTtgt > 1:590859/1‑76 (MQ=255) ggAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTtgtg < 1:822460/76‑1 (MQ=255) agagGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGAc > 1:1246976/1‑76 (MQ=255) agGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACAt < 2:1205706/76‑1 (MQ=255) gCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTc < 1:751410/76‑1 (MQ=255) atcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGGATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGa > 2:885854/1‑76 (MQ=255) tcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGAc > 2:418163/1‑76 (MQ=255) tcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGAc > 1:795798/1‑76 (MQ=255) cgaccCCCATATCCTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACg > 1:1060590/1‑76 (MQ=255) | GTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACG > REL606/649529‑649673 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |