| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 648,086 | T→G | E271A (GAG→GCG) | mrdB ← | cell wall shape‑determining protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 648,086 | 0 | T | G | 97.7% | 177.3 / ‑5.1 | 44 | E271A (GAG→GCG) | mrdB | cell wall shape‑determining protein |
| Reads supporting (aligned to +/- strand): ref base T (0/1); new base G (19/24); total (19/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01 | |||||||||||
TATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACTGAGTGCCGT > REL606/648016‑648158 | tATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTcc < 2:907884/76‑1 (MQ=255)tATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTcc < 2:21113/76‑1 (MQ=255) tCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCAGc > 2:1557391/1‑76 (MQ=255) aacacCCGGCGCAGGATCAGCAGAGTGAAGAGAGCGAGCAGAATCAGAACGCCGACTAATCCCAGCGCTCCCCCCa > 2:1843271/2‑76 (MQ=255) acaGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAg < 2:81052/76‑1 (MQ=255) acaGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAg < 1:1290035/76‑1 (MQ=255) ccGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACcg > 1:372496/1‑76 (MQ=255) ccGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACcg < 2:156906/76‑1 (MQ=255) aTGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGa > 1:528038/1‑76 (MQ=255) aTGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGa < 2:1693397/76‑1 (MQ=255) tGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGAt > 1:1024307/1‑76 (MQ=255) aTCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGCACCGCGAAGATaa > 1:648880/1‑76 (MQ=255) agcagAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGt < 2:510975/76‑1 (MQ=255) gcagAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTc < 1:236236/76‑1 (MQ=255) aaTGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCGGTa < 1:2128310/76‑1 (MQ=255) aaTGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTa < 2:724445/76‑1 (MQ=255) aTGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTAt < 2:142903/76‑1 (MQ=255) tGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATg > 1:1479083/1‑76 (MQ=255) tAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGc > 2:2045344/1‑76 (MQ=255) agagCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGtt < 2:568048/76‑1 (MQ=255) agagCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGtt < 2:1629465/76‑1 (MQ=255) agagCGAGCAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGtt > 2:1660553/1‑76 (MQ=255) gagcgagcAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTc < 1:1096012/76‑1 (MQ=255) gagcgagcAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTc < 2:2143299/76‑1 (MQ=255) gcgagcAGAATCAGAATGCCCACCAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCgg > 1:353968/1‑76 (MQ=255) cgagcACAATCAGAATGCCCACTAATCCCAGCGCTTTCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCggg > 2:1349909/1‑76 (MQ=255) agcAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTa < 2:528038/76‑1 (MQ=255) agcAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTa < 2:1852423/76‑1 (MQ=255) gcAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAg < 2:1893672/76‑1 (MQ=255) gcAGAATCAGAATGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAg < 1:882019/76‑1 (MQ=255) agaatcagaatGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGaa > 1:1735957/1‑76 (MQ=255) agaatGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCaa > 2:1192637/1‑76 (MQ=255) agaatGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCaa > 1:717878/1‑76 (MQ=255) gaatGCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAg < 1:596774/76‑1 (MQ=255) aatGCCCACTAATCCCAGCGCTTCCGCCCGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCCAGc > 2:36169/1‑76 (MQ=255) gCCCACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCtgt > 1:1510377/1‑76 (MQ=255) cccACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCtgtg < 1:283306/76‑1 (MQ=255) cACTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGAc > 2:1126010/1‑76 (MQ=255) aCTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACt < 1:506763/76‑1 (MQ=255) aCTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACt < 1:901246/76‑1 (MQ=255) cTAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACTg < 2:1593353/76‑1 (MQ=255) tAATCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACTGa > 2:2040903/1‑76 (MQ=255) tCCCAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACTGAGTg < 2:1899074/76‑1 (MQ=255) gcccAGCGCCTCCGCCAGTACCGCGAAGATAAAGTCGGCATGGCGTTCGGGTAGAAATTCAAGCTGTGACTGAGTg < 2:658807/75‑1 (MQ=255) cAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACTGAGTGCCg > 2:853954/1‑76 (MQ=255) cAGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACTGAGTGCCg > 1:1762549/1‑76 (MQ=255) aGCGCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACTGAGTGCCGt < 1:1228339/76‑1 (MQ=255) | TATCCACAGCCCGCGCATGATCAGCAGAATGTAGAGAGCGAGCAGAATCAGAATGCCCACTAATCCCAGCTCTTCCGCCAGTACCGCGAAGATAAAGTCGGTATGGCGTTCGGGTAGAAATTCAAGCTGTGACTGAGTGCCGT > REL606/648016‑648158 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |