Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,213,870 | 0 | A | C | 66.7% | 11.9 / 19.9 | 21 | T363P (ACC→CCC) | pdeG | cyclic diguanylate phosphodiesterase |
Reads supporting (aligned to +/- strand): ref base A (7/0); new base C (0/14); total (7/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.60e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTCATATTGGAATTAATTTTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACA > NZ_CP009273/1213786‑1213954 | tttCATATTGGAATTAATTTTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACccgc > 2:641499/1‑90 (MQ=255) tGGAATTAATTTTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGAttt > 1:818085/1‑90 (MQ=255) gaaattatttTTAGCGCTTCGCATATTATTCCGCCGACTTTTGTCGACGAGTTTTTAAA‑TTT‑CCGTGACATTTTTCCCCGCCGCGATTTa < 2:99659/88‑1 (MQ=255) gaaattatttTTAGCGCTTCGATTTTTTTTCCGCCGATTTTTGTCGACGTGTTTTTAAT‑TTT‑CCGTGAAATTTTTCCCCGCCGCGATTTa < 1:866854/88‑1 (MQ=255) taattTTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGATTTAAACCt > 2:3789/1‑90 (MQ=255) atttttAGCGCTTCGCATATTTTTTCGCCGCCGTTTGTCGAGGATTTTTTAAA‑TTC‑CCGTGACATTTTTCCCCGCCGCGATTTAAACCtt < 1:676936/88‑1 (MQ=255) ttAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTc > 1:263397/1‑90 (MQ=255) tCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTc > 1:337335/1‑90 (MQ=255) tttttCGCCGACGTTTGTCGACGAGTTTTTAAA‑TTT‑CCGTGACAGTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCg < 2:241905/88‑1 (MQ=255) ccGCCGACGTTTGTCGAGGAGT‑TTTAAA‑TTTCCCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAg < 2:266137/89‑1 (MQ=255) atgTTTGTGGAGGGGTTTTTAAT‑TTT‑CCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTg < 1:59741/88‑1 (MQ=255) aCGTTTGTCGACGATTTTTTAAT‑TTT‑CCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTg < 1:281588/90‑1 (MQ=255) ttGTCGACGAGTTTTTAAT‑TTT‑CCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATg < 2:797985/90‑1 (MQ=255) gTCGACGAGTGTTAAAA‑TTT‑CGGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGtt < 2:677100/90‑1 (MQ=255) tttAAATTTT‑CCGTGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTc < 2:486046/90‑1 (MQ=255) ttAAATTTT‑CGGTGAAATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCt < 2:589878/90‑1 (MQ=255) tt‑CCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAg > 1:589878/1‑90 (MQ=255) t‑CCGTAACTTTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGc < 1:326571/90‑1 (MQ=255) cccGTGACTTTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGc < 1:343123/89‑1 (MQ=255) tGACATTTTTCCCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGtt < 2:818085/90‑1 (MQ=255) gTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACa > 2:92976/1‑90 (MQ=255) | TTTCATATTGGAATTAATTTTAGCGCTTCGCATATTATTTCGCCGACGTTTGTCGACGAGTGTTTAAA‑TTT‑CCGTGACAGTTTTACCCGCCGCGATTTAAACCTTGTTCTGGAAGTCACCGAGCGTGAGCCATTGAATGTTGATGAAAGTCTGGTTCAGCGGTTGAACA > NZ_CP009273/1213786‑1213954 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |