Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,362,048 | 0 | A | C | 33.3% | 308.8 / 20.5 | 121 | intergenic (‑114/+503) | cadC/pheU | DNA‑binding transcriptional activator CadC/tRNA‑Phe |
Reads supporting (aligned to +/- strand): ref base A (54/26); new base C (0/40); total (54/67) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.48e-14 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
Only 100 of 131 total aligned reads displayed.
TACATATTATTGGTGAATGCAAGACGTTATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACG > NC_000913/4362014‑4362082 | taCATATTATTGGTGAATGCAAGACGTTATTTTTa > 1:5502916/1‑35 (MQ=255) taCATATTATTGGTGAATGCAAGACGTTATTTTTa > 1:6295099/1‑35 (MQ=255) taCATATTATTGGTGAATGCAAGACGTTATTTTTa < 1:7665521/35‑1 (MQ=255) taCATATTATTGGTGAATGCAAGACGTTATTTTTa > 1:7353923/1‑35 (MQ=255) taCATATTATTGGTGAATGCAAGACGTTATTTTTa > 1:4095271/1‑35 (MQ=255) aCATATTATTGGTGAATGCAAGACGTTATTTTTAc > 1:4802552/1‑35 (MQ=255) aCATATTATTGGTGAATGCAAGACGTTATTTTTAc > 1:6968801/1‑35 (MQ=255) aCATATTATTGGTGAATGCAAGACGTTATTTTTAc > 1:957677/1‑35 (MQ=255) cATATTATTGGTGAATGCAAGACGTTATTTTTAcc > 1:1905409/1‑35 (MQ=255) cATATTATTGGTGAATGCAAGACGTTATTTTTAcc > 1:6707647/1‑35 (MQ=255) cATATTATTGGTGAATGCAAGACGTTATTTTTAcc > 1:3947315/1‑35 (MQ=255) cATATTATTGGTGAATGCAAGACGTTATTTTTAcc < 1:6082510/35‑1 (MQ=255) cATATTATTGGTGAATGCAAGACGTTATTTTTAcc < 1:5896742/35‑1 (MQ=255) tattatTGGTGAATGCAAGACGTTATTTTTACCAg > 1:5380610/1‑35 (MQ=255) tattatTGGTGAATGCAAGACGTTATTTTTACCAg < 1:322000/35‑1 (MQ=255) attattGGTGAATGCAAGACGTTATTTTTACCAGc < 1:4542115/35‑1 (MQ=255) attattGGTGAATGCAAGACGTTATTTTTACCAGc < 1:1018107/35‑1 (MQ=255) attattGGTGAATGCAAGACGTTATTTTTACCAGc < 1:5774221/35‑1 (MQ=255) attattGGTGAATGCAAGACGTTATTTTTACCAGc < 1:1125197/35‑1 (MQ=255) ttattGGTGAATGCAAGACGTTATTTTTACCAGcc > 1:5596423/1‑35 (MQ=255) ttattGGTGAATGCAAGACGTTATTTTTACCAGcc > 1:6830326/1‑35 (MQ=255) ttattGGTGAATGCAAGACGTTATTTTTACCAGcc < 1:1743196/35‑1 (MQ=255) ttattGGTGAATGCAAGACGTTATTTTTACCAGcc < 1:3739481/35‑1 (MQ=255) ttattGGTGAATGCAAGACGTTATTTTTACCAGcc > 1:3729590/1‑35 (MQ=255) tattGGTGAATGCAAGACGTTATTTTTACCAGCCa < 1:542065/35‑1 (MQ=255) tattGGTGAATGCAAGACGTTATTTTTACCAGCCa < 1:701883/35‑1 (MQ=255) attGGTGAATGCAAGACGTTATTTTTACCAGCCAt < 1:5357918/35‑1 (MQ=255) attGGTGAATGCAAGACGTTATTTTTACCAGCCAt > 1:7330951/1‑35 (MQ=255) attGGTGAATGCAAGACGTTATTTTTACCAGCCAt > 1:2963440/1‑35 (MQ=255) attGGTGAATGCAAGACGTTATTTTTACCAGCCAt > 1:297130/1‑35 (MQ=255) ttGGTGAATGCAAGACGTTATTTTTCCCAGCCATa < 1:3047210/35‑1 (MQ=255) ttGGTGAATGCAAGACGTTATTTTTACCAGCCATa > 1:4106552/1‑35 (MQ=255) ttGGTGAATGCAAGACGTTATTTTTACCAGCCATa > 1:1115166/1‑35 (MQ=255) ttGGTGAATGCAAGACGTTATTTTTACCAGCCATa > 1:1570507/1‑35 (MQ=255) ggTGAATGCAAGACGTTATTTTTACCAGCCATaaa > 1:3964326/1‑35 (MQ=255) ggTGAATGCAAGACGTTATTTTTACCAGCCATaaa > 1:2440335/1‑35 (MQ=255) gTGAATGCAAGACGTTATTTTTACCAGCCATAAAc < 1:5646173/35‑1 (MQ=255) gTGAATGCAAGACGTTATTTTTACCAGCCATAAAc > 1:2185665/1‑35 (MQ=255) gTGAATGCAAGACGTTATTTTTACCAGCCATAAAc > 1:176331/1‑35 (MQ=255) gAATGCAAGACGTTATTTTTCCCAGCCATAAActt < 1:7541474/35‑1 (MQ=255) aaTGCAAGACGTTATTTTTACCAGCCATAAActtc > 1:2394414/1‑35 (MQ=255) aTGCAAGACGTTATTTTTACCAGCCATAAActtct > 1:329358/1‑35 (MQ=255) aTGCAAGACGTTATTTTTACCAGCCATAAActtct > 1:1586985/1‑35 (MQ=255) tGCAAGACGTTATTTTTACCAGCCATAAActtctt > 1:4921795/1‑35 (MQ=255) gCAAGACGTTATTTTTACCAGCCATAAACTTCTTg > 1:3164104/1‑35 (MQ=255) cAAGACGTTATTTTTCCCAGCCATAAACTTCTTGa < 1:5138476/35‑1 (MQ=255) cAAGACGTTATTTTTACCAGCCATAAACTTCTTGa < 1:3080793/35‑1 (MQ=255) cAAGACGTTATTTTTACCAGCCATAAACTTCTTGa > 1:2530506/1‑35 (MQ=255) cAAGACGTTATTTTTACCAGCCATAAACTTCTTGa > 1:2876002/1‑35 (MQ=255) aaGACGTTTTTTTTCCCAGCCATAAACTTCTTGAt < 1:5803096/35‑1 (MQ=25) aaGACGTTATTTTTCCCAGCCATAAACTTCTTGAt < 1:237392/35‑1 (MQ=255) aaGACGTTATTTTTCCCAGCCATAAACTTCTTGAt < 1:237393/35‑1 (MQ=255) aaGACGTTATTTTTCCCAGCCATAAACTTCTTGAt < 1:433949/35‑1 (MQ=255) aaGACGTTATTTTTACCAGCCATAAACTTCTTGAt > 1:5455664/1‑35 (MQ=255) aaGACGTTATTTTTACCAGCCATAAACTTCTTGAt < 1:5964855/35‑1 (MQ=255) aaGACGTTATTTTTACCAGCCATAAACTTCTTGAt < 1:2161883/35‑1 (MQ=255) aGACGTTATTTTTCCCAGCCATAAACTTCTTGAtt < 1:5143388/35‑1 (MQ=255) aGACGTTATTTTTCCCAGCCATAAACTTCTTGAtt < 1:5292161/35‑1 (MQ=255) aGACGTTATTTTTCCCAGCCATAAACTTCTTGAtt < 1:7131611/35‑1 (MQ=255) aGACGTTATTTTTCCCAGCCATAAACTTCTTGAtt < 1:2597524/35‑1 (MQ=255) aGACGTTATTTTTACCAGCCATAAACTTCTTGAtt > 1:7272180/1‑35 (MQ=255) aGACGTTATTTTTACCAGCCATAAACTTCTTGAtt > 1:2135380/1‑35 (MQ=255) gACGTTATTTTTCCCAGCCATAAACTTCTTGATTa < 1:5019748/35‑1 (MQ=255) gACGTTATTTTTCCCAGCCATAAACTTCTTGATTa < 1:6148048/35‑1 (MQ=255) gACGTTATTTTTCCCAGCCATAAACTTCTTGATTa < 1:3463967/35‑1 (MQ=255) aCGTTATTTTTACCAGCCATAAACTTCTTGATTAc > 1:5009554/1‑35 (MQ=255) aCGTTATTTTTACCAGCCATAAACTTCTTGATTAc > 1:6193621/1‑35 (MQ=255) cGTTATTTTTCCCAGCCATAAACTTCTTGATTACa < 1:5148993/35‑1 (MQ=255) cGTTATTTTTCCCAGCCATAAACTTCTTGATTACa < 1:4888319/35‑1 (MQ=255) ttATTTTTCCCAGCCATAAACTTCTTGATTACATa < 1:5540490/35‑1 (MQ=255) ttATTTTTCCCAGCCATAAACTTCTTGATTACATa < 1:4677926/35‑1 (MQ=255) ttATTTTTACCAGCCATAAACTTCTTGATTACATa > 1:1856596/1‑35 (MQ=255) tATTTTTCCCAGCCATAAACTTCTTGATTACATAg < 1:840049/35‑1 (MQ=255) tATTTTTCCCAGCCATAAACTTCTTGATTACATAg < 1:5986441/35‑1 (MQ=255) tATTTTTACCAGCCATAAACTTCTTGATTACATAg > 1:7537782/1‑35 (MQ=255) tATTTTTACCAGCCATAAACTTCTTGATTACATAg < 1:4208105/35‑1 (MQ=255) tATTTTTACCAGCCATAAACTTCTTGATTACATAg < 1:821441/35‑1 (MQ=255) tATTTTTACCAGCCATAAACTTCTTGATTACATAg < 1:6961025/35‑1 (MQ=255) tATTTTCACCAGCCATAAACTTCTTGATTACATAg < 1:4208146/35‑1 (MQ=255) aTTTTTACCAGCCATAAACTTCTTGATTACATAGt > 1:6356541/1‑35 (MQ=255) aTTTTTACCAGCCATAAACTTCTTGATTACATAGt > 1:6558201/1‑35 (MQ=255) aTTTTTACCAGCCATAAACTTCTTGATTACATAGt > 1:1948653/1‑35 (MQ=255) tttttCCCAGCCATAAACTTCTTGATTACATAGTa < 1:4125149/35‑1 (MQ=255) tttttCCCAGCCATAAACTTCTTGATTACATAGTa < 1:4663140/35‑1 (MQ=255) tttttACCAGCCATAAACTTCTTGATTACATAGTa < 1:7638420/35‑1 (MQ=255) tttttACCAGCCATAAACTTCTTGATTACATAGTa > 1:4655517/1‑35 (MQ=255) tttttACCAGCCATAAACTTCTTGATTACATAGTa > 1:5033653/1‑35 (MQ=255) tttttACCAGCCATAAACTTCTTCATTACATAGTa > 1:2323067/1‑35 (MQ=255) ttttGCCAGCCATAAACTTCTTGATTACATAGTAt < 1:7754112/35‑1 (MQ=255) ttttCCCAGCCATAAACTTCTTGATTACATAGTAt < 1:2016087/35‑1 (MQ=255) ttttCCCAGCCATAAACTTCTTGATTACATAGTAt < 1:2520406/35‑1 (MQ=255) ttttCCCAGCCATAAACTTATTGATTACATAGTAt < 1:3932217/35‑1 (MQ=25) tttCCCAGCCATAAACTTCTTGATTACATAGTAtt < 1:5763360/35‑1 (MQ=255) tttCCCAGCCATAAACTTCTTGATTACATAGTAtt < 1:1612478/35‑1 (MQ=255) tttCCCAGCCATAAACTTCTTGATTACATAGTAtt < 1:5004099/35‑1 (MQ=255) tttACCAGCCATAAACTTCTTGATTACATAGTAtt > 1:7578572/1‑35 (MQ=255) ttACCAGCCATAAACTTCTTGATTACATAGTATTa > 1:2897969/1‑35 (MQ=255) ttACCAGCCATAAACTTCTTGATTACATAGTATTa > 1:6347787/1‑35 (MQ=255) ttACCAGCCATAAACTTCTTGATTACATAGTATTa > 1:5535544/1‑35 (MQ=255) ttACCAGCCATAAACTTCTTGATTACATAGTATTa < 1:1483648/35‑1 (MQ=255) | TACATATTATTGGTGAATGCAAGACGTTATTTTTACCAGCCATAAACTTCTTGATTACATAGTATTACG > NC_000913/4362014‑4362082 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |