Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,352,754 | 0 | A | C | 45.5% | 18.7 / 15.2 | 22 | V91G (GTG→GGG) | arcB | sensory histidine kinase ArcB |
Reads supporting (aligned to +/- strand): ref base A (10/2); new base C (0/10); total (10/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.04e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTA > NC_000913/3352693‑3352815 | attatCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCCTCTCCTCCAGTTTTTGCACCa > 1:977599/1‑65 (MQ=255) attatCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCa > 1:3189512/1‑65 (MQ=255) attatCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCa > 1:733806/1‑65 (MQ=255) ttatCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCCCCTCCAGTTTTTGCCCCAg < 1:683952/65‑1 (MQ=255) tatCTTTTAACTGAACGTTGAGGCTCAAACCGCGCTCGCGCCTCCCCTCCATTTTTTGCCCCAGc < 1:40712/65‑1 (MQ=255) cTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGt > 1:56392/1‑65 (MQ=255) cTTTTAACTGAACGTTGAGGCTCAAACCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGt > 1:3950354/1‑65 (MQ=255) ttttACCTGACCGTGGGGGCTCAAATCGCGCTCGCCCATCTCCTCCAGTTTTTGCCCCAGCCGTg < 1:3183649/65‑1 (MQ=255) aaCTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACac > 1:379348/1‑64 (MQ=255) aCTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGa > 1:1048280/1‑65 (MQ=255) ggAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACg < 1:419474/64‑1 (MQ=255) ggTTGAGCCTCAATTCCCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGt < 1:2454160/64‑1 (MQ=255) cGTTGAGGCTCAATTCGCGCTCGCCCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGt < 1:385061/65‑1 (MQ=255) cGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCATTTTTTGCCCCAGCCGTGACAGACGTTGt < 1:930915/65‑1 (MQ=255) cGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGt < 1:1853231/65‑1 (MQ=255) cGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGt < 1:1169673/65‑1 (MQ=255) gTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTc > 1:1668443/1‑65 (MQ=255) ggCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGAc > 1:3174415/1‑65 (MQ=255) gcgcAAGTCGCGCTCGCGCCTCTCCCCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGAct < 1:2859455/62‑1 (MQ=255) gCTCAACTCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGAct < 1:367895/65‑1 (MQ=255) gCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGAct < 1:1144171/65‑1 (MQ=255) gCTCAAACCGCGCTCGCGCTTCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGAct < 1:1018552/65‑1 (MQ=255) cccAAACCGCCCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGActc < 1:3469752/63‑1 (MQ=255) cTCAAATCGCGCTGGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGActc < 1:2706041/65‑1 (MQ=255) cTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGActc < 1:553459/65‑1 (MQ=255) aaccgcgctcgcgcATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGACtcctcc < 1:15510/62‑1 (MQ=255) ctcgccctcgcgcATCCCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCa < 1:1076530/64‑1 (MQ=255) ctcgcgcTTCTCCGCCATTTTTGGCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCt < 1:1090147/65‑1 (MQ=255) ctcgcgcATCCCCCCCATTTTGCGCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCt < 1:811142/65‑1 (MQ=255) ctcctcCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGacca > 1:3919741/1‑65 (MQ=255) ctcctcCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGAAGacca > 1:3930442/1‑65 (MQ=255) cctccATTTTTTCCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGaccacc < 1:2357597/65‑1 (MQ=255) cctccAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGaccacc > 1:2608574/1‑65 (MQ=255) ccAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGAt < 1:2750627/65‑1 (MQ=255) ccAGTTTTGCCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGAt < 1:2247201/65‑1 (MQ=255) tGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTa > 1:388503/1‑65 (MQ=255) | ATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTA > NC_000913/3352693‑3352815 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |