Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,145,707 | C→T | intergenic (+38/‑157) | dnaG → / → rpoD | DNA primase/RNA polymerase sigma factor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,145,707 | 0 | C | T | 97.4% | 286.7 / ‑0.3 | 77 | intergenic (+38/‑157) | dnaG/rpoD | DNA primase/RNA polymerase sigma factor |
Reads supporting (aligned to +/- strand): ref base C (1/1); new base T (53/22); total (54/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.11e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01 |
TCTGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACA > REL606/3145635‑3145777 | tctGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGcctcc > 2:1980846/1‑75 (MQ=255) tctGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGcctcc < 2:660934/75‑1 (MQ=255) ctGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCg > 1:1365543/1‑75 (MQ=255) ctGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCg > 1:416478/1‑75 (MQ=255) ggACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGAc < 2:417143/75‑1 (MQ=255) ggACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGAc < 1:724735/75‑1 (MQ=255) gACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACa > 1:2040190/1‑75 (MQ=255) aCATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAg > 2:1777145/1‑75 (MQ=255) aCATTAAACCAGGAGCTGGAGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAg < 2:1017862/75‑1 (MQ=255) ttAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCg > 2:1665025/1‑75 (MQ=255) tAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGc > 1:1052783/1‑75 (MQ=255) tAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGc > 2:1263043/1‑75 (MQ=255) tAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGc > 1:259068/1‑75 (MQ=255) aaaCCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCa > 1:369058/1‑75 (MQ=255) aaCCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCAc > 2:1980811/1‑75 (MQ=255) aaCCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCAc > 2:530792/1‑75 (MQ=255) ccAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTg > 1:1659182/1‑75 (MQ=255) cAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTga > 2:535909/1‑75 (MQ=255) cAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTga < 2:1219698/75‑1 (MQ=255) gCTGGCGAAAAAGTGATTTAACGGCTTAAGGGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGc < 2:483172/75‑1 (MQ=255) cTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCa < 2:589168/75‑1 (MQ=255) cTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCa < 1:985845/75‑1 (MQ=255) tGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAg > 1:1284277/1‑75 (MQ=255) tGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAg < 1:701363/75‑1 (MQ=255) ggCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGc > 2:1955407/1‑75 (MQ=255) gCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCg < 1:268974/75‑1 (MQ=255) gAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGc > 1:1503616/1‑75 (MQ=255) aaaaGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCaa > 1:1893611/1‑75 (MQ=255) aaaGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCaaa > 1:1971084/1‑75 (MQ=255) aaaGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCaaa > 2:1330800/1‑75 (MQ=255) aaGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAat > 2:744642/1‑75 (MQ=255) aGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAata > 2:633451/1‑75 (MQ=255) aGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAata > 2:20502/1‑75 (MQ=255) tGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAatata < 1:1558239/75‑1 (MQ=255) tGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAatata < 2:1271820/75‑1 (MQ=255) gATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAatatat > 2:2113792/1‑75 (MQ=255) gATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAatatat > 1:31312/1‑75 (MQ=255) aTTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGGCCCCGACACCCGCACTGAGAGGCAGCGGACAAtatata > 1:495372/1‑75 (MQ=255) tttAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATaa < 2:323424/75‑1 (MQ=255) tttAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATaa > 1:225801/1‑75 (MQ=255) tttAACGGCTTAAGGGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATaa < 2:90033/75‑1 (MQ=255) ttAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGGAAATCTACAAg > 2:37842/1‑75 (MQ=255) tagcGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGt < 2:239577/72‑1 (MQ=255) tAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGt > 2:1508505/1‑75 (MQ=255) aaCGGCTTAAGTGCCGAAGAGGGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTa < 2:1558158/75‑1 (MQ=255) aaCGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTa > 2:1260173/1‑75 (MQ=255) aaCGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTa > 2:2059840/1‑75 (MQ=255) aaCGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTa > 1:832913/1‑75 (MQ=255) ggCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGc > 2:683661/1‑75 (MQ=255) ggCTTAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGc > 2:683641/1‑75 (MQ=255) ttAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCt > 1:1736365/1‑75 (MQ=255) tAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTc < 2:97232/75‑1 (MQ=255) tAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTc > 1:1258916/1‑75 (MQ=255) tAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTc > 1:1720873/1‑75 (MQ=255) tAAGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTc > 1:647740/1‑75 (MQ=255) aaGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCg > 2:1490819/1‑75 (MQ=255) aGTGCCGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGt > 1:329064/1‑75 (MQ=255) ccGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAAtt > 2:431390/1‑75 (MQ=255) ccGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAAtt < 2:259068/75‑1 (MQ=255) ccGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAAtt < 1:892719/75‑1 (MQ=255) cGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTa > 1:542564/1‑75 (MQ=255) cGAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTa < 1:2081346/75‑1 (MQ=255) gAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTAt > 2:2489/1‑75 (MQ=255) gAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTAt > 2:1348703/1‑75 (MQ=255) gAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTAt > 2:1832696/1‑75 (MQ=255) gAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTAt > 2:641708/1‑75 (MQ=255) gAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTAt > 2:2069582/1‑75 (MQ=255) gAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTAt < 1:1298350/75‑1 (MQ=255) gAAGAGCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAACTACGCCCTCGTAATTAt > 2:645443/1‑75 (MQ=255) gagCGATCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGt > 1:1396402/1‑75 (MQ=255) aTCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCg > 2:1366780/1‑75 (MQ=255) aTCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCg > 2:2096076/1‑75 (MQ=255) aTCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCg > 1:1729205/1‑75 (MQ=255) aTCGGGAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCg > 1:1729189/1‑75 (MQ=255) gggAGGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTa < 2:912349/75‑1 (MQ=255) gggAAGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTa > 1:811381/1‑75 (MQ=255) aaGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaac > 1:1336265/1‑75 (MQ=255) aaGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaac < 2:170556/75‑1 (MQ=255) aaGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaac > 2:1264033/1‑75 (MQ=255) aGCCTCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAaaca < 2:85531/75‑1 (MQ=255) | TCTGGACATTAAACCAGGAGCTGGCGAAAAAGTGATTTAACGGCTTAAGTGCCGAAGAGCGATCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACA > REL606/3145635‑3145777 |
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: ‑ |