| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 585,390 | 0 | T | G | 51.8% | 1.9 / 24.5 | 27 | T341P (ACC→CCC) | nfrA | bacteriophage adsorption protein NfrA |
| Reads supporting (aligned to +/- strand): ref base T (0/13); new base G (14/0); total (14/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.99e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-05 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAA > NZ_CP009273/585314‑585479 | cGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGtttt < 1:675796/90‑1 (MQ=255) cAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGAGGCGGttttt > 1:723977/1‑90 (MQ=255) agagCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGG‑GCAGGTTTTgc > 2:369156/1‑90 (MQ=255) agCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTgcgc < 2:1748608/90‑1 (MQ=255) agCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTgcgc < 2:1331913/90‑1 (MQ=255) cTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTg < 1:1636326/90‑1 (MQ=255) cAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGcgt < 1:1682407/90‑1 (MQ=255) cAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGcgt < 1:852828/90‑1 (MQ=255) tGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGG‑GGAGGTTTTTTGCTGCGTCGGACt > 1:1564264/1‑90 (MQ=255) gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGG‑GGAGGTTTTTTGCTGCGTCGTACTGACCCTCtttt > 1:145805/1‑90 (MQ=255) gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGG‑GCAGGTTTTTTGCTGCGTCGTACTGACCTTCtttt > 1:1367752/1‑90 (MQ=255) gCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGCGG‑GGAGGTTTTTGGGTGGGTCGTAGTTACCTTCcttt > 2:1058317/1‑90 (MQ=255) cTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGc < 1:273786/90‑1 (MQ=255) aTAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGG‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTTAAAGCACCGGTAg > 1:1364234/1‑90 (MQ=255) tAACGCTCCTCAAGCATTTCATTGGCGGGGGGGGGGGGGG‑GCAGGTTTTTTGCTGCGGCGTACTTACCTTCTTTTTACAACACCCGtctc > 1:448936/1‑87 (MQ=255) cctcAAGCATTTCATTGGCGGGGGGGGGGGGGG‑GCCGGTTTTTTGCCGCGGCCTTCTTACCCTTTTTTTAAAGCACCGGGAGCGTCgggc > 2:491437/1‑90 (MQ=255) cctcAAGCATTTCATTGGCGGGGGGGGGGGGGA‑GGAGTTTTTTTGCTGCGTCGTACTGACCTTTTTTTTAAAACACCCGGAGCGgcgcgc > 1:236714/1‑90 (MQ=255) aaGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCcaac < 2:238645/90‑1 (MQ=255) aGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCcaaca < 2:924964/90‑1 (MQ=255) cATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCaacaac < 1:315104/90‑1 (MQ=255) aTTTCATTGGCGGGGGGGGGGGGGG‑GCAGGTTTTTTGCCCCGCCGGGCTGTCCCTTTTTTATATGCACCGCTGGCGGCGCGCCAacaaca > 2:843388/1‑90 (MQ=255) ttCATTGGCGGGGGGGGGGGGGG‑GGAGGTTTTTCGCCGCGGCGTACTGACCTTCTTTTTAAAGCACCCGGAGGGTCGCGCCCACAACATa > 2:109661/1‑90 (MQ=255) ttCATTGGCGGGGGGGGGGGGGA‑GCAGTTTTTTCGCTGCGTCGTACTTACCTTCTTTTTACAACCCCGGGAGCGGCGCGCCCacaacaaa > 1:238645/1‑88 (MQ=255) ggCGGGGGGGGGGGGGA‑GCAGTTTTTTCGCTGCGTCGTACTGACCTTCTTTTTACAGCACCGGGAGCGGCGCGCCCACAACAAACTGGCg > 1:31617/1‑90 (MQ=255) gggggggggggggg‑GCGGGTTTTTCGCTGCGTCGTACTGACCTTCTTTTTACAGCACCGGGAGGGGGGCGCCAACAACATACTGGCgggg > 2:610901/1‑88 (MQ=255) gggTGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCgg < 1:1702169/90‑1 (MQ=255) gggTGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCgg < 1:1468224/90‑1 (MQ=255) gggTCGCGA‑GCAGTTTTTGCGCTGCGCCGTAGTGGCCTTCTTTGAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCgg < 2:1553522/90‑1 (MQ=255) gtgtCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCa < 1:27739/90‑1 (MQ=255) tGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCaa < 2:723977/90‑1 (MQ=255) | CGCAGAGCTTCAGCCTTGTTACGGGTCGCCACGCTGACAGCATAACGCTCCTCAAGCATTTCATTGGCGGGGAGGGTGGCGA‑GCAGTTTTTGCGCTGCGTCGTACTGACCTTCTTTTAACAGCACCGGTAGCGTCGCGCCAACAACATACTGGCGGTTGTCGGCAA > NZ_CP009273/585314‑585479 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |