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 | 97.8% | 182.2 / ‑3.8 | 46 | I401V (ATC→GTC) | mrdA | transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2) |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (21/24); total (22/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.78e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01 |
GTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACG > REL606/649529‑649673 | gTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGAcccc < 2:1280535/76‑1 (MQ=255) tcttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCa > 1:1553111/1‑76 (MQ=255) cttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCat < 1:1093402/76‑1 (MQ=255) cttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCat < 2:2137844/76‑1 (MQ=255) ccGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATc < 2:1721604/76‑1 (MQ=255) ccGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATc > 1:461254/1‑76 (MQ=255) cGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCg < 1:492810/76‑1 (MQ=255) ccAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTa > 1:543269/1‑76 (MQ=255) cAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAg > 2:61194/1‑76 (MQ=255) cAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAg < 2:2126707/76‑1 (MQ=255) aGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAgg > 2:596623/1‑76 (MQ=255) gTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGcc < 1:1963246/76‑1 (MQ=255) tCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTATGccg > 1:2074395/1‑75 (MQ=255) tCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCa < 2:1006094/76‑1 (MQ=255) cGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCAc > 2:1947342/1‑76 (MQ=255) tACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTg > 1:1486258/1‑76 (MQ=255) aCCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGa < 1:1930659/76‑1 (MQ=255) aCCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGa > 1:248587/1‑76 (MQ=255) ccGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGAt < 2:876909/76‑1 (MQ=255) gtgtAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATaga > 1:2107312/1‑76 (MQ=255) gtgtAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATaga < 1:197071/76‑1 (MQ=255) tAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATagacga > 1:455899/1‑76 (MQ=255) tAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATagaaga > 1:1067419/1‑76 (MQ=255) aaTGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAgaagaa > 2:1549007/1‑76 (MQ=255) aaTGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAgaagaa > 1:965082/1‑76 (MQ=255) aTGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCCGATAGaagaag > 2:729859/1‑76 (MQ=255) gACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGt < 2:1375493/76‑1 (MQ=255) ccATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTAt < 1:1728701/76‑1 (MQ=255) ccATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTAt < 2:374279/76‑1 (MQ=255) ccATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTAt < 1:1539275/76‑1 (MQ=255) tAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCg < 2:764102/76‑1 (MQ=255) tAACCGAATTTACCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCg < 2:1198512/76‑1 (MQ=255) aCCCATCCATTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTcc > 2:854334/1‑76 (MQ=255) catccatTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGc < 2:1426044/76‑1 (MQ=255) atccatTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCg < 1:578985/76‑1 (MQ=255) catTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATc > 2:2184407/1‑76 (MQ=255) atTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGACGGTATCCGCAGATTCTTCCAGCGATCt > 2:1196926/1‑76 (MQ=255) atTCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCt > 2:163612/1‑76 (MQ=255) ttCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCtt < 2:1825796/76‑1 (MQ=255) tCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTtg < 1:1628642/76‑1 (MQ=255) tCGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTtg < 1:1744246/76‑1 (MQ=255) cGGAGAGGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTtgt > 1:550812/1‑76 (MQ=255) agagGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGAc > 2:605599/1‑76 (MQ=255) gagGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACa < 1:118647/76‑1 (MQ=255) gagGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACa < 2:1775602/76‑1 (MQ=255) gagGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACa < 1:2026177/76‑1 (MQ=255) agGCGATCGACCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACAt > 2:953869/1‑76 (MQ=255) cgatcgacCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCa > 2:1531067/1‑76 (MQ=255) cgatcgacCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCa > 1:719690/1‑76 (MQ=255) cgatcgacCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCa < 1:256874/76‑1 (MQ=255) gatcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAg < 1:1432489/76‑1 (MQ=255) gatcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAg > 1:1342631/1‑76 (MQ=255) gatcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAg > 2:850706/1‑76 (MQ=255) atcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGa < 2:1938707/76‑1 (MQ=255) atcgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGa < 2:1389325/76‑1 (MQ=255) cgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACg < 2:1749323/76‑1 (MQ=255) cgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACg > 2:1600299/1‑76 (MQ=255) cgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACg < 1:1631208/76‑1 (MQ=255) cgaccCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACg > 1:2057728/1‑76 (MQ=255) | GTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTCTTCCAGCGATCTTGTGACATTCAGACG > REL606/649529‑649673 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |