Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,563,690 | C→T | 8.8% | intergenic (+172/+97) | glpD → / ← yzgL | sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding/pseudogene, periplasmic solute binding protein homology |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,563,690 | 0 | C | T | 8.8% | 208.6 / 9.7 | 68 | intergenic (+172/+97) | glpD/yzgL | sn‑glycerol‑3‑phosphate dehydrogenase, aerobic, FAD/NAD(P)‑binding/pseudogene, periplasmic solute binding protein homology |
Reads supporting (aligned to +/- strand): ref base C (32/30); new base T (3/3); total (35/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01 |
CTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAAT > NC_000913/3563592‑3563837 | cTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGTGCGCTTTGttca > 4:928961/1‑108 (MQ=255) cTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGTGCGCTTTGttca < 3:928961/111‑4 (MQ=255) tACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGt > 4:381371/1‑107 (MQ=255) tACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGt < 3:381371/107‑1 (MQ=255) tATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTa > 1:1041614/1‑111 (MQ=255) tATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTa > 4:486829/1‑111 (MQ=255) tATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTa < 3:486829/111‑1 (MQ=255) tATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTa < 2:1041614/111‑1 (MQ=255) tCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACt < 2:357223/95‑1 (MQ=255) tCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACt > 1:357223/1‑95 (MQ=255) cAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACttt < 2:737589/117‑1 (MQ=255) cAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACttt > 1:737589/1‑117 (MQ=255) cAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGc > 3:703574/1‑148 (MQ=255) tACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGt < 2:222269/92‑1 (MQ=255) tACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGt > 1:222269/1‑92 (MQ=255) tACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTc > 1:68856/1‑118 (MQ=255) tACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTc < 2:68856/118‑1 (MQ=255) tGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAAt > 1:939494/1‑112 (MQ=255) tGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAAt < 2:939494/112‑1 (MQ=255) atatTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTc > 4:763617/1‑101 (MQ=255) atatTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTc < 3:763617/101‑1 (MQ=255) ttGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCa < 4:703574/148‑1 (MQ=255) tGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCt < 3:263638/85‑1 (MQ=255) tGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCt > 4:263638/1‑85 (MQ=255) tGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTAc < 1:705327/107‑1 (MQ=255) tGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTAc > 2:705327/1‑107 (MQ=255) tttGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCa < 1:389152/94‑1 (MQ=255) tttGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCa > 2:389152/1‑94 (MQ=255) tGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATAcc < 3:425976/120‑1 (MQ=255) tGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATAcc > 4:425976/1‑120 (MQ=255) cTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCa > 1:798028/1‑84 (MQ=255) cTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCa < 2:798028/84‑1 (MQ=255) ccAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACtt > 3:284921/1‑46 (MQ=255) ccAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACtt < 4:284921/46‑1 (MQ=255) aGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAg > 1:244989/1‑75 (MQ=255) aGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAg < 2:244989/75‑1 (MQ=255) aGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCACCGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCgaga > 4:686248/1‑148 (MQ=255) aaGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAAtt < 4:438698/68‑1 (MQ=255) aaGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAAtt > 3:438698/1‑68 (MQ=255) cGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGt < 1:39026/104‑1 (MQ=255) cGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGt > 2:39026/1‑104 (MQ=255) ttCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGtt > 2:352479/1‑121 (MQ=255) ttCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGtt < 1:352479/121‑1 (MQ=255) aCGTCGCATCCGGCATGAACAAAGCGCACTTTGTGAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTaa < 3:686248/148‑1 (MQ=255) aCGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTccc < 4:585311/121‑1 (MQ=255) aCGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTccc > 3:585311/1‑121 (MQ=255) cGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGc < 1:1006031/77‑1 (MQ=255) cGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGc > 2:1006031/1‑77 (MQ=255) cGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGtt > 3:1013909/1‑117 (MQ=255) cGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGtt < 4:1013909/117‑1 (MQ=255) cGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTcc > 3:685458/1‑119 (MQ=255) cGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTcc < 4:685458/119‑1 (MQ=255) tCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCgaga < 1:692875/132‑1 (MQ=255) tCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCgaga > 2:692875/1‑132 (MQ=255) cGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGcc > 4:897135/1‑81 (MQ=255) cGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGcc < 3:897135/81‑1 (MQ=255) tCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTc < 2:825967/52‑1 (MQ=255) tCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTc > 1:825967/1‑52 (MQ=255) tCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCaa < 2:566982/88‑1 (MQ=255) tCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCaa > 1:566982/1‑88 (MQ=255) cATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACtttt > 2:25372/1‑41 (MQ=255) cATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACtttt < 1:25372/41‑1 (MQ=255) tgaCAAAGTGCGCTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACaaa < 2:415167/87‑1 (MQ=255) tgaCAAAGTGCGCTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACaaa > 1:415167/3‑89 (MQ=255) aCAAAGTGCGCTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATa < 4:549548/73‑1 (MQ=255) aCAAAGTGCGCTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATa > 3:549548/1‑73 (MQ=255) cAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATa > 2:651099/1‑148 (MQ=255) aaaGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATaa > 3:623242/1‑148 (MQ=255) cgcACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAAt > 3:339365/1‑148 (MQ=255) | CTACGCTTATCAGGCCTACGTGGTTTATGCAATATATTGAATTTGCATGGTCTTGTAGGCCAGATAAGACGTTCACGTCGCATCCGGCATGAACAAAGCGCACTTTGTCAAAAATCTAACCTACTTTTTAATTCAGGGAATTACCGCAAAGCCCACGCCATCATACCCAACGTAACAAAACTCAGGCACGTTCCCCACTCGCCCCGAGAAAAATAGCATTAATGCGCCCAGCGCCAGCATAAAAAT > NC_000913/3563592‑3563837 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |