| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 3,442,070 | A→G | 70.3% | 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 | 70.3% | 52.2 / 26.4 | 37 | T497A (ACG→GCG) | thiC | thiamin (pyrimidine moiety) biosynthesis protein |
| Reads supporting (aligned to +/- strand): ref base A (4/7); new base G (13/13); total (17/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.95e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01 | |||||||||||
CGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAT > W3110S.gb/3442011‑3442137 | cGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAagag < 1:2233662/71‑1 (MQ=255) tCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCa < 1:2133535/71‑1 (MQ=255) gATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTgg > 1:2125854/1‑71 (MQ=255) aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCa > 1:1576205/1‑65 (MQ=255) aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATc > 1:1353430/1‑67 (MQ=255) aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATc > 1:2051875/1‑67 (MQ=255) aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTggg > 1:2280549/1‑71 (MQ=255) aTTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTggg < 1:2136000/71‑1 (MQ=255) ttGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTgg < 1:368387/69‑1 (MQ=255) tGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGt > 1:2141194/1‑70 (MQ=255) tGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCa < 1:1075229/60‑1 (MQ=255) gatgatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAAt > 1:918917/1‑70 (MQ=255) gatTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATa > 1:1930075/1‑68 (MQ=255) atTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag > 1:371387/1‑46 (MQ=255) atTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag > 1:996512/1‑46 (MQ=255) tGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAaag > 1:613448/1‑68 (MQ=255) ggCTGGTTTGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAaag > 1:805581/1‑67 (MQ=255) gTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaagaa < 1:1748328/64‑1 (MQ=255) tGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAg < 1:1137740/70‑1 (MQ=255) tGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAg < 1:2250210/70‑1 (MQ=255) tGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAg < 1:808591/70‑1 (MQ=255) tGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAg < 1:573274/70‑1 (MQ=255) tGGCTGCGCGATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGCTAAg > 1:1391108/1‑70 (MQ=255) cgcgATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAgggg < 1:1434470/70‑1 (MQ=255) cgATGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCtt < 1:2182386/71‑1 (MQ=255) cgATGCTCTGTTACGTAACGCAAAAAGAGCATCTGTGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCt < 1:184339/70‑1 (MQ=255) aTGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAagaaga < 1:1493542/53‑1 (MQ=255) tGCTCTGTTACGTAGCGCCAAAAGAGCATCTGGGTCTg > 1:514772/1‑38 (MQ=255) gCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAg > 1:609070/1‑58 (MQ=255) ctctGTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCa > 1:712062/1‑70 (MQ=255) gTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCa < 1:463055/55‑1 (MQ=255) gTTACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCa < 1:418931/55‑1 (MQ=255) ttACGTAACGCCAAAAGAGCATCTTGGTCTGCCCAATAAAGAAGATGTTAAGc < 1:1137194/53‑1 (MQ=255) ttACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGc < 1:812877/53‑1 (MQ=255) ttACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGc < 1:2584777/53‑1 (MQ=255) tACGTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCtat > 1:1022377/1‑69 (MQ=255) gTAGCGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAt < 1:2780753/71‑1 (MQ=255) | CGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAGGGGCTTATCACCTATAAGAT > W3110S.gb/3442011‑3442137 |
| 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: ‑ |