Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 3,591,039 | G→A | 37.8% | 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 | 37.8% | 29.8 / 36.4 | 37 | intergenic (+28/+13) | yihG/yihF | predicted endonuclease/conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (7/16); new base A (6/8); total (13/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.95e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01 |
CTCTCATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTCATCT > W3110S.gb/3590980‑3591101 | ctctcATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTAt > 1:899026/1‑71 (MQ=255) ctcATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTAtt > 1:477650/1‑70 (MQ=255) aTGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTaa > 1:2104192/1‑71 (MQ=255) aTGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTAtttt < 1:480256/69‑1 (MQ=255) aTGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTAttt < 1:878017/68‑1 (MQ=255) aTGTCATCACAGCGGCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTAtttt < 1:4513262/69‑1 (MQ=255) tcatcaCAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGc < 1:1125974/53‑1 (MQ=255) tcatcaCAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGc < 1:3459756/53‑1 (MQ=255) cacaGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACa > 1:870375/1‑70 (MQ=255) acaGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACaa < 1:814973/70‑1 (MQ=255) gTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTc < 1:1043231/71‑1 (MQ=255) tCAGAATAAGTAAGTTTTACTCGCCTCTGAATGTGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGt < 1:2492807/69‑1 (MQ=255) tCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACaa < 1:4193990/64‑1 (MQ=255) gAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACg < 1:528039/65‑1 (MQ=255) aaTAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTaa < 1:4002230/53‑1 (MQ=255) aTAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc > 1:544306/1‑71 (MQ=255) aTAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc < 1:3013447/71‑1 (MQ=255) taagtaagTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAAc > 1:2439206/1‑61 (MQ=255) taagtaagTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc < 1:4364360/70‑1 (MQ=255) taagtaagTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc > 1:3878501/1‑70 (MQ=255) agtTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACAt < 1:55315/71‑1 (MQ=255) ttttACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATg < 1:276667/45‑1 (MQ=255) ttACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATaaa < 1:2929907/70‑1 (MQ=255) tACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAtt > 1:4228644/1‑71 (MQ=255) cTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAttct < 1:4065041/71‑1 (MQ=255) cTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAttct < 1:3566755/71‑1 (MQ=255) cTCGCCTCTGAATGACCAGAGACGAGTGAGTATTTTAATGAACaa < 1:3636750/45‑1 (MQ=255) gCCTCTGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAttctt < 1:1101047/69‑1 (MQ=255) ccTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAttcttctt < 1:1983888/71‑1 (MQ=255) tctGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAAttcttcttc > 1:921021/1‑70 (MQ=255) tGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCt > 1:3975274/1‑71 (MQ=255) tGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAAc > 1:2674149/1‑40 (MQ=255) tGAATGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCg < 1:221287/69‑1 (MQ=255) aaTGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTc < 1:26249/70‑1 (MQ=255) aaTGAGCAGAGACGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTCa > 1:1412125/1‑71 (MQ=255) gAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGGTCAtct > 1:277026/1‑71 (MQ=255) agGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCAc < 1:1190096/38‑1 (MQ=255) | CTCTCATGTCATCACAGCGTCAGAATAAGTAAGTTTTACTCGCCTCTGAATGAGCAGAGGCGAGTGAGTATTTTAATGAACAAAACGTCCGGCACGAGACATAAATTCTTCTTCGCTCATCT > W3110S.gb/3590980‑3591101 |
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: ‑ |