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.4% | 211.5 / NA | 62 | A271T (GCA→ACA) | mrdA | transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2) |
Reads supporting (aligned to +/- strand): ref base C (1/0); major base T (27/33); minor base G (1/0); total (29/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.59e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01 |
GTCATAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGCGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAATATCGTGTCC > REL606/649923‑650063 | gTCATAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTa < 1:1177692/76‑1 (MQ=255) gTCATAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTa > 1:1661851/1‑76 (MQ=255) cATAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTAcc > 1:1025154/1‑76 (MQ=255) aTAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTAccc > 2:1545903/1‑76 (MQ=255) tAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCg < 2:42131/76‑1 (MQ=255) tAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGGGCGGCTAcccc > 2:373346/1‑75 (MQ=255) tAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGCGCGGCTAcctg > 1:435034/1‑74 (MQ=255) ccTAGGCGGGGAACCCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGcc < 2:1544541/75‑1 (MQ=255) aCTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGcc > 1:1059538/1‑76 (MQ=255) aCTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGcc < 1:1567274/76‑1 (MQ=255) cTAGGCGTGGAACCCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCca < 1:1352176/76‑1 (MQ=255) cTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCca < 2:1061225/76‑1 (MQ=255) tAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTCCCACAGCTGTGCGGCTACCCGAcac > 2:647973/1‑75 (MQ=255) tAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCcag < 2:937494/76‑1 (MQ=255) tAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCcag > 1:575439/1‑76 (MQ=255) tAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCcag > 1:108588/1‑76 (MQ=255) tAGGCGAGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCcag < 2:814762/76‑1 (MQ=255) gaGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTCCCACAGCTGTGCGGCTACCCGCcag < 2:749363/75‑1 (MQ=255) ggggAACCCAGCGCCAGCGCCCCGCCTGTACGCGGAGCGGTGACTGCCACAGCTGTGCGGCTACCCGCCAGCAGCg < 2:1455329/74‑1 (MQ=255) gTGGAACCCAGCGCCAGCCCCCCACCTGTCCGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCg < 1:1206210/76‑1 (MQ=255) tGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGt > 2:553483/1‑76 (MQ=255) tGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGt > 1:1309688/1‑76 (MQ=255) ggAACCCAGCGCCAGCCCCCCCCCTGTACGCGGATCGGGGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGtt < 1:1365442/76‑1 (MQ=255) aaCCAGCGCCAGCCCCCCACCGGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCa < 1:1279550/76‑1 (MQ=255) ccggcgcCAGCCCCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAat < 1:476032/73‑1 (MQ=255) ccAGCGCCAGCCCCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAat < 1:1636661/76‑1 (MQ=255) ccAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAat < 2:1668823/76‑1 (MQ=255) cagccccagcccCCCACCTGTCCGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAata < 1:117531/76‑1 (MQ=255) cAGCGCCAGCACCCCACCTGTAGGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAata < 2:471702/76‑1 (MQ=255) cAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAata < 2:835201/76‑1 (MQ=255) cgcCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGCGCGGCTACCCGCCAGCAGCGTTTCAatatat > 2:177456/1‑76 (MQ=255) cgcCAGCACCCCACCTGTACGCGAATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAatatat > 2:109046/1‑76 (MQ=255) cAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGc < 1:1705500/76‑1 (MQ=255) aGCACCCCACCTGTCCGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCt < 1:223408/76‑1 (MQ=255) gCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTg < 2:49204/76‑1 (MQ=255) cACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTgg > 1:1165445/1‑76 (MQ=255) aCCCCACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGa > 2:210573/1‑76 (MQ=255) cccACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGt > 1:383707/1‑76 (MQ=255) cccACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGt > 2:426607/1‑76 (MQ=255) cccACCTGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGt < 1:1489645/76‑1 (MQ=255) cccccTGTCCGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGtt < 2:1671769/73‑1 (MQ=255) tGTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTgaga > 1:1470374/1‑76 (MQ=255) gTACGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGAt > 2:157699/1‑76 (MQ=255) aCGCGGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATcc < 2:1594698/76‑1 (MQ=255) gcgGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAg < 1:737325/76‑1 (MQ=255) gcgGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGGTTGAGATCCAg > 2:590884/1‑76 (MQ=255) cgGATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGc > 1:1218592/1‑76 (MQ=255) ggATCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCg < 1:1557245/76‑1 (MQ=255) gATCGGTGTCTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGt < 1:281702/76‑1 (MQ=255) aTCGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTc < 2:364831/76‑1 (MQ=255) cGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAg < 2:119559/76‑1 (MQ=255) cGGTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAg > 1:1369779/1‑76 (MQ=255) ggTGACTCCCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAgg < 2:152158/76‑1 (MQ=255) gTGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAggg > 1:910662/1‑75 (MQ=255) tGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAggca > 2:1360741/1‑74 (MQ=255) tGACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTa > 2:924547/1‑76 (MQ=255) gACTGCCAGATCTGGGGGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTaa < 2:635006/76‑1 (MQ=255) gACTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTaa < 2:1376536/76‑1 (MQ=255) aCTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTaaa < 2:582855/76‑1 (MQ=255) aCTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTaaa < 2:1003200/76‑1 (MQ=255) cTACCACAGCTGTGCGGCTCCCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAat < 2:1237147/76‑1 (MQ=255) cTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAat > 1:747330/1‑76 (MQ=255) cTACCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAat < 1:1734227/76‑1 (MQ=255) aCCGCAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAatat < 2:934183/76‑1 (MQ=255) aCCACAGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAatat > 1:1255226/1‑76 (MQ=255) cacaGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAATATCg < 2:1268399/76‑1 (MQ=255) aGCTGTGCGGCTCCCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAATATCgtgt < 1:1038402/76‑1 (MQ=255) aGCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAATATCgtgt > 1:805487/1‑76 (MQ=255) gCTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAATATCGTGTc > 2:213742/1‑76 (MQ=255) cTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAATATCGTGTcc > 1:1448770/1‑76 (MQ=255) cTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAATATCGTGTcc < 2:1165445/76‑1 (MQ=255) cTGTGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAATATCGTGTcc > 1:39633/1‑76 (MQ=255) | GTCATAACTAGGCGTGGAAACCAGCGCCAGCACCCCACCTGTACGCGGATCGGTGACTACCACAGCTGCGCGGCTACCCGCCAGCAGCGTTTCAATATATTGCTGGAGTTTGAGATCCAGCGTCAGGTAAATATCGTGTCC > REL606/649923‑650063 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |