| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,352,754 | 0 | A | C | 48.1% | 7.0 / 26.9 | 27 | V91G (GTG→GGG) | arcB | sensory histidine kinase ArcB |
| Reads supporting (aligned to +/- strand): ref base A (12/2); new base C (0/13); total (12/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.98e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAA > NC_000913/3352688‑3352811 | gCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCAc > 1:1609881/1‑68 (MQ=255) atatTATCTTTTAACTGAACGTTGAGGCTCAAAGCGCGCGCGCGCATCTCCTCCAGTTTTTGCCCCAg < 1:3851569/68‑1 (MQ=255) attatCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGcc > 1:207092/1‑68 (MQ=255) attatCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGcc > 1:1391354/1‑68 (MQ=255) atCTTTTAACTAAACGTTGAGGCTAAAACCGCGCTCGCGCTTCTCCTCCAGTTTTTGCCCCAGCCGTg < 1:2360753/68‑1 (MQ=255) tCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCGCCAGTTTTTGCCCCAGCCGTGa < 1:4182290/68‑1 (MQ=255) cTTTTACCTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGAc > 1:2176777/1‑68 (MQ=255) cTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCATCAGCCGTGAc > 1:4055209/1‑68 (MQ=255) cTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGAc > 1:1279112/1‑68 (MQ=255) cTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGAc > 1:3048895/1‑68 (MQ=255) cTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGAc > 1:1645066/1‑68 (MQ=255) tAACTGAACTTTGGGGCACAAATCGCGCTCGCGCATCTCCTCCATTTTTTGCCCCAGCCGTGACAGAc < 1:3947132/68‑1 (MQ=255) tAACTGAACGTTGAGGCTCAAATCGCCCCCGCGCATCTCCCCCAGTTTTTGCCCCAGCCGTGACAGAc < 1:151246/68‑1 (MQ=255) aaCTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACg > 1:1822655/1‑68 (MQ=255) cTGAACGTTGACGCTCAAATCGCGCTCGCGCCTCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGtt > 1:2140960/1‑68 (MQ=255) tGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTg > 1:2015429/1‑68 (MQ=255) gagcGTAGAGGCTCAAACCGCGCGCGATTACCTCCACGAGTTTTTGCCCCAGCCGTGACAGACGTTGt < 1:1392034/65‑1 (MQ=255) gaccGTTGGGGTTCAAATCGCGCTCGCGCATCTCCCCCATTTTTTGCACCAGCCGTGACAGACGTTGt < 1:2159873/65‑1 (MQ=255) gaccGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGt < 1:213040/65‑1 (MQ=255) gaccGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGt < 1:31986/65‑1 (MQ=255) gAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGt < 1:1638655/68‑1 (MQ=255) gAACGTTGAGGCTCAAATCGCCGTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGt < 1:4095350/68‑1 (MQ=255) ggggCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGGCGTTGTCGTGAct < 1:1807111/66‑1 (MQ=255) gAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGAct < 1:3846501/68‑1 (MQ=255) gAGGCTCAAATCGCGCTCGCGCATCCCCCCCATTTTTGGCCCCAGCCGTGACAGACGTTGTCGTGAct < 1:4042498/68‑1 (MQ=255) aGGCTCAAATCGCGCTCGCGCTTCTCCTCCATTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGActc < 1:1483813/68‑1 (MQ=255) ggCTCAAATCGCGCCCGCGCGCCTCCTCCAGTTTTGCCGCCAGCCGTGACAGACGTTGTCGTGActcc < 1:26291/68‑1 (MQ=255) gCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGActcct > 1:997958/1‑68 (MQ=255) gCTAAAACCGCGCTCCCGCATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGActcct < 1:2683289/68‑1 (MQ=255) cTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGActcctc > 1:306169/1‑68 (MQ=255) agtcgcgctcgcccCTCCCCTCCTTTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGt < 1:1780314/66‑1 (MQ=255) aTCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGtt > 1:1316224/1‑68 (MQ=255) gcgctcgcgcATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCt < 1:1210779/68‑1 (MQ=255) gccctcgcgcCTCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCt < 1:1323922/65‑1 (MQ=255) ccgctcgcgcATCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCt < 1:1771593/67‑1 (MQ=255) cgcTTCTCCTCCAGTTTTTGCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCgacgac < 1:2325984/68‑1 (MQ=255) cattttttGCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGaa < 1:2314117/65‑1 (MQ=255) cAGTTTTTCCCCCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGaa < 1:3822095/68‑1 (MQ=255) | GCAATATTATCTTTTAACTGAACGTTGAGGCTCAAATCGCGCTCGCGCATCTCCTCCAGTTTTTGCACCAGCCGTGACAGACGTTGTCGTGACTCCTCCAGTTGCTCGACGACCACCGATAGAA > NC_000913/3352688‑3352811 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |