| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,019,219:1 | +C | 17.2% | intergenic (+168/+58) | ydgA → / ← uidC | conserved hypothetical protein/predicted outer membrane porin protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,019,219 | 1 | . | C | 17.2% | 63.0 / 20.5 | 29 | intergenic (+168/+58) | ydgA/uidC | conserved hypothetical protein/predicted outer membrane porin protein |
| Reads supporting (aligned to +/- strand): ref base . (10/14); new base C (2/3); total (12/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 = 6.98e-01 | |||||||||||
CTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑ACGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCA > minE/1019151‑1019287 | gacgctgaaccgtcttatcatgcctacatagcactgccacgtatGTTTACACCGCATCCGGCataaaaaca < 1:2384417‑M2/27‑3 (MQ=255)cTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑Ac < 1:979563/71‑1 (MQ=255)cTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑Ac < 1:3191438/71‑1 (MQ=255) acgctgaagcgtcttatcatgcctacatagcactgccacgtatGTTTACACCGCATCCGGcataaaaacac < 1:1965463‑M2/28‑4 (MQ=255) cAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑Acgcg > 1:2160539/1‑71 (MQ=255) tgaagcgtcttatcatgcctacatagcactgccacgtatGTTTACACCGCATCCGGCATAAAAACAcgcg > 1:3902961‑M2/40‑70 (MQ=255) aagcgtcttatcatgcctacatagcactgccacgtatGTTTACACCGCATCCGGCATAAAAACACGCGCa < 1:526000‑M2/33‑1 (MQ=255) aagcgtcttatcatgcctacatagcactgccacgtatGTTTACACCGCATCCGGCATAAAAACACGCGCAc > 1:1828998‑M2/38‑71 (MQ=255) cTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑ACGCGCACtt < 1:988757/71‑1 (MQ=255) cTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑ACGCGCACtt < 1:1905708/71‑1 (MQ=255) tttGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑ACGCGCACttt > 1:414176/1‑71 (MQ=255) gtcttatcatgcctacatagcactgccacgtatGTTTACACCGCATCCGGCATAAAAACACGCGCACTTTg < 1:2519438‑M2/38‑1 (MQ=255) gAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑ACGCGCACTTTGt < 1:1289499/70‑1 (MQ=255) aaTTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑ACGCGCACTTTGTc < 1:3356061/70‑1 (MQ=255) gTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑ACGCGCACTTTGTCAACAATCTGa < 1:2413387/71‑1 (MQ=255) tttttGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑ACGCGCACTTTGTca < 1:1616320/61‑1 (MQ=255) ctacatagcactgccacgtatGTTTACACCGCATCCGGCATAAAAACACGCGCACtt < 1:2384928‑M2/36‑1 (MQ=255) tttGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑ACGCGCACTTTGTc > 1:1779547/1‑58 (MQ=255) tGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑ACGCGCACTTTGTCAACAATCTGAAAc < 1:2370010/69‑1 (MQ=255) tGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑ACGCGCACTTTGTCAACAATCTGAAACg < 1:2009320/70‑1 (MQ=255) tGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑ACGCGCACt > 1:604085/1‑44 (MQ=255) aaGGCGTTTACGCCGCATCCGGCATAAAAA‑ACGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTtc > 1:2418765/1‑71 (MQ=255) gTTTACGCCGCATCCGGCATAAACA‑ACGCGCACTTTGTCAGCAATCTGa > 1:3026504/1‑49 (MQ=11) ttACGCCGCATCCGGCATAAAAA‑ACGCGCACTTTGTCAACAATCTTAAACGCCGGAGATTTTTCTCTCCg > 1:3518056/1‑70 (MQ=255) tACGCCGCATCCGGCATAAAAA‑ACGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCgg < 1:3416868/70‑1 (MQ=255) tACGCCGCATCCGGCATAAAAA‑ACGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGc > 1:1989146/1‑71 (MQ=255) cgccgcATCCGGCATAAAAA‑ACGCGCACTTTGTCAACAATc > 1:1567664/1‑41 (MQ=255) gccgcATCCGGCATAAAAA‑ACGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGt > 1:2931215/1‑70 (MQ=255) ggCATAAAAA‑ACGCGCACTTTCTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTAc < 1:343755/71‑1 (MQ=255) cATAAAAA‑ACGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACtt < 1:2617863/71‑1 (MQ=255) tAAAAA‑ACGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCa < 1:459675/71‑1 (MQ=255) aaaa‑aCGCGCACTTTGTCAACAATCTGAAACGCCGGAGAttt < 1:1356792/42‑1 (MQ=255) aaaa‑aCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAg > 1:1383759/1‑70 (MQ=255) aaa‑aCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCa > 1:928691/1‑71 (MQ=255) | CTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAA‑ACGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCA > minE/1019151‑1019287 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |