| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 821,440 | 0 | C | T | 31.6% | 48.7 / 19.4 | 38 | F216F (TTC→TTT) | ybhM | Bax1‑I family protein YbhM |
| Reads supporting (aligned to +/- strand): ref base C (3/23); new base T (12/0); total (15/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.68e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGTGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTT > NC_000913/821376‑821497 | tgtgGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGcttc > 1:3385205/1‑65 (MQ=255)tgtgGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGcttc < 1:2604733/65‑1 (MQ=255) ggAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCttt < 1:951544/65‑1 (MQ=255) aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTATCACTCGATTACTAGATTATTGGCTTCttttg > 1:3495047/1‑63 (MQ=255) aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCtttttttt > 1:1118845/1‑63 (MQ=255) aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCtttttttt > 1:891031/1‑63 (MQ=255) aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCtttttttg > 1:985251/1‑63 (MQ=255) aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCtttttttg > 1:3100004/1‑63 (MQ=255) aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCtttttttg > 1:2009060/1‑63 (MQ=255) aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTgt > 1:987601/1‑64 (MQ=255) aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTgg > 1:190908/1‑65 (MQ=255) aGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTgg < 1:2469393/65‑1 (MQ=255) aGTCCGGTTCGTTGCGCTTTTACAGGCGCATTAACACTCTATTACTATTTTATCGGCtttttttg > 1:2970559/1‑63 (MQ=255) ccGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTtgggt > 1:965023/1‑63 (MQ=255) ccGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTtgggt > 1:3013619/1‑63 (MQ=255) ccGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTtgggt > 1:2373623/1‑63 (MQ=255) ccGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTTGGAt > 1:149291/1‑65 (MQ=255) ggTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATAc < 1:139604/65‑1 (MQ=255) cgcTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGa < 1:967648/65‑1 (MQ=255) cgcTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGa < 1:2189636/65‑1 (MQ=255) ttttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAg < 1:3324954/65‑1 (MQ=255) tttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc < 1:1769538/65‑1 (MQ=255) tttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGc < 1:2346069/65‑1 (MQ=255) ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:3641380/65‑1 (MQ=255) ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:917398/65‑1 (MQ=255) ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:1114905/65‑1 (MQ=255) ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:3456562/65‑1 (MQ=255) ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:2999100/65‑1 (MQ=255) ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:2707462/65‑1 (MQ=255) ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:1768702/65‑1 (MQ=255) ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:1899844/65‑1 (MQ=255) ttACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCt < 1:1272370/65‑1 (MQ=255) tACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCta < 1:2690283/65‑1 (MQ=255) aGGTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTTGGAAACTTGCCCCGGTAGGTATaa > 1:352290/1‑65 (MQ=255) ggTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAAc < 1:1497169/65‑1 (MQ=255) gTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTTGGATACTTGCCGCGAAAGCTATAAcc > 1:9448/1‑64 (MQ=255) gTGCATTAACACTCTATTACTATTTTATCGGCTTTTTTGGGATAATTGCCCCGGTAGCTATAACg > 1:334372/1‑65 (MQ=255) gTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACg < 1:1127152/65‑1 (MQ=255) acacTCTATTACTATTTTATCGGCTTCTTTTGGGTACTTGCCCCGATAGCTATAACGCTTGTCtt > 1:291673/1‑64 (MQ=255) cacTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTgg < 1:2181466/65‑1 (MQ=255) acTCTATTACTATTTTATCGGCTTTTTTTGGATACTTGCCGCGGTAGCTATAACGCTTGTCTGGc > 1:151728/1‑65 (MQ=255) acTCTATTACTATTTTATCGGCTTTTTTGGGGTACTTGCCCCGACAGGTATAACGCGTGTCTGGc > 1:3547257/1‑65 (MQ=255) acTCTATTACTATTTTATCGGCTTTTTTGGGCGACTTGCCCCCGAAACTATAACGCTTTTCTGGc > 1:1343943/1‑65 (MQ=255) acTCTATTACTATTTTATCGGCTTCTTTGGGGTACTTGCCGCGGTAGCTATAACGCTTGTCTggg > 1:3092493/1‑64 (MQ=255) acTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGc < 1:948076/65‑1 (MQ=255) tctATTACTATTTTATCGGCTTTTTTTGGATACTTGCCGCGATAGCTATAAAGCTTGTCTGGCaa > 1:3058848/1‑65 (MQ=255) tATTACTATTTTATCGGCTTTTTTTGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCaaaa > 1:2798593/1‑64 (MQ=255) tATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAg < 1:767457/65‑1 (MQ=255) tATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAg < 1:2164563/65‑1 (MQ=255) aTTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAgg < 1:3344384/65‑1 (MQ=255) aTTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAgg < 1:302165/65‑1 (MQ=255) aCTATTTTATCGGCTTTTTTTGGGTACTTGCTGCGATAGGTATAAAGCTTTTTTGGGAAAGGGAt > 1:1534174/1‑65 (MQ=255) aCTATTTTATCGGCTTTTTTTGGATATTTGGCGCGCTAGCTATAACGCTTGTAGGGGAAAGGGAt > 1:2017331/1‑65 (MQ=255) aCTATTTTATCGGCTTTTTTTGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAgggtg > 1:648521/1‑62 (MQ=255) tATTTTATCGGCTTTTTTTGGGTTGTTGTCGGGCTAGCTATAAAGCTTGTCTTGCAAAGGCATat > 1:3277782/1‑64 (MQ=255) tATTTTATCGGCTTTTTTTGGGTTCTTTCCCCGGTAACTATAAAGCTTTGTTGGCAAAAGAATAc > 1:1400405/1‑65 (MQ=255) tATTTTATCGGCTTTTTTGGGATACTTTCCGCGATAGCCATAAACATTGTCTTGCAAACGCATAc > 1:1434141/1‑65 (MQ=255) ttttATCGGCTTTTTTTGGGTACTTTTCGCGATAGCTATAAACCTTTTCTGGGCAAAGCAAAcgc > 1:1874272/1‑65 (MQ=255) tCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGtttt < 1:3154674/65‑1 (MQ=255) | TGTGGAGTCCGGTTCGTTGCGCTTTTACAGGTGCATTAACACTCTATTACTATTTTATCGGCTTCTTTGGGATACTTGCCGCGATAGCTATAACGCTTGTCTGGCAAAGGCATACGCGTTTT > NC_000913/821376‑821497 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |