Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 703,662 | T→C | 8.9% | K124K (AAA→AAG) | nagB ← | glucosamine‑6‑phosphate deaminase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 703,662 | 0 | T | C | 8.9% | 118.3 / 4.1 | 45 | K124K (AAA→AAG) | nagB | glucosamine‑6‑phosphate deaminase |
Reads supporting (aligned to +/- strand): ref base T (21/20); new base C (2/2); total (23/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01 |
AAATGCAATATGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGTTGATGTTTTCT > W3110S.gb/703596‑703731 | aaaTGCAATATGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGAttttt < 1:2572259/70‑1 (MQ=255) tatGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGAtttt < 1:4243091/61‑1 (MQ=255) tatGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGAtttt < 1:4449367/61‑1 (MQ=255) tatGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATAc < 1:1005232/71‑1 (MQ=255) atGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATAc > 1:50870/1‑70 (MQ=255) tGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTg > 1:1106349/1‑71 (MQ=255) tGACCGTCGTTACCTACACCGCCCATAAACAGATCAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTg < 1:4124940/71‑1 (MQ=255) gACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATCTTttc > 1:2576694/1‑61 (MQ=255) tcgtTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCa < 1:241604/60‑1 (MQ=255) cgtTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCa > 1:1063760/1‑59 (MQ=255) gtTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACt > 1:2313007/1‑71 (MQ=255) tACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCTGCACTc < 1:9345/70‑1 (MQ=255) aCCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGAtttt > 1:2410687/1‑48 (MQ=255) aCCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATa > 1:2109303/1‑57 (MQ=255) aCCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCgg > 1:1516886/1‑71 (MQ=255) aCCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCgg > 1:3597836/1‑71 (MQ=255) aCCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATCTTTTCTTCATACTGGCGGCACTCgg > 1:3552935/1‑71 (MQ=255) cTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACt < 1:1530668/57‑1 (MQ=255) tACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCa < 1:2302147/52‑1 (MQ=255) acacCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGAtttttt > 1:3083516/1‑46 (MQ=255) acacCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTttct > 1:3450918/1‑48 (MQ=255) acacCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATa > 1:3961957/1‑53 (MQ=255) cacCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCa > 1:1055703/1‑50 (MQ=255) cacCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCa > 1:1423521/1‑50 (MQ=255) acCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGa < 1:789977/71‑1 (MQ=255) gCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATAc > 1:864063/1‑49 (MQ=255) ccATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATc < 1:2405670/70‑1 (MQ=255) ccATAAACAGATGAATTTTTCCGTAAGAACGGATCTTTTCTTCATACTGGCGGCACTCGGCGTCGATATc < 1:2629184/70‑1 (MQ=255) tAAACAGATGAATTTTTCCGTAAGAACGGAtttttt < 1:4538330/36‑1 (MQ=255) tAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATc < 1:45351/67‑1 (MQ=255) tAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCggg > 1:2836006/1‑71 (MQ=255) tAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCgg < 1:902098/70‑1 (MQ=255) aaaCAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATa > 1:867047/1‑42 (MQ=255) aaaCAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCgggg > 1:2767919/1‑71 (MQ=255) aaCAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGc > 1:4114768/1‑71 (MQ=255) aaCAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGAGGCACTCGGCGt > 1:1314094/1‑58 (MQ=255) aCAGATGAATTTTTCCGTAAGAACGGATTTTttct < 1:3334032/35‑1 (MQ=255) aGATGAATTTTTCCGTAAGAACGGATTTTttcttc > 1:4501037/1‑35 (MQ=255) aGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGtt < 1:19808/71‑1 (MQ=255) gATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATa < 1:1576736/37‑1 (MQ=255) aTGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGa < 1:1747911/56‑1 (MQ=255) tttCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGata < 1:1328191/51‑1 (MQ=255) tCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGa > 1:1310002/1‑47 (MQ=255) tCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATc < 1:183590/51‑1 (MQ=255) aaGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGtt > 1:521637/1‑70 (MQ=255) gAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGtt < 1:4528665/68‑1 (MQ=255) cGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGTTGATGtt < 1:739123/71‑1 (MQ=255) cGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGTTGATGtt < 1:3779968/71‑1 (MQ=255) cGGATCTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGTTGATGtt < 1:2152143/71‑1 (MQ=255) ggATTTTTTCTTCATACTGGCGGCACTCGGCGTCGata < 1:2533647/38‑1 (MQ=255) aTTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGt > 1:340320/1‑46 (MQ=255) ttttttCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGTTGATGTTTTct > 1:4339549/1‑71 (MQ=255) | AAATGCAATATGACCGTCGTTACCTACACCGCCCATAAACAGATGAATTTTTCCGTAAGAACGGATTTTTTCTTCATACTGGCGGCACTCGGCGTCGATATCCGGGGCGTTGCCGTTGAGAAGGTTGATGTTTTCT > W3110S.gb/703596‑703731 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |