Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 1,262,127 | C→T | intergenic (‑102/+49) | prsA ← / ← ispE | ribose‑phosphate pyrophosphokinase/4‑diphosphocytidyl‑2‑C‑methyl‑D‑erythritol kinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,262,127 | 0 | C | T | 97.2% | 135.2 / ‑3.1 | 36 | intergenic (‑102/+49) | prsA/ispE | ribose‑phosphate pyrophosphokinase/4‑diphosphocytidyl‑2‑C‑methyl‑D‑erythritol kinase |
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (16/19); total (17/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.72e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01 |
GTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGG > REL606/1262054‑1262199 | gtggAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTaa < 2:1253485/76‑1 (MQ=255) ggAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCACGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACg > 2:1649759/1‑76 (MQ=255) aaTCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGa > 2:281614/1‑76 (MQ=255) gCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAg > 2:90262/1‑76 (MQ=255) cGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGt < 2:526786/76‑1 (MQ=255) gAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTg < 1:1020241/76‑1 (MQ=255) gAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTg > 1:123394/1‑76 (MQ=255) gAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTg > 2:679077/1‑76 (MQ=255) aaaaCTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGa > 2:1191497/1‑76 (MQ=255) ccAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGaa < 1:1249404/76‑1 (MQ=255) cAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGaaa < 1:1550928/76‑1 (MQ=255) aGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAAc > 1:1270962/1‑76 (MQ=255) gCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTc < 1:1449320/76‑1 (MQ=255) gCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTc < 2:423418/76‑1 (MQ=255) aTGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAg < 2:1467815/76‑1 (MQ=255) aTCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCtt < 2:118486/76‑1 (MQ=255) aTCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCtt > 2:788385/1‑76 (MQ=255) aTCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCtt > 1:1092849/1‑76 (MQ=255) tCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTg < 1:1838337/76‑1 (MQ=255) tCCTTTCAAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTg > 1:1169075/1‑76 (MQ=255) ccTTGCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGc < 1:30123/76‑1 (MQ=255) tttCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGccc > 2:1005718/1‑76 (MQ=255) tCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCgg > 1:1045249/1‑76 (MQ=255) tCCAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCgg > 1:429192/1‑76 (MQ=255) ccAGGCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGc < 2:541218/76‑1 (MQ=255) ggCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTa < 1:1763905/76‑1 (MQ=255) gCGGTGTATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTaa < 2:1443719/76‑1 (MQ=255) gtgtATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGc < 1:861667/76‑1 (MQ=255) tgtATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCa < 2:1731388/76‑1 (MQ=255) gtATTAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCAt > 1:1812773/1‑76 (MQ=255) ttAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGc > 2:341824/1‑76 (MQ=255) ttAAAGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGc > 2:663800/1‑76 (MQ=255) aaaGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGctc > 2:1461501/1‑76 (MQ=255) aaaGAGCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGctc < 1:88002/76‑1 (MQ=255) agCGCGATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCCTGGCTCtgtg > 1:948410/1‑76 (MQ=255) cgcgATGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCa < 2:923407/76‑1 (MQ=255) gATGTAACGTCTGTAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATg < 2:1660107/76‑1 (MQ=255) aTGTAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATgg < 2:1795630/76‑1 (MQ=255) | GTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAACGTCTGGAACAAGGTGACGTTGTCACCGAAACTCAGCTTGCCCGGCTTAAAGCATGGCTCTGTGCAATGG > REL606/1262054‑1262199 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |