| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 376,160 | G→T | 100% | Q55H (CAG→CAT) | yaiL → | nucleoprotein/polynucleotide‑associated enzyme |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 376,160 | 0 | G | T | 84.2% | 28.1 / 5.0 | 19 | Q55H (CAG→CAT) | yaiL | nucleoprotein/polynucleotide‑associated enzyme |
| Reads supporting (aligned to +/- strand): ref base G (3/0); new base T (10/6); total (13/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.17e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CAGGCGCGTGAAGCTCGGGCGGCGGTAGAAGAAAATAAAAAGGCACAGCTTGAGCGTGATAAACAGCTTAGCGAACAGCAAAAACAAGCGGCGTTGGCGAAAGAATATAAAGCTCAGGTGAAGCAGCTTATTG > W3110S.gb/376095‑376227 | cAGGCGCGTGAAGCTCGGGCGGCGGTAGAAGAAAATAAAAAGGCACAGCTTGAGCGTGATAAACATCTTAg > 1:856396/1‑71 (MQ=255) ggCGCGTGAAGCTCGGGCGGCGGTAGAAGAAAATAAAAAGGCACAGCTTGAGCGTGATAAACATCTTAGCg < 1:375410/71‑1 (MQ=255) gCTCGGGAGGCGGTAGAAGAAAATAAAAAGGCA‑AGCTTGAGCGTGATAAACATCTTa > 1:2795003/1‑57 (MQ=38) gggCGGCGGTAGAAGAAAATAAAAAGGCACAGCTTGAGCGTGATAAACATCTTAGCGAACAGCAAAAACa < 1:1403205/70‑1 (MQ=255) gaagaaAATAAAAAGGCACAGCTTGAGCGTGATAAACATCTTAGCGAACAGCAAAAACAAGCGGCGt < 1:2030135/67‑1 (MQ=255) gaagaaAATAAAAAGGCACAGCTTGAGCGTGATAAACATCTTAGCGAACAGCAAAAACAAGCGGCGt < 1:600273/67‑1 (MQ=255) gaaAATAAAAAGGCACAGCTTGAGCGTGATAAACATCTTAGCGAACAGCAAAAACAAGCGGCGTTGGCg > 1:2340337/1‑69 (MQ=255) aaaTAAAAAGGCACAGCTTGAGCGTGATAAACATCTTAGCGAACAGCAAAAACAAGCGGCGTTGGCGAAAg < 1:2273794/71‑1 (MQ=255) aaTAAAAAGGCACAGCTTGAGCGTGATAAACATCTTAGCGAACAGCAAAAACAAGCGGCGTTGGCGAAAGa < 1:1443768/71‑1 (MQ=255) aTAAAAAGGCACAGCTTGAGCGTGATAAACATCTTAGCGAACAGCAAAAACAAGCGGCGTTg > 1:86593/1‑62 (MQ=255) aTAAAAAGGCACAGCTTGAGCGTGATAAACATCTTAGCGAACAGCAAAAAAAAgcggctt > 1:695002/1‑58 (MQ=255) tAAAAAGGCACAGCTTGAGCGTGATAAACATCTTAGCGaa > 1:2572824/1‑40 (MQ=255) tAAAAAGGCACAGCTTGAGCGTGATAAACATCTTAGCGaa > 1:963736/1‑40 (MQ=255) aTAAACAGCTTAGCGAACAGCAAAAACAAGCGGCGTTg > 1:1174164/1‑38 (MQ=255) aaaCATCTTAGCGAACAGCAAAAACAAGCGGCGTTGGCGAAAGAata > 1:966723/1‑47 (MQ=255) aaaCATCTTAGCGAACAGCAAAAACAAGCGGCGTTGGCGAAAGAATATAAAGCTCAGGTGAAGCAGCTTAt > 1:335469/1‑71 (MQ=255) aaaCATCTTAGCGAACAGCAAAAACAAGCGGCGTTGGCGAAAGAATATAAAGCTCAGGTGAAGCAGCTTAt > 1:2117272/1‑71 (MQ=255) aaaCAGCTTAGCGAACAGCAAAAACAAGCGGCGTTGGCGAAAGaa > 1:592713/1‑45 (MQ=255) cAGCTTAGCGAACAGCAAAAACAAGCGGCGTTGGCGAAAGAATATAAAGCTCAGGTGAAGCAGCTTATTg > 1:280975/1‑70 (MQ=255) | CAGGCGCGTGAAGCTCGGGCGGCGGTAGAAGAAAATAAAAAGGCACAGCTTGAGCGTGATAAACAGCTTAGCGAACAGCAAAAACAAGCGGCGTTGGCGAAAGAATATAAAGCTCAGGTGAAGCAGCTTATTG > W3110S.gb/376095‑376227 |
| 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: ‑ |