| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 899,438 | 0 | A | C | 61.9% | 13.5 / 21.0 | 21 | G162G (GGT→GGG) | ybjP | lipoprotein |
| Reads supporting (aligned to +/- strand): ref base A (7/1); new base C (0/13); total (7/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.88e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATC > NZ_CP009273/899356‑899524 | aTTCCTCACGAAATGCCGGACATTTTACGGGTTTTATTGGTTGATCAGGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTCCccgccg < 1:148854/90‑1 (MQ=255) ttCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGt > 1:446138/1‑90 (MQ=255) gCCGGACATTTTAGGGGGTTTTTTGGTTGATCAAGGCGTTAGCGTTTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTa < 1:646639/90‑1 (MQ=255) gCCGGAAATTTTGGGGGTTTTTTTGGTTAACCAGGGCTTTAGCGTTTCTCGGTGGACAGAGGGAGCGTCCCCGCCGTGGCATGAACGCTa < 2:446138/90‑1 (MQ=255) ggACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCg > 1:628363/1‑90 (MQ=255) attttGCGGGGTTTATTGTTTGACAAGGGGGTTAGCGTTTCCCGAGGGACTGACGGAGCTCCCCCGCCGTGGCATGAACGCTACCGCCTa < 2:461677/88‑1 (MQ=255) tttttttGGTTGATAAGGGCGTTAGCGTTTCTCGATGGACGGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGcac < 2:284837/89‑1 (MQ=255) tttttGGTTAATCAGGGCGTTAGCATTTCCGGATGGACTGACGGGGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACAt < 2:427486/87‑1 (MQ=255) ttATTGGTTGATCAGGGCTTTAGCGTTTCGCGAGGGACTGAGGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACAt < 2:672463/90‑1 (MQ=255) ttATTGGTTGATAAAGGCGTTAGCGATTCTCGAGGGACGGAGGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACAt < 1:394956/90‑1 (MQ=255) ttttGGTTGATCAGGGCTTTAGCGTTCCCCGAGGGACGGCCGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACAtc < 2:574372/88‑1 (MQ=255) tATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACAtc > 1:432280/1‑90 (MQ=255) gTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGa > 2:27018/1‑90 (MQ=255) ttGATAAGGGCTTTAGCGTTTCTCGAGGGACTGAGGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGAt < 1:665037/90‑1 (MQ=255) ttGATAAAGGCGTACGCGTTTCCCGATGGACTGAGGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGAt < 1:468947/90‑1 (MQ=255) tAGCGATTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTgg < 1:18889/90‑1 (MQ=255) gCGTTTCTCGATGGACTGAGGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGcc < 1:222069/90‑1 (MQ=255) tGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAac > 2:194455/1‑90 (MQ=255) ggAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACtt > 1:446514/1‑90 (MQ=255) gCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCAt > 2:386880/1‑90 (MQ=255) cGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATc < 2:446514/90‑1 (MQ=255) | ATTCCTCACGAAATGCCGGACAATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATC > NZ_CP009273/899356‑899524 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |