| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 585,405 | 0 | G | T | 60.0% | 3.8 / 22.5 | 15 | Q336K (CAA→AAA) | nfrA | bacteriophage adsorption protein NfrA |
| Reads supporting (aligned to +/- strand): ref base G (0/6); new base T (9/0); total (9/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.00e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-04 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGT > NZ_CP009273/585319‑585489 | agCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTgcgc < 2:592572/90‑1 (MQ=255)agCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTgcgc < 2:197102/90‑1 (MQ=255) tCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTTGCTGcg > 2:153959/1‑90 (MQ=255) tCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGTTTTTTCGCcgcg > 1:478062/1‑90 (MQ=255) tCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGGTTTTTTGCTGcg > 1:299861/1‑90 (MQ=255) cAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGcgt < 1:543209/90‑1 (MQ=255) cGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACa < 2:299861/90‑1 (MQ=255) aCAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGGGGAGGTTTTTCCCTGGGCTGGCCTTAACTTCTTTTTACAGCaac > 2:330498/1‑88 (MQ=255) aGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGGGGGGGGCAGGTTTTTCGCTGCGGCGTACTGACCTTCTTTTTACACCCCCgg > 1:337625/1‑90 (MQ=255) cctcAAGCATTTCATTGGCGGGGGGGGGGGGGGGCAGGTTTTTTTCTGCGGCGTTCTGACCTTCTTTTTACAACACCGGGAGCGgcgcgc > 2:397343/1‑90 (MQ=255) cctcAAGCATTTCATTGGCGGGGGGGGGGGGGAGCAGGTTTTTTGCTGCGTCCTTCTTGCCTTCTTTTTTAAGGACCCGGAGCGTCgcgc > 2:15772/1‑90 (MQ=255) ctcAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGcc < 2:571972/90‑1 (MQ=255) tcAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGATCTTCTTTTAACAGCACCGGTAGCGTCGCGCca < 2:147239/90‑1 (MQ=255) ttCATTGGCGGGGGGGGGGGGGGGCAGGTTTTTCGCTGCGTCGTACTGACCTTCTTTTTAAAACAACGGGAGCGTCGCGCCCacaacaaa > 1:221442/1‑88 (MQ=255) ggggggggggggggGCAGGTTTTTTGCTGCGGCGTACTGACCTTCTTTTTACAGCACCGGGGGCGGCGCGCCCACCAAAAACTGgggggg > 2:166566/1‑88 (MQ=255) gCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAAc < 1:153959/90‑1 (MQ=255) tttttGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGt < 1:552889/90‑1 (MQ=255) | AGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGAGCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAAACTGTACCGT > NZ_CP009273/585319‑585489 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |