Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 649,991 | C→T | A271T (GCA→ACA) | 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,991 | 0 | C | T | 98.0% | 184.9 / NA | 51 | A271T (GCA→ACA) | mrdA | transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2) |
Reads supporting (aligned to +/- strand): ref base C (0/0); major base T (35/15); minor base G (0/1); total (35/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.14e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01 |
TTGGGTCATAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGCGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAATATCGTGT > REL606/649919‑650061 | ttGGGTCATAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGtgcg < 1:701236/76‑1 (MQ=255) tGGGTCATAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCgg < 2:1295264/76‑1 (MQ=255) tGGGTCATAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCgg > 1:137079/1‑76 (MQ=255) aTAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTAccc > 2:1031943/1‑76 (MQ=255) aaCTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGc > 1:1048420/1‑76 (MQ=255) aaCTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGc > 1:485942/1‑76 (MQ=255) cTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCca < 1:1712532/76‑1 (MQ=255) cTAGGCGGGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGGGCGGCTACCCGCca < 1:129207/76‑1 (MQ=255) tAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCcag > 2:1779480/1‑76 (MQ=255) tAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCcag > 1:17560/1‑76 (MQ=255) tAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCcag > 2:380435/1‑76 (MQ=255) tAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCcag < 2:138881/76‑1 (MQ=255) aGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCcagc > 1:223151/1‑76 (MQ=255) gCGGGGAACCCAGCGCCAGCCCCCCCCCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCcagcag < 2:1324258/76‑1 (MQ=255) cGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCagcagc < 2:235184/76‑1 (MQ=255) gTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCg > 1:861887/1‑76 (MQ=255) gAAACCAGCGCCAGCCCCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGttt < 2:1576290/76‑1 (MQ=255) gAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGttt < 2:21597/76‑1 (MQ=255) gAAACCACCGCCAGCACCCCGCCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGttt < 2:64272/76‑1 (MQ=255) aaaCCAGCGCCAGCCCCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTc < 2:1370858/76‑1 (MQ=255) aaaCCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTc > 2:1825576/1‑76 (MQ=255) aaaCCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTc > 2:1020889/1‑76 (MQ=255) aaaCCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGGTTc > 2:636513/1‑76 (MQ=255) aaCCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCCGCAGCGGTTCa > 1:255911/1‑76 (MQ=255) aaCCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCa > 2:1204554/1‑76 (MQ=255) aCCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCaa > 2:366381/1‑76 (MQ=255) ccAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAat < 1:1861592/76‑1 (MQ=255) cAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGtt < 2:68202/70‑1 (MQ=255) cAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGtt > 1:68202/1‑70 (MQ=255) aGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAatat < 2:1416018/76‑1 (MQ=255) gcCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATAtt > 2:1695785/1‑76 (MQ=255) gccccccccCGGTGCGCGGATCGGTGCCTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATGTATTGCTg < 2:1530102/73‑1 (MQ=255) cACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGCTTCCATATCTTGCTgg > 2:109291/1‑76 (MQ=255) aCCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGa > 2:1185757/1‑76 (MQ=255) aCCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGa > 2:1450031/1‑76 (MQ=255) cACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGttt > 1:412037/1‑76 (MQ=255) ccTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTga > 1:742519/1‑76 (MQ=255) ccTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTga > 2:1776764/1‑76 (MQ=255) tGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGTGTTTgaga > 2:1857924/1‑76 (MQ=255) tGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTgaga < 2:845981/76‑1 (MQ=255) tACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATc > 2:1494761/1‑76 (MQ=255) cgcgGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCa > 2:1622298/1‑76 (MQ=255) tCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCa > 2:673069/1‑76 (MQ=255) cGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAg > 1:630258/1‑76 (MQ=255) ggTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAgg > 2:1433983/1‑76 (MQ=255) tGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTa < 2:137079/76‑1 (MQ=255) tGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTa > 2:1264180/1‑76 (MQ=255) tGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTa > 1:948503/1‑76 (MQ=255) tGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTa > 1:1395625/1‑76 (MQ=255) tGACTACCACAGCTGTGCGGCGCCCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTa < 2:1708205/76‑1 (MQ=255) tACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAata > 1:746418/1‑76 (MQ=255) aCCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAatat > 2:1328904/1‑76 (MQ=255) aCCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAatat < 1:251950/76‑1 (MQ=255) acaGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAATATCgt > 2:711506/1‑76 (MQ=255) aGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGGAAATATCgtgt > 2:1840448/1‑76 (MQ=255) | TTGGGTCATAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGCGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAATATCGTGT > REL606/649919‑650061 |
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: ‑ |