Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 28,902 | 0 | G | A | 66.7% | 47.2 / 33.6 | 39 | A177T (GCG→ACG) | dapB | 4‑hydroxy‑tetrahydrodipicolinate reductase |
Reads supporting (aligned to +/- strand): ref base G (6/7); new base A (11/15); total (17/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CACCGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGTGTG > NC_000913/28823‑28979 | cacCGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGGc < 2:6065/81‑1 (MQ=255) cGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTAGCAATGGGAGAGACGAt < 2:6066/81‑1 (MQ=255) gAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCAc < 2:6067/81‑1 (MQ=255) ttattGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGccc < 1:6068/81‑1 (MQ=255) tattGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGGCGATCGCCCACGCCCt > 1:6130/1‑81 (MQ=255) tattGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCt > 1:6119/1‑81 (MQ=255) attGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCt < 2:6069/80‑1 (MQ=255) gAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGa < 2:6070/80‑1 (MQ=255) aaGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGAt < 1:6072/80‑1 (MQ=255) aaGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGAt < 2:6071/80‑1 (MQ=255) atcatAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATc < 2:6075/81‑1 (MQ=255) atcatAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATc < 2:6074/81‑1 (MQ=255) atcatAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATc < 1:6073/81‑1 (MQ=255) cATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGAt > 2:6122/1‑81 (MQ=255) cATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGAt > 2:6126/1‑81 (MQ=255) aTAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGAtt < 1:6077/81‑1 (MQ=255) aaaGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTgc > 2:6118/1‑81 (MQ=255) gTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGcgcg > 2:6137/1‑81 (MQ=255) tGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGt > 1:6136/1‑81 (MQ=255) gATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTc > 2:6128/1‑81 (MQ=255) aTGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCt < 2:6078/81‑1 (MQ=255) aTGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTc > 2:6145/1‑80 (MQ=255) aTGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTc < 1:6079/80‑1 (MQ=255) gcCGTCAGGCACCGCACTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAg > 1:6131/1‑81 (MQ=255) cGTCAGGCACCGCACTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTc > 2:6146/1‑81 (MQ=255) cGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTc > 1:6125/1‑81 (MQ=255) cGTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTc < 1:6080/81‑1 (MQ=255) gTCAGGCACCGCACTGGCAATGGGAGAGACGCTCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCg > 2:6144/1‑81 (MQ=255) gTCAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCg > 1:6135/1‑81 (MQ=255) cAGGCACCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGt > 1:6123/1‑80 (MQ=255) cACCGCACTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAgg < 1:6081/81‑1 (MQ=255) aCCGCACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGc < 1:6082/81‑1 (MQ=255) cGCACTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCca > 1:6147/1‑81 (MQ=255) cACTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCcaca < 2:6083/81‑1 (MQ=255) cTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACAcc > 1:6134/1‑81 (MQ=255) cTGGCAATGGGAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACAcc < 1:6084/81‑1 (MQ=255) gCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCg < 2:6085/79‑1 (MQ=255) tGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAAc < 1:6086/81‑1 (MQ=255) ggAGAGACGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACg < 1:6087/80‑1 (MQ=255) gagagaCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACgtg > 1:6150/1‑81 (MQ=255) gagaCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGtgtg > 2:6142/1‑81 (MQ=255) | CACCGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGTGTG > NC_000913/28823‑28979 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |