Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,169,380 | T→C | 22.2% | intergenic (+50/+32) | bhsA → / ← ldtC | biofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,169,380 | 0 | T | C | 22.2% | 73.8 / 18.1 | 36 | intergenic (+50/+32) | bhsA/ldtC | biofilm, cell surface and signaling protein/L,D‑transpeptidase linking Lpp to murein |
Reads supporting (aligned to +/- strand): ref base T (15/13); new base C (4/4); total (19/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01 |
GCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATT > NC_000913/1169244‑1169527 | gCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgg < 1:1828420/149‑1 (MQ=255) gAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg < 3:88183/99‑1 (MQ=255) gAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg > 4:320280/1‑99 (MQ=255) gAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg < 3:320280/99‑1 (MQ=255) gAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACACGGCATGTCGCgg > 4:88183/1‑99 (MQ=255) ccTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACa > 4:1054961/1‑149 (MQ=255) cTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgg < 1:804759/92‑1 (MQ=255) cTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgg > 2:804759/1‑92 (MQ=255) agTAATTTATAAATAAGCATTAGCCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCggg < 4:1148976/79‑1 (MQ=255) agTAATTTATAAATAAGCATTAGCCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCggg > 3:1148976/1‑79 (MQ=255) agTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCg > 1:800916/1‑149 (MQ=255) aaTTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGaa < 2:800916/149‑1 (MQ=255) ttATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgg < 3:88600/72‑1 (MQ=255) ttATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCgg > 4:88600/1‑72 (MQ=255) ataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATc > 4:151968/1‑79 (MQ=255) ataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATc < 3:151968/79‑1 (MQ=255) ataaGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcc > 4:349048/1‑149 (MQ=255) aaCCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGt < 3:349048/149‑1 (MQ=255) cTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatga < 2:233910/100‑1 (MQ=255) cTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTtgatga > 1:233910/1‑100 (MQ=255) cTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAc > 2:1247655/1‑132 (MQ=255) cTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAc < 1:1247655/132‑1 (MQ=255) aTTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcc > 3:1037770/1‑133 (MQ=255) aTTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatcc < 4:1037770/133‑1 (MQ=255) ttAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCAc < 3:1054961/149‑1 (MQ=255) ccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCAc < 2:769830/61‑1 (MQ=255) ccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCAc > 1:769830/1‑61 (MQ=255) ccGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGCCGGTGGATGTccgccgtc > 4:922017/1‑84 (MQ=38) ccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccat > 4:1793349/1‑106 (MQ=255) ccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGAcatccat < 3:1793349/106‑1 (MQ=255) ccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt > 2:438990/1‑137 (MQ=255) ccGCGACACGGCATGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTtgt < 1:438990/137‑1 (MQ=255) cgcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaa > 4:1160178/1‑149 (MQ=255) gcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaaa > 3:195994/1‑149 (MQ=255) gcgACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTaaa > 4:826517/1‑149 (MQ=255) aTGCCGTGTCGCGGGGATTTTTTTATCCGGGCGTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGatt < 3:826517/149‑1 (MQ=255) | GCGGATGAGATGGGCGCAAAATCTTTCCGTATTACTTCTGTAACCGGTCCGAATACCCTCCATGGAACAGCAGTAATTTATAAATAAGCATTAACCCTCATTAATGCCTGCTACTGCTGATTTTTTCCCCGCGACATGCCGTGTCGCGGGGATTTTTTTATCCGGGCTTTACAGCGTTTGTGGGCTCACTTGATGACGGCGGACATCCACCGGCATCCCGGAACGGACATCCATCACATGTTGCATCACTTCAGCGTCAGTTTGTGCTGCATCTTTAAATGATT > NC_000913/1169244‑1169527 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |