| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 3,591,039 | G→A | 25.0% | intergenic (+28/+13) | yihG → / ← yihF | predicted endonuclease/conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 3,591,039 | 0 | G | A | 25.0% | 41.5 / 11.0 | 24 | intergenic (+28/+13) | yihG/yihF | predicted endonuclease/conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (7/11); new base A (3/3); total (10/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.65e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CTGACGTCTCTCATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTCATCTC > W3110S.gb/3590973‑3591102 | cTGACGTCTCTCATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCg < 1:1140131/69‑1 (MQ=255) tcATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTAttt > 1:557409/1‑70 (MQ=255) tcatcaCAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTAttt > 1:1126637/1‑65 (MQ=255) tcaCAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCgagtg < 1:1135054/55‑1 (MQ=255) tCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTcc < 1:978343/71‑1 (MQ=255) aGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACgag > 1:1002199/1‑43 (MQ=255) aGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACga > 1:1846326/1‑42 (MQ=255) gAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGc < 1:1111583/39‑1 (MQ=255) gAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGc < 1:460412/39‑1 (MQ=255) aTAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc < 1:2496237/71‑1 (MQ=255) taagtaagTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATg < 1:808110/53‑1 (MQ=255) taagtaagTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGta < 1:1368628/64‑2 (MQ=255) agtaagtTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCg > 1:2625070/1‑64 (MQ=255) agtTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACa > 1:1387617/1‑70 (MQ=255) gtTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATa < 1:2622429/71‑1 (MQ=255) gCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAttcttc < 1:2350113/70‑1 (MQ=255) tctGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAttcttcttc > 1:1225355/1‑70 (MQ=255) tGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAttc > 1:1999899/1‑62 (MQ=255) tGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCt > 1:734186/1‑71 (MQ=255) gAATGAGCAGAGGCGAGTGAGTATTTTAATGAA‑AAAACGTCCGGCACGAGACATaaa > 1:501181/1‑57 (MQ=255) aaTGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCt < 1:1433000/69‑1 (MQ=255) tGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTCAtc < 1:1443035/71‑1 (MQ=255) aGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTCAtctc < 1:2708024/71‑1 (MQ=255) aGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTCAtct < 1:574093/70‑1 (MQ=255) | CTGACGTCTCTCATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTCATCTC > W3110S.gb/3590973‑3591102 |
| 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: ‑ |