Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,020,300 | 0 | A | C | 54.5% | 8.6 / 13.3 | 11 | E391D (GAA→GAC) | guaD | guanine deaminase |
Reads supporting (aligned to +/- strand): ref base A (5/0); new base C (0/6); total (5/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.16e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGA > NZ_CP009273/3020211‑3020370 | gCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGaa > 1:44159/1‑90 (MQ=255) cGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCa > 2:184008/1‑90 (MQ=255) ggcggAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGcc > 2:292571/1‑90 (MQ=255) ttGACGTTTTGATTGGCATTTTTTTCCCTGGAAAAGGGGCTGATTTGGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCt < 1:293158/90‑1 (MQ=255) tGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTa > 1:430790/1‑90 (MQ=255) gATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGAc > 2:513691/1‑90 (MQ=255) tttCCCTGGAAAAGAGGCTTTTTTCGTGGTGATGGACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGt < 2:74865/90‑1 (MQ=255) gaaaaGGGGCTGTTTTCGTGGTGATGAACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAAt < 1:444486/88‑1 (MQ=255) agagGCTGATTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGtt < 1:184008/90‑1 (MQ=255) ggCTGTTTTGGGGGTGAGGAACCCCCCCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGt < 2:216172/90‑1 (MQ=255) gCTGTTTTCGTGGTGATGGACCCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtg < 1:129662/90‑1 (MQ=255) cTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGtga > 2:46796/1‑90 (MQ=255) | GCCACGCTCGGCGGAGCGAAATCTCTGGGCCTTGACGATTTGATTGGCAACTTTTTACCTGGCAAAGAGGCTGATTTCGTGGTGATGGAACCCACCGCCACTCCGCTACAGCAGCTGCGCTATGACAACTCTGTTTCTTTAGTCGACAAATTGTTCGTGA > NZ_CP009273/3020211‑3020370 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |