Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,371,314 | 0 | G | A | 56.0% | 4.4 / 26.8 | 25 | F475F (TTC→TTT) | menD | 2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylic‑acid synthase |
Reads supporting (aligned to +/- strand): ref base G (9/2); new base A (0/14); total (9/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.69e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATC > NZ_CP009273/2371249‑2371397 | aaaGGGACGTTTTGCGGCATCAGATAGAACCCCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTc < 2:1487084/90‑1 (MQ=255) ggACTTTTGGCGGCATAAAAAAGAACCGCCCCCGCTCCTTTTGCGGGGTTGGCAACAGCGAAAAAATTTCCCCGCCGTTGTTGTTcacca < 1:541184/90‑1 (MQ=255) ggACGTTTTGCGGCATCAGATAGAAACGCTCCCGCTCGCTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcacca < 2:1107715/90‑1 (MQ=255) ggACGTTTTGCGGCATCAGAAAAAACCCCCCCCCCTCGTTTTGCGGCGTTGGCAACAGGGAAAAAATTTGCCCGCCGTTGTTGTTcacca < 1:1328272/90‑1 (MQ=255) ggcGTTTTGCGGCATCAAATAAAACCGCCCCCGCTCGCTTTGCGGGGTTGGAAAAAGCGAAAAATTTTGCCCGCCGTTGTTGTTcaccac < 1:705123/88‑1 (MQ=255) gatttttcGCGGCATAAAATAAAACCGCCCCCGCTCGTTTGGCGGCGTTGGCAACAGCGAAAAATTTTGCCCGCCGTTGTTGTTcaccac < 1:1369936/86‑1 (MQ=255) gACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcaccac > 2:946498/1‑90 (MQ=255) gACGTTTTGCGGCATCAGATAAAACCCCCCCCCCTCGTTTGGGGGCGTTGGCAACAGGGAAAAAATTTGCCCGCCGTTGTTGTTcaccac < 1:93013/90‑1 (MQ=255) gACGTTTTGCGGCATCAGATAAAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTcaccac < 2:1050204/90‑1 (MQ=255) ctttttGCGGCATCAAATAAAACCGCTCCCGCTCGCTTTGCGGGGTTGGAAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACaa < 2:1168887/88‑1 (MQ=255) tttGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTa > 1:945914/1‑90 (MQ=255) tttGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTa > 1:1217561/1‑90 (MQ=255) aaCGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGaaa > 1:429511/1‑90 (MQ=255) aaCGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGaaa > 2:786757/1‑90 (MQ=255) acgctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAAc > 1:1483991/1‑90 (MQ=255) ctcacgctcGCTTTGCGGCGTTGGCAACAGCgaga > 1:679362/1‑35 (MQ=255) ctcacgctcGCTTTGCGGCGTTGGCAACAGCgaga < 2:679362/35‑1 (MQ=255) ctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcaccac > 1:459613/1‑63 (MQ=255) ctcacgctcGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTcaccac < 2:459613/63‑1 (MQ=255) cccgctcgttTTGGGGCGTGGGCAACGGCGAAAAAATTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGAc < 2:642212/88‑1 (MQ=255) ccgctcgttTGGCGGCGTGGCAAACCGAGAAAAATTTTCCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACg < 1:473677/89‑1 (MQ=255) ttGGGGTGTTGGCAACAGCGAAAAATTTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGc < 2:944842/90‑1 (MQ=255) gcggcgTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCa > 1:1584442/1‑90 (MQ=255) gcggcgTGGGCAACAGCGAAAAATTTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCa < 2:945914/90‑1 (MQ=255) gCAACAGCGAAAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTga < 2:968580/90‑1 (MQ=255) ccAACAGAGAAAAATTTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTga < 2:803283/89‑1 (MQ=255) cAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATc > 2:478086/1‑90 (MQ=255) | AAATGGACGTTTTGCGGCATCAGATAGAAACGCTCACGCTCGCTTTGCGGCGTTGGCAACAGCGAGAAAATTTGCCCGCCGTTGTTGTTCACCACAATTAATACCAGCGGCGCAGAAACCTGACGCAATAACGCCAGCGCGTTGAGATC > NZ_CP009273/2371249‑2371397 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |