Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 649,391 | T→G | I471L (ATC→CTC) | 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,391 | 0 | T | G | 96.4% | 456.5 / 2.2 | 110 | I471L (ATC→CTC) | mrdA | transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2) |
Reads supporting (aligned to +/- strand): ref base T (4/0); new base G (57/49); total (61/49) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.27e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01 |
Only 100 of 125 total aligned reads displayed.
GCACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGGTGTC > REL606/649319‑649463 | gccccTTTTTGCTGTCTTCGGCGGTGCTCATCAGCAAACGCGGAACCTTCACGATACCGTCATTAATCAGGagcat < 1:2576069/73‑1 (MQ=255) gCACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGagcat > 2:3170885/1‑76 (MQ=255) gCACCTGTTTGCCGTCTGCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGagcat < 1:409239/76‑1 (MQ=255) cACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGagcatc < 2:4768211/76‑1 (MQ=255) cccTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAgcatca < 1:3387175/75‑1 (MQ=255) aCCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAgcatca < 2:2565857/76‑1 (MQ=255) aCCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAgcatca < 1:2868274/76‑1 (MQ=255) aCCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAgcatca < 1:2473519/76‑1 (MQ=255) aCCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAgcatca > 2:3973527/1‑76 (MQ=255) cTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGt < 1:100421/76‑1 (MQ=255) cTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGt > 1:916705/1‑76 (MQ=255) gTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGc > 1:3235590/1‑76 (MQ=255) gTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGc > 2:4364119/1‑76 (MQ=255) gTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGc > 1:4593511/1‑76 (MQ=255) tttGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGcc > 1:291389/1‑76 (MQ=255) tGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCtt > 2:2933462/1‑76 (MQ=255) tGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCtt < 2:4466152/76‑1 (MQ=255) gCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTa < 2:3744283/76‑1 (MQ=255) ccGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTGACGATACCGTCATTAATCAGGAGCATCAGTGCCTTAc < 2:322391/76‑1 (MQ=255) gTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTc > 1:4293199/1‑76 (MQ=255) tCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCa < 1:3275600/76‑1 (MQ=255) cTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCAt > 1:1295516/1‑76 (MQ=255) cTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCAt > 2:2448352/1‑76 (MQ=255) cTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCAt > 2:468176/1‑76 (MQ=255) ttCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATc > 1:1254456/1‑76 (MQ=255) ttCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATc > 2:361252/1‑76 (MQ=255) tCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCt < 2:85781/76‑1 (MQ=255) tCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCt < 1:3711275/76‑1 (MQ=255) ggcggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTgg < 1:3681416/76‑1 (MQ=255) ggcggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTgg > 2:3738864/1‑76 (MQ=255) gcggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGa > 2:4625243/1‑76 (MQ=255) cggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGAt > 2:4512850/1‑76 (MQ=255) ggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGAtt < 2:2538870/76‑1 (MQ=255) tGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTgg > 2:2753300/1‑76 (MQ=255) gCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTggg > 1:1008362/1‑76 (MQ=255) tcatcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGt < 1:637741/76‑1 (MQ=255) catcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTc < 2:4692374/76‑1 (MQ=255) catcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTc > 2:4112189/1‑76 (MQ=255) catcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTc < 2:13048/76‑1 (MQ=255) catcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTc > 2:1754770/1‑76 (MQ=255) catcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTc < 2:2592361/76‑1 (MQ=255) atcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCg > 2:1677070/1‑76 (MQ=255) acaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGc < 2:376357/75‑1 (MQ=255) caGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCGCt > 1:704595/1‑76 (MQ=255) caGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCt < 2:3909470/76‑1 (MQ=255) aGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTg < 1:3418279/76‑1 (MQ=255) cAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTc > 1:621056/1‑76 (MQ=255) cAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTc > 2:1698262/1‑76 (MQ=255) aaaTGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTcc < 2:1173111/76‑1 (MQ=255) aaaTGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTcc < 2:752225/76‑1 (MQ=255) aaaTGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTcc > 2:644673/1‑76 (MQ=255) aaTGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCa > 2:2208282/1‑76 (MQ=255) aaTGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCa < 2:4009558/76‑1 (MQ=255) tGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGt > 1:3204954/1‑76 (MQ=255) tGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGt < 1:1570378/76‑1 (MQ=255) gCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTa < 1:5013671/76‑1 (MQ=255) gCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTa > 1:1476532/1‑76 (MQ=255) gCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTa < 2:2478024/76‑1 (MQ=255) cGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTaa < 1:1534044/76‑1 (MQ=255) aggAACCTCCACGATACCGTCATTAATCAGAAGCATCAGTGCCTTACTCATCTGGATTGGGGTCTCTGTCCAGTaa < 2:1642918/75‑1 (MQ=255) ggAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAAc > 1:1697264/1‑76 (MQ=255) aaCCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAAccc > 1:1055268/1‑76 (MQ=255) aCCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCt < 1:1743696/76‑1 (MQ=255) aCCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCt < 1:4611816/76‑1 (MQ=255) aCCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCt < 1:4618532/76‑1 (MQ=255) ccTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTg < 1:4901635/76‑1 (MQ=255) ccTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTg < 2:1883742/76‑1 (MQ=255) cTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGa > 1:1322486/1‑76 (MQ=255) cTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGa > 1:4012166/1‑76 (MQ=255) cACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTCACCCTGACCg > 2:2882684/1‑76 (MQ=255) aCGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGa > 1:4808882/1‑76 (MQ=255) aCGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGa > 1:3703436/1‑76 (MQ=255) gATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATa < 1:1195731/76‑1 (MQ=255) aTACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATAc > 2:246312/1‑76 (MQ=255) aTACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATAc > 2:3461755/1‑76 (MQ=255) aTACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATAc > 1:4005139/1‑76 (MQ=255) tACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATAcc > 1:1001117/1‑76 (MQ=255) aCCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa > 1:1497489/1‑76 (MQ=255) aCCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa > 1:4327725/1‑76 (MQ=255) aCCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa < 1:3868422/76‑1 (MQ=255) aCCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa < 1:3011245/76‑1 (MQ=255) aCCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa < 1:2358952/76‑1 (MQ=255) aCCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa < 1:4773787/76‑1 (MQ=255) aCCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa < 2:1867920/76‑1 (MQ=255) aCCGTCATTAATCAGGAGCAGCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa < 1:4401289/76‑1 (MQ=255) ccGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCaa > 1:836963/1‑76 (MQ=255) cGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAAc > 1:3995252/1‑76 (MQ=255) tCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCg < 2:2972204/76‑1 (MQ=255) tCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCg > 1:2445147/1‑76 (MQ=255) aTTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGa > 1:2179974/1‑76 (MQ=255) ttAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGaa > 2:1501511/1‑76 (MQ=255) ttAATCAGGAGCATCAGTGCCTTACTCATCCGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGaa > 2:797274/1‑76 (MQ=255) tAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAAt < 1:3078600/76‑1 (MQ=255) tAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAAt > 2:4748754/1‑76 (MQ=255) tAAGCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAAt < 2:1604433/76‑1 (MQ=255) tCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGgt > 1:4671906/1‑76 (MQ=255) cAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGgtg > 2:3095365/1‑76 (MQ=255) aGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGgtgt < 2:1103349/76‑1 (MQ=255) aGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGgtgt < 2:3306768/76‑1 (MQ=255) ggAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGGTGTc < 1:2240570/76‑1 (MQ=255) | GCACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGGTGTC > REL606/649319‑649463 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |