Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 71,085 | 0 | T | C | 35.0% | 10.6 / 11.3 | 20 | S233S (AGT→AGC) | araC | DNA‑binding transcriptional dual regulator AraC |
Reads supporting (aligned to +/- strand): ref base T (8/5); new base C (1/6); total (9/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.03e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGT‑GAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCC > NC_000913/70841‑71323 | cGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGATCCATCCACCGATGGATAATCGGGAACGCGAGGCTTGACAGTACATCAGCGATCAACTGGCAGACAGCAATTTTGATATCGCCACAGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCg > 1:88209/1‑251 (MQ=255) tAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGACGCGATTAACGAGTCGCTCCATCCACCGCTGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCCATTTTGGTCTCGCCCGCTTCTCACCGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTCGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATCAGCCAGGCGaa > 1:102158/1‑250 (MQ=255) tGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGACGCGAGGACCAACGCATCAGCCAGGCGAAGCTGCTTTTGAGCAccac < 1:5100/250‑4 (MQ=255) cgcATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCATTACATCAGCGATCACCTTGCAGACAGCAATGTTGATATCGCCATCGGCGCCCATGATGTTTGCTTTTCGCCGTCGCGGCTGTCACAGCTTGTCCGTCAGCAGTTAGGGATTAGCGTCTTAAGCCGGCGCGAGGACCAACGCATCAGCCAGGCGAAGCTGCtttt < 1:89248/230‑1 (MQ=255) cATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATCAGCCAGGCGAAGCTGCTTTTGAGCACCGCCCGGATGCCTATCg < 1:49976/251‑1 (MQ=255) ngAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTTTCCGTACATCAGCGATCACCTGGCCGACAGCAATTTTGCTATCGACAGCGTCGCACCGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGCTTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGcc > 1:20697/2‑250 (MQ=255) aCCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAgn < 1:33627/251‑2 (MQ=255) gCAATTTTGATATCGCCAGCGTCGCACAGTATGTTTGCTTGTCGCCTTCGCGNCTGNCACANCTNTTCCGCCAGCATTTANGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGNNGNNNCNTCTTTTGa < 2:120527/141‑1 (MQ=255) gCAATTTTGATATCGCCAGCGTCGCACAGTATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGACAGCAGTTAGGGATTAGAGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGa > 1:120527/1‑141 (MQ=255) aGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGT‑ga > 1:76883/1‑248 (MQ=255) gCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGAGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTCTTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGT‑GAAGAAAAAGt > 1:31671/1‑251 (MQ=255) tCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCGGGCGCGAGGACCAACGCATTAGTCGTGCGAATCTGCTTTTGAGCTGTACCCGGATGTCTGGCGCCACCGTCTGTCGCAAGGTTTTTTTTAACGATCGACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCcgag < 1:116007/192‑1 (MQ=255) cTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGT‑GAAGAAAAAGTGAATGATGTCGCCGTCAAGTTGTCATAAtt > 1:156698/1‑251 (MQ=255) cacaTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAgc > 1:41841/1‑75 (MQ=255) caTCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGTGCCAGCCCGAGCGAGTTTCGTGCCGGTTGT‑GAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAAc < 1:159197/251‑1 (MQ=255) tttCCGCCATCAGTTAGGGATTAGCGTCTTAAGCTGGCGCTAGGACCAACGCATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTGGCAATGTTGGTTTTTACGATTAACGCTATTTCTCTCGGGGATTGAAAAAATGCACCGGGGCCAGGCCTAGCGAGTATCG‑GTCCTTTGTTGAAGAAAAAGTGAATGATGGAGCCGTCAAGTTGTCATAATTGGTAACg < 1:159168/247‑1 (MQ=255) gATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCCCCGTCGGTCGCAATGTTGTTTTTGCCGATCAACTCTATTTCTCGCGAGTATTTAAAAAAtn < 1:51721/144‑2 (MQ=255) cgcgAGGACCAACGCATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGGTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGTCCAGCCCGAGCGAGTTCCGTGCCGGTTGT‑GAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTGGACGGAGTAGCATAGGGTTTg < 1:67335/251‑1 (MQ=255) aCCAACGCATTAGTCAGGAGAAACTGCTTTTGAGCACTCCCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCAAGCACGAGCGAGTTTCGTGCCGGTTGT‑GAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACACTTGCCGGCTTGACGGAGTAGCATAGGGTTTGCAGCATc > 1:59582/1‑251 (MQ=255) ccAACGCATCAGCCAGGCGAAGCTGCTTTTGAGCACCACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGGGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGT‑GAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATcn < 1:46904/251‑2 (MQ=255) | CGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCCACCGATGGATAATCGGGTACGCGAGGCTTGTCAGTACATCAGCGATCACCTGGCAGACAGCAATTTTGATATCGCCAGCGTCGCACAGCATGTTTGCTTGTCGCCGTCGCGTCTGTCACATCTTTTCCGCCAGCAGTTAGGGATTAGCGTCTTAAGCTGGCGCGAGGACCAACGCATTAGTCAGGCGAAGCTGCTTTTGAGCACTACCCGGATGCCTATCGCCACCGTCGGTCGCAATGTTGGTTTTGACGATCAACTCTATTTCTCGCGAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGT‑GAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCC > NC_000913/70841‑71323 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |