Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,904,097 | +G | intergenic (+723/+650) | ygcE → / ← queE | putative sugar kinase YgcE/putative 7‑carboxy‑7‑deazaguanine synthase QueE |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,904,093 | 1 | . | G | 92.9% | 86.8 / 0.5 | 28 | intergenic (+719/+654) | ygcE/queE | putative sugar kinase YgcE/putative 7‑carboxy‑7‑deazaguanine synthase QueE |
Reads supporting (aligned to +/- strand): ref base . (0/2); new base G (13/13); total (13/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.84e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01 |
TTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGC‑GGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGG > NC_000913/2903961‑2904218 | ttaCAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGgcgc < 1:1504/133‑1 (MQ=255) gTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGC‑ggg‑gga > 2:102049/1‑131 (MQ=255) cTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGC‑ggg‑ggaaca > 1:120041/1‑128 (MQ=255) cTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCACGGTGGCGC‑ggg‑ggaaca > 1:172830/1‑128 (MQ=255) gCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACAcc > 2:347214/1‑133 (MQ=255) cTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCt > 1:311523/1‑133 (MQ=255) cTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCACCGATGACGC‑GGGCG‑AACACCt > 1:242835/1‑133 (MQ=255) ataAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGtttt > 2:307985/1‑133 (MQ=255) aaaaaCTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCg < 1:53115/133‑1 (MQ=255) aaaaaCTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCg < 2:120041/133‑1 (MQ=255) aaaaaCTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCg < 1:473289/133‑1 (MQ=255) aaaaCTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGC‑GGG‑G‑AACACCTGTTTTCGCa < 1:287308/133‑1 (MQ=255) gagaAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTa < 2:423828/133‑1 (MQ=255) agaAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTAt < 1:307985/133‑1 (MQ=255) gaAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATg < 2:337261/133‑1 (MQ=255) cGGGTTATTTGGTGGCTTGTCCTCGCTGGCCCGGGGAACTTGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTa < 1:155799/133‑1 (MQ=255) tATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCt > 1:376321/1‑133 (MQ=255) tATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCt > 2:332116/1‑133 (MQ=255) ccccGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGgc < 2:311523/133‑1 (MQ=255) cGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGcgcg > 1:274195/1‑133 (MQ=255) gCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCgg > 2:225668/1‑133 (MQ=255) gAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGgc < 2:376321/133‑1 (MQ=255) cTCGACAGAACGGCCTCAGTAGTCTCGTCAGTCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGGTTAACCCGGCTGGCGCAGGGAACACGGGCGCa > 2:238502/1‑133 (MQ=255) aaCGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATAc > 2:335793/1‑133 (MQ=255) gCCTCAGTAGACTCGTCAGGCTCCTGTTTATCGCCGCTGGCACGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAg < 2:172830/133‑1 (MQ=255) agtCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTAt < 1:396692/133‑1 (MQ=255) tCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGCGGAACACGCGCGAAGGGAATACAAAACCGTCCATCTGCTCg > 2:233783/1‑133 (MQ=255) tCAGGCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGAGGAACACGAGCGAACGGAATACAAAACCGTGTATCTGCTCg > 2:208198/1‑133 (MQ=255) gCTCCGGTTTATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCAAGGAATACAAAGACGCGTATCTACTCGGttt > 1:246622/1‑133 (MQ=255) tttATCCCCGCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGG‑TTATCCCCGc < 1:102049/133‑1 (MQ=255) cccGCTGGGGC‑GGG‑GTAACACCGGTTTTCGTAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGgcg < 1:208198/133‑1 (MQ=255) gCTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCgg < 1:332116/133‑1 (MQ=255) cTGGCGCGGGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCggg > 2:330967/1‑133 (MQ=255) | TTACAAGTACTGCTAATATAAAAACTTGAGAAAGAGATAACGGGTTATATGGTGGTTTATCCCCGCTGGCGCGGGGAACTCGACAGAACGGCCTCAGTAGTCTCGTCAGGCTCCGGTTTATCCCCGCTGGCGC‑GGG‑G‑AACACCTGTTTTCGCAAATCTATGGACTATTGCTATTCGGTTTATCCCCGCTGGCGCGGGGAACACGGGCGCACGGAATACAAAGCCGTGTATCTGCTCGGTTTATCCCCGCTGGCGCGGG > NC_000913/2903961‑2904218 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |