| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,233,922 | 0 | G | T | 30.0% | 63.2 / 16.0 | 40 | Q75K (CAA→AAA) | higA | antitoxin/DNA‑binding transcriptional repressor HigA |
| Reads supporting (aligned to +/- strand): ref base G (16/12); new base T (12/0); total (28/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.39e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCA > NC_000913/3233860‑3233983 | aaGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAg < 1:4141098/65‑1 (MQ=255) aGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGc > 1:1762158/1‑65 (MQ=255) aaaCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGc < 1:4071719/65‑1 (MQ=255) aaaCCATATTGATCCATAAGGGTACGAATAACGGCTATACCGCCAGGCATGGCTTGAGCCATGGc < 1:3348024/65‑1 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCAt > 1:897998/1‑65 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCAt > 1:3167333/1‑65 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGATCCCTGGCAt > 1:846366/1‑65 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAt > 1:3128854/1‑65 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAt > 1:4038849/1‑65 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAt > 1:4080845/1‑65 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCACGGCAt > 1:486023/1‑65 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGACATGGcct > 1:1184257/1‑63 (MQ=255) ccATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCACGGCAtt > 1:2077350/1‑65 (MQ=255) aaGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAAtt < 1:4419870/65‑1 (MQ=255) gggTACGAATCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCg > 1:789778/1‑65 (MQ=255) gggTACGAATCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCg > 1:4200036/1‑65 (MQ=255) aCGAATCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGgcg > 1:1053509/1‑65 (MQ=255) aCGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGgcg < 1:4351794/65‑1 (MQ=255) aaTCACGGCTATACCGCCAGGCATGGCTTTAGGCATGGCATTAAATTCCGCCAATTCGGGCGCTg > 1:4060014/1‑65 (MQ=255) aaTCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTg > 1:498037/1‑65 (MQ=255) aaTCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTg > 1:4091909/1‑65 (MQ=255) aaTCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTg > 1:1314482/1‑65 (MQ=255) aTCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGa > 1:4257645/1‑65 (MQ=255) aTCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCCGGCGCGGa > 1:2480805/1‑65 (MQ=255) aTCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGa < 1:4312047/65‑1 (MQ=255) aTCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGa > 1:139095/1‑65 (MQ=255) gCTATCNCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAttcttc > 1:2576640/1‑65 (MQ=255) gCTATACCGCCAGGCATGGCTTTAGCCATGGGATTAAATTCCGCAAATTCGGGCGCTGAttcttc > 1:2212920/1‑65 (MQ=255) gCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAttcttc > 1:181084/1‑65 (MQ=255) gCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAttcttc > 1:640279/1‑65 (MQ=255) cTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTcc < 1:4074662/65‑1 (MQ=255) aCCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcg > 1:966622/1‑65 (MQ=255) aCCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcg > 1:459575/1‑65 (MQ=255) aCCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcg > 1:3482993/1‑65 (MQ=255) aCCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcg > 1:3270297/1‑65 (MQ=255) aCCGCCAGGCATGGCTTTAGCCATGGAATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcg > 1:1350460/1‑65 (MQ=255) aCCGCCAGGCATGGCTTGAGCCCTGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcc > 1:3470419/1‑64 (MQ=255) aCCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCCAATTCGGGCGCTGATTCTTCCCAcg > 1:4173114/1‑65 (MQ=255) aCCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAACTTCGGGCGCTGATTCTTCCCAcg > 1:391300/1‑65 (MQ=255) aCCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcg > 1:104788/1‑65 (MQ=255) aCCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTACCAcg > 1:3388485/1‑65 (MQ=255) aCCGCCAGGCAGGGCTTGNGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcg > 1:3232894/1‑65 (MQ=255) cGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcgcg < 1:326202/65‑1 (MQ=255) cGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcgcg < 1:470591/65‑1 (MQ=255) cGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcgcg > 1:3961571/1‑65 (MQ=255) ccAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGt < 1:3715503/65‑1 (MQ=255) ggCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCTGtctc > 1:786716/1‑62 (MQ=255) gCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTAtt < 1:477481/65‑1 (MQ=255) gCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGc < 1:1458470/65‑1 (MQ=255) cTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGca > 1:2467619/1‑65 (MQ=255) cTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGca > 1:747881/1‑65 (MQ=255) cTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGca > 1:3455955/1‑65 (MQ=255) cTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGca > 1:3268502/1‑65 (MQ=255) cTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGca > 1:245833/1‑65 (MQ=255) cTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGca > 1:861173/1‑65 (MQ=255) cTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTCTCCACGCGGTTATTTTGGca > 1:2474608/1‑65 (MQ=255) | AAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCA > NC_000913/3233860‑3233983 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |