| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 429,850 | T→C | 10.2% | Y8H (TAT→CAT) | yajD → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 429,850 | 0 | T | C | 10.2% | 102.8 / 4.7 | 39 | Y8H (TAT→CAT) | yajD | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base T (19/16); new base C (2/2); total (21/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTT > W3110S.gb/429782‑429919 | tCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACCAt > 1:1279875/1‑71 (MQ=255) cgcgTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGa > 1:615350/1‑69 (MQ=255) cgcgTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGa < 1:3644221/69‑1 (MQ=255) gcgTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGAAAg < 1:287686/71‑1 (MQ=255) cgTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGc < 1:1408397/71‑1 (MQ=255) aCCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGc < 1:910744/70‑1 (MQ=255) aCCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGc < 1:2775210/70‑1 (MQ=255) cccTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTa > 1:208068/1‑71 (MQ=255) cTAGATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGa < 1:2830779/60‑1 (MQ=255) cTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGc > 1:1340125/1‑67 (MQ=255) cTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATc < 1:685337/71‑1 (MQ=255) aCATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTa > 1:186915/1‑46 (MQ=255) cATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGAGCgg > 1:3582688/1‑52 (MQ=255) cAGGTAACTTTTCATGGCTATCATCCCAAAAAACCATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCAt > 1:432862/1‑71 (MQ=255) aGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTAt > 1:3736451/1‑57 (MQ=255) aGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCAtt < 1:2349259/71‑1 (MQ=255) aCTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaaaa < 1:2695387/70‑1 (MQ=255) tttCATGGCTATCATCCCAAAAAACTATGCGCGGTTa > 1:431513/1‑37 (MQ=255) tttCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATc > 1:2501066/1‑70 (MQ=255) ttCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTa > 1:118638/1‑71 (MQ=255) ggCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATc > 1:4367590/1‑64 (MQ=255) gCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGt < 1:3796396/70‑1 (MQ=255) gCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTg > 1:537231/1‑71 (MQ=255) atcatcCCAAAAAACTATGCGCGGTTAGAAAGCGGc < 1:2301180/36‑1 (MQ=255) atcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGaaa > 1:2300896/1‑43 (MQ=255) atcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTaa > 1:730222/1‑52 (MQ=255) tcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTg > 1:1777155/1‑70 (MQ=255) ccAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCg < 1:1265398/70‑1 (MQ=255) cAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAg > 1:497685/1‑42 (MQ=255) cAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCg > 1:1156934/1‑69 (MQ=255) aaaaaaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGt < 1:49880/70‑1 (MQ=255) aaaaaaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCCATCCGTGGGTCTGCGGTc > 1:381118/1‑71 (MQ=255) aaaaaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCg < 1:4382379/71‑1 (MQ=255) aaaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCt > 1:2526193/1‑71 (MQ=255) aaaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCt < 1:2493697/71‑1 (MQ=255) aaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCg > 1:62687/1‑64 (MQ=255) aaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTg < 1:2016122/71‑1 (MQ=255) aCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGt < 1:3384799/71‑1 (MQ=255) cTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGtt < 1:4534586/71‑1 (MQ=255) | TCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTT > W3110S.gb/429782‑429919 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |