Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,274,940 | 0 | A | T | 30.0% | 53.7 / 10.9 | 30 | intergenic (‑33/‑342) | garP/garD | galactarate/D‑glucarate transporter GarP/GarD |
Reads supporting (aligned to +/- strand): ref base A (6/15); new base T (9/0); total (15/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.00e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGATTTGCATCAGCAAGTAGGCCCAGCCAAAAGCGGAGAAGATGTAACCCATCGAAACCGCACTTAACTGCAACTCTTTTGCCACTTCGGTACCAGCAATAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGTAACAATAAAAATAATCAGTAATATTAAATAGCGGGTATGCACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTATTTAAATATTATTTTGAGCATATGCACATAAGGTTGCGCGCTAAAGCACAGATTTGCGCTTTACCTTACCGGGCGGCACTGCAATCCCTGAAATGATTGACATTGATCAC > NC_000913/3274693‑3275172 | gggATTTGCATCAGCAAGTAGGCCCAGCCAAAAGCGGAGAAGATGTAACCCATCGAAACCGCACTTAACTGCAACTCTTTTGCCACTTCGGTACCAGCAATAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGTAACAATAAAAATAATCAGTAATATTAAATAGCGGGTATGCACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTttatta < 1:443086/251‑1 (MQ=255) aTTTGCACCAGCAAGTAGGCCCAGCCAAAAGCGGAGAAGATGTAACCCATCGAAACCGCACTTCACTGCAACTCTTTTGCCACTTCGGTACCAGCACTAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGTAAAAAAAAAAATAATCACATATATTTTATAGCGGGGGTGCCCGCCCTTTTTTTTTTTGTCAAAGGTGTCGCA‑AATCATTTTTTTT‑ACCTCCGGGAATTCTACTGTTTTttattattt > 2:445507/1‑250 (MQ=255) aGTAGGCCCAGCCAAAAGCGGAGAAGATGTAACCCATCGAAACCGCACTTAACTGCAACTCTTTTGCCACTTCGGTACCAGCAATAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGTAACAATAAAAATAATCAGTAATATTAAATAGCGGGTATGCACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGtt < 1:315326/251‑1 (MQ=255) aaaGCGGAGAAGATGTAACCCATCGAAACCGCACTTAACTGCAACTCTTTTGCCACTTCGGTACCAGCAATAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGTAACCATAAAAAAAATCAGGAATATTTAATAGCGGGGGTGCCCGCCTTTCTTTTTTTTGTTAAAGGGGGC‑CCGGATCCTTTTTTTT‑TTCTCGGGTACTTTTTCTGTTTTTTTTTTTTTT‑GGGGGAGGGGTTATTTTTGTGttttt > 2:340883/1‑249 (MQ=255) ggAGAAGATGTAACCCATCGAAACCGCACTTAACTGCAACTCTTTTGCCACTTCGGTACCAGCAATAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGTAAAAAAAAAAATAAACAGTAATATTAAATAGAGGGGATGCACGCCTTTCTTTTTTTTGTCAACAGTGTC‑CAGAAATATTTTTTTT‑ACCTCGGGGACTTTTTTTTATTTTTTTTTTTTT‑TGGGAAGGTGttattttttaggtttatttat > 2:372208/1‑236 (MQ=255) gaagaTGTAACCCATCGAAACCGCACTTAACTGCAACTCTTTTGCCACTTCGGTACCAGCAATAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGTAACAATAAAAATAATCAGTAATATTAAATAGCGGGTATGCACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGcc < 1:340883/251‑1 (MQ=255) aTCGAAACCGCACTTAACTGCAACTCTTTTGCCACTTCGGTACCAGCAATAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGTAACAATAAAAATAATCCGTAATATTAAATAGCGGGTATGCACGCCCTTTTTTTTTTTGTCAACGGTGTC‑CAGAATCATTTTATTT‑TCCCCGGGGACTTATTCTGATTTTTTTTTTTAT‑TGGGAAGGGGTTATTTATGAGTTTTCTTTTTGTCGTAACGACAATg > 2:403250/1‑250 (MQ=255) cGCACTTAACTGCAACTCTTTTGCCACTTCGGTACCAGCAATAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGGAACCATAAAAATAATCAGTAATATTAAATAGCGGGGATGCACGCCATTTTTTTTTTTGTTAACGGTGTC‑CAGGAACATTTTTTTTTACC‑CGGGGACTTTTTCTGCTTTTTTTTTTTTT‑GGGGGAGGGGTTATTTATTAGTTTCTTTTATTGCGTAACGACAATGAACTCgggg > 2:440216/1‑250 (MQ=255) tAACTGCAACTCTTTTGCCACTTCGGTACCAGCAATAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGTAAAAATAAAAAAAAACAGTAATATTTAATAGCGGGGGTGCACGCCCTTCTTTTTTTTGTCAACGGTGTT‑CAGAATCATTTTATTT‑TCCTCCGGTACTTTTTCTGGTTTTTTTTTTTAT‑GGGGGAGGGGTTTTTTTTGAGTTTCATTTATTCCGTAACGCCAATGAACTCGGGGAtttgt > 2:171655/1‑248 (MQ=255) ctctTTTGCCACTTCGGTACCAGCACTAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGTAACAATAAAAATAAAAAGGTATATTTAAAAGCGGGTGTTCAAGCCCTTTTTTTTTTTGTCCACGGGGGC‑CAGAATAATTTTTTTT‑TCCTCGGGGAATTTTTTTGATTTTTTTTTTTTT‑GGGGGAGGTGTTTTTTATTAGTTTCATTTTTTTCGTTAAGACAATGAACTCGGGAATTTGTATaaacg > 2:423004/1‑246 (MQ=255) ttGCCACTTCGGTACCAGCAATAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGTAACAATAAAAATAATCAGTAATATTAAATAGCGGGTATGCACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGtata < 2:147415/240‑1 (MQ=255) ttGCCACTTCGGTACCAGCAATAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGTAACAATAAAAATAATCAGTAATATTAAATAACGGGGATGCACGCCCTTTTTTTTTTCGTCAACGGGGTC‑CAGAATCATTTTTTTT‑TCCCCGGGGACTTTTGCTGGTTTTTTTTTTTTT‑GGGGAAGGGGTTATTTTTGAGTTTCATTTTTTCCGTAACGACAATGAACTCGGGAATTAGtata > 1:147415/1‑240 (MQ=255) gCCACTTCGGTACCAGCAATAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGTAACAATAAAAATAATCAGTAATATTAAATAGCGGGTATGCACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGaa < 1:440216/251‑1 (MQ=255) gTTGCACGATCGGCGTAGTTAACGGCGGTAACAATAAAAATTATTAGTAAAATTTAATAAAGGGGATGCGCGCCTTTTTTTTTTTTGTCAACAGGGGC‑CCGAATCATTTTTTTT‑TACTCGGGTACTTTTGCTGTTTTTTTTTTTTAT‑GGGGGAGGGGTTTTTTTTGAGTTTTTTTTTTTCCGTTACGACAATGAAATCGGGGATTAGTATAAGCAACGCGAGAATAAAAAACATTTTGCCAATGCtaattt > 2:107219/1‑251 (MQ=255) cACGATCGGCGTAGTTAACGGCGGTAACAATAAAAATAATCAGTAATATTAAATAGCGGGTATGCACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAATCAGCGCGAGAATAATAATCATTGTGCAAATGCtaatttaatt < 1:445507/251‑1 (MQ=255) ggCGTAGTTAACGGCGGTTACAATAAAAAACAACACTAATATTTAATAGCGGGTATGCACGCCTTTCTTTTTTTTTTCAACAGTGTC‑AAGAAATATTTTATTT‑AACTCGGGTACTTTTGTTGTTTTTTTTTATTTTTGGGGAGGTTTTTTTTTTTGGGTTT‑ATTTTTTGCGTCACGACAAAGAAATCTGGAATTATTATAAGAAGAGCGAGGATACAAATCATTGTGCAAACTCTAATttaattacatcat > 2:143079/1‑245 (MQ=255) gTAGTTAACGGCGGTAACAATAAAAATAATCAGTAATATTAAATAGCGGGTATGCACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTAttt < 2:68634/251‑1 (MQ=255) aaaaaTAATCAGTAATATTAAATAGCGGGTATGCACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTAttt < 2:496314/230‑1 (MQ=255) aaaaaTAATCAGTAATATTAAATAGCGGGTATGCACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTAttt > 1:496314/1‑230 (MQ=255) aaaaaTAATCAGTAATATTAAATAGCGGGTATGCACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTATTTAAATATTATTTGGAGCatat < 1:423004/250‑1 (MQ=255) gTAATATTAAATAGCGGGTATGCACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTATTTAAATATTATTTTGAGCATATGCACATAAGGtt < 1:372208/251‑1 (MQ=255) aaaTAGCGGGGATGCACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTATTTAAATATTATTTTGAGCATATGCACATAAGGTTGCGCGCt > 1:72696/1‑250 (MQ=255) cACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTATTTAAATATTATTTTGAGCATATGCACATAAGGTTGCGCGCTAAAGCACAGATTTgc > 2:194693/1‑251 (MQ=255) tttCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTATTTAAATATTATTTTGAGCATATGCACATAAGGTTGCGCGCTAAAGCACAGATTTGCGCTttaccttacc < 1:4865/250‑1 (MQ=255) cAACGGTGTC‑CAGAATAATTTAATTT‑ACCTCGGGTACTTATGCTGATTTTTACTATTAT‑GGGGAAGGTCTTATTTATGAGTTTCCTTTATGCCGTAACGACAATGCACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTATTTAAATATTATTTTGAGCATCTGCACATAAGGTTGCGCGCTAAAGCACAGATTTGCGCTTTCCCTTACCGggcggc > 2:528310/1‑251 (MQ=255) aTCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTATTTAAATATTATTTTGAGCATATGCACATAAGGTTGCGCGCTAAAGCACAGATTTGCGCTTTACCTTACCGGGCGGCACTGCAATCCCTg > 1:271276/1‑250 (MQ=255) ttattt‑tCCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTATTTAAATATTATTTTGAGCATATGCACATAAGGTTGCGCGCTAAAGCACAGATTTGCGCTTTACCTTACCGGGCGGCACTGCAATCCCTGAAATGa < 2:6261/250‑1 (MQ=255) aCCTCGGGTACTTATGCTGATTTTTATTATTAT‑TGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTAGAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAAGTAATACTATTTAAATATGATTTTGAGCATATGCACATAAGGTTGCGCGCTAAAGCACAGATTTGCTCTTTACCTTACCGGGCGGCACTGCAATCCCTGAAATGATTGACAt < 1:266368/251‑1 (MQ=255) gggTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTATTTAAATATTATTTTGAGCATATGCACATAAGGTTGCGCGCTAAAGCACAGAGTTGCGCTTTACCTTACCGGGCGGCACTGCAATCCCTGCAATGATTGACATTGATc < 2:72696/251‑1 (MQ=255) ggTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTATTTAAATATTATTTTGAGCATATGCACATAAGGTTGCGCGCTAAAGCACAGATTTGCGCTTTACCTTACCGGGCGGCACTGCAATCCCTg > 2:414919/1‑232 (MQ=255) ggTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTATTTAAATATTATTTTGAGCATATGCACATAAGGTTGCGCGCTAAAGCACAGATTTGCGCTTTACCTTACCGGGCGGCACTGCAATCCCTg < 1:414919/232‑1 (MQ=255) gTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTATTTAAATATTATTTTGAGCATATGCACATAAGGTTGCGCGCTAAAGCACAGATTTGCGCTTTACCTTACCGGGCGGCACTGCAATCCCTGAAATGATGGACATTGATcac < 1:438046/251‑1 (MQ=255) | GGGATTTGCATCAGCAAGTAGGCCCAGCCAAAAGCGGAGAAGATGTAACCCATCGAAACCGCACTTAACTGCAACTCTTTTGCCACTTCGGTACCAGCAATAGACAGCGTTGCACGATCGGCGTAGTTAACGGCGGTAACAATAAAAATAATCAGTAATATTAAATAGCGGGTATGCACGCCTTTCTTTTTTTCGTCAACGGTGTC‑CAGAATCATTTTATTT‑ACCTCGGGTACTTATGCTGATTTTTATTATTAT‑GGGGAAGGTGTTATTTATGAGTTTCATTTATGCCGTAACGACAATGAACTCGGGAATTAGTATAAGCAGCGCGAGAATAATAATCATTGTGCAAATGCTAATTTAATTAATACTATTTAAATATTATTTTGAGCATATGCACATAAGGTTGCGCGCTAAAGCACAGATTTGCGCTTTACCTTACCGGGCGGCACTGCAATCCCTGAAATGATTGACATTGATCAC > NC_000913/3274693‑3275172 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |