Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 649,342 | G→T | T487N (ACC→AAC) | 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,342 | 0 | G | T | 98.2% | 203.2 / ‑2.6 | 54 | T487N (ACC→AAC) | mrdA | transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2) |
Reads supporting (aligned to +/- strand): ref base G (0/1); new base T (28/25); total (28/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.81e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01 |
GTAACCGGAATGAATATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTG > REL606/649269‑649412 | gTAACCGGAATGAATATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTg > 2:1190786/1‑76 (MQ=255) tAACCGGAATGAATATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGc > 2:1301034/1‑76 (MQ=255) aCCGGAATGAATATCGCCGATGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCtc > 1:443992/1‑76 (MQ=255) gcGGAATGAATATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCtca < 1:983917/75‑1 (MQ=255) ccGGAATGAATATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCtca < 1:1190786/76‑1 (MQ=255) ccGGAATGAATATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCtca < 2:1792015/76‑1 (MQ=255) cGGAATGAATATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCtcct > 2:1136117/1‑74 (MQ=255) ggAATGAATATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCtcatc < 2:190704/76‑1 (MQ=255) ggAATGAATATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCtcatc > 2:1757078/1‑76 (MQ=255) atgaatATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGc < 1:63975/76‑1 (MQ=255) tgaatATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCa > 2:1255127/1‑76 (MQ=255) tgaatATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCa < 1:646876/76‑1 (MQ=255) gaatATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCaa < 1:1295803/76‑1 (MQ=255) gaatATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCaa < 2:2969/76‑1 (MQ=255) gaatATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCaa < 1:690088/76‑1 (MQ=255) atatCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAAt < 1:646962/76‑1 (MQ=255) atCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGc > 2:318572/1‑76 (MQ=255) atCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGc > 2:967608/1‑76 (MQ=255) tCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCg > 1:846476/1‑76 (MQ=255) ccGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGaa > 1:1191276/1‑76 (MQ=255) gACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAAcc > 2:1566426/1‑76 (MQ=255) aCGGGCTGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCt < 2:1618497/76‑1 (MQ=255) cGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCtt > 2:1164479/1‑76 (MQ=255) gggCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTc > 1:712971/1‑76 (MQ=255) cGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACg < 1:203777/76‑1 (MQ=255) ggTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGa > 1:1827081/1‑76 (MQ=255) ggTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGa > 1:31413/1‑76 (MQ=255) ggTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGa > 1:1393263/1‑76 (MQ=255) ggTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGa > 1:956080/1‑76 (MQ=255) cATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATAcc > 1:1648262/1‑76 (MQ=255) gCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTc > 1:1588149/1‑76 (MQ=255) ttGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTa < 2:1736379/76‑1 (MQ=255) tGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTaa < 1:1200845/76‑1 (MQ=255) tGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTaa < 2:1417322/76‑1 (MQ=255) tGTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTaa < 2:740675/76‑1 (MQ=255) gTACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAAt < 2:579943/76‑1 (MQ=255) tACCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATc > 2:1594989/1‑76 (MQ=255) aancATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCa > 2:1336586/4‑76 (MQ=255) aCCCATGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCa > 2:744871/1‑76 (MQ=255) cgnaTGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAg > 2:1331702/4‑76 (MQ=255) cccATGGCCCCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAg < 1:1129745/76‑1 (MQ=255) tGGCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGatc < 1:397151/76‑1 (MQ=255) ggCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGatca < 1:641813/76‑1 (MQ=255) gCACCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGatcat > 1:794242/1‑76 (MQ=255) cACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGatcatc < 1:1789962/76‑1 (MQ=255) cACCTATTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGatcatc < 1:1155478/76‑1 (MQ=255) aCCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAtcatca > 2:1249706/1‑76 (MQ=255) aCCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAtcatca < 2:296548/76‑1 (MQ=255) aCCTGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAtcatca > 1:533958/1‑76 (MQ=255) tGTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTg < 1:1589219/76‑1 (MQ=255) gTTTGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGc > 1:971976/1‑76 (MQ=255) tttGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGcc > 1:430517/1‑76 (MQ=255) tGCCGTCTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCtt < 1:1859422/76‑1 (MQ=255) cTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGCTACCGACATTAATCAGTATCATCAGTGCCTTACTCAt > 1:859209/1‑76 (MQ=255) cTTCGGCGTTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCAt < 1:583987/76‑1 (MQ=255) cggcgtTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTg < 1:528145/76‑1 (MQ=255) | GTAACCGGAATGAATATCGCCGACGGGCGGTTCATGCGGTTGTACCCATGGCACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTG > REL606/649269‑649412 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |