Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 536,014 | 0 | T | G | 61.5% | 3.6 / 21.0 | 26 | intergenic (+50/‑7) | allB/ybbY | allantoinase AllB/uracil/xanthine transporter |
Reads supporting (aligned to +/- strand): ref base T (3/7); new base G (16/0); total (19/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.82e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACAC > NZ_CP009273/535928‑536094 | gcCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGgggtg > 2:504059/1‑90 (MQ=255) gAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGCTTAGGGGGGttt > 2:303102/1‑90 (MQ=255) tCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTc < 1:23510/90‑1 (MQ=255) aaTTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCaa < 2:24575/90‑1 (MQ=255) tAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTc < 1:297072/90‑1 (MQ=255) aaaCATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTAAGGGGTGTTTTTGTTCAATTTTTCAGGCa > 2:297072/1‑90 (MQ=255) tCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGGTAAGGGGGGTTTTTTTTTTATTTTTTAGTCAGCcgc > 1:62731/1‑90 (MQ=255) cagcagTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGATTAAGGGGGGTTTTTTTTCAATTTTTCAGTCAGCcgcg > 2:557525/1‑90 (MQ=255) gCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCGGTTTAGGGGGGTTTTTGTTCAATTTTTCAGTCAGCCGCGAAAACCTGTTTTc > 1:233595/1‑90 (MQ=255) gCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTTCAGTCAGCCGCGAAAGCCCGTTATc > 1:317089/1‑90 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTTTCCCGTTAAGGGGTGTTTTTTTTCAATTTTTCAGGCAGCCGCGGAAGCCCGTTATCAggg > 1:365109/1‑89 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTTAGGGGGGTTTTTGTTCAATTTTTCAGTCAGCCCCGAAAGCCTGTTTTTAggg > 2:379324/1‑89 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTTAGGGGGGTTTTTGTTTAATTTTTCAGTCAGCCGCGAAAGCCTGTTTTCAggg > 2:482417/1‑89 (MQ=255) cTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGGGTTTTTGTTTAATTTTTCAGTCAGCCGCGAAAGCCTGTTATCAGGa > 2:382760/1‑90 (MQ=255) aTGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGGGGGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCa > 2:15587/1‑90 (MQ=255) gCCCGTCCTTGCGGCGGGCATTTTCCGGGTTAGGGGTGTTTTTGTTTAATTTTTCAGTCAGCCGCGGAAGCCTGTTTTCAGGGTTTTAGt > 1:37009/1‑90 (MQ=255) gCCCGTCCTTGCGGCGGGCATTCTCCGGGTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGGTTTCAGt > 2:339700/1‑90 (MQ=255) gTCCTTGCGGCGGGCATTCTCCGGGTAAGGGGGGTTTATGTTCAATTTTTCAGTCAGCCGCGAAAACCCGTTTTTAGGATTTTAGTGGtt > 1:115488/1‑90 (MQ=255) tCCTTGCGGCGGGCATTCTCCGGTTAAGGGGGGTTTATGTTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTATCAGGGTTTTAGTGGgtt > 2:361280/1‑90 (MQ=255) gggCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTTTCAGGATTTCAGGGGTTTTTTTTTAttt > 2:297108/1‑90 (MQ=255) gggCATTCTCCGGTTAAGGTGTGTTTATGTTCAATCTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTTTTTAttt > 2:356434/1‑90 (MQ=255) gggCATTCTCCGGTTAAGGGGTGTTTATGTTCAATTTTTCAGTCAGCCGCGAAAGCCTGTTTTTAGGGTTTTTGGGGGTTTTTTTTAttt > 1:515500/1‑90 (MQ=255) gggCATTCTCCGGGTAAGGGGTGTTTATGTTTAATTTTTCAGTCAGCCGCGAAAGCCTGTTATCAGGGTTTTTGTGGGtttttttttttt > 1:149026/1‑90 (MQ=255) ttCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGc < 1:303102/90‑1 (MQ=255) tcCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAac < 1:585102/90‑1 (MQ=255) cGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac < 1:482417/90‑1 (MQ=255) cGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAacac < 2:475261/90‑1 (MQ=255) | GCCGAAAGGTCAATTTATCCTTAAACATCAGCAGTAATCTGGCCCCTGCAATGCCCGTCCTTGCGGCGGGCATTCTCCGGTTAAGGTGTGTTTATGTTCAATTTTGCAGTCAGCCGCGAAAGCCTGTTATCAGGATTTCAGTGGTTTTTCTTTATTTTTTGCAACAC > NZ_CP009273/535928‑536094 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |