| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 630,073 | C→T | 32.4% | G10S (GGC→AGC) | ycbC ← | conserved inner membrane protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 630,073 | 0 | C | T | 32.4% | 41.9 / 28.2 | 37 | G10S (GGC→AGC) | ycbC | conserved inner membrane protein |
| Reads supporting (aligned to +/- strand): ref base C (11/14); new base T (3/9); total (14/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.06e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGAT > minE/630011‑630137 | cagcagCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACtt > 1:2931963/1‑71 (MQ=255) agcagcGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACttt < 1:949606/71‑1 (MQ=255) agcGCCAGGCCAGCGCCAATAATCAACAACATCAGCTGAAGTGGCAGCAACATATTGCCAATCACtttttt < 1:3846338/71‑1 (MQ=255) agcGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACtttttt < 1:1908434/71‑1 (MQ=255) gcgcCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACtttttt > 1:926773/1‑70 (MQ=255) cgcCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACttttt < 1:3516445/68‑1 (MQ=255) ccaggccagCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCTAATCACttt < 1:3482569/64‑1 (MQ=255) aggccagCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCAc > 1:3639124/1‑59 (MQ=255) aataatCAACAACATCAGCGGAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCAt < 1:450117/69‑1 (MQ=255) aatCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGc > 1:351975/1‑64 (MQ=255) atCAACAACATCAGCGGAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCttt < 1:151036/70‑1 (MQ=255) tCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATcc > 1:3731859/1‑66 (MQ=255) aacaacATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTgg > 1:1152850/1‑71 (MQ=255) acaacaTCAGCGGAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGt < 1:3868632/71‑1 (MQ=255) tCAGCGGAAGTGGCAGCAACATATTGCTAATCACt > 1:2011610/1‑35 (MQ=255) aGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCttttt > 1:1555346/1‑71 (MQ=255) aGCGGAAGTGGCAGCAACATAATGCTAATCACTTTTTTCAGTGTAAAAAGCATCCttt > 1:1134721/1‑58 (MQ=255) cGGAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTa < 1:3911514/71‑1 (MQ=255) ggAAGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTc > 1:498505/1‑56 (MQ=255) aaGTGGCAGCAACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAAcc < 1:3661819/71‑1 (MQ=255) ggCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCata > 1:2209454/1‑70 (MQ=255) cagcaACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTc < 1:3880470/48‑1 (MQ=255) agcaACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAg > 1:1951184/1‑70 (MQ=255) agcaACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGc < 1:1480636/71‑1 (MQ=255) gcaACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCt < 1:2324379/47‑1 (MQ=255) gcaACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCa < 1:3902169/71‑1 (MQ=255) caACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCtttttt < 1:1327014/58‑1 (MQ=255) caACATATTGCTAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCtttttt < 1:3056141/58‑1 (MQ=255) aaCATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCa < 1:1260607/69‑1 (MQ=255) aCATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTATGGTTCCTTTTTTAACCATATAGCAAgg > 1:1807958/1‑71 (MQ=255) tatTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGAt < 1:1107602/71‑1 (MQ=255) atTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTcc < 1:116205/46‑1 (MQ=255) atTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTcc < 1:654067/46‑1 (MQ=255) atTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTcc < 1:986185/46‑1 (MQ=255) atTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTcc < 1:3489026/46‑1 (MQ=255) atTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTcc < 1:3456949/46‑1 (MQ=255) gCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTaa > 1:2615746/1‑51 (MQ=255) ccAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCtttttt < 1:3613618/48‑1 (MQ=255) | CAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTGGCAGCAACATATTGCCAATCACTTTTTTCAGTGTAAAAAGCATCCTTTCTGGTTCCTTTTTTAACCATATAGCAAGGGAT > minE/630011‑630137 |
| 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: ‑ |