Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,352,754 | 0 | A | C | 50.0% | 3.0 / 15.4 | 18 | V91G (GTG→GGG) | arcB | sensory histidine kinase ArcB |
Reads supporting (aligned to +/- strand): ref base A (7/2); new base C (0/9); total (7/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.26e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGA > NC_000913/3352693‑3352817 | attatCTTTTAACTGAACGTTGAGGCTCAAATCGCTCTCGCGCATCTCCTCCAGTTTTTGCACCa > 1:658350/1‑65 (MQ=255) attatCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCa > 1:2424366/1‑65 (MQ=255) attatCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCa > 1:2112099/1‑65 (MQ=255) attatCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCa > 1:3457417/1‑65 (MQ=255) cTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGt > 1:3614306/1‑65 (MQ=255) cTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGt > 1:2049604/1‑65 (MQ=255) cTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGt > 1:2081439/1‑65 (MQ=255) cTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGt > 1:2281488/1‑65 (MQ=255) ttttAGCTGAGCGTGGAGGCTCAAGCCGCGCCCGCGTATCTCCTCCAGTTTTTGCACCAGCCGTg < 1:3179562/65‑1 (MQ=255) ttttAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTg < 1:1928123/65‑1 (MQ=255) ttttAACAGAACTTTGAGGCTCAAACCGCGCTCGCGCATCCCCTCCAGTTTTTGCCCCAGCCGTg < 1:2509312/65‑1 (MQ=255) aaCTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAg > 1:303456/1‑65 (MQ=255) cTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGAc < 1:3442555/65‑1 (MQ=255) cGTTGGGGCTCAAACCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGt < 1:1954647/65‑1 (MQ=255) cGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGt < 1:2740609/65‑1 (MQ=255) cGTTGAGGCTCAAATCGCGCTCGCGCATCCCCTCCATTTTTTCCCCCAGCCGTGACAGACGTTGt < 1:2900193/65‑1 (MQ=255) cGTTGAGGCTCAAATCGCCCTCGCGCTTCTCCTCCATTTTTTGCACCAGCCGTGACAGACGTTGt < 1:3796255/65‑1 (MQ=255) cGTTGAGGCTAAAACCGCGCTCGCCCATCTCCGCCAGTTTTGGCCCCAGCCGTGACAGACGTTGt < 1:2609433/65‑1 (MQ=255) cGTTGAGCCACAAATCGCGCTCGCGCATCCCCTCCATTTTTTGCCCCAGCCGTGACAGACGTTGt < 1:1036062/65‑1 (MQ=255) aGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGa < 1:1348979/65‑1 (MQ=255) gCTCAAATCGCGCTCGCGCATCTCCTCCATTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGAct < 1:2501809/65‑1 (MQ=255) gCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGAct < 1:3821145/65‑1 (MQ=255) gCTAAAATTGCCCTCTCGCATCTCCCCCATTTTTTCCACCAGCCGTGACAGACGTTGTCGTGAct < 1:912267/65‑1 (MQ=255) gCTAAAACCGCGCTGGCGCATCTCCTCAAGTTTTTCCACCAGCCGTGACAGACGTTGTCGTGAct < 1:1568914/65‑1 (MQ=255) ccTCAAACCGCGCTCGCGCATCCCCCCCACTTTTTCCCCCAGCCGTGACAGACGTTGTCGTGAct < 1:3643829/64‑1 (MQ=255) cTCAAACCGCGCTCGCCCATCCCCTCCATTTTTTCCCCCAGCCGTGACAGACGTTGTCGTGActc < 1:3224019/65‑1 (MQ=255) tCAATTCGCGCCCGCGCATCTCCTCCAGTTTTTCCACCAGCCGTGACAGACGTTGTCGTGActcc < 1:3178003/65‑1 (MQ=255) tCAACTCGCCCTCGCGCACCCCCCCCTGTTTTGCCCCCAGCCGTGACAGACGTTGTCGTGActcc < 1:1921330/65‑1 (MQ=255) cAAAGCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGActcct < 1:843137/65‑1 (MQ=255) aTCGCACTCGCGCCTCTCCCCCATTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCa < 1:2661051/65‑1 (MQ=255) ccgcTTCTCCCCCAGTTTTTCCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTcga < 1:235618/64‑1 (MQ=255) ctcctcCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGacca > 1:2294676/1‑65 (MQ=255) cttccAGTTTTTGGCCCAGCCGTGACAGACGTTGTCGTGCCTCCTCCAGTTGCTCGACGaccacc < 1:1748005/63‑1 (MQ=255) cctccGTTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGcccacc < 1:2332601/65‑1 (MQ=255) ccAGTTTTTCCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGAt < 1:2219321/65‑1 (MQ=255) cACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGa > 1:900303/1‑65 (MQ=255) | ATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAAAGTAGA > NC_000913/3352693‑3352817 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |