Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 3,442,070 | A→G | 18.7% | T497A (ACG→GCG) | thiC → | thiamin (pyrimidine moiety) biosynthesis protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 3,442,070 | 0 | A | G | 18.7% | 65.7 / 11.4 | 32 | T497A (ACG→GCG) | thiC | thiamin (pyrimidine moiety) biosynthesis protein |
Reads supporting (aligned to +/- strand): ref base A (17/9); new base G (3/3); total (20/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.47e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01 |
TCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATT > W3110S.gb/3442008‑3442138 | tCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGcc > 1:2115002/1‑67 (MQ=255) tCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa < 1:1439559/69‑1 (MQ=255) taCGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGcc > 1:293435/2‑66 (MQ=255) gatgTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCaa < 1:2354276/65‑1 (MQ=255) cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGcc > 1:2498983/1‑66 (MQ=255) ggggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCa < 1:721783/69‑1 (MQ=255) ggggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCa < 1:3627838/69‑1 (MQ=255) gATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGc > 1:2400490/1‑65 (MQ=255) aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTg > 1:909084/1‑69 (MQ=255) ttGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCAt < 1:304969/65‑1 (MQ=255) ttGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTgg > 1:1225045/1‑69 (MQ=255) tGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCAt < 1:2063564/64‑1 (MQ=255) gatgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGc < 1:506604/53‑1 (MQ=255) gatgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGc < 1:3325976/53‑1 (MQ=255) gatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCaa > 1:3690572/1‑66 (MQ=255) atTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaa > 1:2305598/1‑68 (MQ=255) atTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaa > 1:1234553/1‑68 (MQ=255) atTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaa > 1:1419821/1‑68 (MQ=255) atTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaa > 1:1335763/1‑68 (MQ=255) atTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaa > 1:376368/1‑68 (MQ=255) atTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaa > 1:3071428/1‑68 (MQ=255) tGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa < 1:228978/37‑1 (MQ=255) ggTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTggg < 1:2775701/49‑1 (MQ=255) tGGCTGCGCGATTCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaag < 1:1437501/59‑1 (MQ=255) gcgcGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgg > 1:2018475/1‑69 (MQ=255) gcgcGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgg > 1:1502579/1‑69 (MQ=255) tGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCtt < 1:1272152/63‑1 (MQ=255) tGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCAcc > 1:3135282/1‑69 (MQ=255) ttACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCt > 1:1794457/1‑68 (MQ=255) tACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGAt > 1:2031740/1‑45 (MQ=255) tACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCa > 1:3537329/1‑64 (MQ=255) tAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGa > 1:2636603/1‑69 (MQ=255) aCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAtt < 1:3536058/69‑1 (MQ=255) aCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAtt < 1:783755/69‑1 (MQ=255) | TCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGATT > W3110S.gb/3442008‑3442138 |
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: ‑ |