Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,927,126 | 0 | A | T | 21.4% | 54.8 / 11.8 | 28 | intergenic (‑123/‑29) | asnC/asnA | transcriptional activator of asnA; autorepressor/asparagine synthetase A |
Reads supporting (aligned to +/- strand): ref base A (10/12); new base T (2/4); total (12/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.73e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAACGTCTGGGGCTGATCGAAGTCCAGGCGCCGATTCTTAGCCG > NC_000913/3926982‑3927270 | tGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATaa < 5:113680/148‑1 (MQ=255) gATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAAc < 8:100818/149‑1 (MQ=255) gATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAAc < 1:120792/149‑1 (MQ=255) aGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTATGATTAAGCAACaaa < 8:36518/148‑3 (MQ=255) gATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTATGATTAAGCAACaaaaa < 6:67597/149‑5 (MQ=255) ttcTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTTCCAgacaa > 3:114491/1‑149 (MQ=255) tttAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACAGCTTACATTGCCAAACAACGTCa > 5:99192/1‑148 (MQ=255) tAGCCTTCCTGTTTAATGCCTCAAAAGCGAGTTAGGCGTTTTATAGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATACAAAATGAAAACCGCTTCCATTGCCAGACAACGTCaaa < 4:114491/148‑1 (MQ=255) tAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGaaa > 4:51066/1‑149 (MQ=255) gAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACg < 1:85749/126‑1 (MQ=255) gAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACg > 2:85749/1‑126 (MQ=255) ttAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCaa < 3:51066/149‑1 (MQ=255) gAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTtctc > 4:68599/1‑130 (MQ=255) gAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTtctc < 3:68599/130‑1 (MQ=255) aaTGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAACGTCt < 6:99192/148‑1 (MQ=255) gtCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTtctc > 7:142038/1‑111 (MQ=255) gtCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTtctc < 8:142038/111‑1 (MQ=255) tCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGt > 2:13324/1‑93 (MQ=255) tCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGt < 1:13324/93‑1 (MQ=255) tCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTtc < 2:94308/108‑1 (MQ=255) tCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTtc < 2:91532/108‑1 (MQ=255) tCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTtc > 1:94308/1‑108 (MQ=255) tCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTtc > 1:91532/1‑108 (MQ=255) tgttgCTTATGATTAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAAc > 5:38172/1‑117 (MQ=255) tgttgCTTATGATTAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAAc < 6:38172/117‑1 (MQ=255) tgttgCTTATGATTAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAACGTCTg < 5:136830/122‑1 (MQ=255) tgttgCTTATGATTAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAACGTCTg > 6:136830/1‑122 (MQ=255) aaTCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAATCCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGCACGACTGGGGCTGATCGAAGTCCAGGCGCCGATTCTTAGCCg > 7:34533/1‑149 (MQ=255) | TGTCGATCAGATAATTTTCCATAGGATTTTCTTATGCGGATTGATGATTCATTCTATTTTAGCCTTCTTTTTTAATGAATCAAAAGTGAGTTAGGCTTTTTATTGAATGATTATTGCATGTGTGTCGGTTTTTGTTGCTTAATCATAAGCAACAGGACGCAGGAGTATAAAAAATGAAAACCGCTTACATTGCCAAACAACGTCAAATTAGCTTCGTGAAATCTCACTTTTCTCGTCAACTGGAAGAACGTCTGGGGCTGATCGAAGTCCAGGCGCCGATTCTTAGCCG > NC_000913/3926982‑3927270 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |