Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_012967 | 4,169,221 | 1 | . | G | 22.2% | 112.8 / 28.3 | 36 | intergenic (+37/‑176) | rpoC/yjaZ | DNA‑directed RNA polymerase subunit beta'/heat‑shock protein C |
Reads supporting (aligned to +/- strand): ref base . (5/23); new base G (8/0); total (13/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.25e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGAATTTCTGAATACTCA > NC_012967/4169150‑4169291 | tgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTc‑g‑gc < 1:759069/72‑1 (MQ=255) tgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTc‑g‑gc < 1:8653081/72‑1 (MQ=255) tgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTc‑g‑gc < 2:1043845/72‑1 (MQ=255) tgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTc‑g‑gc < 2:7833865/72‑1 (MQ=255) tgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTc‑g‑gc < 2:7217862/72‑1 (MQ=255) tgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTc‑g‑gc < 2:7159700/72‑1 (MQ=255) tgAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTc‑g‑gc > 2:6767454/1‑72 (MQ=255) aaCGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTc‑g‑gc‑‑gg < 2:8125694/72‑1 (MQ=255) aCGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GC‑‑ggg < 2:6240494/72‑1 (MQ=255) cgcgGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GC‑‑GGG‑t < 1:4976981/72‑1 (MQ=255) cGCAGGTCTGGGCGGTTCTTATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GC‑‑GGG‑t < 2:4266892/72‑1 (MQ=255) cGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAACCCCC‑CTCC‑G‑GC‑‑ggg > 2:1797672/1‑72 (MQ=255) cGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GC‑‑GGG‑t < 1:1518471/72‑1 (MQ=255) gaaGGTCTGGGAGGTTCTGATANCGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GC‑‑GGG‑tt < 1:2549504/70‑1 (MQ=255) gCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAA‑CCCC‑CTCC‑G‑GC‑‑gggg > 1:4951444/1‑71 (MQ=255) gCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GC‑‑GGG‑tt < 1:6066304/72‑1 (MQ=255) gCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCC‑CTTC‑G‑GC‑‑ggggg > 1:1981323/1‑70 (MQ=255) caaGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GC‑‑GGG‑tt < 2:1401818/70‑1 (MQ=255) cgggTCTGGGGGGTTCTGATAACGAGTAATCTTTAATCCGAAAATAACGA‑AAAAC‑CCCG‑CTTC‑G‑GC‑‑GGG‑ttt < 1:6298280/70‑1 (MQ=37) gggTCTGGGCGTTTCTAAAAACGAGTAATCTTTAACCCAAAAATAACGA‑AAAAC‑CCCG‑CTTC‑G‑GC‑‑GGG‑tttt < 2:8104007/71‑1 (MQ=19) gTCGGGGCGCTTCTGATAACGAGAACTCTTAAACCCGAAAAAAACGA‑AAAAC‑CCCG‑CTTG‑G‑GC‑‑GGG‑tttttt < 2:4880647/72‑1 (MQ=12) ttCTGATAACGAGTAATCGTTAATCCGCAAATAACGTAAAAAA‑CCCC‑CTCC‑G‑GC‑‑GGG‑GGTTTTTTTtggggg > 2:2547828/1‑72 (MQ=255) ttCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GC‑‑GGG‑GTTTTTTTTggggggg > 2:8699928/1‑71 (MQ=255) ttCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑ACCC‑CCTC‑GTGC‑‑GGG‑GGTTTTTTTtggggg > 1:19371/1‑72 (MQ=255) tCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAg < 2:2318725/72‑1 (MQ=255) tCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCC‑CCTC‑G‑GCGGGGG‑TTTTTTTTTgggggg > 1:6512914/1‑72 (MQ=255) ctaaTAACGAATAGTCGTAAACCCACAAATAACGA‑AAAAC‑CCCG‑CTGG‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGt < 1:6403765/69‑1 (MQ=16) tGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCC‑CCTC‑G‑GCG‑GGG‑GTTTTTTTTggggggggg > 1:4813399/1‑69 (MQ=255) aaaaCAAGTAATCTTTAACCCGCAAAAAACTA‑AAAAC‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTa < 1:4109123/70‑1 (MQ=21) aTAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCC‑CCTC‑G‑GCGGGGG‑TTTTTTTTTGGGGGGGGtt > 1:8683178/1‑72 (MQ=255) aCGAGTAATCGTAAATCCGAAAATAACGT‑AAAAC‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAggg < 2:7334716/72‑1 (MQ=38) cGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GCG‑GGG‑TTTTTTTTTGGGGGGGGTTTTggg > 1:7495163/1‑72 (MQ=25) cGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCC‑CTTC‑G‑GCG‑GGG‑TTTTTTTTTGGGGGGGGTTTTggg > 1:2875456/1‑72 (MQ=25) aGTAATCGTTAATCCGCAAATAACGTAAAAAA‑CCCC‑CTTC‑GCGG‑‑GGG‑TTTTTTTTTGGGGGGGGTTTggggg > 2:312317/1‑71 (MQ=255) aGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GCG‑GGG‑TTTTTTTTTGGGGGGGGTTTGGGGaa > 1:7157952/1‑72 (MQ=25) aaaTCGTAGATGCGAAAAAAACGT‑AAAAC‑CCCG‑ATCC‑G‑GG‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAaga < 2:4581630/71‑1 (MQ=255) attCGTTAACCCGAAAATAACGA‑AAAAC‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAagag < 1:5118942/70‑1 (MQ=25) aaTCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GCG‑GGG‑TTTTTTTTTGGGGGGGGGTTTGGGAAaga > 1:4256162/1‑72 (MQ=25) aTCGTTAATCCGCAAATAACGT‑AAAAACCCCC‑CTTCGG‑GG‑‑GGG‑TTTTTTTTTGGGGGGGTTTTGGGgagaga > 1:107598/1‑72 (MQ=255) ttAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCAttt < 2:1777077/72‑1 (MQ=255) ccGAAAATAACGT‑AAAAC‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAg < 2:6353889/72‑1 (MQ=255) gCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑aa < 1:9152323/72‑1 (MQ=255) aaaaaCAT‑AAACA‑CCCGTCTTC‑G‑GC‑‑GGGTTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑Aat < 2:5730866/69‑1 (MQ=255) aaaaCGT‑AGAAG‑CCCG‑CTTC‑G‑GC‑‑GAG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATAtt < 2:4177847/70‑1 (MQ=38) aaCGT‑AAAAC‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTa < 2:4540793/72‑1 (MQ=255) cGT‑AAAAA‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAg < 2:6983236/72‑1 (MQ=255) gT‑AAAAA‑CCCG‑CTTC‑G‑GCG‑GGG‑TTTTTTTTTGGGGGGGGTTTTGGGAAAAGGAATTTTTC‑G‑AAAATTTTAgg > 1:9032707/1‑72 (MQ=9) aaaaaac‑CCCG‑CTGC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAgg < 1:1900695/70‑1 (MQ=38) aaaaa‑CCCG‑CTTC‑G‑GCG‑GGG‑TTTTTTTTTGGGGGGGGTTTGGGGAAA‑AGCAATTTTTAG‑AATATTTTAGGaa > 2:7164643/1‑72 (MQ=16) aaaaa‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGaa < 2:7619848/72‑1 (MQ=255) aaaa‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGGGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGAAt > 2:6983378/1‑72 (MQ=255) aa‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGAAttt < 1:2386040/72‑1 (MQ=255) a‑CCCG‑CTTC‑G‑GCG‑GGG‑TTTTTTTTGGGGGGGGGTTTTGGGAAA‑AGGATTTTTCAGAAATTTTTAAGGAtttc > 2:9226862/1‑71 (MQ=9) a‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGAATTTc < 1:7004953/72‑1 (MQ=255) ccG‑CTTC‑G‑GCG‑GGG‑TTTTTTTTTGGGGGGGGTTTGGGGAAAAGGAATTTTTTAG‑AAAATTTAATGAAttttt > 2:5191930/1‑70 (MQ=12) g‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGAATTTCTGaa < 2:7477129/72‑1 (MQ=255) cTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGAATTTCTGAAt < 2:9134823/72‑1 (MQ=255) ttC‑G‑GC‑‑GGG‑GTTTTTTTTTGGGGGGGGTTAGGGAAAGAGGCTTTGTCAG‑AAGATTTAAGGAATTTCTGaaaa > 2:6505878/1‑70 (MQ=255) tC‑G‑GG‑‑GGGTTTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGAATTTCTGAATa < 1:1797673/72‑1 (MQ=255) tC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGAATTTCTGAATAc < 1:7433266/72‑1 (MQ=255) tC‑G‑GA‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGAATTTCTGAATAc > 2:8502615/1‑72 (MQ=255) c‑g‑gc‑‑ggG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGAATTTCTGAATACt < 2:7034479/72‑1 (MQ=255) gg‑gc‑‑ggG‑TTTTTCTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGAATTTCTGAATACt < 2:3972363/71‑1 (MQ=255) g‑gc‑‑ggG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGAATTTCTGAATACTc < 1:6641094/72‑1 (MQ=255) g‑gc‑‑ggG‑GTTTTTTTTGGGGGGGGGTTTGGGGAAGAGGATTTTTCAG‑AATATTTTAGGAAATTCTGAATACTc > 2:6230967/1‑72 (MQ=255) gc‑‑ggG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGAATTTCTGAATACTCa > 2:2548160/1‑72 (MQ=255) | TGAACGCAGGTCTGGGCGGTTCTGATAACGAGTAATCGTTAATCCGCAAATAACGT‑AAAAA‑CCCG‑CTTC‑G‑GC‑‑GGG‑TTTTTTTATGGGGGGAGTTTAGGGAAAGAGCATTTGTCAG‑AATATTTAAGGAATTTCTGAATACTCA > NC_012967/4169150‑4169291 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |