Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 551,361 | 2 bp→AA | intergenic (+29/+165) | ybcF → / ← purK | putative carbonate kinase/N5‑carboxyaminoimidazole ribonucleotide synthase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 551,361 | 0 | G | A | 92.9% | 70.7 / ‑1.1 | 28 | intergenic (+29/+166) | ybcF/purK | putative carbonate kinase/N5‑carboxyaminoimidazole ribonucleotide synthase |
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (12/14); total (13/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01 | |||||||||||
* | NC_000913 | 551,362 | 0 | C | A | 92.9% | 80.2 / 1.8 | 28 | intergenic (+30/+165) | ybcF/purK | putative carbonate kinase/N5‑carboxyaminoimidazole ribonucleotide synthase |
Reads supporting (aligned to +/- strand): ref base C (1/1); new base A (12/14); total (13/15) | |||||||||||
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.95e-01 |
TGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGG > NC_000913/551233‑551481 || tGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGa < 2:198256/139‑1 (MQ=255) aGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATc > 1:410191/1‑137 (MQ=255) aGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATc < 2:410191/137‑1 (MQ=255) cgTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACg > 2:704214/1‑105 (MQ=255) cgTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATGCGACg < 1:704214/105‑1 (MQ=255) gATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCAt < 1:96065/139‑1 (MQ=255) ggCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATc < 2:833452/111‑1 (MQ=255) ggCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATc > 1:833452/1‑111 (MQ=255) cgAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGcac > 2:156645/1‑139 (MQ=255) aTTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAAt < 2:575865/139‑1 (MQ=255) aaGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCt < 1:156645/139‑1 (MQ=255) gaCGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatg < 1:856427/139‑1 (MQ=255) cGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTgatgcg > 1:992006/1‑139 (MQ=255) cgAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCtt < 1:513169/135‑1 (MQ=255) cgAAGCGGGGACCTGTATTTCGCTGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCtt > 2:513169/1‑135 (MQ=255) tGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAAt < 1:491186/96‑1 (MQ=38) tGTAATCGTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAAt > 2:491186/1‑96 (MQ=38) gTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAAt < 2:149926/78‑1 (MQ=255) gTAGGCATTAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAAt > 1:149926/1‑78 (MQ=37) tAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCt < 2:297483/85‑1 (MQ=255) tAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCt > 1:297483/1‑85 (MQ=255) tAGACATTTGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGCTGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAgg > 1:135203/1‑139 (MQ=255) ttGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGc > 1:810075/1‑107 (MQ=255) ttGTGCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGc < 2:810075/107‑1 (MQ=255) tgCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt > 2:481655/1‑99 (MQ=255) tgCCTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTAt < 1:481655/99‑1 (MQ=255) cTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGTCGCATCCGGCAATTGGTGcac < 2:669210/65‑1 (MQ=25) cTGATAAGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGcac > 1:669210/1‑65 (MQ=255) || TGGTTATGTCAGAAGCCGTGGTAAACCCGCGTGGATTGGGGCGTTATCGCGAATTGAAGAGACGCTGGCGGGCGAAGCGGGGACCTGTATTTCGCTGTAGTCGTAGGCATTAGACATTTGTGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAACCGGTGCCGCATCCGGCAATTGGTGCACAATGCCTGATGCGATGCTTGACGCATCTTATCAGGCCTACAATGGGTACCGGATCGGTAGG > NC_000913/551233‑551481 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |