| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,695,611 | T→C | 64.2% | intergenic (+104/+62) | nudE → / ← mrcA | ADP‑ribose diphosphatase/fused penicillin‑binding protein 1a murein transglycosylase and murein transpeptidase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,695,611 | 0 | T | C | 64.2% | 1.8 / 19.1 | 39 | intergenic (+104/+62) | nudE/mrcA | ADP‑ribose diphosphatase/fused penicillin‑binding protein 1a murein transglycosylase and murein transpeptidase |
| Reads supporting (aligned to +/- strand): ref base T (6/8); new base C (9/16); total (15/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.40e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGGTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGCAAACTGAAAAGGCGCCGAAGCGCCTTTTTAATCAGAAC > minE/2695551‑2695679 | gggTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGc < 1:1098527/68‑1 (MQ=255)gggTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGc < 1:208468/68‑1 (MQ=255)gggTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGc < 1:1315646/68‑1 (MQ=255)gggTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGc < 1:533532/68‑1 (MQ=255) gTTTCTGCAATTTATTGTGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATa < 1:2572830/69‑1 (MQ=255) gTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATa > 1:2350435/1‑69 (MQ=255) tCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATaa > 1:887537/1‑67 (MQ=255) tCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATaa > 1:2825651/1‑67 (MQ=255) tCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATaa > 1:2390426/1‑67 (MQ=255) tCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATaa > 1:1350806/1‑67 (MQ=255) ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCAC‑CCGCATCTGGCATaa > 1:2081369/2‑66 (MQ=2) ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCAC‑CCGCATCTGGCATaa > 1:778255/2‑66 (MQ=2) ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCAC‑CCGCATCTGGCATaa > 1:2174877/2‑66 (MQ=2) ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCAC‑CCGCATCTGGCATaa > 1:813354/2‑66 (MQ=2) ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCAC‑CCGCATCTGGCATaa > 1:606870/2‑66 (MQ=2) ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCAC‑CCGCATCTGGCATaa > 1:2405166/2‑66 (MQ=2) ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCAC‑CCGCATCTGGCATaa > 1:1046574/2‑66 (MQ=2) ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCAC‑CCGCATCTGGCATaa > 1:2999090/2‑66 (MQ=2) ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCAC‑CCGCATCTGGCATaa > 1:3129678/2‑66 (MQ=2) cACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGCaa < 1:2259720/68‑1 (MQ=255) cACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGCaa < 1:1071653/68‑1 (MQ=255) gATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGCaa < 1:2509534/65‑1 (MQ=255) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:85103/40‑1 (MQ=21) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:2415982/40‑1 (MQ=21) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:2956307/40‑1 (MQ=21) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:3377377/40‑1 (MQ=21) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:407320/40‑1 (MQ=21) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:536930/40‑1 (MQ=21) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:2815865/40‑1 (MQ=21) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:2778066/40‑1 (MQ=21) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:229493/40‑1 (MQ=21) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:1830146/40‑1 (MQ=21) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:1558139/40‑1 (MQ=21) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:1305918/40‑1 (MQ=21) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:1156988/40‑1 (MQ=21) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:1111050/40‑1 (MQ=21) tAACGCGTTCAC‑CCGCATCTGGCATAAACAAGTGCACttt < 1:1054931/40‑1 (MQ=21) tAACGCGTTCA‑ACCGCATCTGGCATAAACAAGTGCACttt < 1:2775946/40‑1 (MQ=19) cACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGAAAACTGAAAAGGCGCCGAAGCGCCTTTTTa > 1:2753902/1‑68 (MQ=255) tATCCGGCATAAACAAGTGCACTTTGTCAGCAAACTGAAAAGGCGCCGAAGCGCCTTTTTAATCAGAAc < 1:1399803/69‑1 (MQ=255) | GGGTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGCAAACTGAAAAGGCGCCGAAGCGCCTTTTTAATCAGAAC > minE/2695551‑2695679 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |