| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,233,922 | 0 | G | T | 44.1% | 16.8 / 33.5 | 34 | Q75K (CAA→AAA) | higA | antitoxin/DNA‑binding transcriptional repressor HigA |
| Reads supporting (aligned to +/- strand): ref base G (7/12); new base T (15/0); total (22/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.40e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACA > NC_000913/3233859‑3233985 | aaaGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGa < 1:1166128/65‑1 (MQ=255) aGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTTAGc > 1:1205064/1‑65 (MQ=255) aGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGc > 1:3836026/1‑65 (MQ=255) aacccATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGc < 1:3539763/62‑1 (MQ=255) aaaCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGc < 1:3212700/65‑1 (MQ=255) aaaCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGc < 1:458738/65‑1 (MQ=255) aaaCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGc < 1:1399936/65‑1 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCAt > 1:2069815/1‑65 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCAt > 1:1324771/1‑65 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAt > 1:313451/1‑65 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAt > 1:3354526/1‑65 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAt > 1:496197/1‑65 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATAGCAt > 1:3035271/1‑65 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCCTGCGCCACGGCAt > 1:387825/1‑65 (MQ=255) aCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCCTGAGCCATggtct > 1:2327382/1‑62 (MQ=255) ccATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCAtt > 1:707782/1‑65 (MQ=255) ccATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAtt > 1:702681/1‑65 (MQ=255) ccATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTACGCCATGGCAtt > 1:3914291/1‑65 (MQ=255) atTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTTGACCAACGCATTAAAt > 1:339641/1‑65 (MQ=255) gATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTcc > 1:1254233/1‑65 (MQ=255) tCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTTAGCCCTGGCATTAAATTCCGc > 1:976887/1‑65 (MQ=255) tCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGc > 1:1232425/1‑65 (MQ=255) tCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGc > 1:1520202/1‑65 (MQ=255) tCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGc > 1:3450961/1‑65 (MQ=255) aaGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAAtt < 1:646472/65‑1 (MQ=255) aGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTTTGCCATGGCATTAAATTCCGCAAATTc > 1:3417634/1‑65 (MQ=255) aGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTc > 1:3833564/1‑65 (MQ=255) gggTACGAATCACGGCTATACCGCCAGGCATGGCTTTAGCCAAGGCATTAAATTCCGGCAATTCg > 1:2491121/1‑65 (MQ=255) aCGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGgcg < 1:928112/65‑1 (MQ=255) aCGAATCACGGCTATACCGCCAGGCATGGCTTGAGACATGGCATTAAATTCCCCAAATTCGGgcg > 1:1408264/1‑65 (MQ=255) aaTCACGGCTATACCGCCAGGCATGGCTTTAGCCCTGGCATTAAATTCCGCAAATTCGGGCGCTg > 1:3034151/1‑65 (MQ=255) aaTCACGGCTATACCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTg > 1:1390507/1‑65 (MQ=255) aaTCACGGCTATACCGCCAGGCATGGCTTGGGCCATGGCATTAAATTCCGCAAATTCGGGCGCTg > 1:1946601/1‑65 (MQ=255) aaTCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTg > 1:3023185/1‑65 (MQ=255) aTCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGa > 1:2685346/1‑65 (MQ=255) gCTATACCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAttcttc > 1:3446082/1‑65 (MQ=255) ataCCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTGCCCa < 1:1545466/65‑1 (MQ=255) aCCGCCAGGCATGGCTTTAGGCCTGGCATTAAATTCCGCAAAATCGGGCGCGGATTCTTTCGAcg > 1:2634374/1‑65 (MQ=255) aCCGCCAGGCATGGCTTTAGCCATGGCATTAACTTCCGCAAATTCGGGCGCTGATTCTTCCCAcg > 1:2133918/1‑65 (MQ=255) aCCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcg > 1:3374907/1‑65 (MQ=255) aCCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcg > 1:1108365/1‑65 (MQ=255) aCCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcg > 1:3978806/1‑65 (MQ=255) aCCGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcg > 1:3941989/1‑65 (MQ=255) aCCGCCAGGCATGGCTTGAGCCATGGCATTCAATTCCGCAAATTCGGGCGCTGATTCTTCCCaag > 1:3871120/1‑63 (MQ=255) cGCCAGGCATGGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcgcg > 1:1161928/1‑65 (MQ=255) cGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcgcg > 1:3976618/1‑65 (MQ=255) gCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCgg > 1:142407/1‑65 (MQ=255) cAGGCATCGCTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGtt > 1:1577356/1‑65 (MQ=255) gCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTAtt < 1:1784459/65‑1 (MQ=255) gCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGc < 1:1201694/65‑1 (MQ=255) gCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGc < 1:2535984/65‑1 (MQ=255) cTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGca > 1:76712/1‑65 (MQ=255) cTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGca > 1:903489/1‑65 (MQ=255) cTTTAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGCTATTTTGGca > 1:1724759/1‑65 (MQ=255) cTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGca > 1:1799589/1‑65 (MQ=255) cTTGAGCCATGGCATTAAATTCCGCAAATTCGGGACCTGATTCTTCCCAACCGGTTATTTTggta > 1:3790956/1‑63 (MQ=255) tGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGcaca < 1:1279018/65‑1 (MQ=255) | AAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACA > NC_000913/3233859‑3233985 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |