Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,825,113 | 0 | G | T | 72.7% | 28.1 / 17.0 | 22 | intergenic (‑176/+27) | ves/spy | HutD family protein Ves/ATP‑independent periplasmic chaperone |
Reads supporting (aligned to +/- strand): ref base G (4/2); new base T (6/10); total (10/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.48e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGA > NC_000913/1825039‑1825191 | aCAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCt > 2:304325/1‑81 (MQ=255) cAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGTCGGTCtt < 2:304303/81‑1 (MQ=255) cAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGTCGGTCt < 1:304302/80‑1 (MQ=255) tttGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGTCGGTCTTAAGt < 1:304304/79‑1 (MQ=255) tGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGTCGGTCTTAAGttttt > 2:304328/1‑81 (MQ=255) agaagagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtttttt < 2:304305/80‑1 (MQ=255) agagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGTCGGTCTTAAGTTTTTTGGCt < 2:304306/81‑1 (MQ=255) aaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGTCGGTCTTAAGTTTTTTGGCTGAAAGattat < 2:304307/81‑1 (MQ=255) aaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGTCGGTCTTAAGTTTTTTGGCTGAAAGattatt > 2:304324/1‑81 (MQ=255) aCCGCCGATCCTGTCCACCGAATTACTGCAAGGTAGTGGACAAGACCGTCGGTCTTAAGTTTTTTGGCTGAAAGATTATTc < 1:304308/81‑1 (MQ=255) gccgATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag > 2:304333/1‑80 (MQ=255) cgATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGTCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagcag > 1:304327/1‑81 (MQ=255) gATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGTCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGt > 1:304332/1‑81 (MQ=255) cACCGCATTACTGCAAGGTAGTGGACAAGACCGTCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAt < 1:304309/81‑1 (MQ=255) gCATTACTGCAAGGTAGTGGACAAGACCGTCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa > 1:304331/1‑81 (MQ=255) gCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGc > 2:304329/1‑81 (MQ=255) cAAGGTAGTGGACAAGACCGTCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGcc < 2:304310/81‑1 (MQ=255) gTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTg > 1:304330/1‑81 (MQ=255) tAGTGGACAAGACCGTCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg < 1:304311/81‑1 (MQ=255) aCAAGACCGTCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTc > 2:304336/1‑81 (MQ=255) aCCGTCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCa < 1:304312/81‑1 (MQ=255) aCCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCa < 1:304313/81‑1 (MQ=255) cggcggTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGa < 2:304314/81‑1 (MQ=255) | ACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGA > NC_000913/1825039‑1825191 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |