| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,499,600 | 0 | T | G | 31.0% | 121.5 / 10.3 | 58 | G51G (GGT→GGG) | yhfL | DUF4223 domain‑containing lipoprotein YhfL |
| Reads supporting (aligned to +/- strand): ref base T (14/26); new base G (18/0); total (32/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.27e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAAT > NC_000913/3499552‑3499649 | ccTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGt < 1:4509757/49‑1 (MQ=255)ccTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTc < 1:2601956/50‑1 (MQ=255) ctgctgCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTcc < 1:4779173/50‑1 (MQ=255) ctgctgCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGGTGGGGGcc > 1:3733439/1‑50 (MQ=255) tgctgcACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTCCGGttc > 1:3386507/1‑47 (MQ=38) gctgcACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGGCCTa > 1:5501591/1‑50 (MQ=255) gctgcACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGGCCTa > 1:5564671/1‑50 (MQ=255) gctgcACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGGCCTa > 1:2707019/1‑50 (MQ=255) ctgcACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTa < 1:3001306/49‑1 (MQ=255) ctgcACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTAc < 1:5054809/50‑1 (MQ=255) tgcACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTTCGGTCCTACt > 1:1983457/1‑50 (MQ=255) tgcACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGGCCTACt > 1:1544649/1‑50 (MQ=255) tgcACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGGCCTACt > 1:1769054/1‑50 (MQ=255) gcACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTg > 1:959625/1‑50 (MQ=255) gcACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGGCCCACTg > 1:1292055/1‑50 (MQ=39) cACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGc < 1:3018832/50‑1 (MQ=255) cccGGCAATTTCTATTTCTAAAATCATTGGCGGTTTCGGGCCTACTGca > 1:1485613/1‑49 (MQ=38) cccGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGcac < 1:3322721/50‑1 (MQ=255) ccGGCAATTTCTATTTCTAAAATCATTGGCGGTTTCGGGCCTACTGca > 1:1764487/1‑48 (MQ=38) cGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGGCCTACTTCACAg > 1:1131610/1‑50 (MQ=39) gCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTa < 1:2905743/50‑1 (MQ=255) cAATTTCTATTTCTAAAATCATTGGCGGTTTCGGTCCTACTGCACAGTa > 1:2880328/1‑49 (MQ=255) cAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTa > 1:4889932/1‑49 (MQ=255) cAATTTCTATTTCTAAAATCATTGGCGGTTGCGGGGCTAATTCACAATa > 1:1422411/1‑49 (MQ=25) cAATTTCTATTTCTAAAATCATTGGCGGTTGCGGGCCTACTTCCCAGTa > 1:2627784/1‑49 (MQ=38) atttctatttctAAAATCATTGGCGGTTTCGGGCCTAATGCCCAGTaaag > 1:2204109/1‑48 (MQ=25) atttctatttctAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGt < 1:2819105/50‑1 (MQ=255) atttctatttctAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGt < 1:1432253/50‑1 (MQ=255) tttctatttctaAAATCATTGGCGGTTTCGGGCCTACTGCACAGTAAAtc > 1:1848291/1‑50 (MQ=38) ttctatttctaAAATCATTGGCGGTTGCGGGCCTACTTCACAGGAAGact > 1:4956193/1‑50 (MQ=37) ttctatttctaAAATCATTGGCGGTTGCGGGCCTACTGCACAGTAAGtct > 1:1569829/1‑50 (MQ=255) ctatttctaAAATCATTGGCGGTTTCGGGCCCACTGCCCAGTAAGTctct > 1:4056444/1‑50 (MQ=37) atttctaAAATCATTGGCGGTTGCGGTCCTACTGCCCAGGAAGTCTCTGa > 1:1656129/1‑50 (MQ=39) tttctaAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTg < 1:4161520/48‑1 (MQ=255) tttctaAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTg < 1:283212/48‑1 (MQ=255) tttctaAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGAc < 1:1271068/50‑1 (MQ=255) tctaAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACa < 1:4842128/49‑1 (MQ=255) aaaaTCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGaaa > 1:3837148/1‑50 (MQ=255) aaaTCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGa < 1:4075183/47‑1 (MQ=255) aaaTCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAAt < 1:1854510/50‑1 (MQ=255) aaaTCATTGGCGGTTGCGGGCCTACTGCACAGTAAGTCTCTGACAGAAAt > 1:1597593/1‑50 (MQ=255) aaTCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATc < 1:2726813/50‑1 (MQ=255) aaTCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATc < 1:1675003/50‑1 (MQ=255) aaTCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATc < 1:468139/50‑1 (MQ=255) aTCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCg < 1:356375/50‑1 (MQ=255) aTCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCg < 1:5349628/50‑1 (MQ=255) aTCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCg < 1:3646306/50‑1 (MQ=255) aTTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGc < 1:1595108/49‑1 (MQ=255) gCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAaca > 1:875972/1‑50 (MQ=255) gCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAac < 1:4320870/49‑1 (MQ=255) cGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAacac > 1:3047670/1‑50 (MQ=255) ggTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACAcc > 1:713111/1‑50 (MQ=255) gTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCg < 1:3526203/50‑1 (MQ=255) gTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCg < 1:864622/50‑1 (MQ=255) ttGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCgg < 1:503313/50‑1 (MQ=255) ttGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCgg > 1:1349304/1‑50 (MQ=255) tGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGt > 1:5379890/1‑50 (MQ=255) gCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGtt > 1:5566239/1‑50 (MQ=255) cGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGtt > 1:5586006/1‑49 (MQ=255) gTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTaa > 1:531743/1‑50 (MQ=255) tCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAAt < 1:1104536/50‑1 (MQ=255) tCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAAt > 1:2711292/1‑50 (MQ=255) tCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAAt > 1:490740/1‑50 (MQ=255) | CCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAAT > NC_000913/3499552‑3499649 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |