| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 752,271 | 0 | A | C | 40.0% | 8.8 / 14.2 | 15 | G175G (GGT→GGG) | ybgD | putative fimbrial protein YbgD |
| Reads supporting (aligned to +/- strand): ref base A (9/0); new base C (0/6); total (9/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.00e-04 | |||||||||||
| 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. | |||||||||||
ATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGT > NC_000913/752208‑752334 | atgatgCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGTTACTTCAc > 1:2098825/1‑65 (MQ=255) cgccAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGt > 1:1095354/1‑65 (MQ=255) caaaaaTAATATTTATTTATAATCCAGGCCGTAGGTGGAGTTAGCGGTTCCTGCCCCTGCCGTGa < 1:485852/63‑1 (MQ=255) ataattCTTATTTATAATCCGGCACGTAGGTTGCGTTAGCGGTTACTCCCCCTGCCGTGACATCg < 1:255706/65‑1 (MQ=255) ataataCTTATTTAAAATCCAGCACGTAGGTTGCGTTGGCGGTTACTCCCCCTGCCGTGACATCg < 1:2632319/65‑1 (MQ=255) aCTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGc > 1:802115/1‑65 (MQ=255) aCTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGc > 1:493372/1‑65 (MQ=255) aCTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGc > 1:1360384/1‑65 (MQ=255) tttatAATCCACCGTGTAGGTTGCTTTAGGGGTTCCTTCCCCTGCCGTGACATCGCCTGCATTAt < 1:467504/65‑1 (MQ=255) tttatAATCCACCCCGTAGGTTGCGTTAGGGTTTACTCCCCCTGCCGTGACATCGACTGCATTAt < 1:370861/65‑1 (MQ=255) tttatAACCCAGCACGTAGTTTGCTTTAGCGTTTCCTTCCCCTGCCGTGACATCGACTGCATTAt < 1:1375373/65‑1 (MQ=255) ataATCCAGCCCTTAGTTTGCTTTAGCGTTTACTCCCCCTGCCGTGACATCGACTGCATTATCaa < 1:303936/65‑1 (MQ=255) ataATCCAGCACGTAGGTTGCGTTAGCGGTTACTCCCCCTGCCGTGACATCGACTGCATTATCaa < 1:1045523/65‑1 (MQ=255) ataACCCAGCACGTGGGTTGCGTTAGCGGTTCCTTCCCCTGCCGTGACATCGACTGCATTATCaa < 1:1500350/65‑1 (MQ=255) aaaaTCCAGCTCGGAGGTTGCTTTAGCGGTTCCTCCCCCTGCCGTGACATCGACTGCATTATCaa < 1:301666/63‑1 (MQ=255) aTCCAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAAttt > 1:3408353/1‑65 (MQ=255) ccAGCGCGTGGGTTGCTTTGGCGTTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGt < 1:126801/65‑1 (MQ=255) ccAGCACGTAGGTTGCGTTAGCGGTTACTTCCCCTGCCGTGACATCGACTGCATTATCAATTTGt < 1:52097/65‑1 (MQ=255) ccAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGt > 1:106348/1‑65 (MQ=255) ccAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGt > 1:2548510/1‑65 (MQ=255) ccAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGt > 1:1723138/1‑65 (MQ=255) ccAGCACGTAGGTTCCGTTAGCGTTTACTGCACCTGCCGTGACATCGACTGCATTATCAATTTGt < 1:2992674/65‑1 (MQ=255) aCGTGGGTGGCGTTAGCGTTTCCTCCCCCTGCCGTGACATCGACTGCATTATCAATTTGTtccat < 1:436254/65‑1 (MQ=255) aCGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTtccat > 1:3061445/1‑65 (MQ=255) aCGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTtccat > 1:3111714/1‑65 (MQ=255) aCGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTtccat > 1:3436992/1‑65 (MQ=255) gTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATTAATTTGTtcaatca > 1:1088497/1‑64 (MQ=255) gTTGCGTTAGCGGTTACTTCACCTGCCGGGACATCGACTGCATTATCAATTTGGTCCATCCAGGc > 1:834807/1‑65 (MQ=255) aGCTGTTCCTCCCCCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGt < 1:998726/65‑1 (MQ=255) aGCGTTTCCTCCCCCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGt < 1:566974/65‑1 (MQ=255) gCGTTTACTCCCCCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGtt < 1:123349/65‑1 (MQ=255) ggTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCa > 1:1141487/1‑65 (MQ=255) ggTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCa > 1:1722630/1‑65 (MQ=255) ttCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCt > 1:526032/1‑65 (MQ=255) cACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGt < 1:547621/65‑1 (MQ=255) cACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGt < 1:249882/65‑1 (MQ=255) | ATGATGCCGCCAAAATAATACTTATTTATAATCCAGCACGTAGGTTGCGTTAGCGGTTACTTCACCTGCCGTGACATCGACTGCATTATCAATTTGTTCCATCCAGGCGAAAAAGTTCAGCGTCTGT > NC_000913/752208‑752334 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |