Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,213,865 | 0 | G | T | 55.0% | 14.2 / 32.6 | 20 | S361I (AGT→ATT) | pdeG | cyclic diguanylate phosphodiesterase |
Reads supporting (aligned to +/- strand): ref base G (9/0); new base T (0/11); total (9/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.95e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGGAATTAATTTTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTT > NZ_CP009273/1213794‑1213949 | tGGAATTAATTTTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGAttt > 1:619417/1‑90 (MQ=255) tGGAATTAATTTTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGAttt > 1:870948/1‑90 (MQ=255) tatttttAGCGTTTCGCATATTTTTTCGCCGACTTTTGTCGACGAGTTTTTAAA‑TTT‑CCGTGACATTTTTACCCGCCGCGATTTAAACCt < 2:681954/87‑1 (MQ=255) ttAGCGTTTCGCATTTTTTTTCGCCGGCTTTTGTCGACGAGT‑TTTAAA‑TTTCCCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTc < 1:745587/90‑1 (MQ=255) cTTCGCATTTTTTTTCGCCGCGTTTTGTCGACGAGTGTTAAAT‑TTT‑CCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAg < 1:682156/90‑1 (MQ=255) cATTTTTTTTCGCCGTCTTTTGTCGACGAGTTTTTAAT‑TTC‑CCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCAcc < 1:1345965/90‑1 (MQ=255) tttCGCCGCCTTTTGTCGACGAGTGTTTAAT‑TTT‑CCGTGACATTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTg < 2:796852/90‑1 (MQ=255) ttCGCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGa > 2:1383125/1‑90 (MQ=255) tCGCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAg > 1:564557/1‑90 (MQ=255) tCGCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAg > 2:1203829/1‑90 (MQ=255) cGCCGGCTTTTGTCGAGGAGTTTTTAAT‑TTT‑CCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGc < 2:870948/90‑1 (MQ=255) cGCCGACGTTTGTCGACGTTTTTTTAAT‑TTC‑CCCAGATATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGc < 1:1418342/90‑1 (MQ=255) gCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGcc > 1:685631/1‑90 (MQ=255) ttttGTCGACGAGTTTTTAAA‑TTT‑CCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGaa < 2:1298/89‑1 (MQ=255) gaGTTTTTAAA‑TTT‑CCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGaa < 2:67458/90‑1 (MQ=255) gatttttAAATTTT‑CCGGGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGaa < 2:869212/86‑1 (MQ=255) tAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTgg > 1:408729/1‑90 (MQ=255) ccGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCg > 1:1187626/1‑90 (MQ=255) ccGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCg > 1:1178534/1‑90 (MQ=255) tGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGtt < 2:361982/90‑1 (MQ=255) | TGGAATTAATTTTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTT > NZ_CP009273/1213794‑1213949 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |